[ZJOI2009]狼和羊的故事【网络流】【最大流(最小割)】

>Description

羊狼圈是一个n*m个矩阵格子，这个矩阵的边缘已经装上了篱笆。在羊狼圈中再加入一些篱笆，将羊狼分开。狼和羊都有属于自己领地，若狼和羊们不能呆在自己的领地，那它们就会变得非常暴躁，不利于他们的成长。 Orez想要添加篱笆的尽可能的短。当然这个篱笆首先得保证不能改变狼羊的所属领地，再就是篱笆必须修筑完整，也就是说必须修建在单位格子的边界上并且不能只修建一部分。

---

>Input
文件的第一行包含两个整数n和m。接下来n行每行m个整数，1表示该格子属于狼的领地，2表示属于羊的领地，0表示该格子不是任何一只动物的领地。

>Output
文件中仅包含一个整数ans，代表篱笆的最短长度。

数据范围
10%的数据 n，m≤3
30%的数据 n，m≤20
100%的数据 n，m≤100

---

>解题思路
最小割，也就是最大流完成

1. $S$向所有狼连一条权值为$inf$的边
2. 所有羊向$T$连一条权值为$inf$的边
3. 所有狼向周围的空地和羊连一条权值为$1$的边
4. 所有空地向周围的其它空地和羊连一条权值为$1$的边

边权为1很容易想到，因为在相邻的两格中间建篱笆的代价为1
这时我们只要求S~T的最小割就可以完成

乱乱的改了好久TT

---

>代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define N 110
#define inf 1 << 30
#define S n * m + 1
#define T n * m + 2
#define int long long
using namespace std;

const int xx[4] = {-1, 0, 0, 1}, yy[4] = {0, -1, 1, 0};
struct line
{
	int to, next, w, op;
} a[240005];
int n, m, ans, t, s[N][N], h[N * N], c[N * N];
bool yd[N * N];

void add (int x, int y, int ww)
{
	a[++t] = (line) {y, h[x], ww, t + 1}; h[x] = t; //顺流
	a[++t] = (line) {x, h[y], 0, t - 1}; h[y] = t; //逆流
}
bool spfa ()
{
	queue<int> st;
	memset (c, 0x7f, sizeof (c));
	memset (yd, 0, sizeof (yd));
	c[S] = 0; yd[S] = 1;
	st.push(S);
	while (!st.empty())
	{
		int u = st.front();
		for (int i = h[u]; i; i = a[i].next)
		  if (a[i].w && c[u] + 1 < c[a[i].to])
		  {
		  	c[a[i].to] = c[u] + 1;
		  	if (a[i].to == T) return 1; //如果S有路径可以到达T，直接返回1
		  	if (!yd[a[i].to])
		  	{
		  		yd[a[i].to] = 1;
		  		st.push(a[i].to);
		  	}
		  }
		yd[u] = 0;
		st.pop();
	}
	return 0;
}
int dfs (int now, int maxf)
{
	if (now == T) return maxf;
	int ret = 0;
	for (int i = h[now]; i; i = a[i].next)
	  if (a[i].w && c[now] + 1 == c[a[i].to])
	  {
	  	int f = dfs (a[i].to, min (maxf - ret, a[i].w)); //求最小割
	  	a[i].w -= f;
	  	a[a[i].op].w += f;
	  	ret += f;
	  	if (ret == maxf) return ret; //优化
	  }
	return ret;
}
int dinic () //dinic算法
{
	int ans = 0;
	while (spfa ()) ans += dfs (S, inf); //如果有增广路就进行增广
	return ans;
}

signed main()
{
	scanf ("%lld%lld", &n, &m);
	for (int i = 1; i <= n; i++)
	  for (int j = 1; j <= m; j++)
	  {
	  	scanf ("%lld", &s[i][j]);
	  	if (s[i][j] == 1) add (S, (i - 1) * m + j, inf);
	  	else if (s[i][j] == 2) add ((i - 1) * m + j, T, inf);
	  }
	for (int i = 1; i <= n; i++)
	  for (int j = 1; j <= m; j++)
	  	if (s[i][j] == 0 || s[i][j] == 1)
	  	  for (int k = 0; k < 4; k++) //四周
	      {
	      	int x = i + xx[k], y = j + yy[k];
			if (x < 1 || x > n || y < 1 || y > m) continue;
	      	if (s[x][y] == 2 || s[x][y] == 0)
	      	   add ((i - 1) * m + j, (x - 1) * m + y, 1);
	      }
	printf ("%lld", dinic ());
	return 0;
}