bzoj 2147 happiness

2127: happiness

Time Limit: 51 Sec  Memory Limit: 259 MB
Submit: 1300  Solved: 623
[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

Sample Output

1210
【样例说明】
两人都选理，则获得100+110+1000的喜悦值。
【数据规模】
对于100%以内的数据，n,m<=100 所有喜悦值均为小于等于5000的非负整数

题解：网络流求最小割

考虑到一个人，文理不可兼得，不妨先建点，然后向源点（文科），汇点（理科）连边，流量（也就是割）是对应喜悦值。（这里的想法是先建个差不多的，有漏洞再拆点啊，建辅助点啊什么的）

    然后再考虑一对朋友之间的共文理喜悦值：

如果都选文，那么需要割掉双方都选理的喜悦值，

如果都选理，那么需要割掉双方都选文的喜悦值，

如果一文一理，那么就都割。

首先对于一个点对，我们考虑到他们之间的关系只在于都选文或者都选理的喜悦值，

而求最小割时我们要割掉这部分权值，那么我们可以建边，然后把边进行处理。

一、

首先我最开始的想法是把边的权值设置为都选文+都选理的和。

这样就成功地满足了一文一理的情况。

但是却无法分开。

二、

这时候我们可以把边建成点，然后两边的人向此点连容量inf的无向边，

然后源点（文科）向此点连都选文的容量，汇点连都选理的。

这样我们就满足了都选文或者都选理的情况。

也就是都选文的话，我们就可以不割“都选文”这种权值，都选理 同理。

但是这是WA的，cheat也不行（实测），因为我们又无法让两人一文一理了。

三、

再往“一”那里考虑，发现我们有时需要把两种权值都割断。

这个时候想到：

a-(len)-->b，当len值一定，把图建成a-->b-->c-->d，只要边长都是len，那么最小割是不变的。

所以我们把“边点”拆成两个点，一个对应S集（文），一个对应T集（理），，

这样我们就可以把原来的inf也改成对应权值，使得割断哪一条边都是最小割值。

这么说可能有些含糊，不妨来张图。

ps:分析转载自http://blog.csdn.net/vmurder/article/details/42609669

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#include<queue>using namespace std;int n,m,i,j,num,tot,sum,numl;int next[600000],v[600000],cost[6000000],point[50003];int deep[50003],cur[50003];const int inf=1e9;void add(int x,int y,int z){  tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; cost[tot]=z;  tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; cost[tot]=0;}bool bfs(int s,int t){ memset(deep,0x7f,sizeof(deep)); for (int i=0;i<=numl;i++)  cur[i]=point[i]; deep[0]=0; queue<int> p; p.push(0); while (!p.empty())  {  int now=p.front(); p.pop();  for (int i=point[now];i!=-1;i=next[i])   {   if(deep[v[i]]>inf&&cost[i])    {    deep[v[i]]=deep[now]+1;    p.push(v[i]);    }   }  } if (deep[t]>inf) return false; else  return true;}int dfs(int now,int t,int limit){if (now==t||!limit)  return limit;int flow=0,f;for (int i=cur[now];i!=-1;i=next[i]) { cur[now]=i; if (deep[v[i]]==deep[now]+1&&(f=dfs(v[i],t,min(limit,cost[i]))))  {  flow+=f; limit-=f;  cost[i]-=f; cost[i^1]+=f;  if (!limit) break;  } }return flow;}int dinic(int s,int t){  int ans=0;  while(bfs(s,t))   ans+=dfs(s,t,inf);  return ans;}int main(){  scanf("%d%d",&n,&m);  tot=-1;  memset(point,-1,sizeof(point));  memset(next,-1,sizeof(next));  num=n*m;  numl=n*m+2*(n-1)*m+2*(m-1)*n+1;  for (i=1;i<=n;i++)   for (j=1;j<=m;j++)    {      int x; scanf("%d",&x); sum+=x;      add(0,(i-1)*m+j,x);    }  for (i=1;i<=n;i++)   for (j=1;j<=m;j++)    {      int x; scanf("%d",&x); sum+=x;      add((i-1)*m+j,numl,x);    }  for (i=1;i<=n-1;i++)   for (j=1;j<=m;j++)    {      int x; scanf("%d",&x); sum+=x;      num++;      add(num,(i-1)*m+j,x);      add(0,num,x);      add(num,i*m+j,x);    }  for (i=1;i<=n-1;i++)   for (j=1;j<=m;j++)    {       int x; scanf("%d",&x); sum+=x;       num++;       add((i-1)*m+j,num,x);       add(i*m+j,num,x);       add(num,numl,x);    }  for (i=1;i<=n;i++)   for (j=1;j<=m-1;j++)    {      int x; scanf("%d",&x); sum+=x;      num++;      add(num,(i-1)*m+j,x);      add(num,(i-1)*m+j+1,x);      add(0,num,x);    }  for (i=1;i<=n;i++)   for (j=1;j<=m-1;j++)    {      int x; scanf("%d",&x); sum+=x;      num++;      add((i-1)*m+j,num,x);      add((i-1)*m+j+1,num,x);      add(num,numl,x);    }  printf("%d",sum-dinic(0,numl));}