BZOJ3894: 文理分科

继续填坑..

经典最小割

对每个点拆三个点a,b,c，分别代表周围是否全文，这个点的文理，周围是否全理
先将所有价值加上
然后每个点的b和周围4个点的a,c有最小割的那个二元关系什么的建边，转一个最小割的模型，价值和-最小割就是最大价值

建边那里，和源点S连代表选文，和汇点T连代表选理
于是连边(S,b,选理的花费)，(b,T,选文的花费)，(S,a,0)，(a,T,周围全文的花费）,c点的连边和a点相反

然后对于二元关系，比如b和一个相邻点的a，不能这个点选理，那个a又选是全文，于是连一条边(b,a,inf)，b和c的类似

然后好像就搞定了

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
#define inf 1e9
using namespace std;

const int dx[] = {0,-1,0,1,0};
const int dy[] = {0,0,1,0,-1};
const int maxn = 310000;
const int maxm = maxn<<2;

int n,m,st,ed;
struct edge{int y,c,nex;}a[maxm]; int len,fir[maxn],fi[maxn];
inline void ins(const int x,const int y,const int c)
{
    a[++len]=(edge){y,c,fir[x]};fir[x]=len;
    a[++len]=(edge){x,0,fir[y]};fir[y]=len;
}

int h[maxn];
queue<int>q;
bool bfs()
{
    for(int i=1;i<=ed;i++) h[i]=0;
    h[st]=1; q.push(st);
    while(!q.empty())
    {
        const int x=q.front(); q.pop();
        for(int k=fir[x],y=a[k].y;k;k=a[k].nex,y=a[k].y) if(a[k].c&&!h[y])
            h[y]=h[x]+1,q.push(y);
    }
    return h[ed]>0;
}
int dfs(int x,int flow)
{
    if(x==ed) return flow;
    int delta=0;
    for(int &k=fi[x],y=a[k].y;k;k=a[k].nex,y=a[k].y)
    {
        if(a[k].c&&h[y]==h[x]+1)
        {
            int minc=dfs(y,min(a[k].c,flow-delta));
            a[k].c-=minc,a[k^1].c+=minc;
            delta+=minc;
        }
        if(delta==flow) return delta;
    }
    return delta;
}
int Flow()
{
    int re=0;
    while(bfs()) 
    {
        for(int i=1;i<=ed;i++) fi[i]=fir[i];
        re+=dfs(st,INT_MAX);
    }
    return re;
}
int ans;

int main()
{
    len=1;

    scanf("%d%d",&n,&m); int N=n*m; st=N*3+1,ed=st+1;
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
    {
        int id=(i-1)*m+j;
        int x; scanf("%d",&x);
        ins(st,id,x); ans+=x;
    }
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
    {
        int id=(i-1)*m+j;
        int x; scanf("%d",&x);
        ins(id,ed,x); ans+=x;
    }
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
    {
        int id=(i-1)*m+j;
        int x; scanf("%d",&x);
        ins(st,id+N,x),ans+=x;
    }
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
    {
        int id=(i-1)*m+j;
        int x; scanf("%d",&x);
        ins(id+2*N,ed,x),ans+=x;
    }

    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
    {
        int idx=(i-1)*m+j;
        for(int k=0;k<5;k++)
        {
            int x=i+dx[k],y=j+dy[k];
            if(x&&y&&x<=n&&y<=m)
            {
                int idy=(x-1)*m+y;
                ins(idx+N,idy,inf);
                ins(idy,idx+2*N,inf);
            }
        }
    }
    printf("%d\n",ans-Flow());

    return 0;
}