#include<bits/stdc++.h>
using namespace std;
const int MOD=998244353;
int main()
{
	int n;
	cin>>n;
	vector<int> a(n+5);
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
	}
	auto poww=[&](long long x,long long y,long long z)
	{
		long long ret=1;
		while(y)
		{
			if(y&1)ret*=x,ret%=z;
			x*=x,x%=z;
			y>>=1;
		}
		return ret;
	};
	vector<int> vis(n+5);
	int maxp2=0;
	map<int,int> mp;
	for(int i=1;i<=n;i++)
	{
		if(not vis[i])
		{
			int cnt=0;
			for(int j=i;not vis[j];j=a[j])
			{
//				cerr<<"? "<<j<<endl;
				vis[j]=1;
				cnt++;
			}
			int cp2=0;
			while(cnt%2==0)
				cnt/=2,cp2++;
			maxp2=max(maxp2,cp2);
			
//			cerr<<"cycle "<<cnt<<endl;
			
			int t=cnt;
			long long phi=1;
			for(int i=2;i*i<=t;i++)
			{
				if(t%i==0)
				{
					phi*=i-1;
					while(t%i==0)t/=i,phi*=i;
					phi/=i;
				}
			}
			if(t!=1)phi*=t-1;
			
//			cerr<<"phi "<<phi<<endl;
			
			t=phi;
			vector<int> pr;
			for(int i=2;i*i<=phi;i++)
			{
				if(phi%i==0)
				{
					pr.push_back(i);
					while(phi%i==0)phi/=i;
				}
			}
			if(phi!=1)pr.push_back(phi);
			
			phi=t;
			for(auto p:pr)
			{
//				cerr<<"trail div "<<p<<endl;
//				cerr<<t/p<<' '<<poww(2,t/p,cnt)<<endl;
				while(t%p==0 and poww(2,t/p,cnt)==1)
				{
					t/=p;
				}
			}
			
			cnt=t;
			for(int x=2;x*x<=cnt;x++)
			{
				if(cnt%x==0)
				{
					int cpk=0;
					while(cnt%x==0)
						cnt/=x,cpk++;
					mp[x]=max(mp[x],cpk);
				}
			}
			if(cnt!=1)mp[cnt]=max(mp[cnt],1);
		}
	}
	long long mul=1;
	for(auto [x,y]:mp)
	{
		for(int i=1;i<=y;i++)
			mul=(mul*x)%MOD;
	}
	cout<<(mul+maxp2)%MOD;
	return 0;
}