[bzoj2127] happiness 最小割

2127: 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

【样例说明】

两人都选理，则获得100+110+1000的喜悦值。

【数据规模】

对于100%以内的数据，n,m<=100 所有喜悦值均为小于等于5000的非负整数

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int INF = 1000000000 + 7 ;
const int N = 100 + 5;
int n,m,x,a[N][N],b[N][N],num[N][N],tot=0,cnt=1,ans;
int head[N*N],h[N*N],q[N*N],cur[N*N];
int T;
struct edge{
    int to,v,next;
}e[N*N*N];
void insert( int u, int v, int w ){ e[++cnt].to = v; e[cnt].v = w; e[cnt].next = head[u]; head[u] = cnt;}
bool bfs(){
    int hed = 0, tail = 1;
    memset(h,-1,sizeof(h));
    q[0] = h[0] = 0;
    while( hed != tail ){
        int u = q[hed++];
        for( int i = head[u]; i; i = e[i].next ){
            int v = e[i].to;
            if( h[v] == -1 && e[i].v ) {
                q[tail++] = v;
                h[v] = h[u] + 1;
            }
        }
    }
    return h[T] != -1;
}
int dfs( int x, int f ){
    if( x == T ) return f;
    int w,used=0;
    for( int i = cur[x]; i; i = e[i].next ){
        int v = e[i].to;
        if( h[v] == h[x] + 1 ){
            w = dfs( v, min( f-used, e[i].v ) );
            used += w; e[i].v -= w; e[i^1].v += w;
            if( e[i].v ) cur[x] = i; 
            if( used == f ) return f;
        }
    }
    if( !used ) h[x] = -1;
    return used;
}
void dinic(){
    while( bfs() ){
        for( int i = 0; i <= T; i++ ) cur[i] = head[i]; 
        ans += dfs( 0, INF ); 
    }
}
int main(){
    scanf("%d%d", &n, &m);
    T = n*m+1;
    for( int i = 1; i <= n; i++ ) for( int j = 1; j <= m; j++ ) scanf("%d", &a[i][j]),tot+=a[i][j],a[i][j]<<=1;
    for( int i = 1; i <= n; i++ ) for( int j = 1; j <= m; j++ ) scanf("%d", &b[i][j]),tot+=b[i][j],b[i][j]<<=1;
    for( int i = 1; i <= n; i++ ) for( int j = 1; j <= m; j++ ) num[i][j] = j+m*(i-1);
    for( int i = 1; i < n; i++ ) for( int j = 1; j <= m; j++ ){
        scanf("%d", &x); a[i][j] += x; a[i+1][j] += x; tot += x;
        insert( num[i][j], num[i+1][j], x ); insert( num[i+1][j], num[i][j], x );
    }
    for( int i = 1; i < n; i++ ) for( int j = 1; j <= m; j++ ){
        scanf("%d", &x); b[i][j] += x; b[i+1][j] += x; tot += x;
        insert( num[i][j], num[i+1][j], x ); insert( num[i+1][j], num[i][j], x );
    }
    for( int i = 1; i <= n; i++ ) for( int j = 1; j < m; j++ ){
        scanf("%d", &x); a[i][j] += x; a[i][j+1] += x; tot += x;
        insert( num[i][j], num[i][j+1], x ); insert( num[i][j+1], num[i][j], x );
    }
    for( int i = 1; i <= n; i++ ) for( int j = 1; j < m; j++ ){
        scanf("%d", &x); b[i][j] += x; b[i][j+1] += x; tot += x;
        insert( num[i][j], num[i][j+1], x ); insert( num[i][j+1], num[i][j], x );
    }
    for( int i = 1; i <= n; i++ ) for( int j = 1; j <= m; j++ ){
        insert( 0, num[i][j], a[i][j] ); insert( num[i][j], 0, 0 );
        insert( num[i][j], T, b[i][j] ); insert( T, num[i][j], 0 );
    }
    dinic();
    printf("%d", tot-(ans>>1) );
    return 0;
}