BZOJ 2127 二元关系

题意：

思路：

先把所有的值加起来

最小割割哪儿 就代表那个地方不选
一减 剩下的就是 最大值了

//By SiriusRen#include <cstdio>#include <cstring>#include <algorithm>using namespace std;const int N=105,M=10005,K=300050,inf=0x3f3f3f3f;int n,m,T,cnt,tot,ans;int A[N][N],B[N][N],C[N][N],D[N][N],E[N][N],F[N][N],id[N][N];int first[M],next[K],v[K],w[K],q[M],vis[M];void Add(int x,int y,int z){w[tot]=z,v[tot]=y,next[tot]=first[x],first[x]=tot++;}void add(int x,int y,int z){Add(x,y,z),Add(y,x,0);}bool tell(){    memset(vis,-1,sizeof(vis)),vis[0]=0;    int head=0,tail=0;    while(head<=tail){        int t=q[head++];        for(int i=first[t];~i;i=next[i])if(!~vis[v[i]]&&w[i])            vis[v[i]]=vis[t]+1,q[++tail]=v[i];    }return ~vis[T];}int zeng(int x,int y){    if(x==T)return y;    int r=0;    for(int i=first[x];~i&&y>r;i=next[i])if(vis[v[i]]==vis[x]+1&&w[i]){        int t=zeng(v[i],min(y-r,w[i]));        w[i]-=t,w[i^1]+=t,r+=t;    }if(!r)vis[x]=-1;    return r;}int main(){    memset(first,-1,sizeof(first));    scanf("%d%d",&n,&m),T=n*m+1;    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)scanf("%d",&A[i][j]),id[i][j]=++cnt;    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)scanf("%d",&B[i][j]);    for(int i=1;i<n;i++)for(int j=1;j<=m;j++)scanf("%d",&C[i][j]);    for(int i=1;i<n;i++)for(int j=1;j<=m;j++)scanf("%d",&D[i][j]);    for(int i=1;i<=n;i++)for(int j=1;j<m;j++)scanf("%d",&E[i][j]);    for(int i=1;i<=n;i++)for(int j=1;j<m;j++)scanf("%d",&F[i][j]);    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)        add(0,id[i][j],A[i][j]*2),add(id[i][j],T,B[i][j]*2),ans+=A[i][j]+B[i][j];    for(int i=1;i<n;i++)for(int j=1;j<=m;j++){        add(0,id[i][j],C[i][j]),add(0,id[i+1][j],C[i][j]);        add(id[i][j],T,D[i][j]),add(id[i+1][j],T,D[i][j]);        add(id[i][j],id[i+1][j],C[i][j]+D[i][j]);        add(id[i+1][j],id[i][j],C[i][j]+D[i][j]);        ans+=C[i][j]+D[i][j];    }    for(int i=1;i<=n;i++)for(int j=1;j<m;j++){        add(0,id[i][j],E[i][j]),add(0,id[i][j+1],E[i][j]);        add(id[i][j],T,F[i][j]),add(id[i][j+1],T,F[i][j]);        add(id[i][j],id[i][j+1],E[i][j]+F[i][j]);        add(id[i][j+1],id[i][j],E[i][j]+F[i][j]);        ans+=E[i][j]+F[i][j];    }ans*=2;    while(tell())while(m=zeng(0,inf))ans-=m;    printf("%d\n",ans/2);}