bzoj 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的非负整数

最小割= =写的有点烦躁。。了

S->B:w1/2，B->T:w2/2，A<-->B:(w1+w2)/2。

所有边乘2做防止小数。

最后再/2就好

#include<cstdio>#include<string>#include<cstring>using namespace std;int head[100001];struct map{ int f; int s,t;     int next;}a[400001];int b[5001];int edge;int p;int q[400001],d[400001];inline void add(int s,int t,int f){     a[edge].next=head[s];     head[s]=edge;     a[edge].s=s;     a[edge].t=t;     a[edge].f=f;}inline bool bfs(){     int l=0,r=0;     memset(q,0,sizeof(q));     r++;     q[r]=0;     memset(d,-1,sizeof(d));     d[0]=0;     int i,k;     while(l<r)     {       l++;          int k=q[l];          for(i=head[k];i!=0;i=a[i].next)          {               if(a[i].f>0&&d[a[i].t]==-1)               {                    d[a[i].t]=d[k]+1;                    r++;                    q[r]=a[i].t;               }          }     }     if(d[p]>=0)          return true;     return false;}inline int dfs(int k,int s){     if(k==p)          return s;     int t=s;     int i;     for(i=head[k];i!=0;i=a[i].next)     {          if(d[a[i].t]==d[k]+1&&a[i].f>0)          {               int xx=dfs(a[i].t,min(s,a[i].f));               a[i].f-=xx;               if(i%2==0)                    a[i-1].f+=xx;               else                    a[i+1].f+=xx;               s-=xx;          }     }     return t-s;}inline int maxflow(){     int s=0;     while(bfs())          s+=dfs(0,2100000000);     return s;}struct happiness{     int a,b,c,d,e,f;}map[101][101];int main(){ //左文右理      int n,m;     scanf("%d%d",&n,&m);     int i,j;     int sum=0;     for(i=1;i<=n;i++)     {          for(j=1;j<=m;j++)          {               scanf("%d",&map[i][j].a);               sum+=map[i][j].a;               edge++;               add(0,(i-1)*m+j,map[i][j].a*2);               edge++;               add((i-1)*m+j,0,0);          }     }     p=n*m+1;     for(i=1;i<=n;i++)     {          for(j=1;j<=m;j++)          {               scanf("%d",&map[i][j].b);               sum+=map[i][j].b;               edge++;               add((i-1)*m+j,p,map[i][j].b*2);               edge++;               add(p,(i-1)*m+j,0);          }     }     for(i=1;i<=n-1;i++)          for(j=1;j<=m;j++)               scanf("%d",&map[i][j].c);     for(i=1;i<=n-1;i++)     {          for(j=1;j<=m;j++)          {               scanf("%d",&map[i][j].d);               sum+=map[i][j].c;               sum+=map[i][j].d;                              edge++;               add(0,(i-1)*m+j,map[i][j].c);               edge++;   add((i-1)*m+j,0,0);               edge++;               add(0,i*m+j,map[i][j].c);               edge++;   add(i*m+j,0,0);      edge++;               add((i-1)*m+j,p,map[i][j].d);               edge++;   add(p,(i-1)*m+j,0);               edge++;               add(i*m+j,p,map[i][j].d);               edge++;   add(p,i*m+j,0);                     edge++;               add(i*m+j,(i-1)*m+j,map[i][j].c+map[i][j].d);               edge++;               add((i-1)*m+j,i*m+j,map[i][j].c+map[i][j].d);          }     }     for(i=1;i<=n;i++)          for(j=1;j<=m-1;j++)               scanf("%d",&map[i][j].e);     for(i=1;i<=n;i++)     {          for(j=1;j<=m-1;j++)          {               scanf("%d",&map[i][j].f);               sum+=map[i][j].e;               sum+=map[i][j].f;                              edge++;               add(0,(i-1)*m+j,map[i][j].e);               edge++;   add((i-1)*m+j,0,0);               edge++;               add(0,(i-1)*m+j+1,map[i][j].e);               edge++;   add((i-1)*m+j+1,0,0);      edge++;               add((i-1)*m+j,p,map[i][j].f);               edge++;   add(p,(i-1)*m+j,0);               edge++;               add((i-1)*m+j+1,p,map[i][j].f);               edge++;   add(p,(i-1)*m+j+1,0);               edge++;                              add((i-1)*m+j,(i-1)*m+j+1,map[i][j].e+map[i][j].f);               edge++;               add((i-1)*m+j+1,(i-1)*m+j,map[i][j].e+map[i][j].f);          }     }     int ans=maxflow();     printf("%d\n",sum-ans/2);     return 0;}