约束 limit.cpp

【一句话题意】现有n个点，n条有向边。要求将n个点分入两个集合A，B，使得对于集合A中的每一个点x，存在一个集合B中的点y，从y到x有一条有向边。问集合A最多能包含几个点。
n<=1e5

【分析】首先需要确定的是，这是一道包着“集合”外壳的图论题。

乍一看，是否与求最大二分图匹配有几分神似？鉴于一句OI名言**“题目是死的，人是活的，出题人是懒的。”**机房一位大佬尝试着使用二分图算法，结果12ms轻松ac。

然而二分图实际上是有可能会被毒瘤数据给卡掉的，所以我们还是有必要学习正解。std使用一种毒瘤的类拓扑排序来做的。之所以是说类，是因为当关于这个图呈$\rho$形(关于$\rho$的小尾巴可能是一条链也可能是棵树)，正常的拓扑排序由于入度不为0，排序会停止。但这个类拓扑排序当节点可以被加入A集合时，有必要强制入队删除朝向他父亲的边。这样就会破坏一部分的$\rho$形。
最后再处理一下剩余的环形。

【code】

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1e5+100;
struct Queue{
	int a[maxn<<1];
	int l,r;
	bool empty(){return l==r;}
	void push(int x){a[r++]=x;r%=maxn<<1;}
	int front(){return a[l];}
	void pop(){l++;l%=maxn<<1;}
}q;
inline void read(int &x){
	x=0;char tmp=getchar();
	while(tmp<'0'||tmp>'9') tmp=getchar();
	while(tmp>='0'&&tmp<='9') x=(x<<1)+(x<<3)+tmp-'0',tmp=getchar();
}
int fa[maxn],cnt[maxn],rd[maxn];
bool vis[maxn];
int n,ans;
int main(){
	cin>>n;
	for(int i=1;i<=n;i++)
		read(fa[i]),rd[fa[i]]++;
	for(int i=1;i<=n;i++) if(!rd[i]) q.push(i),cnt[i]=0,vis[i]=true;
	while(!q.empty()){
		int x=q.front();q.pop();
		if(!cnt[x]&&!vis[fa[x]]) ++ans,cnt[fa[x]]=1,q.push(fa[x]),vis[fa[x]]=true;
		if(!--rd[fa[x]]&&!vis[fa[x]]) q.push(fa[x]);
	}
	memset(vis,0,sizeof(vis));
	for(int i=1;i<=n;i++){
		if(!rd[i]) continue;
		int cnt=0,x=i;
		while(!vis[x]) vis[x]=1,x=fa[x],cnt++;
		ans+=cnt/2;
	}
	cout<<ans<<endl;
	return 0;
}