luogu2472 [SCOI2007]蜥蜴

题目

　　https://www.luogu.org/problem/show?pid=2472

题解

　　这是一道最大流题。

　　最小化剩下的蜥蜴就是最大化出去的蜥蜴，每个位置进行拆点，中间连一条容量为这个位置高度的边。枚举点，如果两个点能到达就连一条容量为∞的边，如果一个点能出去就向超级汇T连一条容量为∞的边。每个有蜥蜴的点由S向它连一条容量为1的边，然后做最大流就行了。答案=蜥蜴数-最大流

代码

//最大流 
#include <cstdio>
#include <algorithm>
#include <cmath>
#define inf 0x3f3f3f3f
#define maxn 100000
using namespace std;
int R, C, D, next[maxn], head[maxn], to[maxn], c[maxn], num[maxn], last[maxn], d[maxn],
	Exit, cnt, S, T, tot=1;
void adde(int a, int b, int v){to[++tot]=b;c[tot]=v;next[tot]=head[a];head[a]=tot;}
void adde2(int a, int b, int v){adde(a,b,v);adde(b,a,0);}
int ISAP(int pos, int in)
{
	int t, flow=0, p;
	if(pos==T)return in;
	for(p=last[pos];p;last[pos]=p=next[p])
		if(c[p] and d[to[p]]+1==d[pos])
		{
			flow+= t=ISAP(to[p],min(in-flow,c[p]));
			c[p]-=t, c[p xor 1]+=t;
			if(Exit or in==flow)return flow;
		}
	Exit = --num[d[pos]]==0;
	++num[++d[pos]];
	last[pos]=head[pos];
	return flow;
}
char getc()
{
	char c=getchar();
	while((c<48 or c>58) and c!='L' and c!='.')c=getchar();
	return c;
}
int table(int x, int y){return C*(x-1)+y;}
void init()
{
	int i, j, ii, jj, t1, t2;
	char ch;
	scanf("%d%d%d",&R,&C,&D);
	S=R*C*2+1, T=S+1;
	for(i=1;i<=R;i++)for(j=1;j<=C;j++)
	{
		ch=getc()-48;
		adde2(table(i,j),table(i,j)+R*C,(int)ch);
		if(abs(i-1)<D or abs(j-1)<D or abs(R-i)<D or abs(C-j)<D)
			adde2(table(i,j)+R*C,T,inf);
	}
	for(i=1;i<=R;i++)for(j=1;j<=C;j++)
	{
		ch=getc();
		if(ch=='L')adde2(S,table(i,j),1),cnt++;
	}
	for(i=1;i<=R;i++)for(j=1;j<=C;j++)for(ii=1;ii<=R;ii++)for(jj=1;jj<=C;jj++)
	{
		if(i==ii and j==jj)continue;
		t1=i-ii, t2=j-jj;
		if(t1*t1+t2*t2<=D*D)adde2(table(i,j)+R*C,table(ii,jj),inf);
	}
}
int main()
{
	int flow=0;
	init();
	while(!Exit)flow+=ISAP(S,inf);
	printf("%d",cnt-flow);
	return 0;
}