最短路之——弗洛伊德算法(floyd)

来源: https://blog.csdn.net/riba2534/article/details/54562440


我们要做的是求出从某一点到达任意一点的最短距离,我们先用邻接矩阵来建图,map[i][j]表示从i点到j点的距离,把自己到自己设为0,把自己到不了的边初始化为无穷大,代码为:

[cpp]  view plain  copy
  1. //初始化  
  2.     for(int i=1; i<=n; i++)  
  3.         for(int j=1; j<=n; j++)  
  4.             if(i==j)  
  5.                 map[i][j]=0;  
  6.             else  
  7.                 map[i][j]=inf;  
  8.     //读入边  
  9.     for(int i=1; i<=m; i++)  
  10.     {  
  11.         scanf("%d%d%d",&t1,&t2,&t3);  
  12.         map[t1][t2]=t3;  
  13.     }  
最后,建好的图可以用表格来表示:

现在,我们来思考,假设我们来找一个中转的点,看他们的路程会不会改变,我们先以1号顶点作为中转点最为例子,制图:


我们发现,图有了变化,我们怎么判断以1号顶点作为中转点图的路程是不是更短呢,我们只需要判断map[i][1]+map[1][j]的路程是不是比map[i][j]的路程更短,就可以判断,

代码为:

[cpp]  view plain  copy
  1. for(int i=1; i<=n; i++)  
  2.     for(int j=1; j<=n; j++)  
  3.         if(map[i][1]+map[1][j]<map[i][j])  
  4.             map[i][j]=map[i][1]+map[1][j];  
现在该怎么办呢,我们接着以2号顶点作为中转点,很简单代码修改一句就就可以:

[cpp]  view plain  copy
  1. for(int i=1; i<=n; i++)  
  2.     for(int j=1; j<=n; j++)  
  3.         if(map[i][2]+map[2][j]<map[i][j])  
  4.             map[i][j]=map[i][2]+map[2][j];  
现在我们是不是发现了一个规律,只要不断的遍历每一个点,并且以每一个点作为中转点看看它的值会不会改变,就可以得到从一个点到任意一个点的最短路径,也就是多源最短路,这就是弗洛伊德算法,代码为:

[cpp]  view plain  copy
  1. for(int k=1; k<=n; k++)  
  2.     for(int i=1; i<=n; i++)  
  3.         for(int j=1; j<=n; j++)  
  4.             if(map[i][k]+map[k][j]<map[i][j])  
  5.                 map[i][j]=map[i][k]+map[k][j];  
这样就可以遍历每个顶点,找出所有的最短路,算法的复杂度为O(n^3).

对于我一开始提出的问题,完整的代码为:

[cpp]  view plain  copy
  1. #include <stdio.h>  
  2. #include <string.h>  
  3. #include <string>  
  4. #include <iostream>  
  5. #include <stack>  
  6. #include <queue>  
  7. #include <vector>  
  8. #include <algorithm>  
  9. #define mem(a,b) memset(a,b,sizeof(a))  
  10. using namespace std;  
  11. const int inf=1<<29;  
  12. int main()  
  13. {  
  14.     int map[10][10],n,m,t1,t2,t3;  
  15.     scanf("%d%d",&n,&m);//n表示顶点个数,m表示边的条数  
  16.     //初始化  
  17.     for(int i=1; i<=n; i++)  
  18.         for(int j=1; j<=n; j++)  
  19.             if(i==j)  
  20.                 map[i][j]=0;  
  21.             else  
  22.                 map[i][j]=inf;  
  23.     //读入边  
  24.     for(int i=1; i<=m; i++)  
  25.     {  
  26.         scanf("%d%d%d",&t1,&t2,&t3);  
  27.         map[t1][t2]=t3;  
  28.     }  
  29.     //弗洛伊德(Floyd)核心语句  
  30.     for(int k=1; k<=n; k++)  
  31.         for(int i=1; i<=n; i++)  
  32.             for(int j=1; j<=n; j++)  
  33.                 if(map[i][k]+map[k][j]<map[i][j])  
  34.                     map[i][j]=map[i][k]+map[k][j];  
  35.     for(int i=1; i<=n; i++)  
  36.     {  
  37.         for(int j=1; j<=n; j++)  
  38.             printf("%10d",map[i][j]);  
  39.         printf("\n");  
  40.     }  
  41.     return 0;  
  42. }  
相关推荐
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页