最大流

Dinic算法是一种用于解决最大流问题的图算法,它基于增广路径的思想。最大流问题是在一个有向图中找到从源点到汇点的最大流量的路径。

复杂度为 \(O(V^2 E)\)

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
constexpr int inf = 1e9;

template <typename T> struct MaxFlow {
struct edge {
int to;
T cap;
edge(int to, T cap) : to(to), cap(cap) {}
};
const int n;
std::vector<edge> E;
std::vector<std::vector<int>> G;
std::vector<int> cur, h;
MaxFlow(int n) : n(n), G(n) {}
bool bfs(int s, int t) {
h.assign(n, -1);
std::queue<int> Q;
h[s] = 0;
Q.push(s);
while (!Q.empty()) {
const int u = Q.front();
Q.pop();
for (int i : G[u]) {
auto [v, c] = E[i];
if (c > 0 && h[v] == -1) {
h[v] = h[u] + 1;
if (v == t) {
return true;
}
Q.push(v);
}
}
}
return false;
}
T dfs(int u, int t, T f) {
if (u == t) {
return f;
}
auto r = f;
for (int &i = cur[u]; i < int(G[u].size()); i++) {
const int j = G[u][i];
auto [v, c] = E[j];
if (c > 0 && h[v] == h[u] + 1) {
auto a = dfs(v, t, std::min(r, c));
E[j].cap -= a;
E[j ^ 1].cap += a;
r -= a;
if (r == 0) {
return f;
}
}
}
return f - r;
}
void add_edge(int u, int v, T cap) {
G[u].push_back(E.size());
E.emplace_back(v, cap);
G[v].push_back(E.size());
E.emplace_back(u, 0);
}
T Dinic(int s, int t) {
T ans = 0;
while (bfs(s, t)) {
cur.assign(n, 0);
ans += dfs(s, t, inf);
}
return ans;
}
};