hdu 2732（网络流）
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2732
思路：如果从某一根柱子能够直接跳到迷宫的外面，那么我们就将这个点连接到汇点。对于哪些不能跳出去但是又有柱子的点，那么 我们就去按照跳跃距离搜寻有没有其他的柱子能够去跳跃，如果能够找到的话，那么连接这两点，并且将容量控制为弧尾节点的柱子数，也正是由于一条弧只能够约 束一个顶点，所以我们需要进行拆点，点内之间流量为本身柱子数。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 using namespace std; 6 #define MAXN 2222 7 #define MAXM 2222222 8 #define inf 1<<30 9 struct Edge{ 10 int v,cap,next; 11 }edge[MAXM]; 12 int head[MAXN]; 13 int pre[MAXN]; 14 int cur[MAXN]; 15 int level[MAXN]; 16 int gap[MAXN]; 17 int n,m,d,vs,vt,NE,NV; 18 19 void Insert(int u,int v,int cap,int cc=0){ 20 edge[NE].v=v;edge[NE].cap=cap; 21 edge[NE].next=head[u];head[u]=NE++; 22 23 edge[NE].v=u;edge[NE].cap=cc; 24 edge[NE].next=head[v];head[v]=NE++; 25 } 26 27 int sap(int vs,int vt){ 28 memset(pre,-1,sizeof(pre)); 29 memset(level,0,sizeof(level)); 30 memset(gap,0,sizeof(gap)); 31 memcpy(cur,head,sizeof(head)); 32 int u=pre[vs]=vs,aug=inf,maxflow=0; 33 gap[0]=NV; 34 while(level[vs]<NV){ 35 loop: 36 for(int &i=cur[u];i!=-1;i=edge[i].next){ 37 int v=edge[i].v; 38 if(edge[i].cap&&level[u]==level[v]+1){ 39 aug=min(aug,edge[i].cap); 40 pre[v]=u; 41 u=v; 42 if(v==vt){ 43 maxflow+=aug; 44 for(u=pre[v];v!=vs;v=u,u=pre[u]){ 45 edge[cur[u]].cap-=aug; 46 edge[cur[u]^1].cap+=aug; 47 } 48 aug=inf; 49 } 50 goto loop; 51 } 52 } 53 int minlevel=NV; 54 for(int i=head[u];i!=-1;i=edge[i].next){ 55 int v=edge[i].v; 56 if(edge[i].cap&&level[v]<minlevel){ 57 minlevel=level[v]; 58 cur[u]=i; 59 } 60 } 61 if(--gap[level[u]]==0)break; 62 level[u]=minlevel+1; 63 gap[level[u]]++; 64 u=pre[u]; 65 } 66 return maxflow; 67 } 68 69 char G1[22][22]; 70 char G2[22][22]; 71 int map[22][22]; 72 int mat[22][22]; 73 74 int main(){ 75 // freopen("1.txt","r",stdin); 76 int _case,t=1,total,sum; 77 scanf("%d",&_case); 78 while(_case--){ 79 scanf("%d%d",&n,&d); 80 memset(map,0,sizeof(map)); 81 memset(mat,0,sizeof(mat)); 82 memset(head,-1,sizeof(head)); 83 for(int i=0;i<n;i++)scanf("%s",G1[i]); 84 for(int i=0;i<n;i++)scanf("%s",G2[i]); 85 m=strlen(G1[0]),NE=0,vs=0,total=0,sum=0; 86 for(int i=0;i<n;i++){ 87 for(int j=0;j<m;j++){ 88 map[i][j]=G1[i][j]-'0'; 89 if(map[i][j]){ mat[i][j]=++total;Insert(2*total-1,2*total,map[i][j]); } 90 } 91 } 92 vt=2*total+1,NV=vt+1; 93 for(int i=0;i<n;i++){ 94 for(int j=0;j<m;j++){ 95 if(G2[i][j]=='L'){ sum++;Insert(vs,2*mat[i][j]-1,1); } 96 } 97 } 98 for(int i=0;i<n;i++){ 99 for(int j=0;j<m;j++)if(mat[i][j]){ 100 for(int x=-d;x<=d;x++){ 101 for(int y=abs(x)-d;y<=d-abs(x);y++){ 102 int ii=i+x,jj=j+y; 103 if(ii<0||ii>=n||jj<0||jj>=m)continue; 104 else if(ii==i&&jj==j)continue; 105 else if(mat[ii][jj]==0)continue; 106 Insert(2*mat[i][j],2*mat[ii][jj]-1,inf); 107 } 108 } 109 if(i<d||j<d||n-i<=d||m-j<=d)Insert(2*mat[i][j],vt,inf); 110 } 111 } 112 int ans=sum-sap(vs,vt); 113 printf("Case #%d: ",t++); 114 if(ans==0){ 115 printf("no lizard was left behind.\n"); 116 }else if(ans==1){ 117 printf("1 lizard was left behind.\n"); 118 }else { 119 printf("%d lizards were left behind.\n",ans); 120 } 121 } 122 return 0; 123 }
