HDU 1569 - 方格取数(2) 二分图最大点权独立集(构图最大流解)

             题意:

                      给你一个m*n的格子的棋盘，每个格子里面有一个非负数。从中取出若干个数，使得任意的两个数所在的格子没有公共边，就是说所取数所在的2个格子不能相邻，并且取出的数的和最大。

             题解:

                      首先抽象成二分图..由于是上下左右四个方向相邻才干扰...可以把所有的点分成奇数点和偶数点(x+y为奇or偶)...把相互冲突的做边..这个模型很像二分图的最大独立集了...二分图的最大独立集=总的点数-二分图的最小覆盖.....但本题是要求没有关系的点的权值最大...也就是二分图的最大点权独立集

                      解法是构图做最大流...把点分成两部分以后..再有一个超级源点和超级汇点..源点向所有左侧的点做边...容量为左侧对应点的权值...左侧点对所有其冲突的右侧点左边.容量为无穷大..右侧点往超级汇点做边..容量为对应右侧点的权值...跑最大流..得到的是最小点权覆盖...那么二分图最大点权独立集=所有点权之和-二分图最小点权覆盖..

Program:

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 3005
#define MAXM 500005
#define oo 1000000007
#define ll long long
using namespace std;
struct Dinic
{
      struct node 
      {
            int x,y,c,next;
      }line[MAXM];   
      int Lnum,_next[MAXN],dis[MAXN];
      void initial(int n) 
      {
             for (int i=0;i<=n;i++) _next[i]=-1;
             Lnum=-1;
      } 
      void addline(int x,int y,int c)
      {
             line[++Lnum].next=_next[x],_next[x]=Lnum;
             line[Lnum].x=x,line[Lnum].y=y,line[Lnum].c=c;
             line[++Lnum].next=_next[y],_next[y]=Lnum;
             line[Lnum].x=y,line[Lnum].y=x,line[Lnum].c=0;             
      }
      bool BFS(int s,int e)
      { 
             queue<int> Q;
             while (!Q.empty()) Q.pop();
             memset(dis,0,sizeof(dis));
             dis[s]=1;
             Q.push(s);
             while (!Q.empty())
             {
                   int h,k;
                   h=Q.front(),Q.pop();
                   if (h==e) return dis[e];
                   for (k=_next[h];k!=-1;k=line[k].next)
                      if (line[k].c && !dis[line[k].y])
                         dis[line[k].y]=dis[h]+1,Q.push(line[k].y);                 
             } 
             return false;
      }
      int dfs(int x,int flow,int e)  
      {     
             if (x==e) return flow;   
             int temp,cost=0;  
             for (int k=_next[x];k!=-1;k=line[k].next)  
             if (line[k].c && dis[line[k].y]==dis[x]+1)  
             {  
                    temp=dfs(line[k].y,min(flow-cost,line[k].c),e);   
                    if (temp)  
                    {  
                           line[k].c-=temp,line[k^1].c+=temp;  
                           cost+=temp;  
                           if (flow==cost) return cost;  
                    }else dis[line[k].y]=-1;  
             }  
             return cost;  
      }  
      int MaxFlow(int s,int e)
      {
             int MaxFlow=0;
             while (BFS(s,e)) 
                MaxFlow+=dfs(s,oo,e); 
             return MaxFlow;
      }
}T;
int main()  
{         
       int R,C,i,j,s,e,v,sum;    
       while (~scanf("%d%d",&R,&C))  
       {  
               T.initial(R*C+5);
               s=R*C+1,e=R*C+2,sum=0;
               for (i=0;i<R;i++)
                   for (j=0;j<C;j++)
                   {
                           scanf("%d",&v),sum+=v;
                           if ((i+j)%2) 
                           {
                                   T.addline(s,i*C+j,v);
                                   if (i) T.addline(i*C+j,(i-1)*C+j,oo);
                                   if (j) T.addline(i*C+j,i*C+j-1,oo);
                                   if (i!=R-1) T.addline(i*C+j,(i+1)*C+j,oo);
                                   if (j!=C-1) T.addline(i*C+j,i*C+j+1,oo);
                           }
                           else 
                                   T.addline(i*C+j,e,v);
                   }  
               printf("%d\n",sum-T.MaxFlow(s,e));
       }  
       return 0;  
}