Loj#2769-「ROI 2017 Day 1」前往大都会【最短路树,斜率优化】

正题

题目链接:https://loj.ac/p/2769


题目大意

给出nnn个点mmm条地铁线路,每条线路是一条路径。

111nnn的最短路且在最短路径的情况下相邻换乘点的距离平方和最大。

1≤n,m,∑si≤1061\leq n,m,\sum s_i\leq 10^61n,m,si106


解题思路

首先肯定是在最短路树上跑,然后考虑怎么求平方和最大。

因为每条可以计算和的路径一定是连续的一段,所以考虑将一条铁路线去掉不在最短路树上的边后分出若干段来。

对于每一段上的转移都是一个fx=max{fy+(dx−dy)2}f_{x}=max\{f_y+(d_x-d_y)^2\}fx=max{fy+(dxdy)2}的形式,可以斜率优化。

但是为了充分转移,需要按照disdisdis从小到大转移。

栈啊什么的开vectorvectorvector就好了。

时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<cctype>
#define mp(x,y) make_pair(x,y)
#define sec second
#define fir first
#define ll long long
using namespace std;
const ll N=1e6+10;
struct node{ll to,next,w;
}a[N];
ll n,m,T,tot,ls[N],dis[N],f[N],top[N],p[N],las[N];
vector<pair<ll,ll> >e[N],E[N],pos[N];
priority_queue<pair<ll,ll> >q;
vector<ll> _x[N],_y[N],s[N];bool v[N];
ll read(){ll x=0,f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=-f;c=getchar();}while(isdigit(c)){x=(x<<1)+(x<<3)+c-48;c=getchar();}return x*f;
}
void addl(ll x,ll y,ll w){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;a[tot].w=w;return;
}
void dij(){memset(dis,0x3f,sizeof(dis));q.push(mp(0,1));dis[1]=0;while(!q.empty()){ll x=q.top().sec;q.pop();if(v[x])continue;v[x]=1;for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(dis[x]+a[i].w<dis[y]){dis[y]=dis[x]+a[i].w;q.push(mp(-dis[y],y));}}}return;
}
bool cmp(ll x,ll y)
{return dis[x]<dis[y];}
double slope(ll x,ll y,ll p)
{return (double)(_y[p][y]-_y[p][x])/(_x[p][y]-_x[p][x]);}
void calc(ll x,ll p,ll t,vector<ll> &s,ll &top){//x点 p线 t站if(!t)return;ll dx=e[p][t].sec;while(top>1&&slope(s[top-1],s[top],p)<=2*dx)top--;if(top){ll y=s[top];f[x]=max(f[x],_y[p][y]-2*dx*_x[p][y]+dx*dx);}
}
void add(ll x,ll p,ll t,vector<ll> &s,ll &top){//x点 p线 t站if(t>=e[p].size()-1)return;ll dx=e[p][t].sec;_y[p][t]=f[x]+dx*dx;_x[p][t]=dx;while(top>1&&slope(s[top-1],t,p)>=slope(s[top-1],s[top],p))top--;s[++top]=t;las[p]=t;return;
}
signed main()
{
//	freopen("city.in","r",stdin);
//	freopen("city.out","w",stdout);n=read();m=read();for(ll i=1;i<=m;i++){ll num,x,t=0;num=read()+1;while(num--){if(E[i].size())t=read();x=read();if(E[i].size()){pair<ll,ll> z=E[i][E[i].size()-1];addl(z.fir,x,t);}//建边 E[i].push_back(mp(x,t));}}dij();for(int i=1;i<=m;i++){for(int j=0;j<E[i].size();j++){int x=E[i][j].first,t=E[i][j].sec,flag=1;if((!j)||dis[E[i][j-1].fir]+E[i][j].sec!=dis[E[i][j].fir]){if(T&&e[T].size()==1){pos[e[T][0].first].pop_back();e[T][0].first=x;flag=0;}else T++,e[T].push_back(mp(x,0));}else e[T].push_back(mp(x,t));if(flag)s[T].push_back(0),_x[T].push_back(0),_y[T].push_back(0);pos[x].push_back(mp(T,e[T].size()-1));//e:线路(点/和前面的长度)//pos:点(线路,对应位置)//s,_x,_y:dp用数组 }}for(ll i=1;i<=T;i++)for(ll j=1;j<e[i].size();j++)e[i][j].sec+=e[i][j-1].sec;for(ll i=1;i<=n;i++)p[i]=i;sort(p+1,p+1+n,cmp);for(ll _i=1;_i<=n;_i++){ll x=p[_i];for(ll _j=0;_j<pos[x].size();_j++)calc(x,pos[x][_j].fir,pos[x][_j].sec,s[pos[x][_j].fir],top[pos[x][_j].fir]);for(ll _j=0;_j<pos[x].size();_j++)add(x,pos[x][_j].fir,pos[x][_j].sec,s[pos[x][_j].fir],top[pos[x][_j].fir]);if(x==n)break;}printf("%lld %lld\n",dis[n],f[n]);return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318655.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

牛客题霸 [两个链表生成相加链表] C++题解/答案

牛客题霸 [两个链表生成相加链表] C题解/答案 题目描述 假设链表中每一个节点的值都在 0 - 9 之间&#xff0c;那么链表整体就可以代表一个整数。 给定两个这种链表&#xff0c;请生成代表两个整数相加值的结果链表。 例如&#xff1a;链表 1 为 9->3->7&#xff0c;链…

2021牛客暑期多校训练营8 J-Tree(思维dp+rmq)

J-Tree Cosmic_Tree题解 赛时队友想到了预处理s→ts\to ts→t路径上的点最远到达的点&#xff0c;后面贪心做的一直wa。正解是dp。 下面的solve(l,r,op)\text{solve}(l,r,\text{op})solve(l,r,op)函数中&#xff0c;当lrlrlr时不满足题意不合法&#xff0c;当op0\text{op}0op…

[翻译] ASP.NET Core 2.2 正式版发布

本文为翻译&#xff0c;原文地址&#xff1a;https://blogs.msdn.microsoft.com/webdev/2018/12/04/asp-net-core-2-2-available-today/我&#xff08;文章作者&#xff09;很高兴地宣布ASP.NET Core 2.2现在作为.NET Core 2.2的一部分提供&#xff01;怎么获取它您可以从 .NET…

理解至上:数位dp(ybtoj-B数计数)

文章目录简要题目描述解析dp定义:试填法代码thanks for reading&#xff01;简要 数位dp&#xff0c;天下第一 最重要的应该有两个&#xff1a; 1.状态转移式的确定 2.试填法不断往后模拟 &#xff08;至今是唯一一道数位dp&#xff0c;究竟重要的是啥我其实也没有太多经验 &am…

P7408-[JOI 2021 Final]ダンジョン 3【贪心,树状数组】

正题 题目链接:https://www.luogu.com.cn/problem/P7408 题目大意 一个有n1n1n1层的地牢&#xff0c;从iii到i1i1i1层要AiA_iAi​点能量&#xff0c;第iii层可以花费BiB_iBi​获得111点能量。 mmm次询问从SiS_iSi​层出发到第TiT_iTi​层在能量上限为UiU_iUi​的情况下至少需…

牛客题霸 [输出二叉树的右视图] C++题解/答案

牛客题霸 [输出二叉树的右视图] C题解/答案 题目描述 请根据二叉树的前序遍历&#xff0c;中序遍历恢复二叉树&#xff0c;并打印出二叉树的右视图 题解&#xff1b; 分两个过程&#xff1a; 先用前序遍历中序遍历恢复二叉树&#xff0c;这个应该都会。。 打印二叉树的后视…

2021“MINIEYE杯”中国大学生算法设计超级联赛(10)Pty loves string(Border+二维数点)

Pty loves string 建立Border树后&#xff0c;发现可以转化成两个子树中相同点的数量&#xff0c;时间戳转化为连续的区间后相当于有两个数组&#xff0c;每次给两个区间&#xff0c;问区间相同点权的数目。 第一个数组作为区间&#xff0c;第二个数组作为权值。将第一个数组建…

微软推出了Cloud Native Application Bundles和开源ONNX Runtime

微软的Microsoft Connect(); 2018年的开发者大会 对Azure和IoT Edge服务进行了大量更新; Windows Presentation Foundation&#xff0c;Windows Forms和Windows UI XAML Library的开源 以及.NET 基金会会员模型的扩展。但那些只是冰山一角。微软还联合Docker发布了Cloud Native…

不止代码:友好城市(动态规划)

解析 先按左端点排序得到一个右端点的新队列&#xff0c;然后就可以发现&#xff1a; 所有合法的方案都是新队列的一个单调递增队列 然后就转化成了最长上升序列的问题 代码 #include<bits/stdc.h> using namespace std; const int N1e6100; int m,n; struct node{int…

二分图匹配(一)

文章目录什么是二分图&#xff1a;例题&#xff1a;NC111768 CF741C题目描述&#xff1a;题解&#xff1a;代码&#xff1a;二分图最大匹配匈牙利算法算法思想&#xff1a;代码&#xff1a;Knig定理二分图最优匹配KM(Kuhn-Munkres)算法算法思路&#xff1a;具体操作代码&#x…

P3980-[NOI2008]志愿者招募【费用流】

正题 题目链接:https://www.luogu.com.cn/problem/P3980 题目大意 nnn天&#xff0c;第iii天需要AiA_iAi​个志愿者。有mmm种志愿者&#xff0c;第iii种从sis_isi​天服务到tit_iti​天&#xff0c;需要cic_ici​元的费用。 求满足条件的最小费用 1≤n≤1000,1≤m≤100001\le…

不止代码:恐狼后卫(ybtoj-区间dp)

文章目录题目描述解析代码thanks for reading!题目描述 一代炉石的眼泪啊 解析 用dp[i][j]表示i与j之间的全部消掉&#xff08;不含两端&#xff09;的最小花费 然后枚举中间最后杀死的狼就行了 本题没有一次AC&#xff0c;因为一开始dp定义成了包含两端&#xff0c;然后因为…

牛客题霸 [二叉树的之字形层序遍历] C++题解/答案

牛客题霸 [二叉树的之字形层序遍历] C题解/答案 题目描述 给定一个二叉树&#xff0c;返回该二叉树的之字形层序遍历&#xff0c;&#xff08;第一层从左向右&#xff0c;下一层从右向左&#xff0c;一直这样交替&#xff09; 例如&#xff1a; 给定的二叉树是{3,9,20,#,#,15…

P4180-[BJWC2010]严格次小生成树【Kruskal,倍增】

正题 题目链接:https://www.luogu.com.cn/problem/P4180 题目大意 nnn个点mmm条边的一张无向图&#xff0c;求它的严格次小生成树。 1≤n≤105,1≤m≤31051\leq n\leq 10^5,1\leq m\leq 3\times 10^51≤n≤105,1≤m≤3105 解题思路 一定存在一种严格次小生成树和最小生成树只…

【.NET Core项目实战-统一认证平台】第九章 授权篇-使用Dapper持久化IdentityServer4...

上篇文章介绍了IdentityServer4的源码分析的内容&#xff0c;让我们知道了IdentityServer4的一些运行原理&#xff0c;这篇将介绍如何使用dapper来持久化Identityserver4&#xff0c;让我们对IdentityServer4理解更透彻&#xff0c;并优化下数据请求&#xff0c;减少不必要的开…

不止代码:生日欢唱(ybtoj-区间dp)

文章目录题目描述解析代码thanks for reading&#xff01;题目描述 解析 这题挺好的 思路&#xff1a;dp[i][j]表示必须把i和j配对&#xff0c;可达到的最大值 首先&#xff1a; dp[i][j]dp[i-1][j-1]a[i]*b[j];然后可以分别尝试把男生或女生往前放弃一段&#xff1a; for(i…

[BZOJ5312]冒险(势能线段树)

[BZOJ5312]冒险 维护一个长度为 n 的序列&#xff0c;支持 m 次操作&#xff0c;操作包括区间按位或一个数&#xff0c;区间按位与一个数&#xff0c;以及查询区间最大值。 线段树每个节点上维护区间与、区间或和区间最大值。 如果一次操作对区间与的影响和对区间或的影响相同&…

牛客题霸 [表达式求值] C++题解/答案

牛客题霸 [表达式求值] C题解/答案 题目描述 请写一个整数计算器&#xff0c;支持加减乘三种运算和括号。 题解&#xff1a; 没有除法emmm 我们从头开始依次判断每个字符 如果是左括号&#xff0c;我们就找右括号&#xff0c;并截取括号内的数字 记录上一次的符号&#xff…

P6085-[JSOI2013]吃货JYY【状压dp,欧拉回路】

正题 题目链接:https://www.luogu.com.cn/problem/P6085 题目大意 nnn个点的一张无向图&#xff0c;有kkk条必走边&#xff0c;mmm条其他边&#xff0c;求从111出发经过必走边后回到起点的最短路径。 2≤n≤13,0≤k≤78,2≤m≤2002\leq n\leq 13,0\leq k\leq 78,2\leq m\leq 2…

打造自己的.NET Core项目模板

前言每个人都有自己习惯的项目结构&#xff0c;有人的喜欢在项目里面建解决方案文件夹&#xff1b;有的人喜欢传统的三层命名&#xff1b;有的人喜欢单一&#xff0c;简单的项目一个csproj就搞定。。反正就是萝卜青菜&#xff0c;各有所爱。可能不同的公司对这些会有特定的要求…