2. 程式人生 > >【啊哈！算法】算法6：只有五行的Floyd最短路算法

【啊哈！算法】算法6：只有五行的Floyd最短路算法

阿新 • • 發佈：2018-01-28
此外 can oom 定義 其中 存在 i++ printf ons

技術分享圖片 暑假，小哼準備去一些城市旅遊。有些城市之間有公路，有些城市之間則沒有，如下圖。為了節省經費以及方便計劃旅程，小哼希望在出發之前知道任意兩個城市之前的最短路程。 技術分享圖片

上圖中有4個城市8條公路，公路上的數字表示這條公路的長短。請註意這些公路是單向的。我們現在需要求任意兩個城市之間的最短路程，也就是求任意兩個點之間的最短路徑。這個問題這也被稱為“多源最短路徑”問題。 現在需要一個數據結構來存儲圖的信息，我們仍然可以用一個4*4的矩陣（二維數組e）來存儲。比如1號城市到2號城市的路程為2，則設e[1][2]的值為2。2號城市無法到達4號城市，則設置e[2][4]的值為∞。另外此處約定一個城市自己是到自己的也是0，例如e[1][1]為0，具體如下。
<ignore_js_op>技術分享圖片 現在回到問題：如何求任意兩點之間最短路徑呢？通過之前的學習我們知道通過深度或廣度優先搜索可以求出兩點之間的最短路徑。所以進行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。通過這個的例子，我們發現每個頂點都有可能使得另外兩個頂點之間的路程變短。好，下面我們將這個問題一般化。
當任意兩點之間不允許經過第三個點時，這些城市之間最短路程就是初始路程，如下。

<ignore_js_op>技術分享圖片 假如現在只允許經過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號頂點的情況下，任意兩點之間的最短路程更新為：

技術分享圖片 通過上圖我們發現：在只通過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];
  6. //經過2號頂點
  7. for(i=1;i<=n;i++)
  8. for(j=1;j<=n;j++)
  9. 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號頂點進行中轉的情況下，求任意兩點之間的最短路程。任意兩點之間的最短路程更新為： 技術分享圖片 最後允許通過所有頂點作為中轉，任意兩點之間最終的最短路程為： 技術分享圖片 整個算法過程雖然說起來很麻煩，但是代碼實現卻非常簡單，核心代碼只有五行：
  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號點的最短路程。其實這是一種“動態規劃”的思想，關於這個思想我們將在《啊哈！算法2——偉大思維閃耀時》在做詳細的討論。下面給出這個算法的完整代碼：
  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);
  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;
  15. //讀入邊
  16. for(i=1;i<=m;i++)
  17. {
  18. scanf("%d %d %d",&t1,&t2,&t3);
  19. e[t1][t2]=t3;
  20. }
  22. //Floyd-Warshall算法核心語句
  23. for(k=1;k<=n;k++)
  24. for(i=1;i<=n;i++)
  25. for(j=1;j<=n;j++)
  26. if(e[i][j]>e[i][k]+e[k][j] )
  27. e[i][j]=e[i][k]+e[k][j];
  29. //輸出最終的結果
  30. for(i=1;i<=n;i++)
  31. {
  32. for(j=1;j<=n;j++)
  33. {
  34. printf("%10d",e[i][j]);
  35. }
  36. printf("\n");
  37. }
  39. return 0;
  40. }
復制代碼

有一點需要註意的是：如何表示正無窮。我們通常將正無窮定義為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表示邊的條數。 接下來m行，每一行有三個數t1、t2 和t3，表示頂點t1到頂點t2的路程是t3。 得到最終結果如下： 技術分享圖片

通過這種方法我們可以求出任意兩個點之間最短路徑。它的時間復雜度是O(N3)。令人很震撼的是它竟然只有五行代碼，實現起來非常容易。正是因為它實現起來非常容易，如果時間復雜度要求不高，使用Floyd-Warshall來求指定兩點之間的最短路或者指定一個點到其余各個頂點的最短路徑也是可行的。當然也有更快的算法，請看下一節：Dijkstra算法。 另外需要註意的是：Floyd-Warshall算法不能解決帶有“負權回路”（或者叫“負權環”）的圖，因為帶有“負權回路”的圖沒有最短路。例如下面這個圖就不存在1號頂點到3號頂點的最短路徑。因為1->2->3->1->2->3->…->1->2->3這樣路徑中，每繞一次1->-2>3這樣的環，最短路就會減少1，永遠找不到最短路。其實如果一個圖中帶有“負權回路”那麽這個圖則沒有最短路。 <ignore_js_op>技術分享圖片 此算法由Robert W. Floyd（羅伯特·弗洛伊德）於1962年發表在“Communications of the ACM”上。同年Stephen Warshall（史蒂芬·沃舍爾）也獨立發表了這個算法。Robert W．Floyd這個牛人是朵奇葩，他原本在芝加哥大學讀的文學，但是因為當時美國經濟不太景氣，找工作比較困難，無奈之下到西屋電氣公司當了一名計算機操作員，在IBM650機房值夜班，並由此開始了他的計算機生涯。此外他還和J.W.J. Williams（威廉姆斯）於1964年共同發明了著名的堆排序算法HEAPSORT。堆排序算法我們將在第七章學習。Robert W．Floyd在1987年獲得了圖靈獎。



碼字不容易啊，轉載麻煩註明出處
【一周一算法】算法6：只有五行的Floyd最短路算法
http://bbs.ahalei.com/thread-4554-1-1.html

【啊哈！算法】算法6：只有五行的Floyd最短路算法

相關推薦

【啊哈！算法】算法6：只有五行的Floyd最短路算法

此外 can oom 定義 其中 存在 i++ printf ons 暑假，小哼準備去一些城市旅遊。有些城市之間有公路，有些城市之間則沒有，如下圖。為了節省經費以及方便計劃旅程，小哼希望在出發之前知道任意兩個城市之前的最短路程。

【啊哈！演算法】系列7：Dijkstra最短路演算法

        上週我們介紹了神奇的只有五行的Floyd最短路演算法，它可以方便的求得任意兩點的最短路徑，這稱為“多源最短路”。本週來來介紹指定一個點（源點）到其餘各個頂點的最短路徑，也叫做“單源最短路徑”。例如求下圖中的1號頂點到2、3、4、5、6號頂點的最短路徑。

【啊哈！演算法】演算法1：最快最簡單的排序——桶排序

《啊哈C》出版之後，很多網友希望能夠有加深的內容，比如資料結構、演算法之類的。今後每週五更新一篇吧。最快最簡單的排序——桶排序　　在我們生活的這個世界中到處都是被排序過的。站隊的時候會按照身高排序，考試的名次需要按照分數排序，網上購物的時候會按照價格排序，電子郵箱中的郵件按照

【啊哈！演算法】演算法2：氣泡排序

　　簡化版的桶排序不僅僅有上一節所遺留的問題，更要命的是：它非常浪費空間！例如需要排序數的範圍是0~2100000000之間，那你則需要申請2100000001個變數，也就是說要寫成int a[2100000001]。因為我們需要用2100000001個“桶”來儲存0~2

【啊哈！演算法】演算法6：只有五行的Floyd最短路演算法

        我們來想一想，根據我們以往的經驗，如果要讓任意兩點（例如從頂點a點到頂點b）之間的路程變短，只能引入第三個點（頂點k），並通過這個頂點k中轉即a->k->b，才可能縮短原來從頂點a點到頂點b的路程。那麼這個中轉的頂點k是1~n中的哪個點呢？甚至有時候不只通過一個點，而是經過兩個點

Floyd最短路算法的解釋.

bsp 先後 結點 優化 ini 空間 計算 方程 適用於 適用於有向/無向圖,本質上是一個動態規劃. D[k][i][j]代表經前k個結點中轉,i到j的距離.可以寫出方程: D[k][i][j]=min{D[k-1][i][j], D[k-1][i][k]+D[k-1][

《啊哈！算法》.啊哈磊.掃描版pdf

link tex 第7章 tom 國際 搜索引擎 微軟亞洲研究院 right 作者 下載地址：網盤下載 內容簡介 · · · · · · 這不過是一本有趣的算法書而已。和別的算法書比較，如果硬要說它有什麽特點的話，那就是你能看懂它。 這是一本

[啊哈！算法].啊哈磊.掃描版高清版pdf免費下載

鏈表 城市 c語言 dot alt 語言 ima 二叉 國際 下載地址：網盤下載 備用地址：網盤下載 內容簡介 · · · · · ·這不過是一本有趣的算法書而已。和別

【圖】實驗題二：實現圖的遍歷算法

malloc for clu %d type font AS arc -- 遞歸深度優先遍歷算法 1 #include<stdio.h> 2 #include<malloc.h> 3 #define MAXV 100 4 #def

震驚，最短路算法！！！！

keyword 信息 dex oom 沒有 true 進行 ted 減少 轉載 http://www.61mon.com/index.php/archives/194/ 文章目錄 一：背景 二：算法過程 三：完整代碼 四：時間復雜度 五：該算法的缺陷

啊哈！演算法】最快最簡單的排序——桶排序

最快最簡單的排序——桶排序 　　在我們生活的這個世界中到處都是被排序過的。站隊的時候會按照身高排序，考試的名次需要按照分數排序，網上購物的時候會按照價格排序，電子郵箱中的郵件按照時間排序……總之很多東西都需要排序，可以說排序是無處不在。現在我們舉個具體的例子來介紹一下

最短路算法

class str word i++ spa ret clas mic func Floyd（hdu 1874） for(int k=0;k<n;k++) for(int i=0;i<n;i++) for(int j=0;j<n;j++) e

Dijkstra最短路算法詳解

logs alt 成了 解決 img eof cnblogs 這就是 時間 想必大家一定會Floyd了吧，Floyd只要暴力的三個for就可以出來，代碼好背，也好理解，但缺點就是時間復雜度高是O(n3)。 於是今天就給大家帶來一種時間復雜度是O(n2)，的算法：Di

（轉）最短路算法 -- Floyd算法

無法 著名 jks href cat mda floyd算法 第七章 blank 轉自：http://blog.51cto.com/ahalei/1383613 暑假，小哼準備去一些城市旅遊。有些城市之間有公路，有些城市之間則沒有，如下圖。為了節省經費以

最短路算法 -- SPFA模板

const 結點 con set spf () 最短路算法 tor acm 一、算法步驟 建立一個隊列，初始時隊列裏只有起始點，再建立一個數組記錄起始點到所有點的最短路徑（該數組的初始值要賦為極大值，該點到它本身的路徑賦為0，下面的模板中該數組為dist[]）。然後執行松

spfa 算法模板 可求帶負權邊的最短路

cto nbsp init ems push name for 如果 class 　　它是隊列優化的Bellman-Ford算法。 　　優化的原理是：下一次松弛操作時被更新dis的點其實與上一次被更新的點有關！如果上一次被更新的點有一條邊指向某點V，那麽在下一次，點V就是可

最短路算法（floyed+Dijkstra+bellman-ford+SPFA）

bool logs pan 動態規劃 優點 push i++ img tails 最短路算法簡單模板 一.floyed算法 首先對於floyed算法來說就是最短路徑的動態規劃解法，時間復雜度為O（n^3） 適用於圖中所有點與點之間的最短路徑的算法，一般適用於點n較小的情

最短路算法 —— Dijkstra算法

math bool see namespace bsp 算法 鄰接矩陣 最短路徑 eat 用途： 解決單源最短路徑問題（已固定一個起點，求它到其他所有點的最短路問題） 算法核心（廣搜）： （1）確定的與起點相鄰的點的最短距離，再根據已確定最短距離的點更新其他與之相鄰的點的最

算法系列-複雜度分析：淺析最好、最壞、平均、均攤時間複雜度

整理自極客時間-資料結構與演算法之美。原文內容更完整具體，且有音訊。購買地址： 上一節，我們講了複雜度的大 O 表示法和幾個分析技巧，還舉了一些常見覆雜度分析的例子，比如 O(1)、O(logn)、O(n)、O(nlogn) 複雜度分析。掌握了這些內容，對於複雜度分析這個知識點，你已經

【ROS-MoveIt！原始碼學習】ROS中機器人模型的構建（Build RobotModel）

1.本文討論的內容和目的 本文從ROS MOVEIT!的入門教程moveit_tutorial中的運動學模型（kinematic_model_tutorial.launch.cpp）入手，目的是想搞明白： （1）ROS是如何匯入機器人模型；（本文內容） （2）如何呼叫逆運動學演算法求