【BZOJ1880】[Sdoi2009]Elaxia的路线 最短路+DP

【BZOJ1880】[Sdoi2009]Elaxia的路线

Description

最近，Elaxia和w**的关系特别好，他们很想整天在一起，但是大学的学习太紧张了，他们 必须合理地安排两个人在一起的时间。Elaxia和w**每天都要奔波于宿舍和实验室之间，他们 希望在节约时间的前提下，一起走的时间尽可能的长。 现在已知的是Elaxia和w**所在的宿舍和实验室的编号以及学校的地图：地图上有N个路 口，M条路，经过每条路都需要一定的时间。 具体地说，就是要求无向图中，两对点间最短路的最长公共路径。

Input

第一行：两个整数N和M（含义如题目描述）。 第二行：四个整数x1、y1、x2、y2（1 ≤ x1 ≤ N，1 ≤ y1 ≤ N，1 ≤ x2 ≤ N，1 ≤ ≤ N），分别表示Elaxia的宿舍和实验室及w**的宿舍和实验室的标号（两对点分别 x1,y1和x2,y2）。 接下来M行：每行三个整数，u、v、l（1 ≤ u ≤ N，1 ≤ v ≤ N，1 ≤ l ≤ 10000），表 u和v之间有一条路，经过这条路所需要的时间为l。 出出出格格格式式式：：： 一行，一个整数，表示每天两人在一起的时间（即最长公共路径的长度）。

Output

一行，一个整数，表示每天两人在一起的时间（即最长公共路径的长度）

Sample Input

9 10
1 6 7 8
1 2 1
2 5 2
2 3 3
3 4 2
3 9 5
4 5 3
4 6 4
4 7 2
5 8 1
7 9 1

Sample Output

3

HINT

对于30%的数据，N ≤ 100；
对于60%的数据，N ≤ 1000；
对于100%的数据，N ≤ 1500，输入数据保证没有重边和自环。

题解：先从x1,y1,x2都跑一边最短路，这样我们就可以快速判断一条边是否在最短路上了。最后从y2跑最短路，然后用f[x]表示从y2走到x时最长公共路径的长度。边最短路边DP即可。

然而好像被hack了。。。还是看大爷的题解吧。

#include <cstring>
#include <cstdio>
#include <iostream>
#include <queue>
#include <utility>
#define mp(A,B) make_pair(A,B)
using namespace std;
typedef pair<int,int> pii;
int n,m,cnt,S1,S2,T1,T2,la,lb;
int to[1000010],next[1000010],val[1000010],head[10010],dis[4][10010],vis[10010],f[10010];
priority_queue<pii> q;
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<'0'||gc>'9')	{if(gc=='-')f=-f;	gc=getchar();}
	while(gc>='0'&&gc<='9')	ret=ret*10+gc-'0',gc=getchar();
	return ret*f;
}
void dijkstra(int S,int t)
{
	int i,u;
	memset(vis,0,sizeof(vis));
	dis[t][S]=0,q.push(mp(0,S));
	while(!q.empty())
	{
		u=q.top().second,q.pop();
		if(vis[u])	continue;
		vis[u]=1;
		for(i=head[u];i!=-1;i=next[i])	if(dis[t][to[i]]>dis[t][u]+val[i])
			dis[t][to[i]]=dis[t][u]+val[i],q.push(mp(-dis[t][to[i]],to[i]));
	}
}
inline void add(int a,int b,int c)
{
	to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
}
int main()
{
	n=rd(),m=rd(),S1=rd(),T1=rd(),S2=rd(),T2=rd();
	int i,a,b,c,u,v;
	memset(head,-1,sizeof(head));
	for(i=1;i<=m;i++)	a=rd(),b=rd(),c=rd(),add(a,b,c),add(b,a,c);
	memset(dis,0x3f,sizeof(dis));
	dijkstra(S1,0),dijkstra(T1,1),dijkstra(S2,2),la=dis[0][T1],lb=dis[2][T2];
	memset(vis,0,sizeof(vis));
	dis[3][T2]=0,q.push(mp(0,T2));
	while(!q.empty())
	{
		u=q.top().second,q.pop();
		if(vis[u])	continue;
		vis[u]=1;
		for(i=head[u];i!=-1;i=next[i])
		{
			v=to[i];
			if(dis[3][v]>dis[3][u]+val[i])
			{
				dis[3][v]=dis[3][u]+val[i],f[v]=0;
				q.push(mp(-dis[3][v],v));
			}
			if(dis[3][v]==dis[3][u]+val[i])
			{
				f[v]=max(f[v],f[u]+val[i]*((dis[1][u]+dis[0][v]+val[i]==la||dis[0][u]+dis[1][v]+val[i]==la)&&(dis[3][u]+dis[2][v]+val[i]==lb||dis[2][u]+dis[3][v]+val[i]==lb)));
			}
		}
	}
	printf("%d",f[S2]);
	return 0;
}//3 2 1 2 1 3 1 2 1 2 3 1