ccz181078

  博客园 :: 首页 :: 博问 :: 闪存 :: 新随笔 :: 联系 :: :: 管理 ::

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

 

posted on 2016-04-24 17:01  nul  阅读(235)  评论(0编辑  收藏  举报