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
1 1
100 110
1
1000
Sample Output
1210
【样例说明】
两人都选理,则获得100+110+1000的喜悦值。
【数据规模】
对于100%以内的数据,n,m<=100 所有喜悦值均为小于等于5000的非负整数
【样例说明】
两人都选理,则获得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;
}