文章目录

  • 1 Dijkstra算法基本原理
  • 2 算法过程图解1(有向图)
  • 3 算法过程图解2(无向图)
  • 4 C++代码
    • 4.1 案例1代码
    • 4.2 案例2邻接矩阵定义
    • 4.3 案例2代码Dijkstra算法

1 Dijkstra算法基本原理

Dijkstra算法是根据贪心算法实现的,首先找出当前点到所有能到达的点之间最短的距离,然后松弛一次继续循环。所谓松弛一次,就是在已经访问过的点中遍历一遍,看看有没有更近的,如果有更近的就更新距离。这样每次找最近的可达点+松弛遍历历史节点的操作,一直重复就能找到最短路径。

算法介绍

迪杰斯特拉(Dijkstra)算法是典型最短路径算法,用于计算一个节点到其他节点的最短路径。
它的主要特点是以起始点为中心向外层层扩展(广度优先搜索思想),直到扩展到终点为止。

基本步骤

通过Dijkstra计算图G中的最短路径时,需要指定起点s(即从顶点s开始计算)。

此外,引进两个集合S和U。S的作用是记录已求出最短路径的顶点(以及相应的最短路径长度),而U则是记录还未求出最短路径的顶点(以及该顶点到起点s的距离)。

初始时,S中只有起点s;U中是除s之外的顶点,并且U中顶点的路径是"起点s到该顶点的路径"。

然后,从U中找出路径最短的顶点,并将其加入到S中;

接着,更新U中的顶点和顶点对应的路径,再从U中找出路径最短的顶点,并将其加入到S中;

接着,更新U中的顶点和顶点对应的路径。

重复该操作,直到遍历完所有顶点。

2 算法过程图解1(有向图)

首先看一个图,目标是找到从1到6的最短距离

首先用一个6*6的二维数组存储距离信息,纵坐标为起点,横坐标为终点,对应的坐标值就是起点到终点的距离,其中,不可达点的距离为∞\infty∞,到自身的距离为0。

接着用一个一维数组来存储路径

从点1开始,遍历二维数组的第一行,发现距离1是最小的,对应点1到点2,选定点2;松弛,没有更近值;

从点2开始,遍历二维数组的第二行,发现距离3是最小的,对应点2到点4,选定点4;松弛,没有更近值;

从点4开始,遍历二维数组的第四行,发现距离4是最小的,对应点4到点3,选定点3;松弛,没有更近值;

从点3开始,遍历二维数组的第三行,发现距离5是最小的,对应点3到点5,选定点5;松弛,没有更近值;

从点5开始,遍历二维数组的第五行,发现距离4是最小的,对应点5到点6,选定点6;松弛,没有更近值;

点6,已经到达目标节点,遍历结束。

3 算法过程图解2(无向图)

注意:B(23)应该为B(13)

初始状态:S是已计算出最短路径的顶点集合,U是未计算除最短路径的顶点的集合!
第1步:将顶点D加入到S中。
此时,S={D(0)}, U={A(∞),B(∞),C(3),E(4),F(∞),G(∞)}。 注:C(3)表示C到起点D的距离是3。

第2步:将顶点C加入到S中。
上一步操作之后,U中顶点C到起点D的距离最短;因此,将C加入到S中,同时更新U中顶点的距离。以顶点F为例,之前F到D的距离为∞;但是将C加入到S之后,F到D的距离为9=(F,C)+(C,D)。
此时,S={D(0),C(3)}, U={A(∞),B(13),E(4),F(9),G(∞)}。

第3步:将顶点E加入到S中。
上一步操作之后,U中顶点E到起点D的距离最短;因此,将E加入到S中,同时更新U中顶点的距离。还是以顶点F为例,之前F到D的距离为9;但是将E加入到S之后,F到D的距离为6=(F,E)+(E,D)。
此时,S={D(0),C(3),E(4)}, U={A(∞),B(13),F(6),G(12)}。

第4步:将顶点F加入到S中。
此时,S={D(0),C(3),E(4),F(6)}, U={A(22),B(13),G(12)}。

第5步:将顶点G加入到S中。
此时,S={D(0),C(3),E(4),F(6),G(12)}, U={A(22),B(13)}。

第6步:将顶点B加入到S中。
此时,S={D(0),C(3),E(4),F(6),G(12),B(13)}, U={A(22)}。

第7步:将顶点A加入到S中。
此时,S={D(0),C(3),E(4),F(6),G(12),B(13),A(22)}。

此时,起点D到各个顶点的最短距离就计算出来了:A(22) B(13) C(3) D(0) E(4) F(6) G(12)

4 C++代码

4.1 案例1代码

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#define Inf 0x3f3f3f3fusing namespace std;int map[1005][1005];int vis[1005],dis[1005];
int n,m;//n个点,m条边void Init ()
{memset(map,Inf,sizeof(map));for(int i=1;i<=n;i++){map[i][i]=0;}
}void Getmap()
{int u,v,w;for(int t=1;t<=m;t++){scanf("%d%d%d",&u,&v,&w);if(map[u][v]>w){map[u][v]=w;map[v][u]=w;}}  }void Dijkstra(int u)
{memset(vis,0,sizeof(vis));for(int t=1;t<=n;t++){dis[t]=map[u][t];}vis[u]=1;for(int t=1;t<n;t++){int minn=Inf,temp;for(int i=1;i<=n;i++){if(!vis[i]&&dis[i]<minn){minn=dis[i];temp=i;}}vis[temp]=1;for(int i=1;i<=n;i++){if(map[temp][i]+dis[temp]<dis[i]){dis[i]=map[temp][i]+dis[temp];}}}}int main()
{scanf("%d%d",&m,&n);Init();Getmap();Dijkstra(n);printf("%d\n",dis[1]);return 0;
}

4.2 案例2邻接矩阵定义

// 邻接矩阵
typedef struct _graph
{char vexs[MAX];       // 顶点集合int vexnum;           // 顶点数int edgnum;           // 边数int matrix[MAX][MAX]; // 邻接矩阵
}Graph, *PGraph;// 边的结构体
typedef struct _EdgeData
{char start; // 边的起点char end;   // 边的终点int weight; // 边的权重
}EData;

Graph是邻接矩阵对应的结构体。

vexs用于保存顶点,vexnum是顶点数,edgnum是边数;matrix则是用于保存矩阵信息的二维数组。例如,matrix[i][j]=1,则表示"顶点i(即vexs[i])"和"顶点j(即vexs[j])"是邻接点;matrix[i][j]=0,则表示它们不是邻接点。

EData是邻接矩阵边对应的结构体。

4.3 案例2代码Dijkstra算法

/** Dijkstra最短路径。* 即,统计图(G)中"顶点vs"到其它各个顶点的最短路径。** 参数说明:*        G -- 图*       vs -- 起始顶点(start vertex)。即计算"顶点vs"到其它顶点的最短路径。*     prev -- 前驱顶点数组。即,prev[i]的值是"顶点vs"到"顶点i"的最短路径所经历的全部顶点中,位于"顶点i"之前的那个顶点。*     dist -- 长度数组。即,dist[i]是"顶点vs"到"顶点i"的最短路径的长度。*/
void dijkstra(Graph G, int vs, int prev[], int dist[])
{int i,j,k;int min;int tmp;int flag[MAX];      // flag[i]=1表示"顶点vs"到"顶点i"的最短路径已成功获取。// 初始化for (i = 0; i < G.vexnum; i++){flag[i] = 0;              // 顶点i的最短路径还没获取到。prev[i] = 0;              // 顶点i的前驱顶点为0。dist[i] = G.matrix[vs][i];// 顶点i的最短路径为"顶点vs"到"顶点i"的权。}// 对"顶点vs"自身进行初始化flag[vs] = 1;dist[vs] = 0;// 遍历G.vexnum-1次;每次找出一个顶点的最短路径。for (i = 1; i < G.vexnum; i++){// 寻找当前最小的路径;// 即,在未获取最短路径的顶点中,找到离vs最近的顶点(k)。min = INF;for (j = 0; j < G.vexnum; j++){if (flag[j]==0 && dist[j]<min){min = dist[j];k = j;}}// 标记"顶点k"为已经获取到最短路径flag[k] = 1;// 修正当前最短路径和前驱顶点// 即,当已经"顶点k的最短路径"之后,更新"未获取最短路径的顶点的最短路径和前驱顶点"。for (j = 0; j < G.vexnum; j++){tmp = (G.matrix[k][j]==INF ? INF : (min + G.matrix[k][j])); // 防止溢出if (flag[j] == 0 && (tmp  < dist[j]) ){dist[j] = tmp;prev[j] = k;}}}// 打印dijkstra最短路径的结果printf("dijkstra(%c): \n", G.vexs[vs]);for (i = 0; i < G.vexnum; i++)printf("  shortest(%c, %c)=%d\n", G.vexs[vs], G.vexs[i], dist[i]);
}

参考文献

Dijkstra算法(一)之C语言详解
https://www.cnblogs.com/skywang12345/p/3711512.html

Dijkstra算法图文详解
https://blog.csdn.net/lbperfect123/article/details/84281300

动画理解Dijkstra算法过程
https://blog.csdn.net/jzj1993/article/details/39157055

Dijkstra算法图文详解和C++代码相关推荐

  1. Dijkstra算法图文详解

    Dijkstra算法算是贪心思想实现的,首先把起点到所有点的距离存下来找个最短的,然后松弛一次再找出最短的,所谓的松弛操作就是,遍历一遍看通过刚刚找到的距离最短的点作为中转站会不会更近,如果更近了就更 ...

  2. ADMM,ISTA,FISTA算法步骤详解,MATLAB代码,求解LASSO优化问题

    ADMM,ISTA,FISTA算法步骤详解,MATLAB代码,求解LASSO优化问题 原创文章!转载需注明来源:©️ Sylvan Ding's Blog ❤️ 实验目的 了解 ADMM, ISTA, ...

  3. Dijkstra 路径规划算法原理详解及 Python 代码实现

    荷兰数学家 E.W.Dijkstra 于 1959 年提出了 Dijkstra 算法,它是一种适用于 非负权值 网络的 单源最短路径算法,同时也是目前求解最短路径问题的理论上最完备.应用最广的经典算法 ...

  4. KMP算法图文详解(为什么是next[0]=-1、next[j]=k和k=next[k])

    文章目录 一:KMP算法解决的问题 二:详解KMP (1)暴力匹配的缺点 (2)最长相同前缀和后缀 (3)究竟怎么回溯 (3)next数组 (4)求解next数组 A:next[0]=-1 B:nex ...

  5. 五一带你找到最省旅游路线【dijkstra算法推导详解】

    前言 五一快到了,小张准备去旅游了! 查了查到各地的机票 因为今年被扣工资扣得很惨,小张手头不是很宽裕,必须精打细算.他想弄清去各个城市的最低开销. [嗯,不用考虑回来的开销.小张准备找警察叔叔说自己 ...

  6. KNN算法原理详解及python代码实现

    KNN算法 算法原理 对数据的要求 算法的优缺点 算法需要注意的点 算法实现(python) 算法原理 计算待测样本与train_data的距离d并保存数组中 对d进行排序,取d最近的k个样本 统计样 ...

  7. pythonsort函数时间复杂度_合并排序算法——时间复杂度详解和python代码实现

    递归形式 递归形式是算法中常用到的一种构造思路.递归允许函数或过程对自身进行调用,是对算法中重复过程的高度概括,从本质层面进行刻画,避免算法书写中过多的嵌套循环和分支语法.因此,是对算法结构很大的简化 ...

  8. 蚁群算法原理详解和matlab代码

    1原理: 蚂蚁在寻找食物源的时候,能在其走过的路径上释放一种叫信息素的激素,使一定范围内的其他蚂蚁能够察觉到.当一些路径上通过的蚂蚁越来越多时,信息素也就越来越多,蚂蚁们选择这条路径的概率也就越高,结 ...

  9. 归并求逆序对算法图文详解

    给定一个长度为 n 的整数数列,请你计算数列中的逆序对的数量. 逆序对的定义如下:对于数列的第 i 个和第 j 个元素,如果满足 i<j 且 a[i]>a[j],则其为一个逆序对:否则不是 ...

最新文章

  1. 翻译BonoboService官网的安装教程
  2. 三星全球首家量产18nm DRAM
  3. P of EAA 总结
  4. php和python哪个学起来简单一点-现在自学php和python那个合适?
  5. 深度学习开发者的AI Party来了!WAVE SUMMIT+2021与你相约双十二
  6. Fedora 安装后需要做的第一件事
  7. 嵌入式IAP开发笔记之一:面向STM32的BootLoader程序
  8. tp5 queue.php,tp5(think-queue)消息队列+supervisor进程管理实现队列常驻进程
  9. 手机MODEM 开发(26)---LTE切换和重选参数中文详解
  10. Lucene的索引不跨平台
  11. Android 一个漂亮的Android日期和时间选择器:DateTimePicker
  12. Servlet获取全路径
  13. python爬取高德poi数据_python3爬虫-高德地图POI数据的爬取
  14. telnet 测试IP和端口命令
  15. 数值分析第二次作业-求解系数矩阵为Hilbert 矩阵的线性方程组
  16. 字节跳动开启员工暴富计划,大手笔回购期权,有人最高能赚175%
  17. 金融经济学二十五讲--徐高
  18. c语言用while语句计算圆周率的近似值,编程计算圆周率的近似值
  19. 小盒子可以在大盒子里面移动
  20. EPICS IOC Shell

热门文章

  1. vc 2010 下打包dll 的问题
  2. Java中遍历文件夹的2种方法
  3. 1759: 多少个元素比它大 二分
  4. MyEclipse中Hibernate+C3P0连接池配置
  5. 随想录(thread类的编写)
  6. delete响应服务器,rest-RESTful-DELETE响应主体应包含什么
  7. 10天智能锁项目实战第1天(了解单片机STM32F401RET6和C语言基础)
  8. 条件测试 linux,Linux的条件测试
  9. jsp简介mysql_jsp+servlet+mysql简单实现用户登陆注册
  10. Windows系统 notepad命令详解,Windows系统打开记事本