hdu 3657 最小割的活用 / 奇偶方格取数类经典题 /最小割

题意：方格取数，如果取了相邻的数，那么要付出一定代价。（代价为2*(X&Y)）（开始用费用流，敲升级版3820，跪。。。）

建图：  对于相邻问题，经典方法：奇偶建立二分图。对于相邻两点连边2*(X&Y)，源->X连边，Y->汇连边，权值w为点权。

ans=总点权-最小割：如果割边是源->X，表示x不要选（是割边，必然价值在路径上最小），若割边是Y-汇点，同理；若割边是X->Y，则表示选Y点且选X点， 割为w( 2*(X&Y) )。

自己的确还没有理解其本质精妙所在。不知何以然也。（开始多敲多了几个else一直跪!）

#include<iostream>
#include<queue>
#include<cstdio>
#include<cstring>
#include<set>
#include<vector>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxv=2550,maxe=20000;
int nume=0;int head[maxv];int e[maxe][3];
void inline adde(int i,int j,int c)
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume++][2]=c;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume++][2]=0;
}
int ss,tt,n,m,k;
int vis[maxv];int lev[maxv];
bool bfs()
{
    for(int i=0;i<maxv;i++)
      vis[i]=lev[i]=0;
    queue<int>q;
    q.push(ss);
    vis[ss]=1;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        for(int i=head[cur];i!=-1;i=e[i][1])
        {
            int v=e[i][0];
            if(!vis[v]&&e[i][2]>0)
            {
                lev[v]=lev[cur]+1;
                vis[v]=1;
                q.push(v);
            }
        }
    }
    return vis[tt];
}
int dfs(int u,int minf)
{
    if(u==tt||minf==0)return minf;
    int sumf=0,f;
    for(int i=head[u];i!=-1&&minf;i=e[i][1])
    {
        int v=e[i][0];
        if(lev[v]==lev[u]+1&&e[i][2]>0)
        {
            f=dfs(v,minf<e[i][2]?minf:e[i][2]);
            e[i][2]-=f;e[i^1][2]+=f;
            sumf+=f;minf-=f;
        }
    }
    if(!sumf) lev[u]=-1;
    return sumf;
}
int dinic()
{
    int sum=0;
    while(bfs())sum+=dfs(ss,inf);
    return sum;
};
int mapp[52][52];int must[maxv];int sums=0;
void read_build()
{
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
      {
        scanf("%d",&mapp[i][j]);
        sums+=mapp[i][j];
       }
     int aa,bb;
    for(int i=0;i<k;i++)
    {
        scanf("%d%d",&aa,&bb);
       must[(aa-1)*m+bb-1]=1;
    }
    for(int i=0;i<n;i++)
      for(int j=0;j<m;j++)
      {
         if((i+j)%2==0)
         {
             if(must[i*m+j])
              adde(ss,i*m+j,inf);
             else
              adde(ss,i*m+j,mapp[i][j]);
             if(i-1>=0)
                 adde(i*m+j,(i-1)*m+j,2*(mapp[i][j]&mapp[i-1][j]));
              if(i+1<n)
                adde(i*m+j,(i+1)*m+j,2*(mapp[i][j]&mapp[i+1][j]));
              if(j-1>=0)
                adde(i*m+j,i*m+j-1,2*(mapp[i][j]&mapp[i][j-1]));
              if(j+1<m)
                adde(i*m+j,i*m+j+1,2*(mapp[i][j]&mapp[i][j+1]));
         }
         else
         {
              if(must[i*m+j])
              adde(i*m+j,tt,inf);
             else
              adde(i*m+j,tt,mapp[i][j]);
         }
      }
  /*  for(int i=0;i<=tt;i++)
      for(int j=head[i];j!=-1;j=e[j][1])
      {
          if(j%2==0)
          printf("%d->%d:%d\n",i,e[j][0],e[j][2]);
      }*/
}
void init()
{
    nume=0;sums=0;
    ss=n*m+2;tt=ss+1;
    for(int i=0;i<=tt;i++)
      {
          head[i]=-1;
          must[i]=0;
      }
}
int main()
{
      while(scanf("%d%d%d",&n,&m,&k)!=EOF)
     {
         init();
        read_build();
       int ans;
       ans=sums-dinic();
      printf("%d\n",ans);
    }
    return 0;
}