【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=0;rg char ch=getchar();
    while(ch<'0'||ch>'9')ch=getchar();
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x;
}
const int maxn=21*21*2+3,S=801,T=802,maxm=(400*2+400*41)<<1;
int fir[maxn],head[maxn],nxt[maxm],dis[maxm],w[maxm],id=1,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,0);
}
int num[21][21];
char yes[21][21],h[21][21];
inline bool BFS() {
    int que[maxn],hd=0,tl=1;
    que[0]=S;bool yes[maxn]={0};
    yes[S]=1,dep[S]=0;
    while(hd^tl) {
    int now=que[hd++];
    for(int i=fir[now];i;i=nxt[i])
        if(w[i]>0&&!yes[t])yes[t]=1,que[tl++]=t,dep[t]=dep[now]+1;
    }
    return yes[T];
}
inline int Dinic(int now,int h) {
    if(now==T)return h;
    int ans=0;
    for(int&i=head[now];i;i=nxt[i])
    if(w[i]>0&&dep[t]==dep[now]+1) {
        int D=Dinic(t,min(h,w[i]));
        w[i]-=D,w[i^1]+=D,ans+=D,h-=D;
        if(h==0)return ans;
    }
    return ans;
}
int main(){
    rg int r=gi(),c=gi(),d=gi(),ans=0;
    rep(i,1,r)scanf("%s",h[i]+1);
    rep(i,1,r)scanf("%s",yes[i]+1);
    rg int Id=0;
    rep(i,1,r)rep(j,1,c)
    if(h[i][j]-'0'){
        num[i][j]=++Id,++Id;add(num[i][j],num[i][j]+1,h[i][j]-'0');
        if(yes[i][j]=='L')add(S,num[i][j],1),ans++;
        if(i<=d||j<=d||r-i+1<=d||c-j+1<=d)add(num[i][j]+1,T,2333);
    }
    rep(i,1,r)rep(j,1,c)if(num[i][j])rep(I,1,r)rep(J,1,c)if(num[I][J]&&(i-I)*(i-I)+(j-J)*(j-J)<=d*d)add(num[i][j]+1,num[I][J],2333);
    while(BFS())memcpy(head,fir,sizeof fir),ans-=Dinic(S,2333);
    printf("%d\n",ans);
    return 0;
}

PS.状态不佳。。。