///////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////
//                       _oo0oo_                         //
//                      o8888888o                        //
//                      88" . "88      ------ hzt1       //
//                      (| -_- |)                        //
//                      0\  =  /0                        //
//                    ___/`---'\___                      //
//                  .' \|     |// '.                     //
//                 / \|||  :  |||// \                    //
//                / _||||| -:- |||||- \                  //
//               |   | \  -  /// |     |                 //
//               | \_|  ''\---/''  |_/ |                 //
//               \  .-\__  '-'  ___/-. /                 //
//             ___'. .'  /--.--\  `. .'___               //
//          ."" '<  `.___\_<|>_/___.' >' "".             //
//         | | :  `- \`.;`\ _ /`;.`/ - ` : | |           //
//         \  \ `_.   \_ __\ /__ _/   .-` /  /           //
//     =====`-.____`.___ \_____/___.-`___.-'=====        //
//                       `=---='                         //
//                                                       //
//                                                       //
//     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~       //
//                                                       //
//                 God-He Bless All.                     //
//           This Code Will Never Explode.               //
//                                                       //
//                                                       //
///////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////
#include<cstdio>
#include<cstring>
#include<cctype>
#include<algorithm>
#include<vector>
#define pb push_back
#define mp make_pair
#define xx first
#define yy second
#define rep(i,a,b) for(int i=(a),i##_end_=(b);i<=i##_end_;i++)
#define dwn(i,a,b) for(int i=(a),i##_end_=(b);i>=i##_end_;i--)
using namespace std;
const int Size=1<<16;
char buffer[Size],*head,*tail;
inline char Getchar() {
    if(head==tail) {
        int l=fread(buffer,1,Size,stdin);
        tail=(head=buffer)+l;
    }
    if(head==tail) return -1;
    return *head++;
}
inline int read() {
    int x=0,f=1;char c=Getchar();
    for(;!isdigit(c);c=Getchar()) if(c=='-') f=-1;
    for(;isdigit(c);c=Getchar()) x=x*10+c-'0';
    return x*f;
}
typedef long long ll;
typedef pair<int,int> pii;
const int maxn=100010;
const int mod=1000000007;
int n,p[maxn],vis[maxn],cir[maxn],deg[maxn],len[maxn],cnt[maxn],f[maxn];
void GG() {puts("0");exit(0);}
int cmp(int a,int b) {
	if(a<b) return 2;
	if(a==b) return 1;
	return 0;
}
int main() {
	n=read();
	rep(i,1,n) deg[p[i]=read()]++;
	rep(i,1,n) if(!vis[i]) {
		int x=i;
		while(!vis[x]) vis[x]=i,x=p[x];
		if(vis[x]!=i) continue;
		while(!cir[x]) cir[x]=1,x=p[x];
	}
	rep(i,1,n) if((!cir[i]&&deg[i]>1)||(cir[i]&&deg[i]>2)) GG();
	rep(i,1,n) if(!deg[i]) {
		int x=i,c=0;
		while(!cir[x]) c++,x=p[x];
		len[x]=c;
	}
	int ans=1;
	rep(i,1,n) if(cir[i]) {
		int x=i,fl=0,fi=0,k=0,la=0;
		while(cir[x]) {
			k++;cir[x]=0;
			if(len[x]) {
				if(!fi) {fi=la=k;fl=len[x];x=p[x];continue;}
				ans=(ll)ans*cmp(len[x],k-la)%mod;la=k;
				if(!ans) GG();
			}
			x=p[x];
		}
		if(!fi) cnt[k]++;
		else ans=(ll)ans*cmp(fl,fi+k-la)%mod;
		if(!ans) GG();
	}
	rep(i,1,n) {
		int tp=1;if(i>1&&i%2==1) tp++;
		f[0]=1;
		rep(j,1,cnt[i]) {
			f[j]=f[j-1]*tp%mod;
			if(j>=2) (f[j]+=(ll)f[j-2]*(j-1)%mod*i%mod)%=mod;
		}
		ans=(ll)ans*f[cnt[i]]%mod;
	}
	printf("%d\n",ans);
	return 0;
}


Measured memory_bytes: 2176