[国家集训队]happiness

一、题目

点此看题

二、解法

UPD2021/1/15：我觉得逻辑是这样的，如果出现要么选这个，要么选那个，就可以考虑用最小割。比如说这道题是要么选文科，要么选理科，选这个等价于割掉那个。

这道题的正解是 总贡献-最小割，可能是题目要求最大值我们就用最小割去减 $--\text{SYDeive}$
考虑建图，模型大概长这样：

• 我们把每个点连$s$，容量为文科的贡献，连$t$，容量为理科的贡献。
• 考虑两人相同的新贡献，我们新建文科点和理科点，把文科点连$s$，容量为文科相同的贡献，把理科点连$t$，容量为理科相同的贡献，文科点和理科点再连接建立关系的两个点，容量为$inf$。

怎么理解这个模型的正确性呢？可以考虑分类讨论，假设文科不同：
可能的组合是$1,2,3$或$1,2,5,6$或$1,3,4,6$，发现只有这几种割的情况，没有与题目冲突，假设理科不同也是同理，所以答案就是$sum-flow$。

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
const int MAXN = 100005;
using namespace std;
int read()
{
    int num=0,flag=1;char c;
    while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
    while(c>='0'&&c<='9')num=(num<<3)+(num<<1)+(c^48),c=getchar();
    return num*flag;
}
int n,m,k,S,T,t,tot=1,ans,sum,f[MAXN],cur[MAXN],dis[MAXN];
int dx[4]={1,-1},dy[4]={0,0,1,-1};
queue<int> q;
struct edge
{
	int v,c,next;
}e[MAXN*10];
int cal(int x,int y)
{
	return (x-1)*m+y;
}
void add_edge(int u,int v,int c)
{
	e[++tot]=edge{v,c,f[u]},f[u]=tot;
	e[++tot]=edge{u,0,f[v]},f[v]=tot;
}
int bfs()
{
	memset(dis,0,sizeof dis);
	dis[S]=1;
	q.push(S);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=f[u];i;i=e[i].next)
		{
			int v=e[i].v;
			if(e[i].c>0 && !dis[v])
			{
				dis[v]=dis[u]+1;
				q.push(v);
			}
		}
	}
	if(!dis[T]) return 0;
	return 1;
}
int dfs(int u,int ept)
{
	if(u==T) return ept;
	int flow=0,tmp=0;
	for(int &i=cur[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(dis[u]+1==dis[v] && e[i].c>0)
		{
			tmp=dfs(v,min(e[i].c,ept));
			if(!tmp) continue;
			ept-=tmp;
			e[i].c-=tmp;
			e[i^1].c+=tmp;
			flow+=tmp;
			if(!ept) break;
		}
	}
	return flow;
}
int main()
{
	n=read();m=read();
	S=0;T=k=cal(n,m)+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			add_edge(S,cal(i,j),t=read());
			sum+=t;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			add_edge(cal(i,j),T,t=read());
			sum+=t;
		}
	for(int i=1;i<=n-1;i++)
		for(int j=1;j<=m;j++)
		{
			int a=cal(i,j),b=cal(i+1,j),c=++k;
			add_edge(S,c,t=read());
			add_edge(c,b,inf);
			add_edge(c,a,inf);
			sum+=t;
		}
	for(int i=1;i<=n-1;i++)
		for(int j=1;j<=m;j++)
		{
			int a=cal(i,j),b=cal(i+1,j),c=++k;
			add_edge(c,T,t=read());
			add_edge(b,c,inf);
			add_edge(a,c,inf);
			sum+=t;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m-1;j++)
		{
			int a=cal(i,j),b=cal(i,j+1),c=++k;
			add_edge(S,c,t=read());
			add_edge(c,b,inf);
			add_edge(c,a,inf);
			sum+=t;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m-1;j++)
		{
			int a=cal(i,j),b=cal(i,j+1),c=++k;
			add_edge(c,T,t=read());
			add_edge(b,c,inf);
			add_edge(a,c,inf);
			sum+=t;
		}
	while(bfs())
	{
		for(int i=0;i<=k;i++)
			cur[i]=f[i];
		ans+=dfs(S,inf);
	}
	printf("%d\n",sum-ans);
}