单源最短路径算法的实现-----迪杰斯特拉(Dijkstra)算法【C/C++】


前言

我们在生活中常常面临对路径选择的决策问题,这就要用到最短路径的算法了。

对于我这种榆木脑袋,显然迪杰斯特拉的这种算法有点高深。主要是我笨。

对于网图来说,最短路径,就是指两个顶点之间经过的边上权值之和最小的路径,并且我们称路径上的第一个顶点就是源点,最后一个顶点式终点。


一、迪杰斯特拉(Dijkstra)算法是什么?

迪杰斯特拉算法是一个按照路径长度递增的次序产生最短路径的算法。

二、实现步骤

1.算法思路

这里先采用邻接表来遍历。

在遍历节点时,找到未遍历节点中权值最小的进行遍历,并且及时更新最短路径长度dist数组[]。

首先设置path[]数组代表路径信息。 dist[] 表示最短路径长度。

int* path = (int*)malloc(sizeof(G.vexnum));
int* dist = (int*)malloc(sizeof(G.vexnum));

2.进入主函数ShortestPath()

1.创建final数组并且初始化path[]、dist[]数组

final数组来表示是否完成对该节点的最短路径求解。final[v]==1表示完成最短路径搜素,反之final[vi]==0表示未完成。

在算法中只有在求得最短路径后才会将final[vi]置为1,也可以简单理解为访问标志数组。

path数组全体初始化为0。

final数组因为最开始并没有完成最短路径求解,故置为0。

dist数组初始化为与vi相连的节点的权值,没连就是INFINITY(65535)。

int* final = (int*)malloc(sizeof(int) * g.vexnum);
	for (int i = 0; i < g.vexnum; i++) {
    
    
		path[i] = 0;
		final[i] = 0;
		dist[i] = INFNITY;
	}
	ArcNode* p = g.vertexlist[vi].firstarc;
	for (p; p != NULL; p = p->nextarc) {
    
    
		dist[p->adjvex] = p->weight;
	}

2.对于节点的初始化

在遍历vi节点时,vi到vi的路径为0,vi到vi之间也不需要求路径,故dist[vi]=0;final[vi]=1;

dist[vi] = 0;
final[vi] = 1;

肯定有人问,那path呢,path代表路径信息,vi时源点自然就是0了,当然初始化时也可以把path全初始化为-1,看个人习惯了。

3.进入主循环

将对刨掉源点的其他节点进行遍历,故外循环次数为g.vexnum-1次。

再在dist数组中找到权值最小并且未完成最短路径搜索的节点,用k来表示该节点下标。

其次找到最小权值k节点后,设置final[k]=1,再对k节点进行遍历,更新dist和path数组。

更新方法:若与k节点相连的节点未完成最短路径搜索并且k节点权值+该节点权值小于dist数组中的源点到该节点的最短路径,那么将更新dist数组中到该节点的最短路径,并且更新path数组,到该节点的前驱为k节点。

	int k;
	for (int v = 1; v < g.vexnum; v++) {
    
    
		int min = INFNITY;
		for (int w = 0; w < g.vexnum; w++) {
    
    
			if (!final[w] && dist[w] < min) {
    
    
				k = w;
				min = dist[w];
			}
		}
		final[k] = 1;
		ArcNode* p = g.vertexlist[k].firstarc;
		while (p != NULL) {
    
    
			if (!final[p->adjvex] && (p->weight + min) < dist[p->adjvex]) {
    
    
				dist[p->adjvex] = min + p->weight;
				path[p->adjvex] = k;
			}
			p = p->nextarc;
		}
	}

三、全部代码(邻接表下)

void ShortestPath(AdjList g, int vi, int* path, int* dist) {
    
    
	int* final = (int*)malloc(sizeof(int) * g.vexnum);
	for (int i = 0; i < g.vexnum; i++) {
    
    
		path[i] = 0;
		final[i] = 0;
		dist[i] = INFNITY;
	}
	ArcNode* p = g.vertexlist[vi].firstarc;
	for (p; p != NULL; p = p->nextarc) {
    
    
		dist[p->adjvex] = p->weight;
	}

	dist[vi] = 0;
	final[vi] = 1;

	int k;
	for (int v = 1; v < g.vexnum; v++) {
    
    
		int min = INFNITY;
		for (int w = 0; w < g.vexnum; w++) {
    
    
			if (!final[w] && dist[w] < min) {
    
    
				k = w;
				min = dist[w];
			}
		}
		final[k] = 1;
		ArcNode* p = g.vertexlist[k].firstarc;
		while (p != NULL) {
    
    
			if (!final[p->adjvex] && (p->weight + min) < dist[p->adjvex]) {
    
    
				dist[p->adjvex] = min + p->weight;
				path[p->adjvex] = k;
			}
			p = p->nextarc;
		}
	}
	free(final);
	return;
}

四、全部代码(邻接矩阵下)

思路大同小异,在初始化时有些不同,其他很相像。

void ShortestPath(AdjMatrix g, int vi, int* path, int* dist) {
    
    
	int* final = (int*)malloc(sizeof(int) * g.vexnum);
	for (int i = 0; i < g.vexnum; i++) {
    
    
		path[i] = 0;
		final[i] = 0;
		dist[i] = g.arc[vi][i];
	}
	dist[vi] = 0;
	final[vi] = 1;

	int k;
	for (int v = 1; v < g.vexnum; v++) {
    
    
		int min = INFNITY;
		for (int w = 0; w < g.vexnum; w++) {
    
    
			if (!final[w] && dist[w] < min) {
    
    
				k = w;
				min = dist[w];
			}
		}
		final[k] = 1;
		ArcNode* p = g.vertexlist[k].firstarc;
		for (int w = 0; w < g.vexnum; w++) {
    
    
			if (!final[w] && (min+g.arc[k][w])<dist[w]) {
    
    
				dist[w]=min+g.arc[k][w];
				path[w]=k;
			}
		}
	}
	free(final);
	return;
}

五、测试代码(邻接表下)

这里就测试一个邻接表下的。

自己花了个图

在这里插入图片描述

因为我的边表建立的时候A是第一个,自然A就是源点。

结果如下

在这里插入图片描述

很完美。

六、简单的BFS算法的求解单源最短路径问题的解决方法

适用于图里面权值为1或者没有权值,只是求步数的情况。
设立一个数组用于保存每一步的层号,利用BFS非递归,把每次的层号保存在数组里再输出即可。

void BFSMINDistance(AdjList g, int v) {
    
    
	int capacity = g.vexnum;
	int d[MAX] = {
    
     0 };
	for (int i = 0; i < g.vexnum; i++) {
    
    
		visited[i] = 0;
	}
	for (int i = 0; i < g.vexnum; i++) {
    
    
		d[i] = -1;
	}
	visited[v] = 1;
	d[v] = 0;
	queue<int>Q;
	Q.push(v);
	while (!Q.empty()) {
    
    
		int front = Q.front();
		Q.pop();
		ArcNode* p = g.vertexlist[front].firstarc;
		for (p; p != NULL; p = p->nextarc) {
    
    
			if (!visited[p->adjvex]) {
    
    
				visited[p->adjvex] = 1;
				d[p->adjvex] = d[front] + 1;
				Q.push(p->adjvex);
			}
		}
	}
	for (int i = 0; i < g.vexnum; i++) {
    
    
		if (i != v) {
    
    
			if (d[i] != -1) {
    
    
				printf("%c到%c的最短路径为%d\n", g.vertexlist[v].data, g.vertexlist[i].data, d[i]);
			}
		}
	}
}

七、简单算法与迪杰斯特拉算法对比测试

对于同一张图
在这里插入图片描述
在这里插入图片描述
测试代码

//迪杰斯特拉的测试代码
int* path = (int*)malloc(sizeof(G.vexnum));
int* dist = (int*)malloc(sizeof(G.vexnum));
printf("输入要查询的节点的最短路径\n");
char c;
scanf(" %c", &c);
printf("=======================迪杰斯特拉最短路径算法=====================\n");
int vi = LocateVertex(&G, c);
ShortestPath(G, 0, path, dist);
printf("================================================================\n");

//BFS法求最短路径的测试代码
printf("输入要查询的节点的最短路径\n");
char c1;
scanf(" %c", &c1);
printf("=======================BFS最短路径最短路径算法=====================\n");
int vi1 = LocateVertex(&G, c1);
BFS_MIN_Distance(G, vi1);
printf("================================================================\n");

首先是带权值图的测试
结果如下
在这里插入图片描述

在这里插入图片描述
显然迪杰斯特拉算法可以处理带权图
接下来是测试权值为1的图
在这里插入图片描述
在这里插入图片描述
显然两种算法结果一致!


总结

很显然这个算法的时间复杂度是O(n²),如果要知道任意顶点到其余所有顶点的最短路径,那么就可以对每一个顶点当作源点进行一次迪杰斯特拉算法。这时候后整个算法的时间复杂度也就成了O(n³)。这个和弗洛伊德算法的时间复杂度一样,但弗洛伊德算法那是相当的优雅。

猜你喜欢

转载自blog.csdn.net/qq_45400167/article/details/125982074