700字范文,内容丰富有趣,生活中的好帮手!
700字范文 > 单源最短路径算法的实现-----迪杰斯特拉(Dijkstra)算法【C/C++】

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

时间:2022-02-14 03:39:33

相关推荐

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

目录

前言一、迪杰斯特拉(Dijkstra)算法是什么?二、实现步骤1.算法思路2.进入主函数ShortestPath()1.创建final数组并且初始化path[]、dist[]数组2.对于节点的初始化3.进入主循环 三、全部代码(邻接表下)四、全部代码(邻接矩阵下)五、测试代码(邻接表下)六、简单的BFS算法的求解单源最短路径问题的解决方法七、简单算法与迪杰斯特拉算法对比测试总结

前言

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

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

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

一、迪杰斯特拉(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³)。这个和弗洛伊德算法的时间复杂度一样,但弗洛伊德算法那是相当的优雅。

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。