bzoj2132: 圈地计划

网络流好题。。
和happiness很想
然后大概想到等效建图。。
于是我又没有弄出来。。
做那题的时候太久了。。已经哈哈

首先将图黑白染色，使(i+j)&1的为白点X，其他为黑点Y，满足白/黑点集之内无相邻关系。
对于白点，连边(S,X,A)(X,T,B)
对于黑点，连边(S,Y,B)(X,Y,A)
两相邻黑白点之间连无向边(CA+CB) 表示A格C与B格C之和。

这类题的构图方法一点要会。。以后还是看看吧。。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int N=105*105;
const int MAX=1<<30;
int n,m;
int a[105][105],b[105][105],c[105][105];
int S,T;
struct qq
{
    int x,y,z,last;
}s[N*20];int num,last[N];
void init (int x,int y,int z)
{
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;
    s[num].last=last[x];
    last[x]=num;
    swap(x,y);z=0;
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;
    s[num].last=last[x];
    last[x]=num;
}
int ans=0;
int fx[4]={0,0,1,-1};
int fy[4]={1,-1,0,0};
int P (int x,int y)
{
    return (x-1)*m+y;
}
void build ()
{
    num=1;memset(last,-1,sizeof(last));
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
            if ((u+i)%2==0)
                swap(a[u][i],b[u][i]);
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
        {
            init(S,P(u,i),a[u][i]);
            init(P(u,i),T,b[u][i]);
            ans=ans+a[u][i]+b[u][i];
            if ((u+i)%2==0)
            {
                for (int j=0;j<4;j++)
                {
                    int nx=u+fx[j],ny=i+fy[j];
                    if(nx>n||nx<1||ny>m||ny<1) continue;
                    init(P(u,i),P(nx,ny),c[u][i]+c[nx][ny]);
                    init(P(nx,ny),P(u,i),c[u][i]+c[nx][ny]);
                    ans=ans+c[u][i]+c[nx][ny];
                }
            }
        }
}
int h[N];
bool bt ()
{
    memset(h,-1,sizeof(h));h[S]=1;
    queue<int> q;
    q.push(S);
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=s[u].last)
        {
            int y=s[u].y;
            if (s[u].z>0&&h[y]==-1)
            {
                h[y]=h[x]+1;
                q.push(y);
            }
        }
    }
    return h[T]!=-1;
}
int mymin (int x,int y){return x<y?x:y;}
int dfs (int x,int f)
{
    if (x==T) return f;
    int s1=0;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (s[u].z>0&&h[y]==h[x]+1&&s1<f)
        {
            int lalal=dfs(y,mymin(s[u].z,f-s1));
            s1+=lalal;
            s[u].z-=lalal;
            s[u^1].z+=lalal;
        }
    }
    if (s1==0) h[x]=-1;
    return s1;
}
int main()
{
    scanf("%d%d",&n,&m);
    S=n*m+1;T=S+1;
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
            scanf("%d",&a[u][i]);
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
            scanf("%d",&b[u][i]);
    for (int u=1;u<=n;u++)
        for (int i=1;i<=m;i++)
            scanf("%d",&c[u][i]);
    build();
     while (bt()==true) ans=ans-dfs(S,MAX);
     printf("%d\n",ans);
    return 0;
}