【bzoj1066】【luogu2472】[SCOI2007]蜥蜴

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec Memory Limit: 162 MB

Description

　　在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃
到边界外。每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不
变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个
石柱上。

Input

　　输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱
，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。

Output

　　输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

HINT

100%的数据满足：1<=r, c<=20, 1<=d<=4

Source

最大流的做法。。。限制在点上，所以拆点，控制流量，每个点流量为高度。S向每个有蜥蜴的点连流量为1的边，表示这里最初有1只蜥蜴。每个能跳出的点向T连inf的边，然后能互相到达的点互相连inf的边。搞定

 // It is made by XZZ

 #include<cstdio>

 #include<algorithm>

 #include<cstring>

 using namespace std;

 #define rep(a,b,c) for(rg int a=b;a<=c;a++)

 #define drep(a,b,c) for(rg int a=b;a>=c;a--)

 #define erep(a,b) for(rg int a=fir[b];a;a=nxt[a])

 #define il inline

 #define rg register

 #define vd void

 #define t (dis[i])

 typedef long long ll;

 il int gi(){

     rg int x=;rg char ch=getchar();

     while(ch<''||ch>'')ch=getchar();

     while(ch>=''&&ch<='')x=x*+ch-'',ch=getchar();

     return x;

 }

 const int maxn=**+,S=,T=,maxm=(*+*)<<;

 int fir[maxn],head[maxn],nxt[maxm],dis[maxm],w[maxm],id=,dep[maxn];

 il vd add(int a,int b,int c){

     nxt[++id]=fir[a],fir[a]=id,dis[id]=b,w[id]=c;

     if(c)add(b,a,);

 }

 int num[][];

 char yes[][],h[][];

 inline bool BFS() {

     int que[maxn],hd=,tl=;

     que[]=S;bool yes[maxn]={};

     yes[S]=,dep[S]=;

     while(hd^tl) {

     int now=que[hd++];

     for(int i=fir[now];i;i=nxt[i])

         if(w[i]>&&!yes[t])yes[t]=,que[tl++]=t,dep[t]=dep[now]+;

     }

     return yes[T];

 }

 inline int Dinic(int now,int h) {

     if(now==T)return h;

     int ans=;

     for(int&i=head[now];i;i=nxt[i])

     if(w[i]>&&dep[t]==dep[now]+) {

         int D=Dinic(t,min(h,w[i]));

         w[i]-=D,w[i^]+=D,ans+=D,h-=D;

         if(h==)return ans;

     }

     return ans;

 }

 int main(){

     rg int r=gi(),c=gi(),d=gi(),ans=;

     rep(i,,r)scanf("%s",h[i]+);

     rep(i,,r)scanf("%s",yes[i]+);

     rg int Id=;

     rep(i,,r)rep(j,,c)

     if(h[i][j]-''){

         num[i][j]=++Id,++Id;add(num[i][j],num[i][j]+,h[i][j]-'');

         if(yes[i][j]=='L')add(S,num[i][j],),ans++;

         if(i<=d||j<=d||r-i+<=d||c-j+<=d)add(num[i][j]+,T,);

     }

     rep(i,,r)rep(j,,c)if(num[i][j])rep(I,,r)rep(J,,c)if(num[I][J]&&(i-I)*(i-I)+(j-J)*(j-J)<=d*d)add(num[i][j]+,num[I][J],);

     while(BFS())memcpy(head,fir,sizeof fir),ans-=Dinic(S,);

     printf("%d\n",ans);

     return ;

 }

PS.状态不佳。。。