bzoj2127: happiness 最小割

Description
高一一班的座位表是个n*m的矩阵，经过一个学期的相处，每个同学和前后左右相邻的同学互相成为了好朋友。这学期要分文理科了，每个同学对于选择文科与理科有着自己的喜悦值，而一对好朋友如果能同时选文科或者理科，那么他们又将收获一些喜悦值。作为计算机竞赛教练的scp大老板，想知道如何分配可以使得全班的喜悦值总和最大。
Input

第一行两个正整数n，m。接下来是六个矩阵

第一个矩阵为n行m列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学选择文科获得的喜悦值。

第二个矩阵为n行m列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学选择理科获得的喜悦值。

第三个矩阵为n-1行m列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学与第i+1行第j列的同学同时选择文科获得的额外喜悦值。

第四个矩阵为n-1行m列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学与第i+1行第j列的同学同时选择理科获得的额外喜悦值。

第五个矩阵为n行m-1列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学与第i行第j+1列的同学同时选择文科获得的额外喜悦值。

第六个矩阵为n行m-1列 

此矩阵的第i行第j列的数字表示座位在第i行第j列的同学与第i行第j+1列的同学同时选择理科获得的额外喜悦值。

Output
输出一个整数，表示喜悦值总和的最大值
Sample Input
1 2
1 1
100 110
1
1000 
Sample Output

1210

第一次写的时候顺利推出公式，然后不知道怎么搞得一直就在wawawaawa tetetetetete

然后我就弃疗了。。。。然后一直在想公式法为什么被叉掉了，壮哉我公式法！！！！！！！！！两个星期后，闲来无事的我在早上一看11点40了，20分钟也没有什么事做，所以就重新推了一下公式。然后就吃饭去了。下午睡醒正当我抱着查一下午的决心去测得时候竟然过了。。。。。再次壮哉我公式法！！！！！人弱没办法，弃疗大法好。。。。

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <queue>
using namespace std;
#define INF 0x3f3f3f3f
#define maxn 444444
int st,ed,en;
struct node
{
    int v,next,c;
}e[maxn];
int dis[maxn],first[maxn];
int n,m;
int q[maxn];
bool bfs()
{
    memset(dis,-1,sizeof(dis));
    dis[st]=0;
    int tail,head;
    tail=head=1;
    tail++;
    q[tail]=st;
    int v,u;
    while(head<tail)
    {
        head++;
        u=q[head];
        for(int i=first[u];i!=-1;i=e[i].next)
        {
            v=e[i].v;
            if(dis[v]!=-1) continue;
            if(e[i].c==0) continue;
            dis[v]=dis[u]+1;
            tail++;
            q[tail]=v;
            if(v==ed) return true;
        }
    }
    return false;
}
int dfs(int x,int mx)
{
    if(x==ed||mx==0) return mx;
    int f,flow=0,v,ret=0;
    for(int i=first[x];i!=-1;i=e[i].next)
    {
        v=e[i].v;
        if(dis[x]+1!=dis[v]) continue;
        if((f=dfs(v,min(mx,e[i].c))))
        {
            e[i].c-=f;
            e[i^1].c+=f;
            flow+=f;
            ret+=f;
            mx-=f;
            if(!mx) break;
        }
    }
    if(ret==0) dis[x]=-1;
    return flow;
}
int dinic()
{
    int tmp=0,maxflow=0;
    while(bfs())
    {
        //for(int i=1;i<=ed;i++) cur[i]=first[i];
        while(tmp=dfs(st,INF)) maxflow+=tmp;
    }
    return maxflow;
}
void add(int a,int b,int c)
{
    //printf("%d %d %d\n",a,b,c);
    e[en].v=b;
    e[en].next=first[a];
    e[en].c=c;
    first[a]=en;
    en++;
 
    e[en].v=a;
    e[en].next=first[b];
    e[en].c=0;
    first[b]=en;
    en++;
}
void init()
{
    en=0;
    memset(first,-1,sizeof(first));
    st=n*m*5+1;
    ed=st+1;
}
int ans;
int getnum(int a,int b,int c)
{
    return (c-1)*(n*m)+m*(a-1)+b;
}
int main()
{
    //freopen("happiness3.in","r",stdin);
    //freopen("happiness.out","w",stdout);
    scanf("%d%d",&n,&m);
    init();
    int a,b,c;
    //xuewen
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(st,getnum(i,j,1),a);
        }
    }
    //xueli
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(getnum(i,j,1),ed,a);
        }
    }
    /
    //tl xuewen
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(st,getnum(i,j,2),a);
            add(getnum(i,j,2),getnum(i,j,1),INF);
            add(getnum(i,j,2),getnum(i+1,j,1),INF);
        }
    }
    //tl xueli
    for(int i=1;i<n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(getnum(i,j,3),ed,a);
            add(getnum(i,j,1),getnum(i,j,3),INF);
            add(getnum(i+1,j,1),getnum(i,j,3),INF);
        }
    }
    //tx xuewen
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(st,getnum(i,j,4),a);
            add(getnum(i,j,4),getnum(i,j,1),INF);
            add(getnum(i,j,4),getnum(i,j+1,1),INF);
        }
    }
    //tx xueli
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<m;j++)
        {
            scanf("%d",&a);
            ans+=a;
            add(getnum(i,j,5),ed,a);
            add(getnum(i,j,1),getnum(i,j,5),INF);
            add(getnum(i,j+1,1),getnum(i,j,5),INF);
        }
    }
    printf("%d\n",ans-dinic());
    return 0;
}
/*
2 2
20 20
20 20
 
10 10
10 10
 
0 0
10 10
0 0
10 10
*/