BZOJ 2127: happiness

2127: happiness

Time Limit: 51 Sec Memory Limit: 259 MB
Submit: 2349 Solved: 1137
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
【样例说明】
两人都选理，则获得100+110+1000的喜悦值。
【数据规模】
对于100%以内的数据，n,m<=100 所有喜悦值均为小于等于5000的非负整数

#include<queue>
#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 50005
#define MAXM 300000
#define P 105
#define INF 0x3f3f3f
struct Edge{ int val,to,next; }e[MAXM];
int head[MAXN],tot;
inline void Add_Edge(int u,int v,int w){
    e[++tot].to=v;e[tot].val=w;
    e[tot].next=head[u];head[u]=tot;
    e[++tot].to=u;e[tot].val=w;
    e[tot].next=head[v];head[v]=tot;
}
int S,T,dep[MAXN];
queue<int> q;
bool BFS(){
    memset(dep,0,sizeof dep );
    while(!q.empty()) q.pop();
    q.push(S);dep[S]=1;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to;
            if(!dep[v]&&e[i].val){
                dep[v]=dep[u]+1;
                if(v==T) return 1;
                q.push(v);
            }
        }
    }
    return 0;
}
int n,m,sum,ans,a[P][P];
int Dinic(int u,int flow){
    if(u==T) return flow;
    int ret=0;
    for(int i=head[u];i&&(flow-ret);i=e[i].next){
        int v=e[i].to;
        if(dep[v]==dep[u]+1&&e[i].val){
            int x=Dinic(v,min(flow-ret,e[i].val));
            if(!x) dep[v]=0;//没有会 T 
            ret+=x; e[i].val-=x; e[i^1].val+=x;
        }
    }
    return ret;
}
void Build(){
    scanf("%d%d",&n,&m); int temp=n*m;;
    S=n*m+2*(n-1)*m+2*n*(m-1)+1,T=n*m+2*(n-1)*m+2*n*(m-1)+2;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) a[i][j]= ++tot;
    tot=1;
    for(int x,i=1;i<=n;i++)//文科
        for(int j=1;j<=m;j++){
            scanf("%d",&x); sum+=x;
            Add_Edge(S,a[i][j],x);
        }
    for(int x,i=1;i<=n;i++)//理科
        for(int j=1;j<=m;j++){
            scanf("%d",&x); sum+=x;
            Add_Edge(a[i][j],T,x);
        }
    for(int x,i=1;i<n;i++)//和下边的同学同时选文科
        for(int j=1;j<=m;j++){
            scanf("%d",&x); sum+=x; ++temp;
            Add_Edge(a[i][j],temp,x),Add_Edge(a[i+1][j],temp,x);
            Add_Edge(S,temp,x);
        }
    for(int x,i=1;i<n;i++)//和下边的同学同时选理科
        for(int j=1;j<=m;j++){
            scanf("%d",&x); sum+=x; ++temp;
            Add_Edge(a[i][j],temp,x);Add_Edge(a[i+1][j],temp,x);
            Add_Edge(temp,T,x);
        }
    for(int x,i=1;i<=n;i++)//和右边的同学同时选文科
        for(int j=1;j<m;j++){
            scanf("%d",&x); sum+=x; ++temp;
            Add_Edge(a[i][j],temp,x);Add_Edge(a[i][j+1],temp,x);
            Add_Edge(S,temp,x);
        }
    for(int x,i=1;i<=n;i++)//和右边的同学同时选li科
        for(int j=1;j<m;j++){
            scanf("%d",&x); sum+=x; ++temp;
            Add_Edge(a[i][j],temp,x);Add_Edge(a[i][j+1],temp,x);
            Add_Edge(temp,T,x);
        }
}
int main(){
    Build();
    while(BFS()) ans+=Dinic(S,INF);
    printf("%d\n",sum-ans);
    return 0;
}

不懂看这—-> http://blog.csdn.net/vmurder/article/details/42609669