HDOJ 3657 - Game 构图最小割

               题意:

                      给了N*M的矩阵..每个单位格子上有数.代表取走这个格子的能获得的分数.但是若相邻的格子取了.就要减去2*(x&y)的分数.x,y分别是这两个格子里的数...现在规定了一些格子必须选.问最后能得到的最大分数..

               题解:

                      该题和HDOJ 3820差不多...只是加了一个有些点必须选..那么把起点到该点 或者 该点到终点的容量设置为无穷大即可..

Program:

#include<iostream>    
#include<algorithm>    
#include<stdio.h>    
#include<string.h>  
#include<time.h> 
#include<map> 
#include<math.h>    
#include<queue>    
#define MAXN 100005 
#define MAXM 2000005   
#define oo 1000000007    
#define ll long long    
using namespace std;   
struct Dinic            
{            
       struct node          
       {           
             int c,u,v,next;          
       }edge[MAXM];          
       int ne,head[MAXN];          
       int cur[MAXN], ps[MAXN], dep[MAXN];          
       void initial()          
       {          
             ne=2;          
             memset(head,0,sizeof(head));           
       }          
       void addedge(int u, int v,int c)          
       {           
             edge[ne].u=u,edge[ne].v=v,edge[ne].c=c,edge[ne].next=head[u];          
             head[u]=ne++;          
             edge[ne].u=v,edge[ne].v=u,edge[ne].c=0,edge[ne].next=head[v];          
             head[v]=ne++;          
       }          
       int MaxFlow(int s,int t)          
       {                               
             int tr, res = 0;          
             int i,j,k,f,r,top;          
             while(1)          
             {          
                    memset(dep, -1, sizeof(dep));          
                    for(f=dep[ps[0]=s]=0,r=1;f!= r;)          
                       for(i=ps[f++],j=head[i];j;j=edge[j].next)          
                         if(edge[j].c&&dep[k=edge[j].v]==-1)          
                         {          
                               dep[k]=dep[i]+1;          
                               ps[r++]=k;          
                               if(k == t){  f=r; break;  }          
                         }          
                    if(dep[t]==-1) break;          
                    memcpy(cur,head,sizeof(cur));          
                    i=s,top=0;          
                    while(1)          
                    {          
                         if(i==t)          
                         {          
                               for(tr=oo,k=0;k<top;k++)          
                                  if(edge[ps[k]].c<tr)          
                                     tr=edge[ps[f=k]].c;          
                               for(k=0;k<top;k++)          
                               {          
                                     edge[ps[k]].c-=tr;          
                                     edge[ps[k]^1].c+=tr;          
                               }          
                               i=edge[ps[top=f]].u;          
                               res+= tr;          
                         }          
                         for(j=cur[i];cur[i];j=cur[i]=edge[cur[i]].next)           
                             if(edge[j].c && dep[i]+1==dep[edge[j].v]) break;           
                         if(cur[i])  ps[top++]=cur[i],i=edge[cur[i]].v;           
                         else          
                         {          
                                 if(!top) break;          
                                 dep[i]=-1;          
                                 i=edge[ps[--top]].u;          
                         }          
                   }          
             }          
             return res;          
      }          
}T;      
int A[55][55];
bool f[55*55];
int main()   
{       
      int n,m,k,s,e,x,h,i,j,sum;      
      while (~scanf("%d%d%d",&n,&m,&k))
      {
               s=n*m+5,e=s+1,T.initial(),sum=0;
               for (i=0;i<n;i++)
                  for (j=0;j<m;j++) 
                      scanf("%d",&A[i][j]),sum+=A[i][j];
               memset(f,false,sizeof(f));
               while (k--) scanf("%d%d",&i,&j),i--,j--,f[i*m+j]=true;
               for (i=0;i<n;i++)
                  for (j=0;j<m;j++)
                  {
                          x=i*m+j;
                          if ((i+j)%2) 
                          {
                                  if (f[x]) T.addedge(s,x,oo);
                                      else  T.addedge(s,x,A[i][j]);
                                  if (i) T.addedge(x,x-m,2*(A[i][j]&A[i-1][j]));
                                  if (j) T.addedge(x,x-1,2*(A[i][j]&A[i][j-1]));
                                  if (i!=n-1) T.addedge(x,x+m,2*(A[i][j]&A[i+1][j])); 
                                  if (j!=m-1) T.addedge(x,x+1,2*(A[i][j]&A[i][j+1]));
                          }else 
                                  if (f[x]) T.addedge(x,e,oo);
                                      else  T.addedge(x,e,A[i][j]);
                  }               
               printf("%d\n",sum-T.MaxFlow(s,e));
      }
      return 0;  
}