参考《大话数据结构-弗洛伊德算法》
题目描述
在带权有向图G中,求G中的任意一对顶点间的最短路径问题,也是十分常见的一种问题。
解决这个问题的一个方法是执行n次迪杰斯特拉算法,这样就可以求出每一对顶点间的最短路径,执行的时间复杂度为O(n3)。
而另一种算法是由弗洛伊德提出的,时间复杂度同样是O(n3),但算法的形式简单很多。
可以将弗洛伊德算法描述如下:
在本题中,读入一个有向图的带权邻接矩阵(即数组表示),建立有向图并按照以上描述中的算法求出每一对顶点间的最短路径长度。
输入格式
输入的第一行包含1个正整数n,表示图中共有n个顶点。其中n不超过50。
以后的n行中每行有n个用空格隔开的整数。对于第i行的第j个整数,如果大于0,则表示第i个顶点有指向第j个顶点的有向边,且权值为对应的整数值;如果这个整数为0,则表示没有i指向j的有向边。当i和j相等的时候,保证对应的整数为0。
输出格式
共有n行,每行有n个整数,表示源点至每一个顶点的最短路径长度。如果不存在从源点至相应顶点的路径,输出-1。对于某个顶点到其本身的最短路径长度,输出0。
请在每个整数后输出一个空格,并请注意行尾输出换行。
样例输入
<span style="color:#333333"><span style="color:#333333">4
0 3 0 1
0 0 4 0
2 0 0 0
0 0 1 0
</span></span>
样例输出
<span style="color:#333333"><span style="color:#333333">0 3 2 1
6 0 4 7
2 5 0 3
3 6 1 0 </span></span>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#define MAXVEX 50
#define MAXEDGE 50//题目要求最大顶点数为50
#define INFINITY_ 65535
typedef struct{//数据结构
int vexs[MAXVEX];
int arc[MAXVEX][MAXVEX];
int numVertexes,numEdges;
}MGraph;
typedef int Patharc[MAXVEX][MAXVEX];//用于存储最短路径下标的数组,参考《大话数据结构-弗洛伊德算法》我认为这里值应该是为后继顶点的下标
typedef int ShortestPathTable[MAXVEX][MAXVEX];//用于存储到个点最短路径的权值和,值表示源点到其余顶点的最短路径长度和
void initial(MGraph*,int);//初始化
void ShortestPath_Floyd(MGraph,Patharc*,ShortestPathTable*);//算法
void print(ShortestPathTable,int);//打印
int main(){
MGraph g;
int n;
scanf("%d",&n);
initial(&g,n);
Patharc P;
ShortestPathTable D;
ShortestPath_Floyd(g,&P,&D);
print(D,n);
}
void initial(MGraph* g,int n){
g->numVertexes=n;
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
scanf("%d",&g->arc[i][j]);
if(g->arc[i][j]==0){//题目要求若两个顶点没有边,则用0表示,这里用65535表示方便算法执行
g->arc[i][j]=INFINITY_;
}
if(i==j){//顶点到自身权值依旧为0
g->arc[i][j]=0;
}
}
}
}
void ShortestPath_Floyd(MGraph g,Patharc* P,ShortestPathTable* D){//不得不说,弗洛伊德算法真的简洁优美
for(int i=0;i<g.numVertexes;i++){
for(int j=0;j<g.numVertexes;j++){
(*D)[i][j]=g.arc[i][j];
(*P)[i][j]=j;
}
}
for(int i=0;i<g.numVertexes;i++){
for(int j=0;j<g.numVertexes;j++){
for(int k=0;k<g.numVertexes;k++){
if((*D)[j][k]>(*D)[j][i]+(*D)[i][k]){
(*D)[j][k]=(*D)[j][i]+(*D)[i][k];
(*P)[j][k]=i;
}
}
}
}
}
void print(ShortestPathTable D,int n){
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
if(D[i][j]==65535)printf("%d ",-1);
else printf("%d ",D[i][j]);
}
printf("\n");
}
}