【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

1

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.状态不佳。。。