[SCOI2007]蜥蜴

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

题解 ：

好久没写网络流了..

对于每个石柱 拆成两个点 称为入点与出点 中间连一条容量为石柱高的边 对于一个石柱能让一个蜥蜴跳出去 就将该石柱的出点与T连一条容量为INF的边 如果一根石柱上有蜥蜴 那就把S与入点连边 根据题意 蜥蜴每一次跳跃 都会使石柱高度减一 也就是出点与入点之间容量减一 这样最大流即是最多能逃离的蜥蜴数 最后将总数减去就行了

代码：

#include<iostream>

#include<algorithm>

#include<cstdio>

#include<cstring>

#include<queue>

#define N 300003

#define inf 1000000000

usingnamespace std;

intdx[10]={0,1,0,-1},dy[10]={1,0,-1,0};

intsz,n,m,vis[23][23],mp[23][23],mp1[23][23],cnt,cur[N],a[23][23];

inttot,point[N],nxt[N],v[N],remain[N],deep[N],nowx,nowy,d,S,T;

boolpd;

voidadd(intx,inty,intz)

{

    tot++; nxt[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z;

    tot++; nxt[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0;

    //cout<<x<<" "<<y<<" "<<z<<endl;

}

boolbfs(ints,intt)

{

    for(inti=1;i<=t;i++) cur[i]=point[i];

    for(inti=1;i<=t;i++) deep[i]=t+1;

    deep[s]=0;

    queue<int> p; p.push(s);

    while(!p.empty()) {

        intnow=p.front(); p.pop();

        for(inti=point[now];i!=-1;i=nxt[i])

         if(deep[v[i]]==t+1&&remain[i])

          deep[v[i]]=deep[now]+1,p.push(v[i]);

    }

    if(deep[t]==t+1) returnfalse;

    returntrue;

}

intdfs(intnow,intt,intlimit)

{

    if(now==t||!limit)returnlimit;

    intf,flow=0;

    for(inti=cur[now];i!=-1;i=nxt[i]) {

        cur[now]=i;

        if(deep[v[i]]==deep[now]+1&&(f=dfs(v[i],t,min(limit,remain[i])))){

            flow+=f; limit-=f;

            remain[i]-=f; remain[i^1]+=f;

            if(!limit) break;

         }

    }

    returnflow;

}

intdinic(ints,intt)

{

    intans=0;

    while(bfs(s,t)) ans+=dfs(s,t,inf);

    returnans;

}

intpow(intx){

    returnx*x;

}

boolcheck(intx1,inty1,intx2,inty2)

{

    if(x1==x2&&y1==y2) returnfalse;

    returnpow(x1-x2)+pow(y1-y2)<=pow(d);

}

intmain()

{

//    freopen("in.txt","r",stdin);

    scanf("%d%d%d",&n,&m,&d);

    S=1; T=n*m*2+2;

    intsz=1; tot=-1;

    memset(point,-1,sizeof(point));

    for(inti=1;i<=n;i++)

     for(intj=1;j<=m;j++)

      mp[i][j]=++sz,mp1[i][j]=++sz;

    for(inti=1;i<=n;i++){

        chars[100]; scanf("%s",s+1);

        for(intj=1;j<=m;j++){

            intt=s[j]-'0'; a[i][j]=t;

            if(t) add(mp[i][j],mp1[i][j],t);

        }

    }

    for(inti=1;i<=n;i++){

        chars[100]; scanf("%s",s+1);

        for(intj=1;j<=m;j++)

         if(s[j]=='L') add(S,mp[i][j],1),cnt++;

    }

     

    for(inti=1;i<=n;i++)

     for(intj=1;j<=m;j++) {

        if(i+d>n||j+d>m||i-d<=0||j-d<=0) add(mp1[i][j],T,inf);

        for(intk=1;k<=n;k++)

         for(intl=1;l<=m;l++)

          if(check(i,j,k,l)&&a[i][j]&&a[k][l]) add(mp1[i][j],mp[k][l],a[i][j]);

     }

//     for(int i = 0; i < tot; i++) printf("%d  %d  %d  %d\n", i, e[i].fr, e[i].to, e[i].v);

    printf("%d\n",cnt-dinic(S,T));

}