[HNOI 2014]道路堵塞
2017-07-27 17:37
260 查看
Description
A国有N座城市,依次标为1到N。同时,在这N座城市间有M条单向道路,每条道路的长度是一个正整数。现在,A国 交通部指定了一条从城市1到城市N的路径,并且保证这条路径的长度是所有从城市1到城市N的路径中最短的。不幸的是,因为从城市1到城市N旅行的人越来越 多,这条由交通部指定的路径经常发生堵塞。现在A国想知道,这条路径中的任意一条道路无法通行时,由城市1到N的最短路径长度是多少。Input
输入文件第一行是三个用空格分开的正整数N、M和L,分别表示城市数目、单向道路数目和交通部指定 的最短路径包含多少条道路。按下来M行,每行三个用空格分开的整数a、b和c,表示存在一条由城市a到城市b的长度为c的单向道路。这M行的行号也是对应 道路的编号,即其中第1行对应的道路编号为1,第2行对应的道路编号为2,...,第M行对应的道路编号为M。最后一行为L个用空格分开的整数 sp(1)...,,sp(L),依次表示从城市1到城市N的由交通部指定的最短路径上的道路的编号。Output
输出文件包含L行,每行为一个整数,第i行(i=1,2...,,L)的整数表示删去编号为sp(i)的道路后从城市1到城市N的最短路径长度。如果去掉后没有从城市1到城市N的路径,则输出一1。Sample Input
4 5 2 1 2 2 1 3 2 3 4 4 3 2 1 2 4 3 1 5
Sample Output
6 6
Hint
100%的数据满足2<N<100000,1<M<200000。所用道路长度大于0小于10000。题目大意
给出一个N个节点,M条边的有向图。给出一条1到N的最短路,问最短路上任意一条边断掉,此时的最短路是多少。题解
考虑不可能每次删掉当前边之后再跑一遍最短路,那必须想办法优化删边和求得最短路过程。我们将最短路上的点从小到大标号。
考虑删边之后,当前的从$1$到$N$的最短路只能是从$1$出发,在最短路径上走一段,再走一段非最短路,最后回到最短路径上。那么如果强制不走当前边,在跑最短路的过程中,只要发现走到了一个最短路径上的点上时(这个点肯定是最短路上标号更大的点),既然在最短路上,就不必要继续松弛了,显然$SPFA$求出的$dis$值+该点到终点的距离$last$就是一个可行的解,更新$ans$。
其实上述过程得到的解还是没有充分利用的。我们会发现,以后的$SPFA$搜到解会和当前这次$SPFA$得到的解重复。这些重复的解的共同特点就是终止的节点的标号都大于这两个的标号。(想象一下,既然第一次搜到了一个标号很大的解,显然它是没有在这个点之前走上最短路的,就是它同时也绕过了后来的一条强制不走的边),所以这个解是可以继续用的。
用堆保存一下走到哪个在最短路上的点导致的更新$ans$的话,就可以在以后反复调用,保证全局最优了。注意每次跑$SPFA$无需清空$dis$,因为最短路从左往右做的时候,$dis$显然递减(想一想最短路的松弛操作),但是到达最短路上的点的标记必须清空。
玄学的时间复杂度...也能过...人生处处是惊喜...难道不是吗?
$SPFA$的时候不要用$STL$中的$queue$,$n$次拓展慢的吓人<代码中沿用,简洁一些(其实我懒,不想改)>。
#include<set> #include<map> #include<queue> #include<stack> #include<cmath> #include<ctime> #include<string> #include<cstdio> #include<vector> #include<cstring> #include<cstdlib> #include<iostream> #include<algorithm> #define LL long long #define RE register #define IL inline using namespace std; const int N=100000; const int M=200000; IL int Read(); int n,m,l,sp[N+5]; int last[N+5],id[N+5]; struct Link { int from,to,cost; }lin[M+5]; struct tt { int to,next,cost,id; }edge[M+5]; int path[N+5],top; IL void Add(int u,int v,int c,int id); struct node { int u,dis; bool operator < (const node &b) const{ return dis>b.dis; } }; priority_queue<node>Q; int dis[N+5]; bool vis[N+5]; IL void SPFA(int donot); int main() { n=Read();m=Read();l=Read(); for (RE int i=1;i<=m;i++) { lin[i].from=Read();lin[i].to=Read();lin[i].cost=Read(); Add(lin[i].from,lin[i].to,lin[i].cost,i); } id[1]=1; for (RE int i=1;i<=l;i++) sp[i]=Read(),id[lin[sp[i]].to]=id[lin[sp[i]].from]+1; for (RE int i=l;i;i--) last[lin[sp[i]].from]=last[lin[sp[i]].to]+lin[sp[i]].cost; memset(dis,127/3,sizeof(dis));dis[1]=0; for (RE int i=1;i<=l;i++) { if (i>1) dis[lin[sp[i-1]].to]=dis[lin[sp[i-1]].from]+lin[sp[i-1]].cost; SPFA(sp[i]); while (!Q.empty()&&id[lin[sp[i]].from]>=id[Q.top().u]) Q.pop(); printf("%d\n",Q.empty() ? -1:Q.top().dis); } return 0; } IL void SPFA(int donot) { queue<int>q; while (!q.empty()) q.pop(); q.push(lin[donot].from); memset(vis,0,sizeof(vis)); while (!q.empty()) { for (RE int i=path[q.front()];i;i=edge[i].next) if (edge[i].id!=donot) { if (dis[edge[i].to]>dis[q.front()]+edge[i].cost) { dis[edge[i].to]=dis[q.front()]+edge[i].cost; if (id[edge[i].to]) Q.push((node){edge[i].to,dis[edge[i].to]+last[edge[i].to]}); else if (!vis[edge[i].to]) { vis[edge[i].to]=1; q.push(edge[i].to); } } } vis[q.front()]=0; q.pop(); } } IL void Add(int u,int v,int c,int id) { edge[++top].to=v; edge[top].cost=c; edge[top].id=id; edge[top].next=path[u]; path[u]=top; } IL int Read() { char c=getchar(); int sum=0; while (c<'0'||c>'9') c=getchar(); while (c>='0'&&c<='9') sum=sum*10+c-'0',c=getchar(); return sum; }
相关文章推荐
- bzoj 3575: [Hnoi2014]道路堵塞
- bzoj3575: [Hnoi2014]道路堵塞
- bzoj 3575: [Hnoi2014]道路堵塞 spfa
- bzoj3575[Hnoi2014]道路堵塞
- BZOJ 3575: [Hnoi2014]道路堵塞
- 3575: [Hnoi2014]道路堵塞
- BZOJ3575 [Hnoi2014]道路堵塞
- bzoj3575 [Hnoi2014]道路堵塞
- 【BZOJ】3575: [Hnoi2014]道路堵塞
- 【bzoj3575】 Hnoi2014—道路堵塞
- bzoj 3575: [Hnoi2014]道路堵塞 最短路
- [HNOI2014]道路堵塞
- BZOJ3575 HNOI2014 道路阻塞
- 【bzoj3573】[HNOI2014]米特运输 树形dp
- noip2014 寻找道路
- bzoj 3571: [Hnoi2014]画框
- 【BZOJ3572】[Hnoi2014]世界树 虚树
- 【bzoj3751】 Hnoi2014—画框
- 【bzoj3576】 Hnoi2014—江南乐
- bzoj3573: [Hnoi2014]米特运输