无负权值最短路径问题的Floyd算法。
二维数组d表示有向图边信息,d[i][j]表示从结点i到结点j的最短路径长度,若i==j,则d[i][j]=0,若i,j之间无直接可达路径,则d[i][j]取值无穷大(为计算方便,取值1000000)
二维数组p表示路径走向信息,若结点i直接可达结点j,则p[i][j]=j,否则p[i][j]=0
计算如下有向图中任意两点的最短路径长度和走向:
define H 1000000
d的初始值:d = {
{0,7,5,12,H,H},
{H,0,H,3,H,H},
{H,H,0,6,H,15},
{12,H,6,0,H,8},
{6,H,13,H,0,H},
{H,4,15,H,9,0}}
p的初始值:p = {
{1,2,3,4,0,0},
{0,2,0,4,0,0},
{0,0,3,4,0,6},
{1,0,3,4,0,6},
{1,0,3,0,5,0},
{0,2,3,0,5,6}};
Floyd算法步骤:
1、边信息矩阵d,路径信息矩阵p
2、i = 1,j = 1,k = 1
3、u = d[i][k]+d[k][j],若u < d[i][j],d[i][j] = u,p[i][j] = p[i][k];否则转到4
4、若j < n,则j = j + 1,转到3;否则转到5
5、若i < n, 则i = i + 1,转到3;否则转到6
6、若k < n, 则k = k + 1,i = 1,j = 1,转到3,否则end
其实就是一个三重循环。代码如下:
#include<iostream>
using namespace std;
#define H 100000
void fun(int d[6][6],int p[6][6])
{
int i,j,k;
for(k=0;k<6;k++)
{
for(i=0;i<6;i++)
{
for(j=0;j<6;j++)
{
if(d[i][k]+d[k][j]<d[i][j])
{
d[i][j]=d[i][k]+d[k][j];
p[i][j]=p[i][k];
}
}
}
}
}
int main()
{
int d[6][6] = {
{0,7,5,12,H,H},
{H,0,H,3,H,H},
{H,H,0,6,H,15},
{12,H,6,0,H,8},
{6,H,13,H,0,H},
{H,4,15,H,9,0}};
int p[6][6] = {
{1,2,3,4,0,0},
{0,2,0,4,0,0},
{0,0,3,4,0,6},
{1,0,3,4,0,6},
{1,0,3,0,5,0},
{0,2,3,0,5,6}};
fun(d,p);
int i,j;
for(i=0;i<6;i++)
{
for(j=0;j<6;j++)
cout<<d[i][j]<<' ';
cout<<endl;
}
cout<<endl;
for(i=0;i<6;i++)
{
for(j=0;j<6;j++)
cout<<p[i][j]<<' ';
cout<<endl;
}
return 0;
}
此时边信息矩阵d中d[i][j]即表示从结点i到结点j的最大路径长度。而路径信息矩阵p记录路径走向。
比如查找结点1到结点5的最短路径走向:
第1行第5列,值为2,表示经过结点2;第2行第5列,值为4,表示经过结点4;第4行第5列,值为6,表示经过结点6,第6行第5列,值为5结束。
得到路径走向:1->2->4->6->5
Floyd算法,可以得到所以结点之间的最短路径长度及走向(若结点之间可达)。