BZOJ3575 HNOI2014 道路阻塞
Posted Fenghr
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了BZOJ3575 HNOI2014 道路阻塞相关的知识，希望对你有一定的参考价值。
3575: [Hnoi2014]道路堵塞
Time Limit: 10 Sec Memory Limit: 128 MB
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
1 2 2
1 3 2
3 4 4
3 2 1
2 4 3
1 5
Sample Output
6
HINT
100%的数据满足2<N<100000，1<M<200000。所用道路长度大于0小于10000。
数据已加强By Vfleaking
你看见我的跳蚤国王了吗？哦它被FFF了？资磁！
这题真是火的冒烟，劲爆无比，谁也想不到正解来自一个很扯淡的证明。
正解：SPFA鬼畜技巧。
因为它已经钦定了最短路，而且要对最短路上每一条边输出一次答案，所以很容易想到正解是枚举BAN哪条边再鬼一鬼。
这里有个十分扯淡的证明：新·最短路一腚是S->之前最短路上某一点p->之前最短路上某一点q->T；
感觉并不知道怎么证明是吧，但如果我说S和p，q和T 可以为同一个点呢？
那么所有的情况都被包含了是的吗？没错是的。您强×把答案放在了一个鬼证明上。
所以按照这个证明的思想，我们枚举删掉最短路上的哪条边，然后把这条边的s加进队列做SPFA。
但是这样不加优化的做法会让你过吗？答案是显然否定的。那么来点更劲的：
先把最短路上的所有边都BAN掉，然后从S到T开始恢复边，并把s丢进队列。
然后每次就可以不清空dis数组了，因为你之前SPFA的结果你都可以调用以前的结果。
然后统计答案就是把走到的钦定最短路上的点的距离加上它到T的距离取个min就好了。这里有一个玄妙的做法就是把答案加上这个点放成结构体扔进堆里，每次弹出到达点在s之前的点因为这是在这条路被BAN的时候SPFA到的，不一定是当前的合法方案，也不会比最优答案更好。所以判完-1后，第一个在s之后的点就不管是不是s更新的点了，却又是一个合法的点，就避免了很多麻烦。
代码还是很好看懂的。
#include <iostream> #include <cstdio> #include <cstdlib> #include <algorithm> #include <vector> #include <cstring> #include <queue> #include <complex> #include <stack> #define LL long long int #define dob double using namespace std; const int N = 200010; struct Node{int from,to,val,next;}E[N]; struct Data{int p,val;bool operator <(const Data &b)const{return val>b.val;}}; int n,m,L,head[N],tot,ban[N],rk[N],In[N],far[N],dis[N],fr[N],s[N]; queue<int>Q;priority_queue<Data>T; int gi() { int x=0,res=1;char ch=getchar(); while(ch>‘9‘||ch<‘0‘){if(ch==‘-‘)res*=-1;ch=getchar();} while(ch<=‘9‘&&ch>=‘0‘)x=x*10+ch-48,ch=getchar(); return x*res; } inline void link(int u,int v,int c) { E[++tot]=(Node){u,v,c,head[u]}; head[u]=tot; } inline void SPFA(int st) { queue<int>Q;Q.push(st); while(!Q.empty()){ int x=Q.front();Q.pop();In[x]=0; for(int e=head[x];e;e=E[e].next){ if(ban[e])continue;int y=E[e].to; if(far[x]+E[e].val<far[y]){ far[y]=far[x]+E[e].val; if(fr[y])T.push((Data){fr[y],far[y]+dis[fr[y]]}); else if(!In[y])Q.push(In[y]=y); } } } } int main() { n=gi();m=gi();L=gi(); for(int i=1;i<=m;++i){ int u=gi(),v=gi(),c=gi(); link(u,v,c); } s[1]=1;fr[1]=1; for(int i=1;i<=L;++i){ ban[rk[i]=gi()]=1; s[i+1]=E[rk[i]].to; fr[s[i+1]]=i+1; } for(int i=L;i;--i) dis[i]=dis[i+1]+E[rk[i]].val; memset(far,127/3,sizeof(far));far[1]=0; SPFA(1); for(int i=1;i<=L;++i){ while(!T.empty() && T.top().p<=i)T.pop(); if(T.empty())printf("-1\n"); else printf("%d\n",T.top().val); far[E[rk[i]].to]=far[s[i]]+E[rk[i]].val; SPFA(s[i+1]); } return 0; }
这还真TMD妙不可言，玄之又玄，精妙绝伦啊(伦boy：嗯？)。
以上是关于BZOJ3575 HNOI2014 道路阻塞的主要内容，如果未能解决你的问题，请参考以下文章