题目大意
已知n只跳跃距离为D的蜥蜴被困在石柱林中,它们可以通过在石柱上跳跃来逃脱,但每个石柱都有一个通过的限制,问最少有多少只蜥蜴最终困在石柱林中。
解题分析
ans可以转换成总蜥蜴数-最大逃脱蜥蜴数。
每个石柱明显说有通过的限制,那么网络流建模,每个石柱拆点,流量限制为通过限制,建立超级源和超级汇,等等… 总之不错的网络流建模题。
解题代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 805
#define maxe 28005
#define INF (((1<<30)-1)<<1)+1
using namespace std;
int n,m,d,k,s,t,num,tot,ans,son[maxe],nxt[maxe],cap[maxe],flow[maxe],lnk[maxn],lst[maxn],dst[maxn],que[maxn];
bool vs[maxn];
void _add(int x,int y,int z){
son[++tot]=y; cap[tot]=z; flow[tot]=0; nxt[tot]=lnk[x]; lnk[x]=tot;
son[++tot]=x; cap[tot]=0; flow[tot]=0; nxt[tot]=lnk[y]; lnk[y]=tot;
}
char getc(){char ch=getchar(); while (('0'>ch||ch>'3')&&ch!='L'&&ch!='.') ch=getchar(); return ch;}
int geti(int id,int x,int y){return id*n*m+(x-1)*m+y;}
void _init(){
freopen("lizard.in","r",stdin);
freopen("lizard.out","w",stdout);
scanf("%d%d%d",&n,&m,&d); tot=1; t=k=2*n*m+2; s=k-1;
memset(lnk,0,sizeof(lnk));
memset(nxt,0,sizeof(nxt));
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++){
int x=getc()-'0';
if (x) _add(geti(0,i,j),geti(1,i,j),x); //石柱的通过限制
}
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (getc()=='L') {num++; _add(s,geti(0,i,j),1);} //超级源向出事蜥蜴位置连边
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
for (int x=max(1,i-d);x<=min(n,i+d);x++)
for (int y=max(1,j-d);y<=min(m,j+d);y++)
if ((i-x)*(i-x)+(j-y)*(j-y)<=d*d&&(i!=x||j!=y))
_add(geti(1,i,j),geti(0,x,y),INF); //相邻石柱连边
for (int i=1;i<=n;i++)
for (int j=1;j<=m;j++)
if (i<=d||j<=d||n-i+1<=d||m-j+1<=d) _add(geti(1,i,j),t,INF); //超级汇
}
bool _bfs(){
memset(vs,0,sizeof(vs));
int hed=0,til=1; vs[t-1]=1; dst[t-1]=0; que[1]=t-1;
while (hed!=til){
int x=que[++hed];
for (int j=lnk[x];j;j=nxt[j])
if (cap[j]>flow[j]&&(!vs[son[j]])){
dst[son[j]]=dst[x]+1;
vs[son[j]]=1; que[++til]=son[j];
if (son[j]==k) return true;
}
}
return false;
}
int _dfs(int x,int t,int now){
if (now==0||x==t) return now; int tem=0;
for (int j=lst[x];j;lst[x]=j=nxt[j])
if (dst[son[j]]==dst[x]+1){
int ew=_dfs(son[j],t,min(now,cap[j]-flow[j]));
if (ew){
flow[j]+=ew; flow[j^1]-=ew; now-=ew; tem+=ew;
if (!now) break;
}
}
return tem;
}
void _solve(){
ans=0;
while (_bfs()){
for (int i=1;i<=k;i++) lst[i]=lnk[i];
ans+=_dfs(s,t,INF);
}
printf("%d",num-ans);
}
int main()
{
_init();
_solve();
return 0;
}