程序设计基础26 弗洛伊德算法

题目描述

在带权有向图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。

请在每个整数后输出一个空格，并请注意行尾输出换行。

样例输入

4
0 3 0 1
0 0 4 0
2 0 0 0
0 0 1 0

样例输出

0 3 2 1 
6 0 4 7 
2 5 0 3 
3 6 1 0 

一，注意事项

    初始化是把不能达到的距离初始化成INF，因为Floyed()子函数中有一项是dis[i][k]+dis[k][j]<dis[i][j]，dis[i][j]是无穷大的话只要有中介点使其连通便可轻易替换。

    还有就是把横纵坐标相等的点初始化为0,

    关于本题还有一项，就是输出的时候若是两点不可达要用-1来替换！当初就是在这WA的。

二，我的代码

#include<cstdio>
#include<algorithm>
using namespace std;
const int INF = 1000000000;
const int max_n = 60;
int N = 0;
int dis[max_n][max_n];
using namespace std;
void Floyed() {
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= N; j++) {
			if (dis[i][j] == 0) {
				dis[i][j] = INF;
			}
		}
	}
	for (int i = 1; i <= N; i++) {
		dis[i][i] = 0;
	}
	for (int k = 1; k <= N; k++) {
		for (int i = 1; i <= N; i++) {
			for (int j = 1; j <= N; j++) {
				if (dis[i][k] != INF && dis[k][j] != INF && dis[i][k] + dis[k][j] < dis[i][j]) {
					dis[i][j] = dis[i][k] + dis[k][j];
				}
			}
		}
	}
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= N; j++) {
			if (dis[i][j] == INF) {
				dis[i][j] = -1;
			}
		}
	}
}
int main() {
	scanf("%d", &N);
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= N; j++) {
			scanf("%d", &dis[i][j]);
		}
	}
	Floyed();
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= N; j++) {
			printf("%d", dis[i][j]);
			if (j != N) {
				printf(" ");
			}
		}
		printf("\n");
	}
	return 0;
}