swust1744: 方格取数问题

链接：http://www.oj.swust.edu.cn/problem/show/1744

题意：中文题。

分析：要求选出的数总和最大。我们将一行的格子黑白染色，然后将所有黑格子与源相连边的容量为他们的值，所有白格子与汇相连边的容量为他们的值，对于每个黑格子与相邻的白格子练一条容量为无穷大的边。这样我们就建成了一个找二分图最小割（最大流）的图，然后我们将所有格子的总和减去这个图跑出的最大流就是我们要的答案啦，我们可以将这个模型看成是求二分图求最大独立集来理解。

代码：

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<sstream>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=205;
const int MAX=1000000100;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=998244353;
const int INF=1000000010;
const double pi=acos(-1.0);
typedef double db;
typedef unsigned long long ull;
int tot,u[N*N],v[4*N*N],cap[4*N*N],flow[4*N*N],pre[4*N*N];
void add(int a,int b,int c,int f) {
    v[tot]=b;cap[tot]=c;flow[tot]=f;pre[tot]=u[a];u[a]=tot++;
}
bool q[N*N];
int d[N*N],dis[N*N],head[N*N];
bool bfs(int a,int b) {
    int i,l=1,r=1;
    memset(q,0,sizeof(q));
    d[1]=a;q[a]=1;dis[a]=0;
    for (;l<=r;l++)
        for (i=u[d[l]];i!=-1;i=pre[i])
        if (!q[v[i]]&&cap[i]>flow[i]) q[v[i]]=1,d[++r]=v[i],dis[v[i]]=dis[d[l]]+1;
    return q[b];
}
int dfs(int a,int b,int en) {
    if (a==en||b==0) return b;
    int ret=0,f;
    for (int& i=u[a];i!=-1;i=pre[i])
    if (dis[v[i]]==dis[a]+1&&(f=dfs(v[i],min(b,cap[i]-flow[i]),en))>0) {
        flow[i]+=f;flow[i^1]-=f;
        ret+=f;b-=f;
        if (b==0) break ;
    }
    return ret;
}
int main()
{
    int a,i,j,n,m,S,T,D,sum=0,ans=0;
    scanf("%d%d", &n, &m);
    S=0;T=n*m+1;
    tot=0;memset(u,-1,sizeof(u));
    for (i=1;i<=n;i++)
        for (j=1;j<=m;j++) {
            scanf("%d", &a);
            D=(i-1)*m+j;ans+=a;
            if (i%2==j%2) {
                add(S,D,a,0),add(D,S,0,0);
                if (i>1) add(D,D-m,INF,0),add(D-m,D,0,0);
                if (j>1) add(D,D-1,INF,0),add(D-1,D,0,0);
                if (i<n) add(D,D+m,INF,0),add(D+m,D,0,0);
                if (j<m) add(D,D+1,INF,0),add(D+1,D,0,0);
            } else add(D,T,a,0),add(T,D,0,0);
        }
    for (i=S;i<=T;i++) head[i]=u[i];
    while (bfs(S,T)) {
        sum+=dfs(S,INF,T);
        for (i=S;i<=T;i++) u[i]=head[i];
    }
    printf("%d\n", ans-sum);
    return 0;
}