bzoj 1726: [Usaco2006 Nov]Roadblocks第二短路

1726: [Usaco2006 Nov]Roadblocks第二短路

Time Limit: 5 Sec   Memory Limit: 64 MB
Submit: 1274   Solved: 605
[ Submit][ Status][ Discuss]

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)

HINT

Source

Gold

求严格次短路

有三种方法啦

1、直接跑最短路时维护次短路

2、跑两遍最短路

一遍从起点出发，另一遍从终点出发

最后枚举每条边（包括反向边）

维护dis[s][from]+dis[to][t]+cost的次小值即可

3、用A*求第k短路（还不会QAQ）

我这里用的是方法1

具体的转移看代码吧

（不是很确定slf优化可不可以加就没加，不过应该是一样的QAQ）

#include<cstdio>
#include<cstring>
#include<queue>
const int M=100007,N=5005,len=1e5;
struct node
{
	int to,next,c;
}e[M*2];
int cnt=0,first[N],disa[N],disb[N];
bool visit[N];
int qu[len+4];
int n,m;
inline void insert(int u,int v,int c)
{
	e[++cnt]=(node){v,first[u],c};first[u]=cnt;
	e[++cnt]=(node){u,first[v],c};first[v]=cnt;
}
void spfa()
{
	memset(disa,127,sizeof(disa));
	memset(disb,127,sizeof(disb));
	
	int head=1,tail=2;
	visit[1]=1;
	disa[1]=0;
	qu[1]=1;
	while(head!=tail)
	{
		int rr=qu[head++];if(head==len)head=0;
		for(int k=first[rr];k;k=e[k].next)
		{
			bool f=0;int to=e[k].to;
			if(disa[rr]+e[k].c<disa[to])
			{	
				disb[to]=std::min(disa[to],disb[rr]+e[k].c);
				disa[to]=disa[rr]+e[k].c,f=1;
			}
			else if(disa[rr]+e[k].c==disa[to]&&disb[rr]+e[k].c<disb[to])//严格小于 
			disb[to]=disb[rr]+e[k].c,f=1;
			else if(disa[rr]+e[k].c>disa[to]&&disa[rr]+e[k].c<disb[to])
			disb[to]=disa[rr]+e[k].c,f=1;
			if(f==1&&visit[to]==0)
			{
				visit[to]=1;
				qu[tail++]=to;if(tail==len)tail=0;
			}
		}
		visit[rr]=0;
	}
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int u,v,q;
		scanf("%d %d %d",&u,&v,&q);
		insert(u,v,q);
	}
	spfa();
	printf("%d\n",disb[n]);
	return 0;
}