最小割
概念
割
对于一个网络流图
割的容量
我们的定义割
最小割
最小割就是求得一个割
证明
最大流最小割定理
定理:
对于任意一个可行流
如果我们求出了最大流
结合前面的不等式,我们可以知道此时
代码
最小割
通过 最大流最小割定理,我们可以直接得到如下代码:
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <queue>
const int N = 1e4 + 5, M = 2e5 + 5;
int n, m, s, t, tot = 1, lnk[N], ter[M], nxt[M], val[M], dep[N], cur[N];
void add(int u, int v, int w) {
ter[++tot] = v, nxt[tot] = lnk[u], lnk[u] = tot, val[tot] = w;
}
void addedge(int u, int v, int w) { add(u, v, w), add(v, u, 0); }
int bfs(int s, int t) {
memset(dep, 0, sizeof(dep));
memcpy(cur, lnk, sizeof(lnk));
std::queue<int> q;
q.push(s), dep[s] = 1;
while (!q.empty()) {
int u = q.front();
q.pop();
for (int i = lnk[u]; i; i = nxt[i]) {
int v = ter[i];
if (val[i] && !dep[v]) q.push(v), dep[v] = dep[u] + 1;
}
}
return dep[t];
}
int dfs(int u, int t, int flow) {
if (u == t) return flow;
int ans = 0;
for (int &i = cur[u]; i && ans < flow; i = nxt[i]) {
int v = ter[i];
if (val[i] && dep[v] == dep[u] + 1) {
int x = dfs(v, t, std::min(val[i], flow - ans));
if (x) val[i] -= x, val[i ^ 1] += x, ans += x;
}
}
if (ans < flow) dep[u] = -1;
return ans;
}
int dinic(int s, int t) {
int ans = 0;
while (bfs(s, t)) {
int x;
while ((x = dfs(s, t, 1 << 30))) ans += x;
}
return ans;
}
int main() {
scanf("%d%d%d%d", &n, &m, &s, &t);
while (m--) {
int u, v, w;
scanf("%d%d%d", &u, &v, &w);
addedge(u, v, w);
}
printf("%d\n", dinic(s, t));
return 0;
}
方案
我们可以通过从源点
void dfs(int u) {
vis[u] = 1;
for (int i = lnk[u]; i; i = nxt[i]) {
int v = ter[i];
if (!vis[v] && val[i]) dfs(v);
}
}
割边数量
如果需要在最小割的前提下最小化割边数量,那么先求出最小割,把没有满流的边容量改成
问题模型 1
有
这是一个经典的 二者选其一 的最小割题目。我们对于每个集合设置源点
注意到当源点和汇点不相连时,代表这些点都选择了其中一个集合。如果将连向
最小割就是最小花费。
问题模型 2
最大权值闭合图,即给定一张有向图,每个点都有一个权值(可以为正或负或
做法:建立超级源点
几个小结论来证明:
- 每一个符合条件的子图都对应流量网络中的一个割。因为每一个割将网络分为两部分,与
相连的那部分满足没有边指向另一部分,于是满足上述条件。这个命题是充要的。 - 最小割所去除的边必须与
和 其中一者相连。因为否则边权是 ,不可能成为最小割。 - 我们所选择的那部分子图,权值和
所有正权值之和 我们未选择的正权值点的权值之和 我们选择的负权值点的权值之和。当我们不选择一个正权值点时,其与 的连边会被断开;当我们选择一个负权值点时,其与 的连边会被断开。断开的边的边权之和即为割的容量。于是上述式子转化为:权值和 所有正权值之和 割的容量。 - 于是得出结论,最大权值和
所有正权值之和 最小割 所有正权值之和 最大流。
习题
贡献者:@Imple@Luo@laocongsc@huaruoji@Ir1d@芊枫@ksyx@Xiaodai@mgt@Henry-ZHR@ouuan@Siyuan
本页面最近更新:2/3/2023, 12:00:00 AM,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用