洛谷P3371 【模板】单源最短路径(弱化版)

P3371 【模板】单源最短路径(弱化版) SPFA 模板。大概思想是把松弛过的点加入队列里面,从队列取出一个点 $i$ 去松弛它的出边 $j$,如果 $j$ 的距离能够被更新,那么把 $j$ 加入到队列里面,它以后可能会更新相邻的点的距离。如果没有可以松弛的点,说明所有点的最短距离都已经得到。 要求不能有负环,但是可以有负权值边。复杂度是 $O(m)$,最坏情况 $O(nm)$,其中 $n$ 是点数,$m$ 是边数。 // Date: Mon Jan 1 18:26:54 2024 #include <climits> #include <cmath> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <functional> #include <iomanip> #include <iostream> #include <map> #include <queue> #include <set> #include <sstream> #include <stack> #include <string> #include <utility> #include <vector> using namespace std; typedef long long ll; typedef unsigned long long ull; typedef vector<int> VI; typedef pair<int, int> PII; template <class T> using pq = priority_queue<T>; template <class T> using pqg = priority_queue<T, vector<T>, greater<T>>; const int INF = 0x3f3f3f3f, MOD = 1e9 + 7; const ll INFL = 0x3f3f3f3f'3f3f3f3f; const double eps = 1e-8; const int dir[8][2] = { {0, 1}, {0, -1}, {1, 0}, {-1, 0}, {1, 1}, {1, -1}, {-1, 1}, {-1, -1}, }; const ull Pr = 131; #define For(i, a, b) for (int i = int(a); i < int(b); ++i) #define Rof(i, a, b) for (int i = int(b) - 1; i >= int(a); --i) #define For1(i, a, b) for (int i = int(a); i <= int(b); ++i) #define Rof1(i, a, b) for (int i = int(b); i >= int(a); --i) #define ForE(i, j) for (int i = h[j]; i != -1; i = ne[i]) #define f1 first #define f2 second #define pb push_back #define has(a, x) (a.find(x) != a.end()) #define nemp(a) (!a.empty()) #define all(a) (a).begin(), (a).end() #define SZ(a) int((a).size()) #define NL cout << '\n'; template <class T> bool ckmin(T &a, const T &b) { return b < a ? a = b, 1 : 0; } template <class T> bool ckmax(T &a, const T &b) { return a < b ? a = b, 1 : 0; } template <typename t> istream &operator>>(istream &in, vector<t> &vec) { for (t &x : vec) in >> x; return in; } template <typename t> ostream &operator<<(ostream &out, vector<t> &vec) { int n = SZ(vec); For(i, 0, n) { out << vec[i]; if (i < n - 1) out << ' '; } return out; } void __print(int x) { cerr << x; } void __print(long x) { cerr << x; } void __print(long long x) { cerr << x; } void __print(unsigned x) { cerr << x; } void __print(unsigned long x) { cerr << x; } void __print(unsigned long long x) { cerr << x; } void __print(float x) { cerr << x; } void __print(double x) { cerr << x; } void __print(long double x) { cerr << x; } void __print(char x) { cerr << '\'' << x << '\''; } void __print(const char *x) { cerr << '\"' << x << '\"'; } void __print(const string &x) { cerr << '\"' << x << '\"'; } void __print(bool x) { cerr << (x ? "true" : "false"); } template <typename T, typename V> void __print(const pair<T, V> &x) { cerr << '{'; __print(x.first); cerr << ", "; __print(x.second); cerr << '}'; } template <typename T> void __print(const T &x) { int f = 0; cerr << '{'; for (auto &i : x) cerr << (f++ ? ", " : ""), __print(i); cerr << "}"; } void _print() { cerr << "]\n"; } template <typename T, typename... V> void _print(T t, V... v) { __print(t); if (sizeof...(v)) cerr << ", "; _print(v...); } #ifdef _DEBUG #define debug1(x) cout << #x " = " << x << endl; #define debug2(x, y) cout << #x " = " << x << " " #y " = " << y << endl; #define debug3(x, y, z) \ cout << #x " = " << x << " " #y " = " << y << " " #z " = " << z << endl; #define dbg(x...) \ cerr << "\e[91m" << __func__ << ":" << __LINE__ << " [" << #x << "] = ["; \ _print(x); \ cerr << "\e[39m" << endl; #else #define debug1 #define debug2 #define debug3 #define dbg(x...) #endif const int N = 10100, M = 500100; int n, m, s, idx, h[N], dis[N]; int e[M], ne[M], w[M]; bool st[N]; void Init() { idx = 0; memset(h, -1, sizeof h); memset(st, false, sizeof st); memset(dis, 0x3f, sizeof dis); } void Add(int a, int b, int c) { e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++; } void spfa(int s) { dis[s] = 0; st[s] = true; queue<int> q({s}); while (nemp(q)) { auto ver = q.front(); q.pop(); st[ver] = false; ForE(i, ver) { int j = e[i]; if (ckmin(dis[j], dis[ver] + w[i]) && !st[j]) { st[j] = true; q.push(j); } } } } void solve() { cin >> n >> m >> s; Init(); For(i, 0, m) { int a, b, c; cin >> a >> b >> c; Add(a, b, c); } spfa(s); For1(i, 1, n) { if (dis[i] == INF) cout << (1LL << 31) - 1 << ' '; else cout << dis[i] << ' '; } cout << '\n'; } int main(void) { #ifdef _DEBUG freopen("3371.in", "r", stdin); #endif std::ios::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); int T = 1; while (T--) { solve(); } return 0; }

January 1, 2024 · 4 min · 847 words

洛谷P4779 【模板】单源最短路径(标准版)

P4779 【模板】单源最短路径(标准版) $Dijkstra$ 模板,求解单源最短路,复杂度 $O(m \log(n))$,其中 $n$ 为点数,$m$ 为边数。 思想是把已经求出最短距离的点加入到集合中,然后把其它点到起点的距离加入到最小堆中,每次从堆顶取出元素 $i$,用它去做松弛操作,如果能够更新相邻点 $j$,那么把 $j$ 的距离加入到最小堆中。可以保证每次从堆顶取出元素之后,这个元素的最短距离就已经被确定,可以直接把它加入到集合中。 算法要求边权值都是正数。 // Date: Fri Dec 29 11:33:23 2023 #include <climits> #include <cmath> #include <cstdio> #include <cstdlib> #include <cstring> #include <algorithm> #include <functional> #include <iomanip> #include <iostream> #include <map> #include <queue> #include <set> #include <sstream> #include <stack> #include <string> #include <utility> #include <vector> using namespace std; const int INF = 0x3f3f3f3f, MOD = 1e9 + 7; const double eps = 1e-8; const int dir[8][2] = { {0, 1}, {0, -1}, {1, 0}, {-1, 0}, {1, 1}, {1, -1}, {-1, 1}, {-1, -1}, }; typedef long long ll; typedef unsigned long long ull; typedef vector<int> VI; typedef pair<int, int> PII; const ull Pr = 131; #define For(i, a, b) for (int i = int(a); i < int(b); ++i) #define Rof(i, a, b) for (int i = int(b) - 1; i >= int(a); --i) #define For1(i, a, b) for (int i = int(a); i <= int(b); ++i) #define Rof1(i, a, b) for (int i = int(b); i >= int(a); --i) #define ForE(i, j) for (int i = h[j]; i != -1; i = ne[i]) #define f1 first #define f2 second #define pb push_back #define has(a, x) (a.find(x) != a.end()) #define nemp(a) (!a.empty()) #define all(a) (a).begin(), (a).end() #define SZ(a) int((a).size()) template <typename t> istream &operator>>(istream &in, vector<t> &vec) { for (t &x : vec) in >> x; return in; } template <typename t> ostream &operator<<(ostream &out, vector<t> &vec) { int n = SZ(vec); For(i, 0, n) { out << vec[i]; if (i < n - 1) out << ' '; } return out; } #ifdef _DEBUG #define debug1(x) cout << #x " = " << x << endl; #define debug2(x, y) cout << #x " = " << x << " " #y " = " << y << endl; #define debug3(x, y, z) \ cout << #x " = " << x << " " #y " = " << y << " " #z " = " << z << endl; #else #define debug1 #define debug2 #define debug3 #endif const int N = 100100, M = 200100; int n, m, s, h[N], st[N], dis[N]; int idx, e[M], ne[M], w[M]; void Init() { idx = 0; memset(h, -1, sizeof h); memset(st, 0, sizeof st); memset(dis, 0x3f, sizeof dis); } void Add(int a, int b, int c) { e[idx] = b, w[idx] = c, ne[idx] = h[a], h[a] = idx++; } void dijkstra(int s) { dis[s] = 0; priority_queue<PII, vector<PII>, greater<PII>> q; q.push({0, s}); while (nemp(q)) { auto t = q.top(); q.pop(); int base = t.f1, ver = t.f2; if (st[ver]) continue; st[ver] = true; ForE(i, ver) { int j = e[i], tmp = base + w[i]; if (tmp < dis[j]) { dis[j] = tmp; q.push({tmp, j}); } } } } int main(void) { #ifdef _DEBUG freopen("4779.in", "r", stdin); #endif std::ios::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL); cin >> n >> m >> s; Init(); For1(i, 1, m) { int a, b, c; cin >> a >> b >> c; Add(a, b, c); } dijkstra(s); For1(i, 1, n) { int tmp = dis[i]; if (tmp == INF) cout << ((1LL << 31) - 1) << ' '; else cout << tmp << ' '; } cout << '\n'; return 0; }

December 29, 2023 · 3 min · 537 words