首页 > 代码库 > bzoj1726[Usaco2006 Nov]Roadblocks第二短路*

bzoj1726[Usaco2006 Nov]Roadblocks第二短路*

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

bzoj1726[Usaco2006 Nov]Roadblocks第二短路*