summaryrefslogtreecommitdiff
path: root/graph/minCostMaxFlow.cpp
diff options
context:
space:
mode:
authorGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
committerGloria Mundi <gloria@gloria-mundi.eu>2024-11-16 01:24:14 +0100
commit98567ec798aa8ca2cfbcb85c774dd470f30e30d4 (patch)
tree5113d5cc24d1ad5f93810b6442ce584a36950dc8 /graph/minCostMaxFlow.cpp
parentad3856a6b766087df0036de0b556f4700a6498c9 (diff)
parent8d11c6c8213f46f0fa19826917c255edd5d43cb1 (diff)
mzuenni tests
Diffstat (limited to 'graph/minCostMaxFlow.cpp')
-rw-r--r--graph/minCostMaxFlow.cpp66
1 files changed, 0 insertions, 66 deletions
diff --git a/graph/minCostMaxFlow.cpp b/graph/minCostMaxFlow.cpp
deleted file mode 100644
index 14a222c..0000000
--- a/graph/minCostMaxFlow.cpp
+++ /dev/null
@@ -1,66 +0,0 @@
-constexpr ll INF = 1LL << 60; // Größer als der maximale Fluss.
-struct MinCostFlow {
- struct edge {
- int to;
- ll f, cost;
- };
- vector<edge> edges;
- vector<vector<int>> adj;
- vector<int> pref, con;
- vector<ll> dist;
- const int s, t;
- ll maxflow, mincost;
-
- MinCostFlow(int n, int source, int target) :
- adj(n), s(source), t(target) {};
-
- void addEdge(int u, int v, ll c, ll cost) {
- adj[u].push_back(sz(edges));
- edges.push_back({v, c, cost});
- adj[v].push_back(sz(edges));
- edges.push_back({u, 0, -cost});
- }
-
- bool SPFA() {
- pref.assign(sz(adj), -1);
- dist.assign(sz(adj), INF);
- vector<bool> inqueue(sz(adj));
- queue<int> queue;
- dist[s] = 0;
- queue.push(s);
- pref[s] = s;
- inqueue[s] = true;
- while (!queue.empty()) {
- int cur = queue.front(); queue.pop();
- inqueue[cur] = false;
- for (int id : adj[cur]) {
- int to = edges[id].to;
- if (edges[id].f > 0 &&
- dist[to] > dist[cur] + edges[id].cost) {
- dist[to] = dist[cur] + edges[id].cost;
- pref[to] = cur;
- con[to] = id;
- if (!inqueue[to]) {
- inqueue[to] = true;
- queue.push(to);
- }}}}
- return pref[t] != -1;
- }
-
- void extend() {
- ll w = INF;
- for (int u = t; pref[u] != u; u = pref[u])
- w = min(w, edges[con[u]].f);
- maxflow += w;
- mincost += dist[t] * w;
- for (int u = t; pref[u] != u; u = pref[u]) {
- edges[con[u]].f -= w;
- edges[con[u] ^ 1].f += w;
- }}
-
- void mincostflow() {
- con.assign(sz(adj), 0);
- maxflow = mincost = 0;
- while (SPFA()) extend();
- }
-};