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[N],enx[N],ev[N],e0[1550],l[1550],ep=2;
bool d1[N],d2[N],d3[N],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;
}