【SCOI2007】bzoj1066 蜥蜴

Description

　　在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃 到边界外。
每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石
柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不
变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个 石柱上。 Input

　　输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱
，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。 Output

　　输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

把每个位置拆成入点和出点，之间连一条容量为石柱高度的边。从源点向有蜥蜴的入点连容量为1的边，从能出去的出点向汇点连容量为无穷大的边，能到达的出点向入点连容量为无穷大的边，然后求最大流。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int oo=0x3f3f3f3f;
int fir[1000],ne[1000000],to[1000000],w[1000000],tot=2,n,r,c,d,
que[1000],f[1000];
/*in:i*c+j out:(i+r)*c+j
from:0 to:1*/
void add(int num,int u,int v,int l)
{
    ne[num]=fir[u];
    fir[u]=num;
    to[num]=v;
    w[num]=l;
}
bool out(int x,int y)
{
    return x<=d||y<=d||x+d>r||y+d>c;
}
bool can(int x1,int y1,int x2,int y2)
{
    return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)<=d*d;
}
void init()
{
    int i,j,x,y;
    char s[100];
    scanf("%d%d%d",&r,&c,&d);
    for (i=1;i<=r;i++)
    {
        scanf("%s",s+1);
        for (j=1;j<=c;j++)
          if (s[j]-'0')
          {
            add(tot++,i*c+j,(i+r)*c+j,s[j]-'0');
            add(tot++,(i+r)*c+j,i*c+j,0);
          }
    }
    for (i=1;i<=r;i++)
    {
        scanf("%s",s+1);
        for (j=1;j<=c;j++)
          if (s[j]=='L')
          {
            n++;
            add(tot++,0,i*c+j,1);
            add(tot++,i*c+j,0,0);
          }
    }
    for (i=1;i<=r;i++)
      for (j=1;j<=c;j++)
        if (out(i,j))
        {
            add(tot++,(i+r)*c+j,1,oo);
            add(tot++,1,(i+r)*c+j,0);
        }
    for (i=1;i<=r;i++)
      for (j=1;j<=c;j++)
        for (x=1;x<=r;x++)
          for (y=1;y<=c;y++)
            if (can(i,j,x,y))
            {
                add(tot++,(i+r)*c+j,x*c+y,oo);
                add(tot++,x*c+y,(i+r)*c+j,0);
            }
}
bool find()
{
    memset(f,0,sizeof(f));
    que[1]=0;
    int hd=1,tl=1,u,v,i;
    f[0]=1;
    while (hd<=tl)
    {
        u=que[hd++];
        for (i=fir[u];i;i=ne[i])
          if (w[i]&&!f[v=to[i]])
          {
            f[v]=f[u]+1;
            que[++tl]=v;
          }
    }
    return f[1];
}
int dfs(int u,int lim)
{
    if (u==1) return lim;
    int ret=0,i,v,x;
    for (i=fir[u];i;i=ne[i])
      if (w[i]&&f[v=to[i]]==f[u]+1)
      {
        x=dfs(v,min(lim-ret,w[i]));
        w[i]-=x;
        w[i^1]+=x;
        ret+=x;
      }
    if (!ret) f[u]=0;
    return ret;
}
int main()
{
    int ans=0,x;
    init();
    while (find())
      while (x=dfs(0,oo))
        ans+=x;
    printf("%d\n",n-ans);
}