#include<deque>
#include<queue>
#include<vector>
#include<algorithm>
#include<iostream>
#include<set>
#include<cmath>
#include<tuple>
#include<string>
#include<chrono>
#include<functional>
#include<iterator>
#include<random>
#include<unordered_set>
#include<unordered_map>
#include<array>
#include<map>
#include<bitset>
#include<iomanip>
#include<list>
#include <numeric>
using namespace std;
typedef unsigned long long int ulint;
typedef long long int llint;
typedef long double lldo;
#define mp make_pair
#define mt make_tuple
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pop_front
#define fir first
#define sec second
#define res resize
#define ins insert
#define era erase
#define RE return 0
//ios::sync_with_stdio(false);
//std::cin.tie(0);
//<< setprecision(20)
const int mod=(int)1e9+7;
const llint big=(llint)(2.19e15+1);
const long double pai=3.141592653589793238462643383279502884197;
const long double ena=2.71828182845904523536;
const long double eps=1e-7;
template <class T,class U>bool mineq(T& a,U b){if(a>b){a=b;return true;}return false;}
template <class T,class U>bool maxeq(T& a,U b){if(a<b){a=b;return true;}return false;}
template <class T> void soun(T& ar)
{sort(ar.begin(),ar.end());ar.erase(unique(ar.begin(),ar.end()),ar.end());}
llint gcd(llint a,llint b){if(a%b==0){return b;}else{return gcd(b,a%b);}}
llint lcm(llint a,llint b){return a/gcd(a,b) *b;}
template<class T,class U> auto LB(T& ve,U in){return lower_bound(ve.begin(),ve.end(),in);}
template<class T,class U> auto UB(T& ve,U in){return upper_bound(ve.begin(),ve.end(),in);}
template<class T,class U> auto LBI(T& ve,U in){return LB(ve,in)-ve.begin();}
template<class T,class U> auto UBI(T& ve,U in){return UB(ve,in)-ve.begin();}
template<class T> void SO(T& ve){sort(ve.begin(),ve.end());}
template<class T> void REV(T& ve){reverse(ve.begin(),ve.end());}
llint rui(llint p,llint n){
	llint ans=1;
	for(int h=0;h<30;h++){
		if(n&(1<<h)){ans*=p;ans%=mod;}
		p*=p;p%=mod;
	}
	//cerr<<"p="<<p<<"ans="<<ans<<endl;
	return ans;
}
int main(void){
	int n,i,j;cin>>n;
	vector<int>a(n);
	vector<int>na(n,-1);
	vector<int>nb(n,-1);
	queue<int>futa;//二つ以上あるものをカウント
	for(i=0;i<n;i++){
		cin>>a[i];a[i]--;
		if(na[a[i]]==-1){na[a[i]]=i;}
		else if(nb[a[i]]==-1){nb[a[i]]=i;futa.push(a[i]);}
		else{cout<<0<<endl;return 0;}
	}
	llint ans=1;
	//aにかぶりがなくなるまで圧縮する
	//同じ数が三つ以上ある->死
	//連続マージは起こらない
	while(futa.size()>0){
		int ter=futa.front();futa.pop();
		int aa=na[ter],bb=nb[ter];
		if(na[aa]!=-1&&nb[bb]!=-1){cout<<0<<endl;return 0;}
		else if(nb[aa]!=-1&&na[bb]!=-1){cout<<0<<endl;return 0;}//合計3以上
		else if(na[aa]!=-1&&na[bb]!=-1){nb[aa]=na[bb];futa.push(aa);}//1 1でのマージ
		else if(nb[aa]==-1&&nb[bb]==-1){ans*=2;ans%=mod;}//2がなく 1 1でもないため、合計は1以下　ans*=2
		a[aa]=-1;a[bb]=-1;
	}
	vector<int>roop(n+1);
	for(int sta=0;sta<n;sta++){
		int gen=a[sta];
		if(gen==-1){continue;}
		if(gen==sta){roop[1]++;continue;}
		int num=1;
		while(-1){
			int og=gen;
			gen=a[gen];
			a[og]=-1;
			num++;
			if(gen==-1){num=0;break;}
			if(gen==sta){break;}
		}
		roop[num]++;
	}
	//roop 対策
	for(i=1;i<=n;i++){
		llint gen=1,now=1,kai=1;
		for(j=roop[i];j>1;j-=2){
			now*=j;now%=mod;
			now*=(j-1);now%=mod;
			now*=i;now%=mod;
			now*=rui(kai*2,mod-2);now%=mod;
			if(i%2==1&&i>1){gen*=4;}
			gen+=now;gen%=mod;
			kai++;
		}
		if(j%2==1&&i%2==1&&i>1){gen*=2;}
		ans*=gen;ans%=mod;
	}cout<<ans<<endl;
	RE;
}

Measured memory_bytes: 1920