【图论】迪杰特斯拉算法

【图论】迪杰特斯拉算法

文章目录

在这里插入图片描述

迪杰特斯拉算法

迪杰特斯拉算法是由荷兰计算机科学家艾兹赫尔·迪杰特斯拉(Edsger W. Dijkstra)在1956年提出的,用于解决单源最短路径问题的经典算法。该算法的目标是从一个起始顶点找到到图中其他顶点的最短路径。

主要特点

  • 适用于带权图,其中权重为非负数。(为什么只适用于非负数,因为迪杰斯特拉的思想是贪心测量,当有负权引入的时候,贪心策略将不再适用)
  • 解决从单个源点到其他所有顶点的最短路径问题。
  • 时间复杂度:当使用优先队列(例如堆)时,复杂度为 O ( E log ⁡ V ) O(E \log V) O(ElogV),其中 V V V 为顶点数量, E E E 为边的数量。

基本思想

Dijkstra算法通过不断探索距离最近的顶点,逐步扩展其最短路径的已知范围,直到找到从源点到所有其他顶点的最短路径。该算法基于贪心策略:每一步选择尚未处理的、距离源点最近的顶点进行扩展。

算法步骤

  1. 初始化
    • 将起始顶点的距离设为0,其余所有顶点的距离设为∞(表示不可达)。
    • 使用一个优先队列(或最小堆)来存储顶点及其当前的最短距离。
  2. 取距离源点最近的顶点,并标记为已处理。
  3. 对于该顶点的每个邻接顶点,更新其最短距离(如果通过当前顶点可以获得更短的路径,则更新距离)。
  4. 重复步骤2和3,直到所有顶点都被处理过,或优先队列为空。

示例

在这里插入图片描述

实现迪杰斯特拉算法

基本步骤

首先假设我们有如下的一个图:

在这里插入图片描述


灰色的点代表起点,我们需要从起点开始算每个点到起点的最短路径。
第一步按照迪杰斯特拉的表述应该将起点到起点的最短路径初始化为0,起点到另外其他点的距离初始化为正无穷。

在这里插入图片描述


这里我们更新完s点之后,应该更新和s点相连的点的最短路径了,由于之前s到t点和到y点的最短路径是正无穷,由于st和sy都小于两个最短路径,所以更新两个最短路径。

在这里插入图片描述


根据贪心策略,更新出来的两个最短路径,sy更小,所以我们应该更新y相连的路径。

在这里插入图片描述


所以接下来我们应该更新y连接的点的最短路径。

在这里插入图片描述


由于原本s到t的最短路径是10,但是现在的最短路径更新了之后是8,所以更新结果,由于之前s到x的最短路径是正无穷,所以现在将最短路径更新为14,s到z 也是相同的,因为之前的最短路径是正无穷,所以现在将最短路径更新为7.
在这三条最短路径中选择最短的那条:

在这里插入图片描述


这里就应该以z为新的起点

在这里插入图片描述


更新z连接的顶点,z一共有两条边,一条是zs,一条是zx,由于s到s是最近的0,所以这里不需要更新,由于之前s到x的距离是14,所以这里更新s到x的距离。
接下来再从剩下的边中,选择最小的路径。

在这里插入图片描述


从t点开始只有一条路径,所以这里t到x,tx是1,由于之前的st的最短路径是8,所以此时的最短路径是9,9<13所以这里应该更新最短路径为9

在这里插入图片描述


这里最短路径算是更新完了。

算法思路

由于这里我们涉及到最短路径,所以应该开辟一个dist数组,我们来想一想一个dist数组是否能解决问题,很显然是不能的,我们还需要一个数组,记录当前最短路径的前一个顶点的下标,在遍历的时候我们可以索引每一个顶点了。
代码展示:

voidDijkstra(const V& src, vector<W>& dist, vector<int>& pPath){//获取起点的下标size_t srci =GetVertexIndex(src);//确定节点的数量size_t n = _vertexs.size(); dist.resize(n, MAX_W); pPath.resize(n, MAX_W);//自己到自己的距离是0 dist[srci]=0;//自己的前一个是自己 pPath[srci]= srci;//已经确定最短路径的顶点的集合 vector<bool>S(n, false);for(size_t j=0;j<n;j++){//选择最短路径顶点且不在s中更新其他路径int u =0;//最小的那个点 W min = MAX_W;//最小权值for(size_t i =0;i < n;i++){if(S[i]== false && dist[i]< min){//选出最小的点 u = i;//更新最小权值 min = dist[i];}}//u被选出来 S[u]= true;//松弛更新u连接出去的顶点v srci->u u->vfor(size_t v =0;v < n;v++){//确定u连接出去的所有边if(S[v]== false && _matrix[u][v]!= MAX_W &&(dist[u]+ _matrix[u][v])< dist[v]){ dist[v]= dist[u]+ _matrix[u][v];//将v的父亲更新为u pPath[v]= u;}}}}

打印路径节点和最短路径:

//打印最短路径voidPrinrtShotPath(const V& src, vector<W>& dist, vector<int>& pPath){//获取起点的下标size_t srci =GetVertexIndex(src);//确定节点的数量size_t n = _vertexs.size();for(size_t i =0;i < n;i++){if(i != srci){// 先找出i顶点的路径 vector<int> path;size_t parenti = i;//先将自己存进去(自己不是原点先push进去)while(parenti != srci){ path.push_back(parenti); parenti = pPath[parenti];} path.push_back(srci);//将路径逆置reverse(path.begin(), path.end());//打印出路径for(auto e : path){ cout << _vertexs[e]<<"->";}//打印最短路径 cout << dist[i]; cout << endl;}}}

因为根据最后一个节点去推上一个节点,推完之后数组中的节点会是一个反着的路径,所以在打印的时候,应该先把数组逆置,逆置之后再进行打印。

总结

在本文中,我们深入探讨了迪杰斯特拉算法的原理与应用。作为一种经典的最短路径算法,迪杰斯特拉算法通过优先队列有效地解决了从单一源点到其他所有节点的最短路径问题。我们分析了其时间复杂度和空间复杂度,了解了在不同图形结构下的性能表现。

通过示例和实现,我们不仅掌握了算法的基本步骤,还体验了其在实际应用中的重要性。无论是在交通导航、网络路由还是各种优化问题中,迪杰斯特拉算法都发挥着不可或缺的作用。

希望本文能够帮助你更好地理解迪杰斯特拉算法,并为你在图论和算法领域的进一步学习打下坚实的基础。如果你有任何疑问或想法,欢迎在评论区与我交流!

Read more

第十六届蓝桥杯省赛(软件类真题)C/C++ 大学A组

第十六届蓝桥杯省赛(软件类真题)C/C++ 大学A组

大纲: A.寻找质数 B:黑白棋 题目&解析&代码 A题 题目解析 本题的目标是枚举质数并计数,直到数到第2025个。由于2025不算太大,第2025个质数大约在17000~18000之间,完全可以在合理时间内通过简单枚举得到。 解题步骤: 从2开始遍历每个整数,判断它是否是质数。 质数判断采用试除法:对于一个数n,只需检查从2到√n的所有整数是否能整除n。若存在能整除的数,则n不是质数;否则是质数。 每找到一个质数,计数器加1。 当计数器达到2025时,输出当前的质数并结束。 优化点: 除了2以外,偶数不可能是质数,因此可以跳过偶数判断(直接步进2)。 在isPrime函数中,可以先处理特殊情况(n<2返回false),然后单独判断偶数,再对奇数进行试除,步进也可以设为2。 C++ 参考代码 以下代码实现了上述算法,并输出第2025个质数。 cpp

By Ne0inhk
初学二叉搜索树踩坑多?C++ 从原理到代码,搞定增删查全流程

初学二叉搜索树踩坑多?C++ 从原理到代码,搞定增删查全流程

🎬 个人主页:Vect个人主页 🎬 GitHub:Vect的代码仓库 🔥 个人专栏: 《数据结构与算法》《C++学习之旅》《计算机基础》 ⛺️Per aspera ad astra. 文章目录 * 1. 二叉搜索树相关概念 * 2. 二叉搜索树的操作 * 2.1. 查找节点 * 2.2. 插入节点 * 2.3. 删除节点 * 3. 二叉搜索树的实现 * 4. 二叉搜索树的应用 * 4.1. K模型 * 4.2. KV模型 1. 二叉搜索树相关概念 如下图所示,二叉搜索树(binary search tree)满足下列条件: 1. 对于根节点,左子树中所有节点的值<根节点的值&

By Ne0inhk

第25章-C++初级实战案例(20个)

案例1:温度转换器 案例描述 实现摄氏度与华氏度之间的相互转换。 知识点 * 基本输入输出 * 数学运算 * 函数封装 完整代码 #include<iostream>#include<iomanip>usingnamespace std;// 摄氏度转华氏度doublecelsiusToFahrenheit(double celsius){return celsius *9.0/5.0+32.0;}// 华氏度转摄氏度doublefahrenheitToCelsius(double fahrenheit){return(fahrenheit -32.0)*5.0/9.0;}intmain(){int choice;double temp, result; cout <<

By Ne0inhk