2127: happiness
Time Limit: 51 Sec Memory Limit: 259 MB
Submit: 1815 Solved: 878
[Submit][Status][Discuss]
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的非负整数
本题和上一题很相似，可以用相同的做法
但是本题可以优化
相邻两个人具有经典的棋盘黑白染色的特征，虽然本题貌似与二分图无关，这是句废话
能不能在一人一点的基础上直接处理出神秘加成
放弃文同时也是放弃了 神秘加成_文/2
两个人都放弃文就成功的放弃了整个 神秘加成_文
理同理
但是一人选文一人选理两个神秘加成都不能要啊，但这样每个神秘加成只放弃了1/2
那就再加两条边，链接相邻两个人，容量为(神秘加成_文+神秘加成_理)/2，割的时候一文一理时必定割去这条两条边之一
提示：
那两条边可以直接加无向边
注意：
煞笔Candy?数组开成[10005][10005] bzoj不是T就是M
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=1e4+,M=1e5+,INF=1e9;
inline int read(){
char c=getchar();int x=,f=;
while(c<''||c>''){if(c=='-')f=-; c=getchar();}
while(c>=''&&c<=''){x=x*+c-''; c=getchar();}
return x*f;
} int n,m,num,s,t,a[][],b[][],c,f[][],g[][],sum;
struct edge{
int v,c,f,ne;
}e[M<<];
int cnt,h[N];
inline void ins(int u,int v,int c){
cnt++;
e[cnt].v=v;e[cnt].c=c;e[cnt].f=;e[cnt].ne=h[u];h[u]=cnt;
cnt++;
e[cnt].v=u;e[cnt].c=;e[cnt].f=;e[cnt].ne=h[v];h[v]=cnt;
}
inline void ins2(int u,int v,int c){
cnt++;
e[cnt].v=v;e[cnt].c=c;e[cnt].f=;e[cnt].ne=h[u];h[u]=cnt;
cnt++;
e[cnt].v=u;e[cnt].c=c;e[cnt].f=;e[cnt].ne=h[v];h[v]=cnt;
}
int q[N],head,tail,vis[N],d[N];
bool bfs(){
memset(vis,,sizeof(vis));
memset(d,,sizeof(d));
head=tail=;
d[s]=;vis[s]=;
q[tail++]=s;
while(head!=tail){
int u=q[head++];
for(int i=h[u];i;i=e[i].ne){
int v=e[i].v;
if(!vis[v]&&e[i].c>e[i].f){
vis[v]=;
d[v]=d[u]+;
q[tail++]=v;
if(v==t) return true;
}
}
}
return false;
}
int cur[N];
int dfs(int u,int a){
if(u==t||a==) return a;
int flow=,f;
for(int &i=cur[u];i;i=e[i].ne){
int v=e[i].v;
if(d[v]==d[u]+&&(f=dfs(v,min(a,e[i].c-e[i].f)))>){
flow+=f;
e[i].f+=f;
e[((i-)^)+].f-=f;
a-=f;
if(a==) break;
}
}
return flow;
}
int dinic(){
int flow=;
while(bfs()){
for(int i=s;i<=t;i++) cur[i]=h[i];
flow+=dfs(s,INF);
}
return flow;
}
inline int id(int i,int j){return (i-)*m+j;}
int main(){
//freopen("in.txt","r",stdin);
n=read();m=read();
s=;t=n*m+;
for(int i=;i<=n;i++) for(int j=;j<=m;j++) a[i][j]=read()<<,sum+=a[i][j]>>;
for(int i=;i<=n;i++) for(int j=;j<=m;j++) b[i][j]=read()<<,sum+=b[i][j]>>;
for(int i=;i<=n-;i++) for(int j=;j<=m;j++){
c=read();a[i][j]+=c;a[i+][j]+=c;sum+=c;
f[i][j]+=c;
//ins2(id(i,j),id(i+1,j),c);
}
for(int i=;i<=n-;i++) for(int j=;j<=m;j++){
c=read();b[i][j]+=c;b[i+][j]+=c;sum+=c;
f[i][j]+=c;
//ins2(id(i,j),id(i+1,j),c);
}
for(int i=;i<=n;i++) for(int j=;j<=m-;j++){
c=read();a[i][j]+=c;a[i][j+]+=c;sum+=c;
g[i][j]+=c;
//ins2(id(i,j),id(i,j+1),c);
}
for(int i=;i<=n;i++) for(int j=;j<=m-;j++){
c=read();b[i][j]+=c;b[i][j+]+=c;sum+=c;
g[i][j]+=c;
//ins2(id(i,j),id(i,j+1),c);
}
for(int i=;i<=n;i++) for(int j=;j<=m;j++){
ins(s,id(i,j),a[i][j]);
ins(id(i,j),t,b[i][j]);
ins2(id(i,j),id(i+,j),f[i][j]);
ins2(id(i,j),id(i,j+),g[i][j]);
}
int ans=dinic();
printf("%d",sum-(ans>>));
}