bzoj 1412: [ZJOI2009]狼和羊的故事（最小割）

1412: [ZJOI2009]狼和羊的故事

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 2399   Solved: 1229
[ Submit][ Status][ Discuss]

Description

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

Input

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

Output

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

Sample Input

2 2
2 2
1 1

Sample Output

2

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

HINT

Source

[ Submit][ Status][ Discuss] 

题解：最小割等于最大流。

所有的狼向源点连一条容量为INF 的边，所有的羊向汇点连一条容量为INF的边，然后矩阵中的每个节点向他相邻的节点连边，容量为一，注意因为矩阵中走时没有方向所以一条边要正反加两次。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define N 200003
#define inf 1000000000
using namespace std;
int n,m;
int next[N],point[N],v[N],remain[N],a[103][103];
int deep[N],cur[N],num[N],tot,laste[N];
void add(int x,int y,int z)
{
	tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;
	tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0;
}
int addflow(int s,int t)
{
	int now=t; int ans=inf;
	while (now!=s)
	 {
	 	ans=min(ans,remain[laste[now]]);
	 	now=v[laste[now]^1];
	 }
	now=t;
	while (now!=s)
	 {
	 	remain[laste[now]]-=ans;
	 	remain[laste[now]^1]+=ans;
	 	now=v[laste[now]^1];
	 }
	return ans;
}
void bfs(int s,int t)
{
	for (int i=s;i<=t;i++)
	 deep[i]=t;
	deep[t]=0;
	queue<int> p;
	p.push(t);
	while (!p.empty())
	{
		int now=p.front(); p.pop();
		for (int i=point[now];i!=-1;i=next[i])
		 if (deep[v[i]]==t&&remain[i^1])
		  {
		  	deep[v[i]]=deep[now]+1;
		  	p.push(v[i]);
		  }
	}
}
int isap(int s,int t)
{
	bfs(s,t);
	int ans=0;
	for (int i=s;i<=t;i++)
	 cur[i]=point[i];
	for (int i=s;i<=t;i++)
	 num[deep[i]]++;
	int now=s;
	while (deep[s]<t)
	{
		if (now==t)
		 {
		 	ans+=addflow(s,t);
		 	now=s;
		 }
		bool f=false;
		for (int i=cur[now];i!=-1;i=next[i])
		 if (deep[v[i]]+1==deep[now]&&remain[i])
		  {
		  	f=true;
		  	laste[v[i]]=i;
		  	cur[now]=i;
		  	now=v[i];
		  	break;
		  }
		if (!f)
		{
			int minn=t;
			for (int i=point[now];i!=-1;i=next[i])
			 if (remain[i]) minn=min(deep[v[i]],minn);
			if (!--num[deep[now]]) break;
			deep[now]=minn+1;
			num[deep[now]]++;
			cur[now]=point[now];
			if (now!=s)
			 now=v[laste[now]^1];
		}
	}
    return ans;
}
int main()
{
	freopen("a.in","r",stdin);
	freopen("my.out","w",stdout);
	scanf("%d%d",&n,&m);
	int num=n*m+2;
	memset(next,-1,sizeof(next));
	memset(point,-1,sizeof(point));
	tot=-1;
	for (int i=1;i<=n;i++)
	 for (int j=1;j<=m;j++)
	  {
	  	int x; scanf("%d",&x);
	  	int t=(i-1)*m+j+1;
	  	if (x==2) 
	  	 add(1,t,inf);
	  	if (x==1)
	  	 add(t,num,inf);
	  	if (j!=m)
	  	 add(t,(i-1)*m+j+2,1),add(t+1,t,1);
	  	if (i!=n)
	  	 add(t,i*m+j+1,1),add(i*m+j+1,t,1);
	  }
	int t=isap(1,num);
	printf("%d\n",t);
}