BZOJ2127 happiness-最小割

题目链接：右转进入题目

题目大意：请自行参考原题

解：

最小割经典题，有一种计算边权的方法并不会。

直接建图：

新建S表示文科，T表示理科。

对于每个同学X，建立(S,X,X同学学习文科的喜悦度),(X,T,X同学学习理科的喜悦度)

对于每对相邻的同学X和Y，其共同学习某文科的喜悦度为C，建立(S,新加一个节点,C),(新加的那个节点,X,无穷大),(新加的那个节点,Y,无穷大)

共同理科同理，换成T换个顺序即可。

统计非无穷大的边权之和，减去最小割就是答案。

代码：

//BZOJ2127
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<climits>
#define INF INT_MAX
#define P(x,y) (~-(x)*m+(y))
#define MAXN 60000
#define MAXM 300000
using namespace std;
struct edges{
	int to,resf,pre;
}e[MAXM];int etop;
int h[MAXN],cur[MAXN],lev[MAXN];
bool vis[MAXN];queue<int> q;
bool bfs(int s,int t)
{
	memset(lev,0,sizeof(lev));
	memset(vis,false,sizeof(vis));
	while(!q.empty()) q.pop();
	q.push(s);lev[s]=0;vis[s]=true;
	while(!q.empty())
	{
		int x=q.front();q.pop();
		for(int i=h[x];i;i=e[i].pre)
			if(!vis[e[i].to]&&e[i].resf)
			{
				vis[e[i].to]=true;
				lev[e[i].to]=lev[x]+1;
				q.push(e[i].to);
			}
	}
	return vis[t];
}
int dfs(int s,int t,int a)
{
	if(s==t||!a) return a;
	int flow=0,f;
	for(int &i=cur[s];i;i=e[i].pre)
		if(lev[e[i].to]==lev[s]+1&&(f=dfs(e[i].to,t,min(a,e[i].resf)))>0)
		{
			a-=f;flow+=f;
			e[i].resf-=f;
			e[((i-1)^1)+1].resf+=f;
			if(!a) break;
		}
	return flow;
}
int add_edge(int u,int v,int resc)
{
	etop++;
	e[etop].to=v;
	e[etop].resf=resc;
	e[etop].pre=h[u];
	h[u]=etop;
	return 0;
}
int main()
{
	int n,m;scanf("%d%d",&n,&m);
	int s=0,t=5*m*n+1,totans=0;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			int w;scanf("%d",&w);totans+=w;
			add_edge(s,P(i,j),w);add_edge(P(i,j),s,0);
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			int l;scanf("%d",&l);totans+=l;
			add_edge(P(i,j),t,l);add_edge(t,P(i,j),0);
		}
	for(int i=1;i<n;i++)
		for(int j=1;j<=m;j++)
		{
			int gw;scanf("%d",&gw);totans+=gw;
			int A=m*n+P(i,j),S1=P(i,j),S2=P(i+1,j);
			add_edge(s,A,gw);add_edge(A,s,0);
			add_edge(A,S1,INF);add_edge(S1,A,0);
			add_edge(A,S2,INF);add_edge(S2,A,0);
		}
	for(int i=1;i<n;i++)
		for(int j=1;j<=m;j++)
		{
			int gl;scanf("%d",&gl);totans+=gl;
			int B=2*m*n+P(i,j),S1=P(i,j),S2=P(i+1,j);
			add_edge(B,t,gl);add_edge(t,B,0);
			add_edge(S1,B,INF);add_edge(B,S1,0);
			add_edge(S2,B,INF);add_edge(B,S2,0);
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<m;j++)
		{
			int gw;scanf("%d",&gw);totans+=gw;
			int A=3*m*n+P(i,j),S1=P(i,j),S2=P(i,j+1);
			add_edge(s,A,gw);add_edge(A,s,0);
			add_edge(A,S1,INF);add_edge(S1,A,0);
			add_edge(A,S2,INF);add_edge(S2,A,0);
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<m;j++)
		{
			int gl;scanf("%d",&gl);totans+=gl;
			int B=4*m*n+P(i,j),S1=P(i,j),S2=P(i,j+1);
			add_edge(B,t,gl);add_edge(t,B,0);
			add_edge(S1,B,INF);add_edge(B,S1,0);
			add_edge(S2,B,INF);add_edge(B,S2,0);
		}
	int flow=0;
	while(bfs(s,t))
	{
		for(int i=s;i<=t;i++)
			cur[i]=h[i];
		flow+=dfs(s,t,INF);
	}
	printf("%d\n",totans-flow);return 0;
}