TLE on7 玄关

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=2e5+10,mod=1e9+7,inf=1e17;
int n,a[N],head[N],idx,cnt;
struct edge{
	int v,next;
}e[N*2];
void con(int u,int v){
	e[++idx].v=v;
	e[idx].next=head[u];
	head[u]=idx;
}
int ksm(int a,int b){
	int ans=1;
	while(b){
		if(b&1)ans*=a;
		a*=a;
		a%=mod;
		ans%=mod;
		b>>=1;
	}
	return ans;
}
int vis[N],nxt[N],fl;
void fr(int u,int fa){
	for(int i=head[u];i;i=e[i].next){
		int v=e[i].v;
		if(v==fa) continue;
		if(!vis[v]){
			nxt[v]=u;
			vis[v]=1;
			fr(v,u);
		}
		else{
			if(fl) continue;
			fl=1;
			int p=u;
			while(p!=v&&p){
				cnt++;
				p=nxt[p];
			}
		}
	}
}
inline int read(){
	int w=1,z=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') {
		if(ch=='-') w=-1;
		ch=getchar();
	}
	while(ch>='0'&&ch<='9') {
		z=z*10+ch-'0';
		ch=getchar();
	}
	return w*z;
}
signed main(){
	n=read();
	int ans=1;
	for(int i=1;i<=n;i++){
		int x=read();
		con(i,x),con(x,i);
	}	
	int sum=0;
	for(int i=1;i<=n;i++){
		if(vis[i]) continue;
		cnt=0;
		fl=0;
		fr(i,0);
		cnt++;
		ans=(ans*((ksm(2,cnt)-2+mod)%mod)%mod);
		sum+=cnt;
	}
	ans=(ans*ksm(2,n-sum))%mod;
	printf("%lld",ans);
	
	return 0;
}