1. 程式人生 > >spfa求次短路

spfa求次短路

read define struct digi 路徑 測試 can 不存在 記得


思路:
先算出每個點到1的最短路d1[i],記錄下路徑,然後枚舉最短路上的邊

刪掉之後再求一遍最短路,那麽這時的最短路就可能是答案。

既然這樣為甚麽不用A*求次短路呢?
因為A*求次短路處理不了無向圖啊,他會來回的走。
ε=(′ο`*)))唉

上菜:
集合位置

題目描述

每次有大的活動,大家都要在一起“聚一聚”,不管是去好樂迪,還是避風塘,或者湯姆熊,大家都要玩的痛快。還記得心語和花兒在跳舞機上的激情與釋放,還記得草草的投籃技藝是如此的高超,還記得狗狗的槍法永遠是‘S‘……還有不能忘了,胖子的歌聲永遠是讓我們驚叫的!!

今天是野貓的生日,所以想到這些也正常,只是因為是上學日,沒法一起去玩了。但回憶一下那時的甜蜜總是一種幸福嘛。。。

但是每次集合的時候都會出現問題!野貓是公認的“路盲”,野貓自己心裏也很清楚,每次都提前出門,但還是經常遲到,這點讓大家很是無奈。後來,野貓在每次出門前,都會向花兒咨詢一下路徑,根據已知的路徑中,總算能按時到了。

現在提出這樣的一個問題:給出n個點的坐標,其中第一個為野貓的出發位置,最後一個為大家的集合位置,並給出哪些位置點是相連的。野貓從出發點到達集合點,總會挑一條最近的路走,如果野貓沒找到最近的路,他就會走第二近的路。請幫野貓求一下這條第二最短路徑長度。

輸入輸出格式

輸入格式:
第一行是兩個整數n(1<=n<=200)和m,表示一共有n個點和m條路,以下n行每行兩個數xi,yi,(-500<=xi,yi<=500),代表第i個點的坐標,再往下的m行每行兩個整數pj,qj,(1<=pj,qj<=n),表示兩個點相通。

輸出格式:
只有一行包含一個數,為第二最短路線的距離(保留兩位小數),如果存在多條第一短路徑,則答案就是第一最短路徑的長度;如果不存在第二最短路徑,輸出-1。

輸入輸出樣例

輸入樣例#1:
3 3
0 0
1 1
0 2
1 2
1 3
2 3
輸出樣例#1:
2.83
說明

各個測試點1s

裸的次短路啦。

上代碼。註意無解情況判定。

技術分享圖片
 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cmath>
 5 #include<queue>
 6 #include<cstring>
 7
#include<string> 8 #include<map> 9 #define ll long long 10 #define DB double 11 #define eps 1e-3 12 #define inf 2147483647 13 using namespace std; 14 inline int read() 15 { 16 int x=0,w=1;char ch=getchar(); 17 while(!isdigit(ch)){if(ch==-) w=-1;ch=getchar();} 18 while(isdigit(ch)) x=(x<<3)+(x<<1)+ch-0,ch=getchar(); 19 return x*w; 20 } 21 const int N=1e6+10; 22 struct node{ 23 int u,v;DB c; 24 int ne; 25 }e[N]; 26 int h[N],tot,n,m,pre[N]; 27 DB x[N],y[N],d[N]; 28 void add(int u,int v,DB c) 29 { 30 tot++;e[tot]=(node){u,v,c,h[u]};h[u]=tot; 31 } 32 DB dis(int u,int v) 33 { 34 return sqrt((x[u]-x[v])*(x[u]-x[v])+(y[u]-y[v])*(y[u]-y[v])); 35 } 36 bool v[N]; 37 queue<int>q; 38 void spfa() 39 { 40 for(int i=1;i<=n;++i) d[i]=inf,v[i]=0; 41 q.push(1);d[1]=0; 42 while(!q.empty()) 43 { 44 int ff=q.front();q.pop();v[ff]=0; 45 for(int i=h[ff];i;i=e[i].ne) 46 { 47 int rr=e[i].v; 48 if(d[rr]>d[ff]+e[i].c) 49 { 50 pre[rr]=i; 51 d[rr]=d[ff]+e[i].c; 52 if(!v[rr]) v[rr]=1,q.push(rr); 53 } 54 } 55 } 56 } 57 DB SPFA(int k) 58 { 59 for(int i=1;i<=n;++i) d[i]=inf,v[i]=0; 60 q.push(1);d[1]=0; 61 while(!q.empty()) 62 { 63 int ff=q.front();q.pop();v[ff]=0; 64 for(int i=h[ff];i;i=e[i].ne) 65 { 66 int rr=e[i].v; 67 if(i==k) continue; 68 if(d[rr]>d[ff]+e[i].c) 69 { 70 d[rr]=d[ff]+e[i].c; 71 if(!v[rr]) v[rr]=1,q.push(rr); 72 } 73 } 74 } 75 return d[n]; 76 } 77 int main() 78 { 79 n=read();m=read(); 80 for(int i=1;i<=n;++i) scanf("%lf%lf",&x[i],&y[i]); 81 tot=1; 82 for(int u,v,i=1;i<=m;++i) 83 { 84 u=read();v=read(); 85 add(u,v,dis(u,v)); 86 add(v,u,dis(v,u)); 87 } 88 spfa(); 89 int nw=n; 90 DB ans=inf; 91 while(nw!=1) 92 { 93 ans=min(ans,SPFA(pre[nw])); 94 nw=e[pre[nw]].u; 95 } 96 if(ans==inf) printf("-1"); 97 else printf("%.2lf",ans); 98 return 0; 99 }
View Code

(?′?‵?)I L???????

spfa求次短路