LG·狼和羊的故事【最小割】

Description–

“狼爱上羊啊爱的疯狂，谁让他们真爱了一场；狼爱上羊啊并不荒唐，他们说有爱就有方向…”
$Youth$ 听到这首歌，心想：狼和羊如此和谐，为什么不尝试羊狼合养呢？
说干就干！ $Youth$ 的羊狼圈可以看作一个 $n\ast m$ 个矩阵格子，这个矩阵的边缘已经装上了篱笆。
可是 $Drunk$ 很快发现狼再怎么也是狼，它们总是对羊垂涎三尺，那首歌只不过是一个动人的传说而已。
所以 $Youth$ 决定在羊狼圈中再加入一些篱笆，还是要将羊狼分开来养。 通过仔细观察，$Youth$ 发现狼和羊都有属于自己领地，若狼和羊们不能呆在自己的领地，那它们就会变得非常暴躁，不利于他们的成长。
$Youth$ 想要添加篱笆的尽可能的短。当然这个篱笆首先得保证不能改变狼羊的所属领地，再就是篱笆必须修筑完整，也就是说必须修建在单位格子的边界上并且不能只修建一部分。

---

Input–

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

Output–

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

---

Sample Input–

2 2
2 2 
1 1 

Sample Output–

2

---

说明–

$10\%$ 的数据 $n，m\le 3$

$30\%$ 的数据 $n，m\le 20$

$100\%$ 的数据 $n，m\le 100$

---

解题思路–

• 原点向所有狼连流量 $inf$ 的边

• 所有羊向汇点连流量 $inf$ 的边

• 所有点向四周连流量为 $1$ 的边。

然后求最小割即可

---

代码–

#include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
#define INF 1 << 30

using namespace std;

const int fx[4][2] = {{-1, 0}, {0, -1}, {0, 1}, {1, 0}};

int n, m, s, t, o, tt = 1, a[110][110], ls[11000], dis[11000];
bool pd[11000];

struct ooo
{
	int to, x, next;
}f[250000];

int num(int aa, int bb)
{
	return (aa - 1) * m + bb;
}

void add(int aa, int bb, int cc)
{
	f[++tt] = (ooo){bb, cc, ls[aa]}, ls[aa] = tt;
	f[++tt] = (ooo){aa, 0, ls[bb]}, ls[bb] = tt;
}

void connect()
{
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= m; ++j)
      {
      	  scanf("%d", &a[i][j]);
      	  if (a[i][j] == 1) add(s, num(i, j), INF);
      	  if (a[i][j] == 2) add(num(i, j), t, INF);
	  }
	for (int i = 1; i <= n; ++i)
	  for (int j = 1; j <= m; ++j)
	    for (int k = 0; k < 4; ++k)
	    {
	    	int xx = i + fx[k][0];
	    	int yy = j + fx[k][1];
	    	if (xx > 0 && xx <= n && yy > 0 && yy <= m)
	    	  add(num(i, j), num(xx, yy), 1);
		}
}

bool spfa()
{
	queue <int> q; q.push(s);
	memset(pd, 0, sizeof(pd));
	memset(dis, 0x7f, sizeof(dis));
	pd[s] = 1, dis[s] = 0;
	while (!q.empty())
	{
		int u = q.front(); q.pop();
		for (int i = ls[u]; i; i = f[i].next)
		  if (f[i].x && dis[o = f[i].to] > dis[u] + 1)
		  {
		  	  dis[o] = dis[u] + 1;
		  	  if (o == t) return 1;
		  	  if (!pd[o])
		  	  {
		  	  	  pd[o] = 1;
				  q.push(o); 
			  }
		  }
		pd[u] = 0;
	}
	return 0;
}

int cfy(int w, int minn)
{
	if (w == t) return minn;
	int re = 0;
	for (int i = ls[w]; i; i = f[i].next)
	  if (f[i].x && dis[o = f[i].to] == dis[w] + 1)
	  {
	  	  int ff = cfy(o, min(minn - re, f[i].x));
	  	  f[i].x -= ff, f[i ^ 1].x += ff, re += ff;
	  	  if (re == minn) return re;
	  }
	return re;
}

int dinic()
{
	int ans = 0;
	while (spfa()) ans += cfy(s, INF);
	return ans;
}

int main()
{
	scanf("%d%d", &n, &m);
	s = n * m + 1, t = n * m + 2;
	connect();
	printf("%d", dinic());
	
	return 0;
}