题目
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;
}