【一句话题意】现有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;
}