文理分科

题目

传送门 to luogu

题目中的 “看到这个题目的人肯定都没有纠结过” 好真实啊……

思路

遇事不决，量子力学 考虑 $\mathtt{d}\mathtt{p}$ 。但是我们很快就吃瘪了，因为一个人的选择会造成深远的影响，估计只有状态压缩能做。贪心就更不行了。

那么考虑建立某种模型。根据题目标签可知，我们要建立网络流图。那么 $science$ 和 $art$ 是很好得到的。然而 $same$ 怎么算？

这时候我又想起 最大权闭合子图 来了。别说，这东西还真好用。$same$ 完全可以看成一个独立于每个人之外的物品！不过它的选择需要一些先决条件！

$same\_art$ 需要 $5$ 个人选择文科。那么就用一个点代表这个贡献，$S$ 向它连容量为价值的边，割掉它相当于不选择这玩意儿。不割掉的话，它向这五个人连 $+\infty$ 边，这五个人向 $T$（代表理科）的连边就必然被割掉，于是必然选文科。$same\_science$ 类似，它向 $T$ 连边嘛。

这个题好像比较基础。我借此来理解一下网络流：

• 一般都是考虑最小割。只是求解的时候利用最大流罢了。
• 如果一个主体进行两个冲突的选择，就 $S\to x\to T$ ，容量分别为权值。
• 如果一个选择将会带来别的选择，可以 $S\to x\to y\to T$ ，表示留 $x$ 则 $y$ 选 $S$ ；或者 $S\to y\to x\to T$ ，表示不留 $x$ 则 $y$ 选 $T$ 。

如果以后还有别的姿势，当然是不会放到这个博客里的啦 见到再说吧。

复杂度 $$\mathcal{O}(n^3{m}^3),\Omega (knm)$$

我这辈子犯过的最大错误：试图计算网络流复杂度。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
typedef long long int_;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}

const int MaxM = 50005;
struct Edge{
	int to, nxt, val;
	Edge(){ /* nothing */ ; }
	Edge(int T,int N,int V){
		to = T, nxt = N, val = V;
	}
};
Edge e[MaxM*10];
int head[MaxM], cntEdge;
void addEdge(int a,int b,int c){
	e[cntEdge] = Edge(b,head[a],c);
	head[a] = cntEdge ++;
	e[cntEdge] = Edge(a,head[b],0);
	head[b] = cntEdge ++;
}

queue< int > q; int dis[MaxM];
bool bfs(int s,int t,int n){
	while(!q.empty()) q.pop();
	memset(dis,-1,n<<2);
	q.push(s), dis[s] = 0;
	while(!q.empty()){
		s = q.front(); q.pop();
		for(int i=head[s]; ~i; i=e[i].nxt)
			if(e[i].val && !(~dis[e[i].to])){
				dis[e[i].to] = dis[s]+1;
				q.push(e[i].to);
			}
	}
	return dis[t] != -1;
}
int cur[MaxM]; // current edge optimization
int dfs(int x,int inFlow,const int &t){
	int sum = 0; if(x == t) return inFlow;
	for(int &i=cur[x]; ~i; i=e[i].nxt)
		if(dis[e[i].to] == dis[x]+1 && e[i].val){
			int d = min(inFlow-sum,e[i].val);
			d = dfs(e[i].to,d,t);
			e[i].val -= d, e[i^1].val += d;
			if((sum += d) == inFlow) break;
		}
	if(!sum) dis[x] = -1; // eat s**t
	return sum;
}
const int infty = (1<<30)-1;
int dinic(int s,int t,int n){
	int res = 0;
	while(bfs(s,t,n)){
		memcpy(cur,head,n<<2);
		res += dfs(s,infty,t);
	}
	return res;
}

int main(){
	memset(head,-1,MaxM<<2);
	int n = readint(), m = readint();
	int ans = 0; // start calc now!
	for(int i=0; i<n; ++i)
	for(int j=0; j<m; ++j){
		int v = readint(); ans += v;
		addEdge(0,i*m+j+1,v);
	}
	for(int i=0; i<n; ++i)
	for(int j=0; j<m; ++j){
		int v = readint(); ans += v;
		addEdge(i*m+j+1,n*m+1,v);
	}
	int id = n*m+1; // allocate id
	for(int i=0; i<n; ++i)
	for(int j=0; j<m; ++j){
		int v = readint(); ans += v;
		addEdge(0,++ id,v);
		for(int x=i-1; x<=i+1; ++x)
		for(int y=j-1; y<=j+1; ++y)
		if((i-x)*(i-x)+(y-j)*(y-j) <= 1
		&& 0 <= x && x < n && 0 <= y && y < m)
			addEdge(id,x*m+y+1,infty);
	}
	for(int i=0; i<n; ++i)
	for(int j=0; j<m; ++j){
		int v = readint(); ans += v;
		addEdge(++ id,n*m+1,v);
		for(int x=i-1; x<=i+1; ++x)
		for(int y=j-1; y<=j+1; ++y)
		if((i-x)*(i-x)+(y-j)*(y-j) <= 1
		&& 0 <= x && x < n && 0 <= y && y < m)
			addEdge(x*m+y+1,id,infty);
	}
	printf("%d\n",ans-dinic(0,n*m+1,id+1));
	return 0;
}