Johnson 全源最短路

最新推荐文章于 2024-03-10 23:27:48 发布

Ctrl AC

最新推荐文章于 2024-03-10 23:27:48 发布

阅读量355

点赞数

CC 4.0 BY-SA版权

分类专栏： # 算法小结文章标签：算法

本文链接：https://blog.csdn.net/m0_62434776/article/details/124496104

算法小结专栏收录该内容

40 篇文章

订阅专栏

这篇博客详细介绍了如何运用Johnson算法来找到加权有向图中的最短路径。首先，通过添加一个源节点s并连接所有其他节点，使得新图中所有路径的起点统一。然后，利用Bellman-Ford或SPFA算法更新s到所有点的最短距离。接着，调整原图的边权重，并移除s节点。最后，从每个节点出发重新计算最短路径。代码示例展示了完整的Johnson算法实现过程，适用于解决大规模图的最短路径问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

“ Ctrl AC！一起 AC！”

步骤：

· 对于原图，新增一个结点s，使该点与所有点之间建立一条边，权重为零。

· 对新图使用bellman或spfa求得新点s到其他点的最短距离，h[0]....h[V-1]

· 对原图进行边权更新，w`(u,v)=w(u,v)+(h[u]-h[v])

· 移除新的顶点s，以每个顶点为起点求出最短路

代码：

配套题目：Johnson

#include <cstring>
#include <iostream>
#include <queue>
#define INF 1e9
using namespace std;
struct edge {
    int v, w, next;
} e[10005];
struct node {
    int dis, id;
    bool operator<(const node& a) const { return dis > a.dis; }
    node(int d, int x) { dis = d, id = x; }
};
int head[5005], vis[5005], t[5005];
int cnt, n, m;
long long h[5005], dis[5005];
void addedge(int u, int v, int w) {
    e[++cnt].v = v;
    e[cnt].w = w;
    e[cnt].next = head[u];
    head[u] = cnt;
}
bool spfa(int s) {
    queue<int> q;
    memset(h, 63, sizeof(h));
    h[s] = 0, vis[s] = 1;
    q.push(s);
    while (!q.empty()) {
        int u = q.front();
        q.pop();
        vis[u] = 0;
        for (int i = head[u]; i; i = e[i].next) {
            int v = e[i].v;
            if (h[v] > h[u] + e[i].w) {
                h[v] = h[u] + e[i].w;
                if (!vis[v]) {
                    vis[v] = 1;
                    q.push(v);
                    t[v]++;
                    if (t[v] == n + 1) return false;
                }
            }
        }
    }
    return true;
}
void dijkstra(int s) {
    priority_queue<node> q;
    for (int i = 1; i <= n; i++) dis[i] = INF;
    memset(vis, 0, sizeof(vis));
    dis[s] = 0;
    q.push(node(0,s));
    while (!q.empty()) {
        int u = q.top().id;
        q.pop();
        if (vis[u]) continue;
        vis[u] = 1;
        for (int i = head[u]; i; i = e[i].next) {
            int v = e[i].v;
            if (dis[v] > dis[u] + e[i].w) {
                dis[v] = dis[u] + e[i].w;
                if (!vis[v]) q.push(node(dis[v],v));
            }
        }
    }
    return;
}
int main() {
    cin >> n >> m;
    for (int i = 1; i <= m; i++) {
        int u, v, w;
        cin >> u >> v >> w;
        addedge(u, v, w);
    }
    for (int i = 1; i <= n; i++) addedge(0, i, 0);
    if (!spfa(0)) {
        cout << -1 << endl;
        return 0;
    }
    for (int u = 1; u <= n; u++)
        for (int i = head[u]; i; i = e[i].next) e[i].w += h[u] - h[e[i].v];
    for (int i = 1; i <= n; i++) {
        dijkstra(i);
        long long ans = 0;
        for (int j = 1; j <= n; j++) {
            if (dis[j] == INF)
                ans += j * INF;
            else
                ans += j * (dis[j] + h[j] - h[i]);
        }
        cout << ans << endl;
    }
    return 0;
}

感谢阅读！！！

“ Ctrl AC！一起 AC！”