漫画:如何求图的最短路径? | 技术头条

640?wx_fmt=gif

640?wx_fmt=jpeg

作者 | 程序员小灰

本文经授权转载自程序员小灰(ID:chengxuyuanxiaohui)

640?wx_fmt=jpeg

640?wx_fmt=jpeg

—————  第二天  —————

640?wx_fmt=jpeg

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

640?wx_fmt=png

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

 

小灰的思路如下:

第一步,利用迪杰斯特拉算法的距离表,求出从顶点A出发,到其他各个顶点的最短距离:

640?wx_fmt=png

第二步,继续使用迪杰斯特拉算法,求出从顶点B出发,到其他各个顶点的最短距离。

第三步,从顶点C出发,到各个顶点的最短距离。

第四步,从顶点D出发......

.......

就像这样,一直遍历到顶点G。

这个思路的时间复杂度是多少呢?

假如图中有n个顶点,如果不考虑堆优化,一次迪杰斯特拉算法的时间复杂度是O(n^2)。所以,把每一个顶点都计算一遍,总的时间复杂度是O(n^3)。

640?wx_fmt=jpeg

640?wx_fmt=jpeg

640?wx_fmt=jpeg

640?wx_fmt=jpeg

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

 

640?wx_fmt=jpeg

举一个栗子:

640?wx_fmt=png

上图的顶点A和顶点C没有直接相连的边,它们之间的直接距离是无穷大。

如果以B作为“中继顶点”,此时A到C的最短路径就是A-B-C,最短距离是3+2=5。

640?wx_fmt=jpeg

再举一个栗子:

640?wx_fmt=png

上图的顶点A和顶点C直接相连,距离是6。但是存在一条“迂回”路径A-B-C,距离是3+2=5<6。

所以,经过中继顶点B,从A到C的最短距离可以是5。

640?wx_fmt=jpeg

下面我们来看一看Floyd算法的详细步骤。

1.要实现Floyd算法,首先需要构建带权图的邻接矩阵:

640?wx_fmt=png

在邻接矩阵当中,每一个数字代表着从某个顶点到另一个顶点的直接距离,这个距离是没有涉及到任何中继顶点的。

2.此时假定只允许以顶点A作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

B和C之间的距离原本是无穷大,此时以A为中继,距离缩短为AB距离+AC距离=

5+2=7。

更新对应矩阵元素(橙色区域代表顶点A到其他顶点的临时距离):

640?wx_fmt=png

3.接下来以顶点A、B作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

A和D之间的距离原本是无穷大,此时以B为中继,距离缩短为AB距离+BD距离=5+1=6。

A和E之间的距离原本是无穷大,此时以B为中继,距离缩短为AB距离+BE距离=5+6=11。

更新对应矩阵元素(橙色区域代表顶点B到其他顶点的临时距离):

640?wx_fmt=png

4.接下来以顶点A、B、C作为中继顶点,那么各顶点之间的距离会变成什么样子呢?

A和F之间的距离原本是无穷大,此时以C为中继,距离缩短为AC距离+CF距离=2+8=10。

更新对应矩阵元素(橙色区域代表顶点C到其他顶点的临时距离):

640?wx_fmt=png

以此类推,我们不断引入新的中继顶点,不断刷新矩阵中的临时距离。

最终,当所有顶点都可以作为中继顶点时,我们的距离矩阵更新如下:

640?wx_fmt=png

此时,矩阵中每一个元素,都对应着某顶点到另一个顶点的最短距离。

640?wx_fmt=jpeg

640?wx_fmt=jpeg

为什么这么说呢?让我们回顾一下动态规划的两大要素:

问题的初始状态
问题的状态转移方程式

对于寻找图的所有顶点之间距离的问题,初始状态就是顶点之间的直接距离,也就是邻接矩阵。

而问题的状态转移方程式又是什么呢?

假设新引入的中继顶点是n,那么:

顶点i 到 顶点j 的新距离 = Min(顶点i 到 顶点j 的旧距离,顶点i 到 顶点n 的距离+顶点n 到 顶点j 的距离)

640?wx_fmt=jpeg

 
  1.  

 
  1. final static int INF = Integer.MAX_VALUE;

  2.  

  3. public static void floyd(int[][] matrix){

  4. //循环更新矩阵的值

  5. for(int k=0; k<matrix.length; k++){

  6. for(int i=0; i<matrix.length; i++){

  7. for(int j=0; j<matrix.length; j++){

  8. if(matrix[i][k] == INF || matrix[k][j] == INF) {

  9. continue;

  10. }

  11. matrix[i][j] = Math.min(matrix[i][j], matrix[i][k] + matrix[k][j]);

  12. }

  13. }

  14. }

  15. // 打印floyd最短路径的结果

  16. System.out.printf("最短路径矩阵: \n");

  17. for (int i = 0; i < matrix.length; i++) {

  18. for (int j = 0; j < matrix.length; j++)

  19. System.out.printf("%3d ", matrix[i][j]);

  20. System.out.printf("\n");

  21. }

  22. }

  23.  

 

  1.  

  2.  

  3. public static void main(String[] args) {

  4. int[][] matrix = {

  5. {0, 5, 2, INF, INF, INF, INF},

  6. {5, 0, INF, 1, 6, INF, INF},

  7. {2, INF, 0, 6, INF, 8, INF},

  8. {INF, 1, 6, 0, 1, 2, INF},

  9. {INF, 6, INF, 1, 0, INF, 7},

  10. {INF, INF, 8, 2, INF, 0, 3},

  11. {INF, INF, INF, INF, 7, 3, 0}

  12. };

  13. floyd(matrix);

  14.  

    }

     

640?wx_fmt=jpeg

640?wx_fmt=jpeg

2019我为什么建议你学Python?

https://edu.csdn.net/topic/python115?utm_source=csdn_bw

 

【END】

640?wx_fmt=jpeg

作为码一代,想教码二代却无从下手:

听说少儿编程很火,可它有哪些好处呢?

孩子多大开始学习比较好呢?又该如何学习呢?

最新的编程教育政策又有哪些呢?

下面给大家介绍CSDN新成员:极客宝宝(ID:geek_baby)

戳他了解更多↓↓↓

640?wx_fmt=jpeg

 热 文 推 荐 

☞ 手机基带芯片激荡 30 年!

 硬核追星!杨超越杯决赛是些啥编程项目?

 前端代码的整洁之道 | 技术头条

准码农自述:我是如何从二本到日本留学,从 C 语言到人工智能的?

☞ 阿里半跪过,任正非差点跳楼,京东被骗光钱:成年人的生活哪有容易二字?

☞ 回报率29%! 大神用情感分析创建一个比特币交易算法, 原来交易玩的是心理战

☞ 硬核粉丝 | 清华双胞胎“YCY Dance Now”杀进超越杯编程大赛决赛

☞ 异构计算=未来?一文带你秒懂3大主流异构

☞ 她说:为啥程序员都特想要机械键盘?这答案我服!

 

System.out.println("点个在看吧!");
console.log("点个在看吧!");
print("点个在看吧!");
printf("点个在看吧!\n");
cout << "点个在看吧!" << endl;
Console.WriteLine("点个在看吧!");
Response.Write("点个在看吧!");
alert("点个在看吧!")
echo "点个在看吧!"

 
你点的每个“在看”,我都认真当成了喜欢
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值