最短路径Floyd算法【图文详解】
Floyd算法
1.定义概览
Floyd-Warshall算法(Floyd-Warshall algorithm)是解决任意两点间的最短路径的一种算法,可以正确处理有向图或负权的最短路径问题,同时也被用于计算有向图的传递闭包。Floyd-Warshall算法的时间复杂度为O(N3),空间复杂度为O(N2)。
2.算法描述
1)算法思想原理:
Floyd算法是一个经典的动态规划算法。用通俗的语言来描述的话,首先我们的目标是寻找从点i到点j的最短路径。从动态规划的角度看问题,我们需要为这个目标重新做一个诠释(这个诠释正是动态规划最富创造力的精华所在)
从任意节点i到任意节点j的最短路径不外乎2种可能,1是直接从i到j,2是从i经过若干个节点k到j。所以,我们假设Dis(i,j)为节点u到节点v的最短路径的距离,对于每一个节点k,我们检查Dis(i,k) + Dis(k,j) < Dis(i,j)是否成立,如果成立,证明从i到k再到j的路径比i直接到j的路径短,我们便设置Dis(i,j) = Dis(i,k) + Dis(k,j),这样一来,当我们遍历完所有节点k,Dis(i,j)中记录的便是i到j的最短路径的距离。
2).算法描述:
a.从任意一条单边路径开始。所有两点之间的距离是边的权,如果两点之间没有边相连,则权为无穷大。
b.对于每一对顶点 u 和 v,看看是否存在一个顶点 w 使得从 u 到 w 再到 v 比己知的路径更短。如果是更新它。
3).Floyd算法过程矩阵的计算----十字交叉法
方法:两条线,从左上角开始计算一直到右下角 如下所示
给出矩阵,其中矩阵A是邻接矩阵,而矩阵Path记录u,v两点之间最短路径所必须经过的点
相应计算方法如下:
最后A3即为所求结果
3.算法代码实现
1 typedef struct 2 { 3 char vertex[VertexNum]; //顶点表 4 int edges[VertexNum][VertexNum]; //邻接矩阵,可看做边表 5 int n,e; //图中当前的顶点数和边数 6 }MGraph; 7 8 void Floyd(MGraph g) 9 { 10 int A[MAXV][MAXV]; 11 int path[MAXV][MAXV]; 12 int i,j,k,n=g.n; 13 for(i=0;i<n;i++) 14 for(j=0;j<n;j++) 15 { 16 A[i][j]=g.edges[i][j]; 17 path[i][j]=-1; 18 } 19 for(k=0;k<n;k++) 20 { 21 for(i=0;i<n;i++) 22 for(j=0;j<n;j++) 23 if(A[i][j]>(A[i][k]+A[k][j])) 24 { 25 A[i][j]=A[i][k]+A[k][j]; 26 path[i][j]=k; 27 } 28 } 29 }
算法时间复杂度:O(n3)
对上面代码还可以进一步优化?答案是肯定的,不信你看下面这个小故事告诉你诀窍:
![081028xjgvimgz7882qdu7.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081028xjgvimgz7882qdu7.png)
上图中有4个城市8条公路,公路上的数字表示这条公路的长短。请注意这些公路是单向的。我们现在需要求任意两个城市之间的最短路程,也就是求任意两个点之间的最短路径。这个问题这也被称为“多源最短路径”问题。
![081028o2n5ebn8hdeh9e5l.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081028o2n5ebn8hdeh9e5l.png)
现在回到问题:如何求任意两点之间最短路径呢?通过之前的学习我们知道通过深度或广度优先搜索可以求出两点之间的最短路径。所以进行n2遍深度或广度优先搜索,即对每两个点都进行一次深度或广度优先搜索,便可以求得任意两点之间的最短路径。可是还有没有别的方法呢?
![081029zdxxq919ttqt8tu8.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081029zdxxq919ttqt8tu8.png)
如现在只允许经过1号顶点,求任意两点之间的最短路程,应该如何求呢?只需判断e[i][1]+e[1][j]是否比e[i][j]要小即可。e[i][j]表示的是从i号顶点到j号顶点之间的路程。e[i][1]+e[1][j]表示的是从i号顶点先到1号顶点,再从1号顶点到j号顶点的路程之和。其中i是1~n循环,j也是1~n循环,代码实现如下。
1 for(i=1;i<=n;i++) 2 { 3 for(j=1;j<=n;j++) 4 { 5 if ( e[i][j] > e[i][1]+e[1][j] ) 6 e[i][j] = e[i][1]+e[1][j]; 7 } 8 }
在只允许经过1号顶点的情况下,任意两点之间的最短路程更新为:
![081029itl7z7m4l9qqg56d.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081029itl7z7m4l9qqg56d.png)
通过上图我们发现:在只通过1号顶点中转的情况下,3号顶点到2号顶点(e[3][2])、4号顶点到2号顶点(e[4][2])以及4号顶点到3号顶点(e[4][3])的路程都变短了。
接下来继续求在只允许经过1和2号两个顶点的情况下任意两点之间的最短路程。如何做呢?我们需要在只允许经过1号顶点时任意两点的最短路程的结果下,再判断如果经过2号顶点是否可以使得i号顶点到j号顶点之间的路程变得更短。即判断e[i][2]+e[2][j]是否比e[i][j]要小,代码实现为如下。
1 //经过1号顶点 2 for(i=1;i<=n;i++) 3 for(j=1;j<=n;j++) 4 if (e[i][j] > e[i][1]+e[1][j]) e[i][j]=e[i][1]+e[1][j]; 5 //经过2号顶点 6 for(i=1;i<=n;i++) 7 for(j=1;j<=n;j++) 8 if (e[i][j] > e[i][2]+e[2][j]) e[i][j]=e[i][2]+e[2][j];
在只允许经过1和2号顶点的情况下,任意两点之间的最短路程更新为:
![081029e7gjlaaul4zk7z4n.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081029e7gjlaaul4zk7z4n.png)
![081029pd747o8o87o07o7l.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081029pd747o8o87o07o7l.png)
![081030h7tmht7cs2h7qftu.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081030h7tmht7cs2h7qftu.png)
整个算法过程虽然说起来很麻烦,但是代码实现却非常简单,核心代码只有五行:
1 for(k=1;k<=n;k++) 2 for(i=1;i<=n;i++) 3 for(j=1;j<=n;j++) 4 if(e[i][j]>e[i][k]+e[k][j]) 5 e[i][j]=e[i][k]+e[k][j];
这段代码的基本思想就是:最开始只允许经过1号顶点进行中转,接下来只允许经过1和2号顶点进行中转……允许经过1~n号所有顶点进行中转,求任意两点之间的最短路程。用一句话概括就是:从i号顶点到j号顶点只经过前k号点的最短路程。其实这是一种“动态规划”的思想!
动态转移方程:dp[k][i][j]=min(dp[k-1][i][j],dp[k-1][i][k]+dp[k-1][k][j]);
下面给出这个算法的完整代码:
1 #include <stdio.h> 2 int main() 3 { 4 int e[10][10],k,i,j,n,m,t1,t2,t3; 5 int inf=99999999; //用inf(infinity的缩写)存储一个我们认为的正无穷值 6 //读入n和m,n表示顶点个数,m表示边的条数 7 scanf("%d %d",&n,&m); 8 9 //初始化 10 for(i=1;i<=n;i++) 11 for(j=1;j<=n;j++) 12 if(i==j) e[i][j]=0; 13 else e[i][j]=inf; 14 //读入边 15 for(i=1;i<=m;i++) 16 { 17 scanf("%d %d %d",&t1,&t2,&t3); 18 e[t1][t2]=t3; 19 } 20 21 //Floyd-Warshall算法核心语句 22 for(k=1;k<=n;k++) 23 for(i=1;i<=n;i++) 24 for(j=1;j<=n;j++) 25 if(e[i][j]>e[i][k]+e[k][j] ) 26 e[i][j]=e[i][k]+e[k][j]; 27 28 //输出最终的结果 29 for(i=1;i<=n;i++) 30 { 31 for(j=1;j<=n;j++) 32 { 33 printf("%10d",e[i][j]); 34 } 35 printf("\n"); 36 } 37 38 return 0; 39 }
有一点需要注意的是:如何表示正无穷。我们通常将正无穷定义为99999999,因为这样即使两个正无穷相加,其和仍然不超过int类型的范围(C语言int类型可以存储的最大正整数是2147483647)。在实际应用中最好估计一下最短路径的上限,只需要设置比它大一点既可以。例如有100条边,每条边不超过100的话,只需将正无穷设置为10001即可。如果你认为正无穷和其它值相加得到一个大于正无穷的数是不被允许的话,我们只需在比较的时候加两个判断条件就可以了,请注意下面代码中带有下划线的语句。
1 //Floyd-Warshall算法核心语句 2 for(k=1;k<=n;k++) 3 for(i=1;i<=n;i++) 4 for(j=1;j<=n;j++) 5 if(e[i][k]<inf && e[k][j]<inf && e[i][j]>e[i][k]+e[k][j]) 6 e[i][j]=e[i][k]+e[k][j];
上面代码的输入数据样式为:
1 4 8 2 1 2 2 3 1 3 6 4 1 4 4 5 2 3 3 6 3 1 7 7 3 4 1 8 4 1 5 9 4 3 12
第一行两个数为n和m,n表示顶点个数,m表示边的条数。
![081030is22w3mmnz3r33m3.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081030is22w3mmnz3r33m3.png)
![081030elthvel6et6k886y.png](http://bbs.ahalei.com/data/attachment/forum/201403/25/081030elthvel6et6k886y.png)
此算法由Robert W. Floyd(罗伯特·弗洛伊德)于1962年发表在“Communications of the ACM”上。同年Stephen Warshall(史蒂芬·沃舍尔)也独立发表了这个算法。Robert W.Floyd这个牛人是朵奇葩,他原本在芝加哥大学读的文学,但是因为当时美国经济不太景气,找工作比较困难,无奈之下到西屋电气公司当了一名计算机操作员,在IBM650机房值夜班,并由此开始了他的计算机生涯。此外他还和J.W.J. Williams(威廉姆斯)于1964年共同发明了著名的堆排序算法HEAPSORT。堆排序算法我们将在第七章学习。Robert W.Floyd在1978年获得了图灵奖。
作 者:Angel_Kitty
出 处:https://www.cnblogs.com/ECJTUACM-873284962/
关于作者:阿里云ACE,目前主要研究方向是Web安全漏洞以及反序列化。如有问题或建议,请多多赐教!
版权声明:本文版权归作者和博客园共有,欢迎转载,但未经作者同意必须保留此段声明,且在文章页面明显位置给出原文链接。
特此声明:所有评论和私信都会在第一时间回复。也欢迎园子的大大们指正错误,共同进步。或者直接私信我
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角【推荐】一下。您的鼓励是作者坚持原创和持续写作的最大动力!
欢迎大家关注我的微信公众号IT老实人(IThonest),如果您觉得文章对您有很大的帮助,您可以考虑赏博主一杯咖啡以资鼓励,您的肯定将是我最大的动力。thx.
![](https://images.cnblogs.com/cnblogs_com/ECJTUACM-873284962/1318325/o_donate.png)
我的公众号是IT老实人(IThonest),一个有故事的公众号,欢迎大家来这里讨论,共同进步,不断学习才能不断进步。扫下面的二维码或者收藏下面的二维码关注吧(长按下面的二维码图片、并选择识别图中的二维码),个人QQ和微信的二维码也已给出,扫描下面👇的二维码一起来讨论吧!!!
![](https://images.cnblogs.com/cnblogs_com/ECJTUACM-873284962/1318325/o_chat.png)
欢迎大家关注我的Github,一些文章的备份和平常做的一些项目会存放在这里。