[BZOJ1412][ZJOI2009]狼和羊的故事

题面戳我
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

题解

所有可能的相邻都需要被分割开,可能的相邻包括“羊->空地->狼”以及“羊->狼”,所以再规定一下方向就好了:
新建源汇,源点连向所有羊,容量为(inf)
所有的羊/空地,连向相邻的空地/狼,容量为(1)
所有狼连向汇点,容量为(inf)
最小割即可

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int N = 10005;
const int inf = 1e9;
struct edge{int to,next,w;}a[N<<4];
int n,m,s,t,mp[105][105],head[N],cnt=1,dep[N],cur[N];
int dx[4]={1,0,-1,0},dy[4]={0,1,0,-1};
queue<int>Q;
int gi()
{
	int x=0,w=1;char ch=getchar();
	while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
	if (ch=='-') w=0,ch=getchar();
	while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
	return w?x:-x;
}
void link(int u,int v,int w)
{
	a[++cnt]=(edge){v,head[u],w};head[u]=cnt;
	a[++cnt]=(edge){u,head[v],0};head[v]=cnt;
}
bool bfs()
{
	memset(dep,0,sizeof(dep));
	dep[s]=1;Q.push(s);
	while (!Q.empty())
	{
		int u=Q.front();Q.pop();
		for (int e=head[u];e;e=a[e].next)
			if (a[e].w&&!dep[a[e].to])
				dep[a[e].to]=dep[u]+1,Q.push(a[e].to);
	}
	return dep[t];
}
int dfs(int u,int flow)
{
	if (u==t)
		return flow;
	int ret=0;
	for (int &e=cur[u];e;e=a[e].next)
		if (a[e].w&&dep[a[e].to]==dep[u]+1)
		{
			int temp=dfs(a[e].to,min(flow,a[e].w));
			flow-=temp;ret+=temp;
			a[e].w-=temp;a[e^1].w+=temp;
		}
	return ret;
}
int Dinic()
{
	int res=0;
	while (bfs())
	{
		for (int i=t;i;i--) cur[i]=head[i];
		while (int temp=dfs(s,inf)) res+=temp;
	}
	return res;
}
int main()
{
	n=gi();m=gi();s=n*m+1;t=s+1;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
		{
			mp[i][j]=gi();
			if (mp[i][j]==2) link(s,(i-1)*m+j,inf);
			if (mp[i][j]==1) link((i-1)*m+j,t,inf);
		}
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++)
			if (mp[i][j]!=1)
				for (int k=0;k<4;k++)
				{
					int x=i+dx[k],y=j+dy[k];
					if (x<1||x>n||y<1||y>m||mp[x][y]==2) continue;
					link((i-1)*m+j,(x-1)*m+y,1);
				}
	printf("%d
",Dinic());
	return 0;
}
原文地址:https://www.cnblogs.com/zhoushuyu/p/8146104.html