算法学习专栏——Dijkstra(二)+bellman-ford
前言
我们现在介绍的所有这些算法都是前面的一些怪物想出来,所以你实在想不出来,也是非常正常的,每一个算法都可以写成一本大书,然后我们必须在短短的数天之内把它们完全搞懂,这就已经是一件非常具有挑战性的事情了,而我现在就是带领着大家做这么一件事情,按照我原本自己的学习周期,这些算法至少是花费了我将近5个月甚至更长的时间,但是现在我必须在2个月内将你们教会还需要熟练运用,这何尝不是一件难事,所以呢,如果你现在开始焦虑,请放下心来。痛苦的日子(我说的基础部分,高阶的另说)终将过去,开窍的日子指日可待!
一、Dijkstra求最短路 II(中等-)
给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环,所有边权均为非负值。
请你求出 1 号点到 n 号点的最短距离,如果无法从 1 号点走到 n 号点,则输出 −1。
输入格式
第一行包含整数 n 和 m。
接下来 m 行每行包含三个整数 x,y,z,表示存在一条从点 x 到点 y 的有向边,边长为 z。
输出格式
输出一个整数,表示 1 号点到 n 号点的最短距离。
如果路径不存在,则输出 −1。
数据范围
1≤n,m≤1.5×10^5^
图中涉及边长均不小于 0,且不超过 10000。
数据保证:如果最短路存在,则最短路的长度不超过 10^9^。
输入样例:
1 |
|
输出样例:
1 |
|
代码
答案(请自己先思考一下再参考)
#include < iostream> #include < queue> #include < algorithm> #include < cstring> using namespace std; const int N = 150000; typedef pair PII; int n, m; int e[N], ne[N], h[N], w[N], idx; int dist[N]; bool st[N]; void add(int x,int y,int z) { e[idx] = y,ne[idx] = h[x],w[idx] = z,h[x] = idx++; } int dijkstra() { memset(dist, 0x3f, sizeof dist); priority_queue< PII, vector< PII>, greater< PII>> q; dist[1] = 0; q.push({0, 1}); while (q.size()) { auto t = q.top(); q.pop(); int distance = t.first, k = t.second; if (st[k]) continue; st[k] = true; for (int i = h[k]; i != -1; i = ne[i]) { int j = e[i]; if (!st[j] && dist[j] > distance + w[i]) { dist[j] = distance + w[i]; q.push({dist[j], j}); } } } if (dist[n] == 0x3f3f3f3f) return -1; return dist[n]; } int main() { ios::sync_with_stdio(false),cin.tie(0),cout.tie(0); cin >> n >> m; memset(h, -1, sizeof h); while (m--) { int x, y, z; cin >> x >> y >>z; add(x, y, z); } int t = dijkstra(); cout << t; return 0; }
二、有边数限制的最短路(中等-)
给定一个 n 个点 m 条边的有向图,图中可能存在重边和自环, 边权可能为负数。
请你求出从 11 号点到 n 号点的最多经过 k 条边的最短距离,如果无法从 1 号点走到 n 号点,输出 impossible
。
注意:图中可能 存在负权回路 。
输入格式
第一行包含三个整数 n,m,k。
接下来 m 行,每行包含三个整数 x,y,z,表示存在一条从点 x 到点 y 的有向边,边长为 z。
点的编号为 1∼n。
输出格式
输出一个整数,表示从 1 号点到 n 号点的最多经过 k 条边的最短距离。
如果不存在满足条件的路径,则输出 impossible
。
数据范围
1≤n,k≤500
1≤m≤10000
1≤x,y≤n
任意边长的绝对值不超过 10000。
输入样例:
1 |
|
输出样例:
1 |
|
思路
代码
答案(请自己先思考一下再参考)
#include < iostream> #include < cstring> #include < algorithm> using namespace std; const int N = 510, M = 1e4 + 10; int dist[N], backup[N]; int n, m, k; struct Edge{ int a, b, w; }edges[M]; void bellman_ford() { memset(dist, 0x3f, sizeof dist); dist[1] = 0; for (int i = 0; i < k; i++) { memcpy(backup, dist, sizeof dist); for (int j = 1; j <= m; j++) { int a = edges[j].a, b = edges[j].b, w = edges[j].w; dist[b] = min(dist[b], backup[a] + w); } } if (dist[n] > 0x3f3f3f3f / 2) { puts("impossible"); return; } cout << dist[n]; return; } int main() { ios::sync_with_stdio(false),cin.tie(0),cout.tie(0); cin >> n >> m >> k; for (int i = 1; i <= m; i++) { int a, b, w; cin >> a >> b >> w; edges[i] = {a, b, w}; } bellman_ford(); return 0; }