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

一行，一个整数，表示每天两人在一起的时间（即最长公共路径的长度）
求出x1-y1，x2-y2的最短路图的交集并求最长链即为答案，由于题意不明还要把x2,y2交换再计算一次

#include<cstdio>
#include<queue>
int n,m,x1,y1,x2,y2,ans=0;
const int N=1550*1550;
int es
,enx
,ev
,e0[1550],l[1550],ep=2;
bool d1
,d2
,d3
,in[1550],a1[1550],a2[1550],ed1[1550],ed2[1550];
int ins[1550];
void sp(int s,int t){
for(int i=1;i<=n;i++)l[i]=2147483647;
l[s]=0;
std::queue<int>q;
q.push(s);
while(!q.empty()){
int w=q.front();q.pop();
in[w]=0;
for(int i=e0[w];i;i=enx[i]){
int u=es[i];
if(l[w]+ev[i]<l[u]){
l[u]=l[w]+ev[i];
if(!in[u])in[u]=1,q.push(u);
}
}
}
}
bool dfs1(int w){
if(w==y1)return 1;
if(ed1[w])return a1[w];
bool t=0;
for(int i=e0[w];i;i=enx[i]){
int u=es[i];
if(l[w]+ev[i]==l[u])
t|=(d1[i]=dfs1(u));
}
a1[w]=t;ed1[w]=1;
return t;
}
bool dfs2(int w){
if(w==y2)return 1;
if(ed2[w])return a2[w];
bool t=0;
for(int i=e0[w];i;i=enx[i]){
int u=es[i];
if(l[w]+ev[i]==l[u])
t|=(d2[i]=d3[i^1]=dfs2(u));
}
a2[w]=t;ed2[w]=1;
return t;
}
void chk(bool*d1,bool*d2){
std::queue<int>q;
for(int i=1;i<=n;i++)if(!ins[i])q.push(i);
while(!q.empty()){
int w=q.front();q.pop();
for(int i=e0[w];i;i=enx[i]){
if(!d1[i]||!d2[i])continue;
int u=es[i];
if(l[w]+ev[i]>l[u])l[u]=l[w]+ev[i];
if(!--ins[u])q.push(u);
}
}
for(int i=1;i<=n;i++){
if(l[i]>ans)ans=l[i];
l[i]=0;
}
}
int main(){
scanf("%d%d",&n,&m);
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
for(int i=0,a,b,c;i<m;i++){
scanf("%d%d%d",&a,&b,&c);
es[ep]=b;enx[ep]=e0[a];ev[ep]=c;e0[a]=ep++;
es[ep]=a;enx[ep]=e0[b];ev[ep]=c;e0[b]=ep++;
}
sp(x1,y1);
dfs1(x1);
sp(x2,y2);
dfs2(x2);
for(int i=1;i<=n;i++){
l[i]=0;
for(int j=e0[i];j;j=enx[j]){
if(d1[j]&d2[j])++ins[es[j]];
}
}
chk(d1,d2);
for(int i=1;i<=n;i++){
l[i]=0;
for(int j=e0[i];j;j=enx[j]){
if(d1[j]&d3[j])++ins[es[j]];
}
}
chk(d1,d3);
printf("%d",ans);
return 0;
}