bzoj1066 [SCOI2007]蜥蜴 网络流复制点模型

比较基本的复制点一个接收一个出边的网络流模型，

用自己向自己的边限制流出个数

注意细节：

1、判断是<1不是<0

2、n*m+

3、ij、kl

码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
using namespace std;
queue<int>q;
#define N 10000000
#define inf 1000000009
int tot=-1,xia[N],cnt,zhong[N],v[N],hou[N],yuan[N],d[N],s,t,i,j,k,l,n,m,dd;
char lin,ch[205][205];
void jia(int a,int b,int c)
{
	++tot,hou[tot]=yuan[a],yuan[a]=tot,zhong[tot]=b,v[tot]=c;	
}
void jian(int a,int b,int c)
{
jia(a,b,c);
jia(b,a,0);	
} 
bool bfs()
{
	int i;
	memset(d,0x7f,sizeof(d));
	for(i=1;i<=t;i++)
	xia[i]=yuan[i];
	q.push(s);
	d[s]=0;
	while(!q.empty())
	{
		int st=q.front();
		q.pop();
		for(i=xia[st];i!=-1;i=hou[i])
		{
			int nd=zhong[i];
			if(!v[i]||d[nd]<inf)continue;
			d[nd]=d[st]+1;
			q.push(nd);	
		}		
	}
	return d[t]<inf;
 
}
int dfs(int o,int t,int limit)
{
	if(o==t||limit==0)return limit;
	int i,flow=0,f;
	for(i=xia[o];i!=-1;i=hou[i])
	{
		xia[o]=i;
		int nd=zhong[i];
		if(d[nd]==d[o]+1&&(f=dfs(nd,t,min(limit,v[i]))))
		{
			flow+=f;
			limit-=f;
			v[i]-=f;
			v[i^1]+=f;			
		if(limit==0)break;
		}	
	}
	return flow;
}
int dinic()
{
	int daan=0;
	while(bfs())
	{
		daan+=dfs(s,t,inf);
	}	
	return daan;
}
int main()
{
	memset(yuan,-1,sizeof(yuan));
scanf("%d%d%d",&n,&m,&dd);
s=n*m*2+1;
t=n*m*2+2;
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
{
scanf("%c",&ch[i][j]);
while(ch[i][j]<'0'||ch[i][j]>'3')scanf("%c",&ch[i][j]);	
}	
for(i=1;i<=n;i++)
for(j=1;j<=m;j++)
{
scanf("%c",&lin);
while(lin!='.'&&lin!='L')scanf("%c",&lin);
	if(lin=='L')
	{
		++cnt;
	jian(s,(i-1)*m+j,1);	
	}
	
}	
	
	for(i=1;i<=n;i++)
	for(j=1;j<=m;j++)
    {
    	jian((i-1)*m+j,n*m+(i-1)*m+j,ch[i][j]-'0');
    	for(k=max(i-dd,1);k<=min(i+dd,n);k++)
    	for(l=max(j-dd,1);l<=min(j+dd,m);l++)	
    	{
    		if(fabs(k-i)+fabs(l-j)<=dd)
    		jian(n*m+(i-1)*m+j,(k-1)*m+l,inf);
		}    	
		if(i-dd<1||i+dd>n||j-dd<1||j+dd>m)jian(n*m+(i-1)*m+j,t,inf);		
		
	}
	printf("%d",cnt-dinic());	
}