Dijkstra求最短路与次短路

花了一个晚上加上午两节课的时间来思考这个问题,一开始进了个误区,后来发现不对,找到了正确思路,现在来做个总结

最短路:

假设有如下 无向图:


每条边有权,要求从A到G的最短路,设数组d[i]用来记录每个点到A的最短路,D[][]用来记录权值d[0]=0。先来说一下我一开始的误区,在看了一些介绍这个算法的文章后我大概知道了是要一步一步来更新每个点的最小值。于是我是这样想的从A到B遍历,A到A的距离是0,然后是B,d[1]=D[0][1]=2,然后是C,d[2]=min(D[0][2], d[1]+D[1][2])=5,接下来是D,d[3]=min(d[1]+D[1][3], d[2]+D[2][3]),我用图来解释这个过程:







最后结果是的d[6]=17,但实际上问题在E和F上已经出现了,由于我们先算的是E,在这时数组里F的距离还是无穷大,所以E只能从B得到一个2,再加上BE=10就是12。但算过F后发现,E的最短距离应该是F的8加上EF=3,结果是11,于是G的最短距离应该是11+5=16。问题就在于我们按ABCDEFG的顺序依次简单地求当前点的最短距离时,每个点的距离只能更新一次,无法从后续计算中得到该点更近的距离。这种方法中”当前点“是未求值的点,而正确的做法应该把”当前点“定为已求过的点。也就是说,我们要做的是从该点出发求临近点的最小值,而不是站在该点,用已求过值的临近点来求该点的值。


正确的过程应该是这样,计算每个点时都要用这个转移方程 d[i]=min(d[i], d[j]+D[i][j])其中j是当前点,i是与j相连的点:

从A开始,A与B,C相连,更新B,C,d[1]=2, d[2]=5;

转到B,B与CDE相连,d[2]=min(d[2], d[1]+D[2][1])=5,同理,d[3]=7, d[4]=12

…………

到F时,F与DEG相连,所以在更新E的时候就会有d[4]=min(d[4], d[5]+D[5][4])=11;

这样做就可以用在E之后的F来更新E的值。


这一过程可以直接将d[i]放到循环里实现,也可以用优先队列,将(d[i], i)入列。当从某个点j求得的d[i]大于数组中已有的d[i]时,就不入列。这个优先队列的过程在纸上画一画就很清楚了。


忘了说了,在转移到下一个点时必须转移到d[i]最小的点,这个条件非常重要!


次短路:


接下来是次短路。次短路的做法和最短路几乎一模一样。在过程中同时记录每个点的最短路和次短路,然后用同样的方法更新每个点的值。要注意的就是计算与当前点相邻点的次短路时需要同时用到该点的最短路和次短路,也就是要求三个值的最小值:min(d2[i], d[j]+D[i][j], d2[j]+D[i][j]),其中数组d2[i]就是每点的次短路。要注意的是要先求最短路再求次短路,保证最短路小于次短路。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Dijkstra算法单源最短路径的经典算法,其基本思想是通过逐步扩展生成最短路径集合,最终得到源点到所有其它点的最短路径。 以下是C++实现: ```c++ #include <iostream> #include <vector> #include <queue> #include <cstring> using namespace std; const int INF = 0x3f3f3f3f; // 定义正无穷 struct Edge { int to, w; Edge(int to, int w) : to(to), w(w) {} }; vector<Edge> G[100010]; // 邻接表存图 int dist[100010]; // 存储最短路径长度 bool vis[100010]; // 标记是否已经确定最短路径 void dijkstra(int s) { memset(dist, INF, sizeof(dist)); // 初始化距离为正无穷 memset(vis, false, sizeof(vis)); // 初始化标记为未确定最短路径 dist[s] = 0; // 源点到自己的距离为0 priority_queue<pair<int, int>, vector<pair<int, int>>, greater<pair<int, int>>> q; // 小根堆 q.push(make_pair(0, s)); // 将源点入队 while(!q.empty()) { int u = q.top().second; // 取出当前距离最小的点 q.pop(); if(vis[u]) continue; // 如果已经确定最短路径,直接跳过 vis[u] = true; // 标记为已确定最短路径 for(auto e : G[u]) { // 遍历所有相邻的点 int v = e.to; int w = e.w; if(dist[v] > dist[u] + w) { // 如果当前路径更优 dist[v] = dist[u] + w; // 更新最短路径距离 q.push(make_pair(dist[v], v)); // 将该点加入小根堆 } } } } int main() { int n, m, s; cin >> n >> m >> s; for(int i = 0; i < m; i++) { int u, v, w; cin >> u >> v >> w; G[u].push_back(Edge(v, w)); } dijkstra(s); for(int i = 1; i <= n; i++) { if(dist[i] == INF) cout << "INF" << endl; // 如果不连通,输出INF else cout << dist[i] << endl; } return 0; } ``` 输入格式:第一行输入三个整数n,m,s,表示图的点数、边数和源点编号。接下来m行每行三个整数u,v,w,表示一条从u到v的有向边,边权为w。 输出格式:输出n行,每行一个整数,表示源点到每个点的最短路径长度。若不连通,则输出INF。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值