bzoj1412[狼和羊的故事]最小割

http://www.lydsy.com/JudgeOnline/problem.php?id=1412

很容易看出题目是要求把狼和羊分成两个独立的集合，并使代价最小。

很明显是最小割嘛。

把狼和羊分别与起点和终点相连。

一开始我没考虑空地的情况，后来才知道空地之间要连一条边。

然后巴拉巴拉上ISAP:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
const int inf=0x3f3f3f3f;
int n,m;
int ss[105][105];
struct edge
{
	int v,next,cap;
}e[200000];
int head[10020];
int S,T;
int dis[10020];
int k;
void add(int u,int v,int cap)
{
	e[k].v=v;
	e[k].cap=cap;
	e[k].next=head[u];
	head[u]=k++;
	
	e[k].v=u;
	e[k].cap=0;
	e[k].next=head[v];
	head[v]=k++;
}
int get(int i,int j)
{
	return m*(i-1)+j;
}
int dx[]={0,0,1,-1};
int dy[]={1,-1,0,0};
void work(int i,int j)
{
	for(int l=0;l<4;l++)
	{
		int nx=i+dx[l];
		int ny=j+dy[l];
		if(nx>=1&&nx<=n&&ny>=1&&ny<=m)
		{
			/*if(ss[nx][ny]==1&&ss[i][j]==2)
			{
				add(get(i,j),get(nx,ny),1);
			}
			if(ss[nx][ny]==0)
			{
				if(ss[i][j]==2)add(get(i,j),get(nx,ny),1);
				if(ss[i][j]==1)add(get(nx,ny),get(i,j),1);
				if(ss[i][j]==0)add(get(i,j),get(nx,ny),1);
			}*/
			add(get(i,j),get(nx,ny),1);
		}
	}
	
}
bool bfs()
{
	queue<int >q;
	memset(dis,-1,sizeof(dis));
	q.push(S);
	dis[S]=0;
	
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=head[u];~i;i=e[i].next)
		{
			int v=e[i].v;
			if(dis[v]==-1&&e[i].cap>0)
			{
				dis[v]=dis[u]+1;
				if(v==T)return true;
				q.push(v);
			}
		}
	}
	return false;
}
int dfs(int u,int low)
{
	if(u==T||low==0)return low;
	int ret=0;
	for(int i=head[u];~i;i=e[i].next)
	{
		int v=e[i].v;
		if(dis[v]==dis[u]+1)
		{
			int f=dfs(v,min(low,e[i].cap));
			e[i].cap-=f;
			e[i^1].cap+=f;
			ret+=f;
			low-=f;
			if(low==0)return ret;
		}
	}
	dis[u]=-1;
	return ret;
}
void ISAP()
{
	int ans=0;
	while(bfs())
	{
		int x;
		if(x=dfs(S,inf))ans+=x;
	}
	printf("%d\n",ans);
}
int main()
{
	memset(head,-1,sizeof(head));
	k=0;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			scanf("%d",&ss[i][j]);
		}
	}
	S=0,T=n*m+1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(ss[i][j]==2)
			{
				add(S,get(i,j),inf);
			}
			else if(ss[i][j]==1)
			{
				add(get(i,j),T,inf);				
			}
			work(i,j);
		}
	}
	ISAP();
	return 0;
}