#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<P,int> PP;
int n;

ll gcd(ll a,ll b){
	if(b==0LL)return a;
	return gcd(b,a%b);
}

ll lcm(ll a,ll b){
	return a/gcd(a,b)*b;
}

void f(vector<int>& a,vector<int>& c,vector<int>& d,vector<int>& s){
	for(int i=0;i<n;i++){
		scanf("%d",&a[i]);
		a[i]--;
	}
	int cnt=0;
	for(int i=0;i<n;i++){
		if(c[i]==-1){
			int j=i,num=0;
			do{
				c[j]=cnt;
				d[j]=num;
				num++;
				j=a[j];
			}while(j!=i);
			s.push_back(num);
			cnt++;
		}
	}
} 

vector<int> a,b,a_gr,a_ord,a_size,b_gr,b_ord,b_size;
map<PP,vector<P> > mp;

int main(void){
	scanf("%d",&n);
	a=b=a_ord=b_ord=vector<int>(n,0);
	a_gr=b_gr=vector<int>(n,-1);
	f(a,a_gr,a_ord,a_size);
	f(b,b_gr,b_ord,b_size);
	ll ans=0;
	for(int i=0;i<n;i++){
		int g=gcd(a_size[a_gr[i]],b_size[b_gr[i]]);
		mp[PP(P(a_gr[i],b_gr[i]),(a_ord[i]%g-b_ord[i]%g+g)%g)].push_back(P(a_ord[i],b_ord[i]));
	}
	map<PP,vector<P> >::iterator it=mp.begin();
	for(;it!=mp.end();it++){
		P tmp=it->first.first;
		vector<P> vp=it->second;
		int A=a_size[tmp.first],B=b_size[tmp.second],S=vp.size();
		ll l=lcm(A,B),g=gcd(A,B),P=vp[0].first,Q=vp[0].second;
		vector<ll> u(S+1),w(B/g);
		u[S]=l;
		int j=Q/g;
		for(int i=0;i<B/g;i++){
			(j+=A/g)%=B/g;
			w[j]=B/g-i-1;
		}
		for(int i=1;i<S;i++){
			ll X=vp[i].first,Y=vp[i].second;
			ll num=(P+A-X)%A;
			(Y+=num)%=B;
			num+=w[Y/g]*A;
			u[i]=num;
		}
		sort(u.begin(),u.end());
		for(int i=0;i<S;i++){
			ll sa=(u[i+1]-u[i])%MOD;
			(ans+=sa*(sa-1LL)/2LL)%=MOD;
		}
	}
	printf("%lld\n",ans);
	return 0;
}