题目
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;
}