Floyd算法详解——包括解题步骤与编程

最新推荐文章于 2023-06-27 14:45:43 发布

HOLD ON!

最新推荐文章于 2023-06-27 14:45:43 发布

阅读量1.3k

点赞数

分类专栏：算法

原文链接：https://blog.youkuaiyun.com/SweeNeil/article/details/88713677

版权

算法专栏收录该内容

46 篇文章

订阅专栏

Floyd算法详解——包括解题步骤与编程

SweeNeil
展开
一、Floyd算法原理
Floyd算法是一个经典的动态规划算法，它又被称为插点法。该算法名称以创始人之一、1978年图灵奖获得者、斯坦福大学计算机科学系教授罗伯特·弗洛伊德命名。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的最短路径的距离。

二、Floyd算法内容步骤
（一）Floyd算法内容

（二）Floyd算法步骤

三、Floyd算法解题步骤
以下纯手工操作~

上面最后有个问题，更正一下，应该是

eg：d13 = 45，R13 = 5，R15 = 5（或者）R51 = 1，也就是13之间经过5，然后再看15之间有没有经过其他点，如果经过了继续查表，由于通过R15或R51可以看出，他们之间没有经过其他点，所以最短路径为1-》5-》3，最短距离为45.

四、Floyd算法编程
编程求解几座城市之间的最短距离，以及最短距离所经过的城市。

floyd.c：

#include <string.h>
#include <stdio.h>

#define NUMS 12
#define INF 65535

typedef struct
{
   char vertex[NUMS];
   int edges[NUMS][NUMS];
   int n,e;
}Graph;

void Dispath(int A[][NUMS],int path[][NUMS],int n);

void ReadGraph(Graph *G)
{
   int i,j;
   FILE * fp = fopen("floyd.txt","rw");
   G->n = NUMS;
   G->e = NUMS * NUMS;
   for(i=0; i<NUMS; i++)
   {
       for(j=0; j<NUMS; j++)
       {
           fscanf(fp,"%d",&(G->edges[i][j]));
           printf("%d \t",G->edges[i][j]);
       }
       printf("\n");
   }
}

void Floyd(Graph G)
{
   int A[NUMS][NUMS],path[NUMS][NUMS];
   int i,j,k;
   for (i=0;i<G.n;i++)
   {
       for (j=0;j<G.n;j++)
       {
           A[i][j]=G.edges[i][j];
           path[i][j]=-1;
       }
   }
   for (k=0;k<G.n;k++)
   {
       for (i=0;i<G.n;i++)
       {
           for (j=0;j<G.n;j++)
           {
               if (A[i][j]>A[i][k]+A[k][j])
               {
                   A[i][j]=A[i][k]+A[k][j];
                   path[i][j]=k;
               }
           }
       }
   }
   Dispath(A,path,G.n);
}

void Ppath(int path[][NUMS],int i,int j)
{
   int k;
   k=path[i][j];
   if (k==-1)
   {
       return;
   }
   Ppath(path,i,k);
   printf("%d,",k + 1);
   Ppath(path,k,j);
}

void Dispath(int A[][NUMS],int path[][NUMS],int n)
{
   int i,j;
   for (i=0;i<n;i++)
   {
       for (j=0;j<n;j++)
       {
           if (A[i][j]==INF)
           {
               if (i!=j)
               {
                   printf("从%d到%d没有路径\n",i+1,j+1);
               }
           }
           else
           {
               printf(" 从%d 到 %d => 最短路径长度为 : %d , 路径站点为 :",i+1,j+1,A[i][j]);
               printf("%d,",i + 1);
               Ppath(path,i,j);
               printf("%d\n",j + 1);
           }
       }
   }
}

int main()
{
   Graph G;
   ReadGraph(&G);
   Floyd(G);
   return 0;
}
城市数据 floyd.txt 如下：

0   20   65535   35   30   65535   65535   65535   65535   65535   65535   65535
20   0   25   65535   65535   65535   65535   65535   65535   65535   65535   65535
65535   25   0   5   65535   22   65535   65535   60   65535   75   65535
35   65535   5   0   65535   20   65535   65535   40   65535   65535   65535
30   65535   65535   65535   0   25   20   30   65535   65535   65535   65535
65535   65535   22   20   25   65535   12   40   65535   65535   65535   65535
65535   65535   65535   65535   20   65535   0   65535   65535   40   65535   65535
65535   65535   65535   65535   30   12   65535   0   35   65535   65535   90
65535   65535   65535   65535   65535   40   65535   35   0   30   35   45
65535   65535   65535   65535   65535   65535   40   65535   30   0   10   20
65535   65535   65535   65535   65535   65535   65535   65535   35   10   0   15
65535   65535   65535   65535   65535   65535   65535   90   45   20   15   0
其中 65535 表示不可直达，运行结果如下：

————————————————
版权声明：本文为优快云博主「SweeNeil」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.youkuaiyun.com/SweeNeil/article/details/88713677