[学习笔记]网络流·最小割题目选讲

• 已经 2 周没写博客了
• 这也是 2019 的第一篇博客，一个新的开始

引入

• 最小割的概念：给定一个图，含源点 $S$ 和汇点 $T$ ，每条有向边有容量
• 选出一些边删掉，使得没有从 $S$ 到 $T$ 的路径
• 求选出的边的容量和的最小值
• 最大流和最小割都是线性规划问题中的一类
• 最小割可以看成是最大流的对偶问题
• 定理：最大流 $=$ 最小割
• 证明
• 显然最小割不能小于最大流
• 我们知道，跑完最大流之后，残余网络中没有从 $S$ 到 $T$ 的路径
• 图被分为两个集合
• 根据流量平衡条件得到，流出 $S$ 所在集合的流量等于流入 $T$ 所在集合的流量
• 而横跨 $S$ 所在集合和 $T$ 所在集合的所有边构成一个割，且等于最大流
• 于是这个割就是最小割
• 证毕
• 最小割在题目中最常见的模型就是互斥选择模型（有这个模型更好的定义请在评论区指出）
• 例如： $n$ 个变量，每个变量取 $0$ 有一定代价，取 $1$ 有一定代价，外加 $m$ 个约束 $(i,j,x)$ ，表示 $i$ 和 $j$ 取不同值有额外代价 $x$ ，求最小代价和
• 下面我在 BZOJ 上选出了 $10$ 道最小割题
• 这 $10$ 道题都是互斥选择的模型

（1）[BZOJ1943][Shoi2007]Vote 善意的投票

题意

• 即为引入中举的栗子
• $n$ 个变量，第 $i$ 个变量为 $x_i$
• 当 $x_i$ 取值等于 $a_i$ 时不产生代价，否则产生 $1$ 的代价
• 有 $m$ 个二元组 $(i,j)$
• 表示 $(i,j)$ 取值不同则会产生 $1$ 的代价
• 求最小代价和
• $2\le n\le 300$

做法

• 源点 $S$ 汇点 $T$
• 对于所有的 $1\le i\le n$ ，对第 $i$ 个变量 $x_i$ 建点 $i$
• 如果 $a_i=0$
• 则由 $S$ 向 $i$ 连容量为 $1$ 的边
• 由 $i$ 向 $T$ 连容量为 $0$ 的边
• 否则由 $S$ 向 $i$ 连容量为 $0$ 的边
• $i$ 向 $T$ 连容量为 $1$ 的边
• 这时候边 $<S,i>$ 和 $<i,T>$ 显然有且仅有一条边被割掉
• 可以看成如果保留边 $<S,i>$ 则 $x_i$ 取 $0$
• 否则 $x_i$ 取 $1$
• 考虑处理变量间的影响
• 对于二元组 $(i,j)$ ，连边 $<i,j>$ 和 $<j,i>$ ，容量都为 $1$
• 这样的含义是，如果保留了边 $<S,i>$ 和 边 $<j,T>$ ，那么还需要额外割掉边 $<i,j>$ ，否则存在一条 $S$ 到 $T$ 的路径 $S->i->j->T$
• $<j,i>$ 同理
• 实际上容量为 $0$ 的边没必要建，在上面只是为了便于理解而建出来
• 求一遍最小割即为答案
• 此外，如果二元组 $(i,j)$ 是限制 $x_i$ 必须等于 $x_j$ ，则只需要把 $<i,j>$ 和 $<j,i>$ 的容量设为 $\infty$ （表示该边割不掉）即可

代码

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define For(i, a, b) for (i = a; i <= b; i++)
#define Edge(u) for (int e = adj[u], v = go[e]; e; e = nxt[e], v = go[e])
#define NF(u) for (int &e = cur[u], v = go[e]; e; e = nxt[e], v = go[e])

inline int read()
{
   
	int res = 0; bool bo = 0; char c;
	while (((c = getchar()) < '0' || c > '9') && c != '-');
	if (c == '-') bo = 1; else res = c - 48;
	while ((c = getchar()) >= '0' && c <= '9')
		res = (res << 3) + (res << 1) + (c - 48);
	return bo ? ~res + 1 : res;
}

template <class T>
inline T Min(const T &a, const T &b) {
   return a < b ? a : b;}

const int N = 610, M = 6e5 + 5, INF = 0x3f3f3f3f;

int n, m, ecnt = 1, nxt[M], adj[N], go[M], cap[M], S, T, lev[N],
len, que[N], cur[N], ans;

void add_edge(int u, int v, int w)
{
   
	nxt[++ecnt] = adj[u]; adj[u] = ecnt; go[ecnt] = v; cap[ecnt] = w;
	nxt[++ecnt] = adj[v]; adj[v] = ecnt; go[ecnt] = u; cap[ecnt] = 0;
}

bool bfs()
{
   
	int i;
	For (i, S, T) cur[i] = adj[i], lev[i] = -1;
	lev[que[len = 1] = S] = 0;
	For (i, 1, len)
	{
   
		int u = que[i];
		Edge(u) if (cap[e] && lev[v] == -1)
		{
   
			lev[que[++len] = v] = lev[u] + 1;
			if (v == T) return 1;
		}
	}
	return 0;
}

int dinic(int u, int flow)
{
   
	if (u == T) return flow;
	int res = 0, delta;
	NF(u) if (cap[e] && lev[u] < lev[v])
	{
   
		delta = dinic(v, Min(cap[e], flow - res));
		if (delta)
		{
   
			cap[e] -= delta; cap[e ^ 1] += delta;
			res += delta; if (res == flow) break;
		}
	}
	if (res != flow) lev[u] = -1;
	return res;
}

int main()
{
   
	int i, x, y;
	n = read(); m = read();
	S = 1; T = n + 2;
	For (i, 1, n)
	{
   
		x = read();
		if (x) add_edge(1 + i, T, 1);
		else add_edge(S, 1 + i, 1);
	}
	while (m--) x = read(), y = read(),
		add_edge(x + 1, y + 1, 1), add_edge(y + 1, x + 1, 1);
	while (bfs()) ans += dinic(S, INF);
	std::cout << ans << std::endl;
	return 0;
}

（2）[BZOJ2127]happiness

题意

• $n\times m$ 的变量矩阵 $x_{i,j}$
• 对于任意的 $1\le i\le n,1\le j\le m$
• 如果 $x_{i,j}=0$ 则有 $a_{i,j}$ 的收益
• 如果 $x_{i,j}=1$ 则有 $b_{i,j}$ 的收益
• 对于每个相邻（有且仅有一条公共边）的格子 $(i,j)$ $(u,v)$
• 如果 $x_{i,j}=x_{u,v}=0$ 则有 $c_{i,j,u,v}$ 的收益
• 如果 $x_{i,j}=x_{u,v}=1$ 则有 $d_{i,j,u,v}$ 的收益
• $n,m\le 100$ ， $0\le a_{i,j},b_{i,j},c_{i,j,u,v},d_{i,j,u,v}\le 5000$

做法

• 转化问题
• 如果 $x_{i,j}=0$ 则有 $b_{i,j}$ 的代价
• 如果 $x_{i,j}=1$ 则有 $a_{i,j}$ 的代价
• 对于每个相邻（有且仅有一条公共边）的格子 $(i,j)$ $(u,v)$
• 如果 $(i,j)=0$ 或 $(u,v)=0$ 则有 $d_{i,j,u,v}$ 的代价
• 如果 $(i,j)=1$ 或 $(u,v)=1$ 则有 $c_{i,j,u,v}$ 的代价
• 求
• $$\sum a+\sum b+\sum c+\sum d-最小代价和$$
• 发现与上一题不同的地方就在于 c 和 d 不相等
• 所以我们考虑设
• $$p_{i,j}=a_{i,j}+\sum _{i\le u,j\le v}{c}_{i,j,u,v}$$
• $$q_{i,j}=b_{i,j}+\sum _{i\le u,j\le v}{d}_{i,j,u,v}$$
• 由 $S$ 向点 $(i,j)$ 连容量为 $p_{i,j}$ 的边
• 点 $(i,j)$ 向 $T$ 连容量为 $q_{i,j}$ 的边
• 对于每个 $(i,j,u,v)$ （ $i\le u,j\le v$ ）
• 由点 $(i,j)$ 向点 $(u,v)$ 连容量为 $c_{i,j,u,v}$ 的边
• 由点 $(u,v)$ 向点 $(i,j$