bzoj1726[Usaco2006 Nov]Roadblocks第二短路
题意:
求无向图点1到n的次短路(长度严格小于最短路)。点数≤5000,边数≤100000。
题解:
求源点为1的单源最短路和源点为n的单源最短路。然后枚举每个点,如果某点到点1和点n的距离和不等于1到n的最短路距离且最小则答案为它。
代码:
1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #include <queue> 5 #define inc(i,j,k) for(int i=j;i<=k;i++) 6 #define maxn 5010 7 #define INF 0x3fffffff 8 using namespace std; 9 10 inline int read(){ 11 char ch=getchar(); int f=1,x=0; 12 while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();} 13 while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar(); 14 return f*x; 15 } 16 struct e{int f,t,w,n;}es[maxn*40]; int g[maxn],ess,d[2][maxn],n,m; bool inq[maxn]; queue<int>q; 17 void pe(int f,int t,int w){es[++ess]=(e){f,t,w,g[f]}; g[f]=ess; es[++ess]=(e){t,f,w,g[t]}; g[t]=ess;} 18 void spfa(int s,bool a){ 19 while(!q.empty())q.pop(); memset(inq,0,sizeof(inq)); inc(i,1,n)d[a][i]=INF; 20 q.push(s); inq[s]=1; d[a][s]=0; 21 while(!q.empty()){ 22 int x=q.front(); q.pop(); inq[x]=0; 23 for(int i=g[x];i;i=es[i].n)if(d[a][es[i].t]>d[a][x]+es[i].w){ 24 d[a][es[i].t]=d[a][x]+es[i].w; 25 if(!inq[es[i].t])q.push(es[i].t),inq[es[i].t]=1; 26 } 27 } 28 } 29 int main(){ 30 n=read(); m=read(); inc(i,1,m){int a=read(),b=read(),c=read(); pe(a,b,c);} 31 spfa(1,0); spfa(n,1); int mn=INF; 32 inc(i,1,ess){ 33 if(d[0][es[i].f]+es[i].w+d[1][es[i].t]<mn&&d[0][es[i].f]+es[i].w+d[1][es[i].t]!=d[0][n]) 34 mn=d[0][es[i].f]+es[i].w+d[1][es[i].t]; 35 } 36 printf("%d",mn); return 0; 37 }
20160907