luogu2711小行星

题目

　　https://www.luogu.org/problem/show?pid=2711

题解

　　我说了自己的想法之后他们都不信，还说我智障，然而我AC了。

　　其实是最小割，一共4排点，x y y z。S向每个x点连c=1的边代表限流，y和y之间连边c=1，z和T之间连边c=1，如果一个小行星是(x,y,z)，就把对应的点之间连起来。现在考虑，如果图中存在一个流从S流经x流经y流经z注入T，说明这颗小行星依然存在，对应的你应该割掉一条边使流消失，也就是这颗行星被去掉了。（割的边的位置对应你去掉的是哪个面），所以就是最小割了。

　　最大流=最小割，跑出最大流就是答案。

　　这道题目可以扩展到N维，方法一样。

代码

//最大流
#include <cstdio>
#include <algorithm>
#include <cstring>
#define maxn 1000000
#define inf 0x3f3f3f3f
using namespace std;
int S, T, head[maxn], to[maxn], c[maxn], d[maxn], num[maxn], last[maxn], Exit, tot=1,
	vis[3][600][600], next[maxn];
void adde(int a, int b, int cc){to[++tot]=b;c[tot]=cc;next[tot]=head[a];head[a]=tot;}
void adde2(int a, int b, int cc){adde(a,b,cc),adde(b,a,0);}
int isap(int pos, int in)
{
	int flow=0, t;
	if(pos==T)return in;
	for(int &p=last[pos];p;p=next[p])
		if(c[p] and d[to[p]]+1==d[pos])
		{
			flow+= t=isap(to[p],min(in-flow,c[p]));
			c[p]-=t, c[p xor 1]+=t;
			if(Exit or in==flow)return flow;
		}
	Exit=--num[d[pos]]==0;
	++num[++d[pos]];
	last[pos]=head[pos];
	return flow;
}
void build()
{
	int n, x, y, z, t=500, i;
	scanf("%d",&n);
	S=2001, T=S+1;
	for(i=1;i<=t;i++)adde2(S,i,1),adde2(t+t+t+i,T,1),adde2(t+i,t+t+i,1);
	while(n--)
	{
		scanf("%d%d%d",&x,&y,&z);
		if(vis[1][x][y]==0)vis[1][x][y]=1,adde2(x,t+y,1);
		if(vis[2][y][z]==0)vis[2][y][z]=1,adde2(t+t+y,t+t+t+z,1);
	}
}
int main()
{
	int ans;
	build();
	for(ans=0;!Exit;ans+=isap(S,inf));
	printf("%d",ans);
	return 0;
}