BZOJ 1066 最大流

    这道题应该能够感觉到用网络流来做，将一个石柱与其他蜥蜴能够跳到的石柱连一条流量为该石柱高度的边，（网上大多认为题目中的距离是欧几里得距离，我是用曼哈顿距离做的，也过了），表示这个石柱最多能够供这么多只蜥蜴跳，由于一个点可能连出去很多条边，而我们需要保证的是从这些边流走的流量和要小于H，因此我们需要拆点限流，即对于每个节点建立一个虚拟节点，并将虚拟节点与该节点连一条流量为H的边，当有其他点要连这个点时，就连到它对应的虚拟节点上，当它要连其他点时，就用它本身和那些点的虚拟节点连边。我们将源点与所有有蜥蜴的石柱连一条流量为1的边，将所有能跳出去的石柱与汇点连一条流量为正无穷的边，然后跑一遍dinic就好。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 100005
#define INF 1000000000
char s[50][50];
int pre[maxn],last[maxn],other[maxn],cap[maxn];
int n,m,r,d[1005],S,T,l,q[maxn+5],tot,ans;

void connect(int x,int y,int z)
{
	//printf("%d %d %d\n",x,y,z);
	pre[l]=last[x];
	last[x]=l;
	other[l]=y;
	cap[l]=z;
	l++;
	swap(x,y);
	pre[l]=last[x];
	last[x]=l;
	other[l]=y;
	cap[l]=0;
	l++;	
}

bool bfs(void)
{
	memset(d,0,sizeof d);
	int h=0,t=1;
	q[1]=S;d[S]=1;
	while (h!=t)
	{
		h=h%maxn+1;
		int u=q[h];
		for (int p=last[u];p!=-1;p=pre[p]) 
		{
			int v=other[p];
			if (d[v]||cap[p]==0) continue;
			d[v]=d[u]+1;
		//	printf("%d %d %d %d\n",u,v,d[v],cap[p]);
			if (v==T) return 1;
			t=t%maxn+1;
			q[t]=v;	
		}
	}
	return 0;
}

int dinic(int u,int flow)
{
	int rest=flow;
	if (u==T) return flow;
	for (int p=last[u];p!=-1;p=pre[p]) 
	{
		int v=other[p];
		if (rest>0&&cap[p]>0&&d[v]==d[u]+1) 
		{
			int temp=dinic(v,min(rest,cap[p]));
			rest-=temp;
			cap[p]-=temp;
			cap[p^1]+=temp;	
		}
	}
	return flow-rest;
}

int main()
{
	scanf("%d%d%d",&n,&m,&r);
	memset(last,-1,sizeof last);
	S=0;T=n*m*2+1;
	for (int i=1;i<=n;i++)	
		scanf("%s",s[i]);
	for (int i=1;i<=n;i++)
		for (int j=0;j<m;j++)
		{
			if (s[i][j]!='0') 
			{
				int a=(i-1)*m+j+1;
				int c=s[i][j]-'0';
				connect(a+n*m,a,c);
				int up=(i-r)<1?1:i-r;
				int down=(i+r)>n?n:i+r;
				int L=(j+1-r)<1?1:j+1-r;
				int R=(j+1+r)>n?n:j+1+r;
				for (int k=up;k<=down;k++)
					for (int w=L;w<=R;w++)
					{
						if ((k==i&&w==j+1)||abs(k-i)+abs(w-(j+1))>r||s[k][w-1]=='0') continue;
						connect(a,n*m+(k-1)*m+w,c);	
					}
				if (i-r<1||i+r>n||j-r<0||j+r>=m) connect(a,T,INF);
			}
		}
	for (int i=1;i<=n;i++) 
		scanf("%s",s[i]);
	for (int i=1;i<=n;i++)
		for (int j=0;j<m;j++)
			if (s[i][j]=='L') 
			{
				connect(S,n*m+(i-1)*m+j+1,1);
				tot++;
			}
	while (bfs())
		ans+=dinic(S,INF);	
	printf("%d\n",tot-ans);
	return 0;	
}