Floyd算法详解(附matlab代码)

什么是Floyd算法

Floyd算法又称为插点法,是一种利用动态规划的思想寻找给定的加权图中多源点之间最短路径的算法,与Dijkstra算法类似。该算法名称以创始人之一、1978年图灵奖获得者、斯坦福大学计算机科学系教授罗伯特·弗洛伊德命名。

Floyd算法可以给出网络中任意两个节点之间的最短路径,因此它是比Dijkstra更一般的算法。Floyd算法的思想是将n个节点的网络表示为n行n列的矩阵,而矩阵中的元素(i,j)表示从节点i到节点j的距离dij,如果两点直接没有边相连,则相应的元素就是无穷(∞).

算法的基本步骤

  • 定义n×n的方阵序列D-1, D0 , … Dn-1,

  • 初始化: D-1=C

    • D-1[i][j]=边<i,j>的长度,表示初始的从i到j的最短路径长度,即它是从i到j的中间不经过其他中间点的最短路径。
  • 迭代:设Dk-1已求出,如何得到Dk(0≤k≤n-1)?

    • Dk-1[i][j]表示从i到j的中间点不大于k-1的最短路径p:i…j,
    • 考虑将顶点k加入路径p得到顶点序列q:i…k…j,
    • 若q不是路径,则当前的最短路径仍是上一步结果:Dk[i][j]= Dk-1[i][j];
    • 否则若q的长度小于p的长度,则用q取代p作为从i到j的最短路径
  • 因为q的两条子路径i…k和k…j皆是中间点不大于k-1的最短路径,所以从i到j中间点不大于k的最短路径长度为:
    Dk[i][j]=min{ Dk-1[i][j], Dk-1[i][k] +Dk-1[k][j] }

实例讲解

 

 注:Floyd算法可以存在负权重

 matlab代码

function Floyd(w,router_direction,MAX)

%w为此图的距离矩阵

%router_direction为路由类型:0为前向路由;非0为回溯路由

%MAX是数据输入时的∞的实际值

len=length(w);

flag=zeros(1,len);

%根据路由类型初始化路由表

R=zeros(len,len);

for i=1:len

if router_direction==0%前向路由

R(:,i)=ones(len,1)*i;

else %回溯路由

R(i,:)=ones(len,1)*i;

end

R(i,i)=0;

end

disp('');

disp('w(0)');

dispit(w,0);

disp('R(0)');

dispit(R,1);

%处理端点有权的问题

for i=1:len

tmp=w(i,i)/2;

if tmp~=0

w(i,:)=w(i,:)+tmp;

w(:,i)=w(:,i)+tmp;

flag(i)=1;

w(i,i)=0;

end

end

%Floyd算法具体实现过程

for i=1:len

for j=1:len

if j==i || w(j,i)==MAX

continue;

end

for k=1:len

if k==i || w(j,i)==MAX

continue;

end

if w(j,i)+w(i,k)<w(j,k) %Floyd算法核心代码

w(j,k)=w(j,i)+w(i,k);

if router_direction==0%前向路由

R(j,k)=R(j,i);

else %回溯路由

R(j,k)=R(i,k);

end

end

end

end

%显示每次的计算结果

disp(['w(',num2str(i),')'])

dispit(w,0);

disp(['R(',num2str(i),')'])

dispit(R,1);

end

%中心和中点的确定

[Center,index]=min(max(w'));

disp(['中心是V',num2str(index)]);

[Middle,index]=min(sum(w'));

disp(['中点是V',num2str(index)]);

end

function dispit(x,flag)

%x:需要显示的矩阵

%flag:为0时表示显示w矩阵,非0时表示显示R矩阵

len=length(x);

s=[];

for j=1:len

if flag==0

s=[s sprintf('%5.2f\t',x(j,:))];

else

s=[s sprintf('%d\t',x(j,:))];

end

s=[s sprintf('\n')];

end

disp(s);

disp('---------------------------------------------------');

end

  • 4
    点赞
  • 25
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
根据提供的引用内容,没有提到Floyd pthread。但是可以根据Floyd算法的特点,使用多线程(pthread)来加速Floyd算法的执行。 Floyd算法的时间复杂度为O(N^3),对于大规模的图来说,计算时间会非常长。因此,可以使用多线程来加速Floyd算法的执行。具体来说,可以将图分成若干个子图,每个子图分配一个线程来计算,最后将结果合并即可。 下面是一个使用pthread实现的Floyd算法的示例代码: ```c #include <stdio.h> #include <stdlib.h> #include <pthread.h> #define MAXN 1000 #define MAX_THREAD_NUM 4 int n, m; int d[MAXN][MAXN]; int thread_num; void *floyd(void *arg) { int k = *(int *)arg; int i, j, l; for (k = k * n / thread_num; k < (k + 1) * n / thread_num; k++) { for (i = 0; i < n; i++) { for (j = 0; j < n; j++) { if (d[i][k] != -1 && d[k][j] != -1) { if (d[i][j] == -1 || d[i][j] > d[i][k] + d[k][j]) { d[i][j] = d[i][k] + d[k][j]; } } } } } return NULL; } void floyd_parallel() { pthread_t threads[MAX_THREAD_NUM]; int i; for (i = 0; i < thread_num; i++) { int *arg = malloc(sizeof(int)); *arg = i; pthread_create(&threads[i], NULL, floyd, arg); } for (i = 0; i < thread_num; i++) { pthread_join(threads[i], NULL); } } int main() { scanf("%d%d", &n, &m); int i, j; for (i = 0; i < n; i++) { for (j = 0; j < n; j++) { if (i == j) { d[i][j] = 0; } else { d[i][j] = -1; } } } for (i = 0; i < m; i++) { int u, v, w; scanf("%d%d%d", &u, &v, &w); d[u][v] = w; } scanf("%d", &thread_num); floyd_parallel(); for (i = 0; i < n; i++) { for (j = 0; j < n; j++) { printf("%d ", d[i][j]); } printf("\n"); } return 0; } ```
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

行走的参考文献

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值