最短路径问题

最短路径问题

问题

如果从图中某一顶点（源点）到达另一顶点（终点）的路径可能不止一条，如何找到一条路径使得沿此路径上各边的权值总和（称为路径长度）达到最小。

实例

用Floyd算法求解下图各个顶点的最短距离。写出Floyd算法的伪代码和给出距离矩阵（顶点之间的最短距离矩阵），按实验报告模板编写算法。

Floyd（弗洛伊德）算法

算法思想：
从任意节点i到任意节点j的最短路径不外乎2种可能：1)直接从节点i到节点j，2)从节点i经过若干个节点k到节点j。所以，我们假设arcs(i,j)为节点i到节点j的最短路径的距离，对于每一个节点k，我们检查arcs(i,k) + arcs(k,j) < arcs(i,j)是否成立，如果成立，证明从节点i到节点k再到节点j的路径比节点i直接到节点j的路径短，我们便设置arcs(i,j) = arcs(i,k) + arcs(k,j)，这样一来，当我们遍历完所有节点k，arcs(i,j)中记录的便是节点i到节点j的最短路径的距离。

实现

void floyd() {
	FILE *fin;
	int graph[MAX_SIZE][MAX_SIZE] = { 0 };//邻接矩阵
	int vn, en;//点数、边数
	int startV, endV, w;//起始点、终点、权值
	int i, j, k;
	fin = fopen("input.txt", "r");
	fscanf(fin, "%d %d", &vn, &en);
	//初始化图
	for (i = 0; i < vn; i++) {
		for (j = 0; j < vn; j++) {
			if (i != j) {
				graph[i][j] = MAX_WEIGHT;
			}
			else {
				graph[i][j] = 0;
			}
		}
	}
	//读图
	for (i = 0; i < en; i++) {
		fscanf(fin, "%d %d %d", &startV, &endV, &w);
		graph[startV - 1][endV - 1] = w;
	}
	//Floyd算法
	for (k = 0; k < vn; k++) {
		//k为过度点
		for (i = 0; i < vn; i++) {
			for (j = 0; j < vn; j++) {
				//当通过过渡点的路径小于当前最短路径时，更新最短路径
				if (graph[i][k] + graph[k][j] < graph[i][j]) {
					graph[i][j] = graph[i][k] + graph[k][j];
				}
			}
		}
		//第k次更新结束

		printf("并入点%d后\n", k + 1);

		for (i = 0; i < vn; i++) {
			for (j = 0; j < vn; j++) {
				printf("%5d", graph[i][j]);
			}
			printf("\n");
		}

	}
	printf("注：%d表示不可达\n",MAX_WEIGHT);
}

Dijkstra（迪克斯特拉）算法

算法思想：
设G=(V,E)是一个带权有向图，把图中顶点集合V分成两组，第一组为已求出最短路径的顶点集合（用S表示，初始时S中只有一个源点，以后每求得一条最短路径 , 就将加入到集合S中，直到全部顶点都加入到S中，算法就结束了），第二组为其余未确定最短路径的顶点集合（用U表示），按最短路径长度的递增次序依次把第二组的顶点加入S中。在加入的过程中，总保持从源点v到S中各顶点的最短路径长度不大于从源点v到U中任何顶点的最短路径长度。此外，每个顶点对应一个距离，S中的顶点的距离就是从v到此顶点的最短路径长度，U中的顶点的距离，是从v到此顶点只包括S中的顶点为中间顶点的当前最短路径长度。

实现

void dijkstra() {
	FILE *fin;
	int graph[MAX_SIZE][MAX_SIZE] = { 0 };//邻接矩阵
	int vn, en;//点数、边数
	int startV, endV, w;//起始点、终点、权值
	int i, j, k;
	int d[10];//保存最短路径长度
	int final[10];//若final[i] = 1则说明 顶点vi已在集合S中
	int origin;//起始点
	fin = fopen("input2.txt", "r");
	fscanf(fin, "%d %d", &vn, &en);
	//初始化图
	for (i = 0; i < vn; i++) {
		for (j = 0; j < vn; j++) {
			if (i != j) {
				graph[i][j] = MAX_WEIGHT;
			}
			else {
				graph[i][j] = 0;
			}
		}
	}
	//读图
	for (i = 0; i < en; i++) {
		fscanf(fin, "%d %d %d", &startV, &endV, &w);
		graph[startV - 1][endV - 1] = w;
	}
	printf("请输入起始点:");
	scanf("%d",&origin);
	origin = origin - 1;
	//Dijkstra算法
	//初始化建立可达关系表
	for (i = 0; i < vn; i++) {
		final[i] = 0;
		d[i] = graph[origin][i];
	}
	//初始化 起始点属于集合S
	d[origin] = 0;
	final[origin] = 1;
	//开始主循环 每次求得起始点到某个顶点v的最短路径 并加v到集合S中
	for (i = 1; i < vn; i++) {
		int min = MAX_WEIGHT;
		//选择当前出离起始点最近且不在集合S中的顶点v
		for (j = 0; j < vn; j++) {
			if (!final[j]) {
				if (d[j]<min) {
					k = j;
					min = d[j];
				}
			}
		}
		final[k] = 1;//选出该点后加入到合集S中
		//更新当前最短路径和距离
		for (j = 0; j < vn; j++) {
			if (!final[j] && (min + graph[k][j] < d[j])&&graph[k][j]!=0) {
				d[j] = min + graph[k][j];
			}
		}
	}
	for (i = 0; i < vn; i++) {
		printf("从点%d 到点%d 的最短路径为 %d\n", origin+1,i+1, d[i]);
	}
	printf("注：%d表示不可达\n", MAX_WEIGHT);
}

github源码地址