[网络流24题]方格取数问题（最小割）

方格取数问题

这题和骑士共存问题差不多

1、对这张图进行黑白染色，黑色的放在一个集合，白色放在一个集合，不能共同取的一定是颜色不同的，即当前取的点是黑色，那么它上下左右的白色不能取

2、让黑色点作为二分图的左部，白色点作为二分图的右部，将每个黑色点与它上下左右的白色点建边，即有取了该点，就不能取的点建边，题中是点权，我们把它化为边权，即连的边的容量为该黑色点的点权。

3、对于这个二分图来说，问题等价于最少删除多少边，能让最终的流量为0。所以对有冲突的点进行建边，这样如果一个点要选，那么和他相连的点就要割开，当前建立的图是对冲突进行建立，这不就是最小割吗，而最小割等于最大流。

4、建立一个源点0，将0与所有黑点连边，边权为黑色点的点权，建立一个汇点n * m + 1，将所有白色点与汇点连边，边权为白色点的点权。

#pragma GCC optimize(3)
#include<bits/stdc++.h>
using namespace std;
const int N=55*55;
struct E{
    int to,next,v;
}e[4*N];
int dx[]={1,-1,0,0},dy[]={0,0,-1,1};
int h[N],tot;
int n,m,k;
int ans;
int deep[N],cur[N];
void add(int a,int b,int v){
    e[tot]={b,h[a],v},h[a]=tot++;
    e[tot]={a,h[b],0},h[b]=tot++;
}
int bfs(int s,int t){
      memset(deep,0,sizeof deep);
      queue<int> que;
      deep[s]=1;
      que.push(s);
      while(!que.empty()){
          int x=que.front();
          que.pop();
          for(int i=h[x];~i;i=e[i].next){
              int to=e[i].to,v=e[i].v;
              if(v && !deep[to]){
                  deep[to]=deep[x]+1;
                  que.push(to);
              }
          }
      }
      return deep[t];
}
int dfs(int s,int t,int flow){
    if(s==t) return flow;
    int sum=0;
    for(int &i=cur[s];~i;i=e[i].next){
        int to=e[i].to,v=e[i].v;
        if(v && deep[to]==deep[s]+1){
            int Next=dfs(to,t,min(flow,v));
            e[i].v-=Next,e[i^1].v+=Next;
            flow-=Next,sum+=Next;
        }
    }
    if(!sum) deep[s]=-2;
    return sum;
}
void Dinic(int s,int t){
    int INF=1e9;
    ans=0;
    while(bfs(s,t)) {
        for(int i=0;i<=n*m+1;i++) cur[i]=h[i];
        ans+=dfs(s,t,INF);
    }
}
int main(){
    
    memset(h,-1,sizeof h);
    tot=0;
    int sum=0;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
            int q;scanf("%d",&q);
            sum+=q;
            if((i+j)&1){
                add(0,(i-1)*m+j,q);
                for(int k=0;k<4;k++){
                    int x=i+dx[k],y=j+dy[k];
                    if(x>=1 && x<=n && y>=1 &&y<=m){
                        add((i-1)*m+j,(x-1)*m+y,q);
                    }
                }
            }
            else add((i-1)*m+j,n*m+1,q);
        }
	}  
	
    Dinic(0,n*m+1);
    printf("%d\n",sum-ans);
		
    return 0;
}