hdu 2732（网络流）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2732

思路：如果从某一根柱子能够直接跳到迷宫的外面，那么我们就将这个点连接到汇点。对于哪些不能跳出去但是又有柱子的点，那么 我们就去按照跳跃距离搜寻有没有其他的柱子能够去跳跃，如果能够找到的话，那么连接这两点，并且将容量控制为弧尾节点的柱子数，也正是由于一条弧只能够约 束一个顶点，所以我们需要进行拆点，点内之间流量为本身柱子数。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define MAXN 2222
#define MAXM 2222222
#define inf 1<<30
struct Edge{
   int v,cap,next;
}edge[MAXM];
int head[MAXN];
int pre[MAXN];
int cur[MAXN];
int level[MAXN];
int gap[MAXN];
int n,m,d,vs,vt,NE,NV;

void Insert(int u,int v,int cap,int cc=0){
   edge[NE].v=v;edge[NE].cap=cap;
   edge[NE].next=head[u];head[u]=NE++;

   edge[NE].v=u;edge[NE].cap=cc;
   edge[NE].next=head[v];head[v]=NE++;
}

int sap(int vs,int vt){
   memset(pre,-1,sizeof(pre));
   memset(level,0,sizeof(level));
   memset(gap,0,sizeof(gap));
   memcpy(cur,head,sizeof(head));
   int u=pre[vs]=vs,aug=inf,maxflow=0;
   gap[0]=NV;
   while(level[vs]<NV){
      loop:
      for(int &i=cur[u];i!=-1;i=edge[i].next){
         int v=edge[i].v;
         if(edge[i].cap&&level[u]==level[v]+1){
            aug=min(aug,edge[i].cap);
            pre[v]=u;
            u=v;
            if(v==vt){
               maxflow+=aug;
               for(u=pre[v];v!=vs;v=u,u=pre[u]){
                  edge[cur[u]].cap-=aug;
                  edge[cur[u]^1].cap+=aug;
               }
               aug=inf;
            }
            goto loop;
         }
      }
      int minlevel=NV;
      for(int i=head[u];i!=-1;i=edge[i].next){
         int v=edge[i].v;
         if(edge[i].cap&&level[v]<minlevel){
            minlevel=level[v];
            cur[u]=i;
         }
      }
      if(--gap[level[u]]==0)break;
      level[u]=minlevel+1;
      gap[level[u]]++;
      u=pre[u];
   }
   return maxflow;
}

char G1[22][22];
char G2[22][22];
int map[22][22];
int mat[22][22];

int main(){
 //  freopen("1.txt","r",stdin);
   int _case,t=1,total,sum;
   scanf("%d",&_case);
   while(_case--){
      scanf("%d%d",&n,&d);
      memset(map,0,sizeof(map));
      memset(mat,0,sizeof(mat));
      memset(head,-1,sizeof(head));
      for(int i=0;i<n;i++)scanf("%s",G1[i]);
      for(int i=0;i<n;i++)scanf("%s",G2[i]);
      m=strlen(G1[0]),NE=0,vs=0,total=0,sum=0;
      for(int i=0;i<n;i++){
         for(int j=0;j<m;j++){
            map[i][j]=G1[i][j]-'0';
            if(map[i][j]){ mat[i][j]=++total;Insert(2*total-1,2*total,map[i][j]); }
         }
      }
      vt=2*total+1,NV=vt+1;
      for(int i=0;i<n;i++){
         for(int j=0;j<m;j++){
            if(G2[i][j]=='L'){ sum++;Insert(vs,2*mat[i][j]-1,1); }
         }
      }
      for(int i=0;i<n;i++){
         for(int j=0;j<m;j++)if(mat[i][j]){
            for(int x=-d;x<=d;x++){
               for(int y=abs(x)-d;y<=d-abs(x);y++){
                  int ii=i+x,jj=j+y;
                  if(ii<0||ii>=n||jj<0||jj>=m)continue;
                  else if(ii==i&&jj==j)continue;
                  else if(mat[ii][jj]==0)continue;
                  Insert(2*mat[i][j],2*mat[ii][jj]-1,inf);
               }
            }
            if(i<d||j<d||n-i<=d||m-j<=d)Insert(2*mat[i][j],vt,inf);
         }
      }
      int ans=sum-sap(vs,vt);
      printf("Case #%d: ",t++);
      if(ans==0){
         printf("no lizard was left behind.\n");
      }else if(ans==1){
         printf("1 lizard was left behind.\n");
      }else {
         printf("%d lizards were left behind.\n",ans);
      }
   }
   return 0;
}