[SCOI2007]蜥蜴

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 4086 Solved: 2065

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

Pku 2711 Leapin' Lizards

思路

最大流；

s--->1--->石柱--->柱高--->石柱--->t;

代码实现

 1 #include<cstdio>
 2 #include<cstring>
 3 const int maxn=1e3;
 4 const int maxm=1e6;
 5 const int inf=2147483647;
 6 inline int min_(int x,int y){return x<y?x:y;}
 7 inline int max_(int x,int y){return x>y?x:y;}
 8 int n,m,l,s,t,ans;
 9 int a;
10 char ch[30];
11 int h[maxn],hs=1;
12 int et[maxm],en[maxm],ew[maxm];
13 void add(int u,int v,int w){
14     ++hs,et[hs]=v,ew[hs]=w,en[hs]=h[u],h[u]=hs;
15     ++hs,et[hs]=u,ew[hs]=0,en[hs]=h[v],h[v]=hs;
16 }
17 int d[maxn],q[maxn],head,tail;
18 void dfs(){
19     memset(d,0,sizeof(d));
20     head=tail=0;
21     d[s]=1,q[head++]=s;
22     while(head>tail){
23         a=q[tail++];
24         for(int i=h[a];i;i=en[i])
25         if(!d[et[i]]&&ew[i]){
26             d[et[i]]=d[a]+1;
27             if(et[i]==t) return;
28             q[head++]=et[i];
29         }
30     }
31 }
32 int ap(int k,int w){
33     if(k==t) return w;
34     int uw=w;
35     for(int i=h[k];i&&uw;i=en[i])
36     if(d[et[i]]==d[k]+1&&ew[i]){
37         int nw=ap(et[i],min_(uw,ew[i]));
38         if(nw) ew[i]-=nw,ew[i^1]+=nw,uw-=nw;
39         else d[et[i]]=0;
40     }
41     return w-uw;
42 }
43 void Dinic(){while(dfs(),d[t]) ans-=ap(s,inf);}
44 int main(){
45     scanf("%d%d%d",&n,&m,&l);
46     s=0,t=1;
47     for(int i=1;i<=n;i++){
48         scanf("%s",ch+1);
49         for(int j=1;j<=m;j++){
50             a=ch[j]-'0';
51             int qp=i*m+j,zp=qp+500;
52             if(a){
53                 add(qp,zp,a);
54                 for(int ii=max_(1,i-l);ii<=min_(i+l,n);ii++)
55                 for(int jj=max_(1,j-l);jj<=min_(j+l,m);jj++)
56                 if((ii-i)*(ii-i)+(jj-j)*(jj-j)<=l*l) add(zp,ii*m+jj,inf);
57                 if(i-l<1||i+l>n||j-l<1||j+l>m) add(zp,t,inf);
58             }
59         }
60     }
61     for(int i=1;i<=n;i++){
62         scanf("%s",ch+1);
63         for(int j=1;j<=m;j++){
64             int qp=i*m+j;
65             if(ch[j]=='L') ans++,add(s,qp,1);
66         }
67     }
68     Dinic();
69     printf("%d
",ans);
70     return 0;
71 }