【概述】

对于单源最短路来说，有时间复杂度为 O(E+VlogV) 要求权值非负的 Dijkstra，时间复杂度为 O(VE) 适用于带负权值的 Bellman Ford

对于全源最短路来说，除了时间复杂度为 O(V*V*V) 利用动态规划思想的 Floyd 算法外，可以认为是单源最短路径的推广，即分别以每个顶点为源点求其至其他顶点的最短距离

对于每个顶点利用 Ford 算法，时间复杂度为 O(V*V*E)，由于图中顶点都是连通的，边的数量可能会比点多，这个时间并没有比 Floyd 更优；而对于每个顶点利用 Dijkstra 算法，时间复杂度为 O(V*E+V*V*logV)，时间复杂度更优，但问题在于 Dijkstra 要求图中边权非负，不适合通用情况

针对 Dijkstra 不能适用于存在负边权的情况，Donald B. Johnson 提出了对所有边的权值进行重赋值的算法，使得边的权值非负，从而可以利用 Dijkstra 进行最短路的计算，即 Johnson 算法

【算法原理】

Johnson 算法的关键在于引入了一个势能函数 h[]，其作为边的映射关系

假设有一个新的点 n+1，其指向其他所有点的边权均为 0，计算这个新点到其他所有点的最短路径数组 h[]：对于边 (u,v)，其权值修改为 dis(u,v)=dis(u,v)+h[u]-h[v]

新增点 n+1 到 u、v 两点的最短路径分别为 h[u]、h[v]，这两点边的权重为 dis(u,v)，那么满足不等式：h[v]<=h[u]+dis(u,v)，那么可以保证新的权重为非负值

对所有边进行重赋值后，再进行 n 次 Dijkstra 来解决 n 个单源最短路问题

设 s 到 t 的最短路经过了 v1,v2,...,vk，那么边的累计和 dis'(s,t)=dis(s,t)+h[s]-h[t]，原来的最短路为：dis(s,t)=dis'(s,t)+h[v]-h[u]

综上，Johnson 算法的描述如下：

对于给定图 G=(V,E)，新增一顶点 S，对 S 到图中所有点都建一条边，得到新图 G'
对图 G' 中点 S 使用 Ford 算法计算单源最短路，得到势能函数 h[]
对原图 G 中所有的边进行重赋值：对于每条边 (u,v)，其新的权值为 dis(u,v)+h[u]-h(v)
对原图 G 的每个顶点运行 Dijkstra，求得全源最短路径

【模版】

vector<Pair> edge[N];
int dis[N][N];
int h[N];
bool vis[N];
void SPFA(int n) {memset(h, INF, sizeof(h));memset(vis, false, sizeof(vis));h[n + 1] = false;queue<int> Q;Q.push(n + 1);while (!Q.empty()) {int u = Q.front();Q.pop();for (int i = 0; i < edge[u].size(); i++) {int v = edge[u][i].first;int w = edge[u][i].second;if (h[v] > h[u] + w) {h[v] = h[u] + w;if (!vis[v]) {vis[v] = true;Q.push(v);}}}}
}
void Dijkstra(int S) {memset(dis, INF, sizeof(dis));dis[S][S] = 0;priority_queue<Pair> Q;Q.push(make_pair(0, S));while (!Q.empty()) {Pair u = Q.top();Q.pop();if (dis[S][u.second] < u.first)continue;for (int i = 0; i < edge[u.second].size(); i++) {int v = edge[u.second][i].first;int w = edge[u.second][i].second;if (dis[S][v] > dis[S][u.second] + w) {dis[S][v] = dis[S][u.second] + w;Q.push(make_pair(dis[S][v], v));}}}
}
void Johnson(int n) {SPFA(n); //计算n+1号点到其他点的距离for (int i = 1; i <= n; i++) //对所有边重新赋值for (int j = 0; j < edge[i].size(); j++)edge[i][j].second += h[i] - h[edge[i][j].first];for (int i = 1; i <= n; i++) //对所有点跑一次DijkDijkstra(i);
}
int main() {int n, m;scanf("%d%d", &n, &m);for (int i = 1; i <= m; i++) {int x, y, w;scanf("%d%d%d", &x, &y, &w);edge[x].push_back(make_pair(y, w));}for (int i = 1; i <= n; i++) //新点到其他点建边,边权为0edge[n + 1].push_back(make_pair(i, 0));Johnson(n);for (int i = 1; i <= n; i++) {for (int j = 1; j <= n; j++) {if (dis[i][j] == INF)printf("-1 ");else printf("%d ", dis[i][j] + h[j] - h[i]);printf("\n");}}return 0;
}

图论 —— 最短路 —— Johnson 算法相关推荐

图论-最短路Dijkstra算法详解超详有图解
整体来看dij就是从起点开始扩散致整个图的过程,为什么说他稳定呢,是因为他每次迭代,都能得到至少一个结点的最短路.(不像SPFA,玄学复杂度) 但是他的缺点就是不能处理带负权值的边,和代码量稍稍复杂. ...
图论 —— 最短路 —— Floyd 算法
[概述] Floyd 算法又称为插点法,是一种用于寻找给定的加权图中多源点之间最短路径的算法. 其最大特点是可以计算出现负边权时的最短路,实际应用中,很多题目不是问如何用 Floyd 求最短路,而是用 ...
图论 —— 最短路 —— Dijkstra 算法
[概述] Dijkstra 算法是单源最短路径算法,即计算起点只有一个的情况到其他点的最短路径,其无法处理存在负边权的情况. 其时间复杂度是:O(E+VlogV) [算法分析] 将点分为两类,一类是已 ...
图论 —— 最短路 —— Bellman-Ford 算法与 SPFA
[概述] Bellman-Ford算法适用于计算单源最短路径,即:只能计算起点只有一个的情况. 其最大特点是可以处理存在负边权的情况,但无法处理存在负权回路的情况. 其时间复杂度为:O(V*E),其中 ...
【图论】Johnson算法(多源最短路优化)
样例题的话随便找个Floyd算法自己测一下吧,下面上封装后的Johnson板子: (好啦-既然你们懒那我给你们找个吧..AcWing.854) 链接:https://www.acwing.com/pr ...
图论最短路：Bellman-Ford与其优化SPFA算法的一点理解
文章目录前言一.对Bellman-Ford的深入理解 1. Bellman-Ford有什么用? 2. 什么是松弛操作? 3. Bellman-Ford的k次迭代意义? 4. 一个重要定理 5. 对 ...
【阿良的算法之路】图论最短路算法模板
图论: [阿良的算法之路]图论最短路算法模板 [模板]dirjkstra单源最短路径 [模板]Bellman-Ford多源最短路 [模板]Spfa求最短路 [模板]Spfa判断负环 [模板]Floya ...
超详细！图论最短路算法与极简c++代码（配题目）
最短路算法知识结构图 /* 最短路问题{单源最短路{1.所有边权都是正数{1.朴素Dijkstra算法 O(n^2) 适用于稠密图2.堆优化的Dijkstra O(mlog(n)) 适用于稀疏图}2. ...
图论最短路及生成树(Prim,Djikstra,Spfa,Bellan-ford,kruskal,topsort)
图论在算法中具有举足轻重的地位,只有学好图才能游刃有余.本文章将介绍图论中一些基础算法,可以说总结的十分全面,文章结尾也会分析各算法的差异,清晰易懂.并附上代码模板. 图论(最短路.生成树) 一.拓扑 ...
【挑战程序设计】- 2.5 图论(最短路、最小生成树)
2.5 图论(最短路.最小生成树) 文章目录 2.5 图论(最短路.最小生成树) 2.5.1 定义们 2.5.2 图的表示 2.5.3 图的搜索 2.5.4 最短路问题单源1:bellman-for ...

图论 —— 最短路 —— Johnson 算法

【概述】

【算法原理】

【模版】

图论 —— 最短路 —— Johnson 算法相关推荐

最新文章

热门文章