bzoj3894 文理分科

Description 文理分科是一件很纠结的事情！（虽然看到这个题目的人肯定都没有纠 结过）
小P所在的班级要进行文理分科。他的班级可以用一个n*m的矩阵进行 描述，每个格子代表一个同学的座位。每位同学必须从文科和理科中选择
一科。同学们在选择科目的时候会获得一个满意值。满意值按如下的方式 得到：
1．如果第i行第秒J的同学选择了文科，则他将获得art[i][j]的满意值，如 果选择理科，将得到science[i][j]的满意值。
2．如果第i行第J列的同学选择了文科，并且他相邻（两个格子相邻当且 仅当它们拥有一条相同的边）的同学全部选择了文科，则他会更开
心，所以会增加same_art[i][j]的满意值。 3．如果第i行第j列的同学选择了理科，并且他相邻的同学全部选择了理
科，则增加same_science[i]j[]的满意值。 小P想知道，大家应该如何选择，才能使所有人的满意值之和最大。请
告诉他这个最大值。

Input 第一行为两个正整数：n，m 接下来n术m个整数，表示art[i][j]； 接下来n术m个整数．表示science[i][j]；
接下来n术m个整数，表示same_art[i][j]；

Output 输出为一个整数，表示最大的满意值之和

先把所有收益累加，每个点往左右连边代表放弃的代价。每五个点单独拉出来一个点代表放弃全选。求最小割。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int s=40005,t=40006,oo=0x3f3f3f3f; 
int fir[40010],ne[2000010],to[2000010],w[2000010],que[40010],dep[40010],
a[110][110],b[110][110],id[110][110],
xx[4]={1,-1,0,0},yy[4]={0,0,1,-1},n,m,num;
void add(int u,int v,int x)
{
    num++;
    ne[num*2]=fir[u];
    fir[u]=num*2;
    to[num*2]=v;
    w[num*2]=x;
    ne[num*2+1]=fir[v];
    fir[v]=num*2+1;
    to[num*2+1]=u;
    w[num*2+1]=0;
}
bool ok(int x,int y)
{
    return x>=1&&y>=1&&x<=n&&y<=m;
}
bool bfs()
{
    int hd=1,tl=1,u,v;
    memset(dep,0,sizeof(dep));
    que[1]=s;
    dep[s]=1;
    while (hd<=tl)
    {
        u=que[hd++];
        for (int i=fir[u];i;i=ne[i])
            if (!dep[v=to[i]]&&w[i])
            {
                dep[v]=dep[u]+1;
                que[++tl]=v;
            }
    }
    return dep[t];
}
int dfs(int u,int lim)
{
    if (u==t) return lim;
    int ret=0,x,v;
    for (int i=fir[u];i&&ret<lim;i=ne[i])
        if (w[i]&&dep[v=to[i]]==dep[u]+1)
        {
            x=dfs(v,min(lim-ret,w[i]));
            w[i]-=x;
            w[i^1]+=x;
            ret+=x;
        }
    if (!ret) dep[u]=0;
    return ret;
}
void check()
{
    for (int i=1;i<=t;i++)
        for (int j=fir[i];j;j=ne[j])
            if (w[j])
                printf("%d->%d:%d\n",i,to[j],w[j]);
}
int main()
{
    int ans=0,x1,y1,x=0,y;
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
            id[i][j]=++x;
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&x);
            ans+=x;
            add(id[i][j],t,x);
        }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&x);
            ans+=x;
            add(s,id[i][j],x);
        }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&x);
            ans+=x;
            add(id[i][j]+n*m,t,x);
            add(id[i][j],id[i][j]+n*m,oo);
            for (int k=0;k<4;k++)
                if (ok(x1=i+xx[k],y1=j+yy[k]))
                    add(id[x1][y1],id[i][j]+n*m,oo);
        }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=m;j++)
        {
            scanf("%d",&x);
            ans+=x;
            add(s,id[i][j]+2*n*m,x);
            add(id[i][j]+2*n*m,id[i][j],oo);
            for (int k=0;k<4;k++)
                if (ok(x1=i+xx[k],y1=j+yy[k]))
                    add(id[i][j]+2*n*m,id[x1][y1],oo);
        }
    //check();
    while (bfs()) ans-=dfs(s,oo);
    printf("%d\n",ans);
}