首页 > 代码库 > [SCOI2007]蜥蜴

[SCOI2007]蜥蜴

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 4086  Solved: 2065
[Submit][Status][Discuss]

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

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\n",ans);70     return 0;71 }

 

[SCOI2007]蜥蜴