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

bzoj1726 [Usaco2006 Nov]Roadblocks第二短路

1726: [Usaco2006 Nov]Roadblocks第二短路

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 826  Solved: 394
[Submit][Status]

Description

贝茜把家搬到了一个小农场,但她常常回到FJ的农场去拜访她的朋友。贝茜很喜欢路边的风景,不想那么快地结束她的旅途,于是她每次回农场,都会选择第二短的路径,而不象我们所习惯的那样,选择最短路。 贝茜所在的乡村有R(1<=R<=100,000)条双向道路,每条路都联结了所有的N(1<=N<=5000)个农场中的某两个。贝茜居住在农场1,她的朋友们居住在农场N(即贝茜每次旅行的目的地)。 贝茜选择的第二短的路径中,可以包含任何一条在最短路中出现的道路,并且,一条路可以重复走多次。当然咯,第二短路的长度必须严格大于最短路(可能有多条)的长度,但它的长度必须不大于所有除最短路外的路径的长度。

Input

* 第1行: 两个整数,N和R,用空格隔开

* 第2..R+1行: 每行包含三个用空格隔开的整数A、B和D,表示存在一条长度为 D(1 <= D <= 5000)的路连接农场A和农场B

Output

* 第1行: 输出一个整数,即从农场1到农场N的第二短路的长度

Sample Input

4 4
1 2 100
2 4 200
2 3 250
3 4 100


Sample Output

450

输出说明:

最短路:1 -> 2 -> 4 (长度为100+200=300)
第二短路:1 -> 2 -> 3 -> 4 (长度为100+250+100=450)

 

次短路的做法好多吧……

可以在spfa的时候顺便维护个次短的长度

也可以正反各一遍最短路,然后枚举每一条边的dis[1][l]+dis[r][n]+v更新次短路(我就是这么做的)

还可以像k短路一样倒着跑最短路然后正着A*

#include<cstdio>#include<iostream>#include<cstring>using namespace std;struct ed{int x,y,z;}b[300010];struct edge{int to,next,v;}e[300010];int head[10001];int dis1[10001];int dis2[10001];int q[1000010];bool mrk[10001];int n,m,cnt,mx,ans=2147483647;inline void ins(int u,int v,int w){	e[++cnt].to=v;	e[cnt].next=head[u];	e[cnt].v=w;	head[u]=cnt;}inline void insert(int u,int v,int w){	ins(u,v,w);	ins(v,u,w);}inline void spfa(int S,int *dis){	memset(mrk,0,sizeof(mrk));	memset(q,0,sizeof(q));	dis[S]=0;q[1]=S;mrk[S]=1;	int t=0,w=1;	while (t<w)	{		int now=q[++t];		for (int i=head[now];i;i=e[i].next)		  if (dis[e[i].to]>dis[now]+e[i].v)			{				dis[e[i].to]=dis[now]+e[i].v;				if (!mrk[e[i].to])				{					mrk[e[i].to]=1;					q[++w]=e[i].to;				}			}		mrk[now]=0;	}}int main(){	scanf("%d%d",&n,&m);	for (int i=1;i<=m;i++)	{		scanf("%d%d%d",&b[i].x,&b[i].y,&b[i].z);		insert(b[i].x,b[i].y,b[i].z);	}	memset(dis1,127/3,sizeof(dis1));	memset(dis2,127/3,sizeof(dis2));	spfa(1,dis1);	spfa(n,dis2);	mx=dis1[n];	for (int i=1;i<=m;i++)	{		int now=dis1[b[i].x]+dis2[b[i].y]+b[i].z;		now=min(now,dis1[b[i].y]+dis2[b[i].x]+b[i].z);		if (now!=mx&&now<ans)ans=now;	}	for (int i=1;i<=n;i++)	  for(int j=head[i];j;j=e[j].next)	    ans=min(ans,dis1[i]+dis2[i]+e[j].v*2);	printf("%d\n",ans);}

  

bzoj1726 [Usaco2006 Nov]Roadblocks第二短路