題目描述

最近，Elaxia和w**的關系特別好，他們很想整天在一起，但是大學的學習太緊張了，他們 必須合理地安排兩個人在一起的時間。

Elaxia和w**每天都要奔波於宿舍和實驗室之間，他們 希望在節約時間的前提下，一起走的時間盡可能的長。

現在已知的是Elaxia和w**所在的宿舍和實驗室的編號以及學校的地圖：地圖上有N個路 口，M條路，經過每條路都需要一定的時間。 具體地說，就是要求無向圖中，兩對點間最短路的最長公共路徑。

輸入輸出格式

輸入格式：

第一行：兩個整數N和M（含義如題目描述）。

第二行：四個整數x1、y1、x2、y2（1 ≤ x1 ≤ N，1 ≤ y1 ≤ N，1 ≤ x2 ≤ N，1 ≤ y2 ≤ N），分別表示Elaxia的宿舍和實驗室及w**的宿舍和實驗室的標號（兩對點分別 x1,y1和x2,y2）。

接下來M行：每行三個整數，u、v、l（1 ≤ u ≤ N，1 ≤ v ≤ N，1 ≤ l ≤ 10000），表 u和v之間有一條路，經過這條路所需要的時間為l。

輸出格式：

一行，一個整數，表示每天兩人在一起的時間（即最長公共路徑的長度）

輸入輸出樣例

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

3

說明

對於30%的數據，N ≤ 100；

對於60%的數據，N ≤ 1000；

對於100%的數據，N ≤ 1500，輸入數據保證沒有重邊和自環。

題解

　　不會，抄題解的

　　就這裏->這個

 1 //minamoto
 2 #include<iostream>
 3 #include<cstdio>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<queue>
 7 #define inf 0x3f3f3f3f
 8 using namespace std;
 9 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
10
 
 char buf[1<<21],*p1=buf,*p2=buf;
11 template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
12 inline int read(){
13     #define num ch-‘0‘
14     char ch;bool flag=0;int res;
15     while(!isdigit(ch=getc()))
16     (ch==‘-‘)&&(flag=true);
17     for(res=num;isdigit(ch=getc());res=res*10+num);
18     (flag)&&(res=-res);
19     #undef num
20     return res;
21 }
22 const int N=1505,M=2500005;
23 int n,m,X1,Y1,X2,Y2;
24 int ver[M],Next[M],head[N],edge[M],from[M],tot;
25 int ver0[M],Next0[M],head0[N],edge0[M],from0[M],tot0,im[M];
26 int dis[6][N],sol[N],cnt[N];bool vis[N];
27 queue<int> q;
28 inline void add(int u,int v,int e){
29     ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e,from[tot]=u;
30 }
31 inline void add0(int u,int v,int e,bool x){
32     ver0[++tot0]=v,Next0[tot0]=head0[u],head0[u]=tot0,edge0[tot0]=e,im[tot0]=x;
33 }
34 void spfa(int s,int I){
35     memset(dis[I],0x3f,sizeof(dis[I]));
36     q.push(s),dis[I][s]=0;
37     while(!q.empty()){
38         int u=q.front();q.pop();vis[u]=0;
39         for(int i=head[u];i;i=Next[i]){
40             int v=ver[i];
41             if(dis[I][v]>dis[I][u]+edge[i]){
42                 dis[I][v]=dis[I][u]+edge[i];
43                 if(!vis[v]) vis[v]=1,q.push(v);
44             }
45         }
46     }
47 }
48 void init(){
49     for(int i=1;i<=tot;++i){
50         int u=from[i],v=ver[i];
51         if(dis[1][u]+edge[i]+dis[2][v]==dis[1][Y1]){
52             if(dis[3][u]+edge[i]+dis[4][v]==dis[3][Y2]
53             ||dis[4][u]+edge[i]+dis[3][v]==dis[3][Y2])
54             add0(u,v,edge[i],1);
55             else add0(u,v,edge[i],0);
56             ++cnt[v];
57         }
58     }
59 }
60 void solve(){
61     q.push(X1);
62     while(!q.empty()){
63         int u=q.front();q.pop();
64         for(int i=head0[u];i;i=Next0[i]){
65             int v=ver0[i];
66             if(!(--cnt[v])) q.push(v);
67             cmax(sol[v],sol[u]+edge0[i]*im[i]);
68         }
69     }
70 }
71 int main(){
72     n=read(),m=read();
73     X1=read(),Y1=read(),X2=read(),Y2=read();
74     for(int i=1;i<=m;++i){
75         int u=read(),v=read(),e=read();
76         add(u,v,e),add(v,u,e);
77     }
78     spfa(X1,1),spfa(Y1,2),spfa(X2,3),spfa(Y2,4);
79     init(),solve();
80     printf("%d\n",sol[Y1]);
81     return 0;
82 }

BZOJ1880: [Sdoi2009]Elaxia的路線（最短路）