最短路径——Floyd-Warshall算法

最新推荐文章于 2020-09-15 12:01:54 发布

weixin_30535043

最新推荐文章于 2020-09-15 12:01:54 发布

阅读量120

点赞数

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/orange1438/p/4054649.html

版权

Floyd-Warshall算法，简称Floyd算法，用于求解任意两点间的最短距离，时间复杂度为O(n^3)。

我们平时所见的Floyd算法的一般形式如下：

1 void Floyd()
2 {
3      int i,j,k;
4      for(k=1;k<=n;k++)
5          for(i=1;i<=n;i++)
6              for(j=1;j<=n;j++)
7                  if(dist[i][k]+dist[k][j]<dist[i][j])
8                      dist[i][j]=dist[i][k]+dist[k][j];
9 ｝

注意下第6行这个地方，如果dist[i][k]或者dist[k][j]不存在，程序中用一个很大的数代替。最好写成if(dist[i][k]!=INF && dist[k][j]!=INF && dist[i][k]+dist[k][j]<dist[i][j])，从而防止溢出所造成的错误。
上面这个形式的算法其实是Floyd算法的精简版，而真正的Floyd算法是一种基于DP(Dynamic Programming)的最短路径算法。

例题分析：

设图G中n 个顶点的编号为1到n。令c [i, j, k]表示从i 到j 的最短路径的长度，其中k 表示该路径中的最大顶点，也就是说c[i,j,k]这条最短路径所通过的中间顶点最大不超过k。因此，如果G中包含边<i, j>，则c[i, j, 0] =边<i, j> 的长度；若i= j ，则c[i,j,0]=0；如果G中不包含边<i, j>，则c (i, j, 0)= +∞。c[i, j, n] 则是从i 到j 的最短路径的长度。对于任意的k>0，通过分析可以得到：中间顶点不超过k 的i 到j 的最短路径有两种可能：该路径含或不含中间顶点k。若不含，则该路径长度应为c[i, j, k-1]，否则长度为 c[i, k, k-1] +c [k, j, k-1]。c[i, j, k]可取两者中的最小值。状态转移方程：c[i, j, k]=min{c[i, j, k-1], c [i, k, k-1]+c [k, j, k-1]}，k＞0。这样，问题便具有了最优子结构性质，可以用动态规划方法来求解。

为了进一步理解，观察上面这个有向图：若k=0, 1, 2, 3，则c[1,3,k]= +∞；c[1,3,4]= 28；若k = 5, 6, 7，则c [1,3,k] = 10；若k=8, 9, 10，则c[1,3,k] = 9。因此1到3的最短路径长度为9。下面通过程序来分析这一DP过程，对应上面给出的有向图：

 1 #include <iostream>
 2 using namespace std;
 3 
 4 const int INF = 100000;
 5 int n=10,map[11][11],dist[11][11][11];
 6 void init()
 7 {
 8     int i,j;
 9     for(i=1;i<=n;i++)
10         for(j=1;j<=n;j++)
11             map[i][j]=(i==j)?0:INF;
12     map[1][2]=2,map[1][4]=20,map[2][5]=1;
13     map[3][1]=3,map[4][3]=8,map[4][6]=6;
14     map[4][7]=4,map[5][3]=7,map[5][8]=3;
15     map[6][3]=1,map[7][8]=1,map[8][6]=2;
16     map[8][10]=2,map[9][7]=2,map[10][9]=1;
17 }
18 void floyd_dp()
19 {
20     int i,j,k;
21     for(i=1;i<=n;i++)
22         for(j=1;j<=n;j++)
23             dist[i][j][0]=map[i][j];
24     for(k=1;k<=n;k++)
25         for(i=1;i<=n;i++)
26             for(j=1;j<=n;j++){
27                 dist[i][j][k]=dist[i][j][k-1];
28                 if(dist[i][k][k-1]+dist[k][j][k-1]<dist[i][j][k])
29                     dist[i][j][k]=dist[i][k][k-1]+dist[k][j][k-1];
30             }
31 }
32 int main()
33 {
34     int k,u,v;
35     init();
36     floyd_dp();
37     while(cin>>u>>v,u||v)
38     {
39         for(k=0;k<=n;k++)
40         {
41             if(dist[u][v][k]==INF) cout<<"+∞"<<endl;
42              else cout<<dist[u][v][k]<<endl;
43         }
44     }
45     return 0;
46  }

Floyd-Warshall算法不仅能求出任意2点间的最短路径，还可以保存最短路径上经过的节点。下面用精简版的Floyd算法实现这一过程，程序中的图依然对应上面的有向图。

 1 #include <iostream>
 2 using namespace std;
 3 
 4 const int INF = 100000;
 5 int n=10,path[11][11],dist[11][11],map[11][11];
 6 void init(){
 7     int i,j;
 8     for(i=1;i<=n;i++)
 9         for(j=1;j<=n;j++)
10             map[i][j]=(i==j)?0:INF;
11     map[1][2]=2,map[1][4]=20,map[2][5]=1;
12     map[3][1]=3,map[4][3]=8,map[4][6]=6;
13     map[4][7]=4,map[5][3]=7,map[5][8]=3;
14     map[6][3]=1,map[7][8]=1,map[8][6]=2;
15     map[8][10]=2,map[9][7]=2,map[10][9]=1;
16 }
17 void floyd(){
18     int i,j,k;
19     for(i=1;i<=n;i++)
20         for(j=1;j<=n;j++)
21             dist[i][j]=map[i][j],path[i][j]=0;
22     for(k=1;k<=n;k++)
23         for(i=1;i<=n;i++)
24             for(j=1;j<=n;j++)
25                 if(dist[i][k]+dist[k][j]<dist[i][j])
26                     dist[i][j]=dist[i][k]+dist[k][j],path[i][j]=k;
27 }
28 void output(int i,int j){
29     if(i==j) return;
30     if(path[i][j]==0) cout<<j<<' ';
31     else{
32         output(i,path[i][j]);
33         output(path[i][j],j);
34     }
35 }
36 int main(){
37     int u,v;
38     init();
39     floyd();
40     while(cin>>u>>v,u||v){
41         if(dist[u][v]==INF) cout<<"No path"<<endl;
42         else{
43             cout<<u<<' ';
44             output(u,v);
45             cout<<endl;
46         }
47     }
48     return 0;
49 }