1. 程式人生 > >[SDOI2009]Elaxia的路線 SPFA+Topo

[SDOI2009]Elaxia的路線 SPFA+Topo

emp pop back blog etc bbf pro ret dig

P2149 [SDOI2009]Elaxia的路線

題目描述

最近,Elaxia和w的關系特別好,他們很想整天在一起,但是大學的學習太緊張了,他們 必須合理地安排兩個人在一起的時間。Elaxia和w每天都要奔波於宿舍和實驗室之間,他們 希望在節約時間的前提下,一起走的時間盡可能的長。 現在已知的是Elaxia和w**所在的宿舍和實驗室的編號以及學校的地圖:地圖上有N個路 口,M條路,經過每條路都需要一定的時間。 具體地說,就是要求無向圖中,兩對點間最短路的最長公共路徑。

輸入輸出格式

輸入格式:

第一行:兩個整數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。

輸出格式:

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

輸入輸出樣例

輸入樣例#1:
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
輸出樣例#1:
3

說明

對於30%的數據,N <= 100;

對於60%的數據,N <= 1000;

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

從 四個點 分別跑一次SPFA 記錄dis值

然後找出重復路徑 跑Topo 排序 找重復路徑的最大值

不過由於路徑的方向問題 有可能有重復的路徑無法被記錄

只需要把 x2 y2 倒過來跑一遍就好了

技術分享
  1 #include <queue>
  2 #include <cstdio>
  3 #include <cctype>
  4 #include <vector>
  5 
  6 const int MAXN=2010;
  7 const int INF=0x7fffffff;
  8 
  9 int n,m,x1,y1,x2,y2,ans;
 10 
 11 int dis[6][MAXN],num[MAXN],in
[MAXN],q[MAXN*1000]; 12 13 bool vis[MAXN]; 14 15 inline void read(int&x) { 16 int f=1;register char c=getchar(); 17 for(x=0;!isdigit(c);c==-&&(f=-1),c=getchar()); 18 for(;isdigit(c);x=x*10+c-48,c=getchar()); 19 x=x*f; 20 } 21 22 struct node { 23 int to,val; 24 node() {} 25 node(int to,int val):to(to),val(val) {} 26 }; 27 28 std::vector<node> Graph[MAXN],New[MAXN]; 29 30 void SPFA(int S,int flag) { 31 for(register int i=1;i<=n;++i) dis[flag][i]=INF,vis[i]=false; 32 dis[flag][S]=0; 33 std::queue<int> Q; 34 Q.push(S); 35 while(!Q.empty()) { 36 int u=Q.front(); 37 Q.pop(); 38 vis[u]=false; 39 for(register int i=0;i<Graph[u].size();++i) { 40 node p=Graph[u][i]; 41 if(dis[flag][p.to]>dis[flag][u]+p.val) { 42 dis[flag][p.to]=dis[flag][u]+p.val; 43 if(!vis[p.to]) Q.push(p.to),vis[p.to]=true; 44 } 45 } 46 } 47 } 48 49 void topo() { 50 int head=0,tail=0; 51 for(int i=1;i<=n;++i) { 52 if(!in[i]) q[++tail]=i; 53 dis[5][i]=0; 54 } 55 while(head<tail) { 56 int u=q[++head]; 57 for(int i=0;i<New[u].size();++i) { 58 node v=New[u][i]; 59 dis[5][v.to]=dis[5][u]+v.val; 60 if(!--in[v.to]) q[++tail]=v.to; 61 } 62 } 63 for(int i=1;i<=n;++i) 64 ans=ans<dis[5][i]?dis[5][i]:ans; 65 return; 66 } 67 68 int hh() { 69 read(n);read(m); 70 read(x1);read(y1);read(x2);read(y2); 71 for(register int x,y,z,i=1;i<=m;++i) { 72 read(x);read(y);read(z); 73 Graph[x].push_back(node(y,z)); 74 Graph[y].push_back(node(x,z)); 75 } 76 SPFA(x1,1); 77 SPFA(x2,2); 78 SPFA(y1,3); 79 SPFA(y2,4); 80 for(int i=1;i<=n;++i) 81 for(int j=0;j<Graph[i].size();++j) { 82 node v=Graph[i][j]; 83 if(dis[1][i]+dis[3][v.to]+v.val==dis[1][y1]&&dis[2][i]+dis[4][v.to]+v.val==dis[2][y2]) 84 New[i].push_back(node(v.to,v.val)),++in[v.to]; 85 } 86 topo(); 87 SPFA(y2,2);SPFA(x2,4); 88 for(int i=1;i<=n;++i) New[i].clear(),in[i]=0; 89 for(int i=1;i<=n;++i) 90 for(int j=0;j<Graph[i].size();++j) { 91 node v=Graph[i][j]; 92 if(dis[1][i]+dis[3][v.to]+v.val==dis[1][y1]&&dis[2][i]+dis[4][v.to]+v.val==dis[2][x2]) 93 New[i].push_back(node(v.to,v.val)),++in[v.to]; 94 } 95 topo(); 96 printf("%d\n",ans); 97 return 0; 98 } 99 100 int sb=hh(); 101 int main(int argc,char**argv) {;}
代碼

[SDOI2009]Elaxia的路線 SPFA+Topo