【网络流第四弹】最大点权独立集 ——HDU 1565 方格取数（1）

题目：点击打开链接

网络流学习的下一步，最大点权独立集。多和最小点权覆盖集放到一起使用。分别求图中不相邻的点权最大/或者边中至少有一点在集合中，求最小总权值的问题。

公式：

最大点权独立集=总权值-最小点权覆盖集。

最小点权覆盖集=图的最小割值=最大流。和普通网络流一样，可以使用ISAP,EK,DINIC等多种方法。

建图：

多采用奇偶建立二分图的方法。建边时脑子一定要清楚。相邻点要染不同的颜色，源点为0，汇点为最后的那个点。染色方法基本上是判断if(边+边)%2==0源点与其连接，反之汇点与其连接。如果两点相邻的话，连接值为INF，无穷大。

脑子一定要清醒。。wa了若干次，用的是最快的sap.

#include <iostream>#include <cstring>#include <string>#include <iomanip>using namespace std;typedef  struct {int v,next,val;} edge;const int INF=0x3F3F3F3F; const int MAXN=2000; const int MAXM=50000; edge e[MAXM]; int p[MAXN],eid; inline void init(){memset(p,-1,sizeof(p));eid=0;}  //有向 inline void insert1(int from,int to,int val) {     e[eid].v=to;     e[eid].val=val;     e[eid].next=p[from];     p[from]=eid++;     swap(from,to);     e[eid].v=to;     e[eid].val=0;     e[eid].next=p[from];     p[from]=eid++; }  //无向 inline void insert2(int from,int to,int val) {     e[eid].v=to;     e[eid].val=val;     e[eid].next=p[from];     p[from]=eid++;     swap(from,to);     e[eid].v=to;     e[eid].val=val;     e[eid].next=p[from];     p[from]=eid++; } int n,m;//n为点数 m为边数 int h[MAXN]; int gap[MAXN]; int source,sink; inline int dfs(int pos,int cost) {     if (pos==sink)     {         return cost;     }     int j,minh=n-1,lv=cost,d;     for (j=p[pos];j!=-1;j=e[j].next)     {         int v=e[j].v,val=e[j].val;         if(val>0)         {             if (h[v]+1==h[pos])             {                 if (lv<e[j].val) d=lv;                 else d=e[j].val;                                  d=dfs(v,d);                 e[j].val-=d;                 e[j^1].val+=d;                 lv-=d;                 if (h[source]>=n) return cost-lv;                 if (lv==0) break;             }             if (h[v]<minh)    minh=h[v];         }     }     if (lv==cost)     {         --gap[h[pos]];         if (gap[h[pos]]==0) h[source]=n;         h[pos]=minh+1;         ++gap[h[pos]];     }     return cost-lv; } int sap(int st,int ed) {     source=st;     sink=ed;     int ret=0;     memset(gap,0,sizeof(gap));     memset(h,0,sizeof(h));     gap[st]=n;     while (h[st]<n)     {         ret+=dfs(st,INF);     }     return ret; } int main() {      int pt,ma;     while(cin>>pt)     {         int tmp;          init();          n=2000;           m=10000;           int sum=0; //最大点权独立集=总权值-最小点权覆盖集=总权值-最小割          int st=0;  //start point          int ed=(pt*pt)+1;          for(int i=1;i<=pt;i++)          {              for(int j=1;j<=pt;j++)              {                  scanf("%d",&tmp);                  sum+=tmp;                  if((i+j)%2==0)  //偶数序号点                {                    insert1(st,(i-1)*pt+j,tmp);                    if(i>1)                    insert1((i-1)*pt+j,(i-2)*pt+j,INF);                    if(j>1)                    insert1((i-1)*pt+j,(i-1)*pt+j-1,INF);  //顺序一定不能错                     if(i<pt)                    insert1((i-1)*pt+j,(i)*pt+j,INF);                    if(j<pt)                    insert1((i-1)*pt+j,(i-1)*pt+j+1,INF);                }                else                    insert1((i-1)*pt+j,ed,tmp);                              }                      }          int res=sap(st,ed);          cout<<sum-res<<endl;      }      return 0; }