【Floyd最短路径算法】

【Floyd最短路径算法】

求解每对顶点间的最短路径问题
允许图中有带负权值的边，但不允许有包含带负权值的边组成的回路。
时间复杂度O(|V|³),空间复杂度O(|V|²)。
主要步骤为计算递推方阵(开始为邻接矩阵)
A ( k ) [ i ] [ j ] = M i n { A ( k − 1 ) [ i ] [ j ] , A ( k − 1 ) [ i ] [ k ] + A ( k − 1 ) [ k ] [ j ] } ,     k = 0 , 1 , … , n − 1 A^{(k)}[i][j]=Min\{ A^{(k-1)}[i][j],A^{(k-1)}[i][k]+A^{(k-1)}[k][j] \},\ \ \ k=0,1,\dots,n-1 A(k)[i][j]=Min{A(k−1)[i][j],A(k−1)[i][k]+A(k−1)[k][j]},   k=0,1,…,n−1
其中A^(k)[ i ] [ j ]是从v_i到v_j、中间顶点序号不大于k的最短路径的长度

void Floyd(MGraph G) {	// Floyd算法 求解每对顶点间的最短路径问题
	EdgeType arcs[MaxVertexNum][MaxVertexNum];
	int i, j, k;
	// 初始为邻接矩阵
	for (i = 0; i < G.vexnum; i++){
		for (j = 0; j < G.vexnum; j++) {
			arcs[i][j] = G.Edge[i][j];
			printf("%d ", arcs[i][j]);
		}
		printf("\n");
	}
	// 迭代
	for (k = 0; k < G.vexnum; k++) { // 迭代n次
		for(i=0;i<G.vexnum;i++)	// 对矩阵每个元素
			for (j = 0; j < G.vexnum; j++) {
				if (arcs[i][j] > arcs[i][k] + arcs[k][j] && i!=j)  // 对角线不更新
					arcs[i][j] = arcs[i][k] + arcs[k][j];
			}
	}
	// 打印迭代结果
	for (i = 0; i < G.vexnum; i++) {
		for (j = 0; j < G.vexnum; j++) {
			printf("%d ", arcs[i][j]);
		}
		printf("\n");
	}
}