Dijkstra堆优化版
题目描述
小明是一位科学家,他需要参加一场重要的国际科学大会,以展示自己的最新研究成果。
小明的起点是第一个车站,终点是最后一个车站。然而,途中的各个车站之间的道路状况、交通拥堵程度以及可能的自然因素(如天气变化)等不同,这些因素都会影响每条路径的通行时间。
小明希望能选择一条花费时间最少的路线,以确保他能够尽快到达目的地。
题目链接:https://kamacoder.com/problempage.php?pid=1047
文章讲解:https://programmercarl.com/kamacoder/0047.参会dijkstra朴素.html
思考
堆优化 dijkstra。
朴素 Dijkstra 的时间复杂度分析
- 总体复杂度:
- 外层循环:遍历所有节点,确保每个节点都被访问 -
- 内层循环:在未访问节点中寻找距离源点最近的节点 -
关键问题:什么情况下朴素版本效率低下?
稀疏图(边数相对于节点数较少)
为什么稀疏图效率低?
在稀疏图中,内层 for 循环会遍历大量节点,但实际上很多节点之间没有直接连接,导致大量无效的遍历操作。
Dijkstra 三部曲的耗时操作
Dijkstra 算法的三个核心步骤:
第一步:选择距离源点最近且未访问的节点
当前做法:遍历 minDist 数组,查找最小值 -
问题:为什么每次都要重新遍历整个数组?
第二步:标记该节点为已访问
耗时:
第三步:更新相邻节点的距离
当前做法:遍历所有节点,更新 minDist 数组 -
问题:为什么要遍历所有节点,而不是只遍历实际相邻的节点?
优化
我们需要的是快速获取当前距离源点最近的未访问节点,而小顶堆可以快速获取最小值。
针对第一步的优化
原来:线性搜索 minDist 数组找最小值
现在:使用小顶堆,堆顶就是最小值。
时间复杂度:从
针对第三步的优化
原来:遍历所有节点更新距离
现在:只遍历当前节点的实际邻居节点
问题:如何高效存储和访问邻居节点?使用邻接表而不是邻接矩阵
数据结构
使用小顶堆:存储 (节点编号,源点到该节点的权值) 对。C++使用优先级队列实现小顶堆。
使用邻接表:存储图的结构,便于快速访问邻居。除了邻居编号,还要存储权值,使用 pair 不太方便区分,可以使用自定义的结构体。
注意! 构建邻接表是需要预先分配大小,否则后面添加会导致越界。执行 grah[p1].push_back(Edge(p2,val))
时,由于 grah
是空的vector
,访问 grah[p1]
会导致越界访问
算法流程
- 初始化:将源点加入堆
- 主循环:
- 从堆顶取出距离最小的节点
- 如果该节点已访问,跳过(处理重复入堆问题)
- 标记为已访问
- 遍历其邻居,更新距离并加入堆
复杂度分析
时间复杂度:
- V:节点数
- E:边数
- 堆操作:每条边最多导致一次入堆操作,每次操作
- 总操作数:最多 E 次入堆 + V 次出堆
空间复杂度:
- 邻接表:
- 堆:最多
个元素
代码实现
#include <iostream>
#include <vector>
#include <list>
#include <queue>
#include <climits>
using namespace std;
// 小顶堆
class mycomparison {
public:
bool operator()(const pair<int, int>& lhs, const pair<int, int>& rhs) {
return lhs.second > rhs.second;
}
};
// 定义一个结构体来表示带权重的边
struct Edge {
int to; // 邻接顶点
int val; // 边的权重
Edge(int t, int w): to(t), val(w) {} // 构造函数
};
int main() {
int n, m, p1, p2, val;
cin >> n >> m;
vector<list<Edge>> grid(n + 1);
for(int i = 0; i < m; i++){
cin >> p1 >> p2 >> val;
// p1 指向 p2,权值为 val
grid[p1].push_back(Edge(p2, val));
}
int start = 1; // 起点
int end = n; // 终点
// 存储从源点到每个节点的最短距离
std::vector<int> minDist(n + 1, INT_MAX);
// 记录顶点是否被访问过
std::vector<bool> visited(n + 1, false);
// 优先队列中存放 pair<节点,源点到该节点的权值>
priority_queue<pair<int, int>, vector<pair<int, int>>, mycomparison> pq;
// 初始化队列,源点到源点的距离为0,所以初始为0
pq.push(pair<int, int>(start, 0));
minDist[start] = 0; // 起始点到自身的距离为0
while (!pq.empty()) {
// 1. 第一步,选源点到哪个节点近且该节点未被访问过 (通过优先级队列来实现)
// <节点, 源点到该节点的距离>
pair<int, int> cur = pq.top(); pq.pop();
if (visited[cur.first]) continue;
// 2. 第二步,该最近节点被标记访问过
visited[cur.first] = true;
// 3. 第三步,更新非访问节点到源点的距离(即更新minDist数组)
for (Edge edge : grid[cur.first]) { // 遍历 cur指向的节点,cur指向的节点为 edge
// cur指向的节点edge.to,这条边的权值为 edge.val
if (!visited[edge.to] && minDist[cur.first] + edge.val < minDist[edge.to]) { // 更新minDist
minDist[edge.to] = minDist[cur.first] + edge.val;
pq.push(pair<int, int>(edge.to, minDist[edge.to]));
}
}
}
if (minDist[end] == INT_MAX) cout << -1 << endl; // 不能到达终点
else cout << minDist[end] << endl; // 到达终点最短路径
}