bzoj1880 [Sdoi2009]Elaxia的路线

#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;
}