暑假，小哼准备去一些城市旅游。有些城市之间有公路，有些城市之间则没有，如下图。为了节省经费以及方便计划旅程，小哼希望在出发之前知道任意两个城市之前的最短路程。
        上图中有4个城市8条公路，公路上的数字表示这条公路的长短。请注意这些公路是单向的。我们现在需要求任意两个城市之间的最短路程，也就是求任意两个点之间的最短路径。这个问题这也被称为“多源最短路径”问题。
        现在需要一个数据结构来存储图的信息，我们仍然可以用一个4*4的矩阵（二维数组e）来存储。比如1号城市到2号城市的路程为2，则设e[1][2]的值为2。2号城市无法到达4号城市，则设置e[2][4]的值为∞。另外此处约定一个城市自己是到自己的也是0，例如e[1][1]为0，具体如下。
        现在回到问题：如何求任意两点之间最短路径呢？通过之前的学习我们知道通过深度或广度优先搜索可以求出两点之间的最短路径。所以进行n2遍深度或广度优先搜索，即对每两个点都进行一次深度或广度优先搜索，便可以求得任意两点之间的最短路径。可是还有没有别的方法呢？
        我们来想一想，根据我们以往的经验，如果要让任意两点（例如从顶点a点到顶点b）之间的路程变短，只能引入第三个点（顶点k），并通过这个顶点k中转即a->k->b，才可能缩短原来从顶点a点到顶点b的路程。那么这个中转的顶点k是1~n中的哪个点呢？甚至有时候不只通过一个点，而是经过两个点或者更多点中转会更短，即a->k1->k2b->或者a->k1->k2…->k->i…->b。比如上图中从4号城市到3号城市（4->3）的路程e[4][3]原本是12。如果只通过1号城市中转（4->1->3），路程将缩短为11（e[4][1]+e[1][3]=5+6=11）。其实1号城市到3号城市也可以通过2号城市中转，使得1号到3号城市的路程缩短为5（e[1][2]+e[2][3]=2+3=5）。所以如果同时经过1号和2号两个城市中转的话，从4号城市到3号城市的路程会进一步缩短为10。通过这个的例子，我们发现每个顶点都有可能使得另外两个顶点之间的路程变短。好，下面我们将这个问题一般化。
        当任意两点之间不允许经过第三个点时，这些城市之间最短路程就是初始路程，如下。
 

假如现在只允许经过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循环，代码实现如下。

for(i=1;i<=n;i++)

[cpp] view plaincopy
  1. {
  2. for(j=1;j<=n;j++)
  3. {
  4. if ( e[i][j] > e[i][1]+e[1][j] )
  5. e[i][j] = e[i][1]+e[1][j];
  6. }
  7. }

在只允许经过1号顶点的情况下，任意两点之间的最短路程更新为：

        通过上图我们发现：在只通过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]要小，代码实现为如下。
[cpp] view plaincopy
  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号顶点的情况下，任意两点之间的最短路程更新为：

        通过上图得知，在相比只允许通过1号顶点进行中转的情况下，这里允许通过1和2号顶点进行中转，使得e[1][3]和e[4][3]的路程变得更短了。
        同理，继续在只允许经过1、2和3号顶点进行中转的情况下，求任意两点之间的最短路程。任意两点之间的最短路程更新为：
        最后允许通过所有顶点作为中转，任意两点之间最终的最短路程为：
        整个算法过程虽然说起来很麻烦，但是代码实现却非常简单，核心代码只有五行：
[cpp] view plaincopy
  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号点的最短路程。

[cpp] view plaincopy
  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. for(i=1;i<=n;i++)
  10. for(j=1;j<=n;j++)
  11. if(i==j) e[i][j]=0;
  12. else e[i][j]=inf;
  13. //读入边
  14. for(i=1;i<=m;i++)
  15. {
  16. scanf("%d %d %d",&t1,&t2,&t3);
  17. e[t1][t2]=t3;
  18. }
  19. //Floyd-Warshall算法核心语句
  20. for(k=1;k<=n;k++)
  21. for(i=1;i<=n;i++)
  22. for(j=1;j<=n;j++)
  23. if(e[i][j]>e[i][k]+e[k][j] )
  24. e[i][j]=e[i][k]+e[k][j];
  25. //输出最终的结果
  26. for(i=1;i<=n;i++)
  27. {
  28. for(j=1;j<=n;j++)
  29. {
  30. printf("%10d",e[i][j]);
  31. }
  32. printf("\n");
  33. }
  34. return 0;
  35. }
        另外需要注意的是：Floyd-Warshall算法不能解决带有“负权回路”（或者叫“负权环”）的图，因为带有“负权回路”的图没有最短路。例如下面这个图就不存在1号顶点到3号顶点的最短路径。因为1->2->3->1->2->3->…->1->2->3这样路径中，每绕一次1->-2>3这样的环，最短路就会减少1，永远找不到最短路。其实如果一个图中带有“负权回路”那么这个图则没有最短路。

傻子也能看懂的弗洛伊德算法（转）相关推荐

  1. Floyd-傻子也能看懂的弗洛伊德算法

    暑假,小哼准备去一些城市旅游.有些城市之间有公路,有些城市之间则没有,如下图.为了节省经费以及方便计划旅程,小哼希望在出发之前知道任意两个城市之前的最短路程. 上图中有4个城市8条公路,公路上的数字表 ...

  2. Floyd-傻子也能看懂的弗洛伊德算法（转）

    暑假,小哼准备去一些城市旅游.有些城市之间有公路,有些城市之间则没有,如下图.为了节省经费以及方便计划旅程,小哼希望在出发之前知道任意两个城市之前的最短路程. 上图中有4个城市8条公路,公路上的数字表 ...

  3. 最短路径算法 Floyd-傻子也能看懂的弗洛伊德算法（转）

    讲的清楚的1 https://www.cnblogs.com/wangyuliang/p/9216365.html 讲的清楚的2 ahttps://zhuanlan.zhihu.com/p/13911 ...

  4. 人人都能看懂的EM算法推导

    作者丨August@知乎(已授权) 来源丨https://zhuanlan.zhihu.com/p/36331115 编辑丨极市平台 估计有很多入门机器学习的同学在看到EM算法的时候会有种种疑惑:EM ...

  5. em算法 实例 正态分布_人人都能看懂的EM算法推导

    ↑ 点击蓝字 关注极市平台作者丨August@知乎(已授权)来源丨https://zhuanlan.zhihu.com/p/36331115编辑丨极市平台 极市导读 EM算法到底是什么,公式推导怎么去 ...

  6. 争取能让大家都能看懂的 DFA 算法

    为什么要学习这个算法 我们公司一直都有的一个敏感词检测服务,前一段时间遇到了瓶颈,因为词库太多了导致会有一些速度过慢,而且一个正则表达式已经放不下了,需要进行拆分正则才可以. 正好我以前看过有关 df ...

  7. 我这种傻子都能看懂的 java oi流

    此文章较为适合我这种初学者 废话不多说直接进入正题 观看如下代码 public class 铜中 {public static void main(String []args) {File yf=ne ...

  8. 傻子也能看懂的迪杰斯特拉算法（转）

    本周来来介绍指定一个点(源点)到其余各个顶点的最短路径,也叫做"单源最短路径".例如求下图中的1号顶点到2.3.4.5.6号顶点的最短路径.           与Floyd-Wa ...

  9. 傻子都能看懂的马拉车Manacher

    Manacher's Algorithm 马拉车算法操作及原理 package advanced_001;public class Code_Manacher {public static char[ ...

  10. 匈牙利算法——你一定可以看懂的图论算法

    第一次讲解一个算法,也是最近才学到的一个算法,我的语言表达能力可能不是很强,排版也不是特别好,但是真的有认真在做一篇文章,如果有好的意见可以告诉我,非常感谢. 匈牙利算法是一种在多项式时间内求解任务分 ...

最新文章

  1. linux下使用extundelete恢复ext3/ext4分区数据
  2. C/C++位域结构深入解析
  3. React中props与state的区别
  4. MySQL数据存储目录查找
  5. 激战服务器位置,如何选服务器 《激战2》服务器设置讲解
  6. .net框架读书笔记---CLR内存管理\垃圾收集（二）
  7. 信息学奥赛一本通C++语言——1128：图像模糊处理
  8. mysql age字段类型_mysql 列类型
  9. sencha touch Container控件 setRecord 与 setData的区别
  10. BZOJ - 2819 Nim
  11. LibCef中的一些坑
  12. 卧槽，物色了一款隐秘拍摄神器，别乱用！
  13. 地图比例尺的计算（转）
  14. SpringBoot系列之Spring Data MongoDB教程
  15. 雷电2接口_厉害！跟Type-C长得一毛一样的雷电3，原来区别这么大，别搞混了
  16. Python xlwings Excel单元格部分内容添加下划线
  17. 苹果手机点击事件无效
  18. 记忆力很差怎么学计算机,九个方法增强记忆力 学生记忆力很差怎么办
  19. Android实现收款成功金额的语音播报功能(Nice tone)
  20. 大脑，字节，背部，Bun头-程序员的优先事项

热门文章

  1. 多数投票算法(Boyer-Moore Algorithm)
  2. Windows rundll32的用法-脚本命令管理打印机
  3. Javascript中大于和小于
  4. openwrt定时任务 2021-10-08
  5. 刻录光驱不能访问出现函数不正确问题解决
  6. 智者务其实，愚者争虚名
  7. 以太网未识别的网络win10_win10系统遇到以太网无法识别网络如何解决
  8. 带你学会C++文字页面类项目——3.答题器制作
  9. potoshope cs5 序列号
  10. 许远东受邀上海管理科技论坛做《LTD数字化经营方法论》分享