// #pragma GCC target("avx2")  // CPU 処理並列化
// #pragma GCC optimize("O3")  // CPU 処理並列化
// #pragma GCC optimize("unroll-loops")  // 条件処理の呼び出しを減らす
#include<stdio.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<deque>
#include<stack>
#include<string>
#include<string.h>
#include<vector>
#include<set>
#include<map>
#include<bitset>
#include<stdlib.h>
#include<cassert>
#include<time.h>
#include<bitset>
#include<numeric>
#include<unordered_set>
#include<complex>
using namespace std;
const long long mod=1000000007;
const long long inf=mod*mod;
const long long d2=(mod+1)/2;
const double EPS=1e-10;
const double INF=1e+10;
const double PI=acos(-1.0);
const int C_SIZE = 5121000;
namespace{
	long long fact[C_SIZE];
	long long finv[C_SIZE];
	long long inv[C_SIZE];
	long long Comb(int a,int b){
	 	if(a<b||b<0)return 0;
	 	return fact[a]*finv[b]%mod*finv[a-b]%mod;
	}
	void init_C(int n){
		fact[0]=finv[0]=inv[1]=1;
		for(int i=2;i<n;i++){
			inv[i]=(mod-(mod/i)*inv[mod%i]%mod)%mod;
		}
		for(int i=1;i<n;i++){
			fact[i]=fact[i-1]*i%mod;
			finv[i]=finv[i-1]*inv[i]%mod;
		}
	}
	long long pw(long long a,long long b){
		if(a<0LL)return 0;
		if(b<0LL)return 0;
		long long ret=1;
		while(b){
			if(b%2)ret=ret*a%mod;
			a=a*a%mod;
			b/=2;
		}
		return ret;
	}
	long long pw_mod(long long a,long long b,long long M){
		if(a<0LL)return 0;
		if(b<0LL)return 0;
		long long ret=1;
		while(b){
			if(b%2)ret=ret*a%M;
			a=a*a%M;
			b/=2;
		}
		return ret;
	}
	int pw_mod_int(int a,int b,int M){
		if(a<0)return 0;
		if(b<0)return 0;
		int ret=1;
		while(b){
			if(b%2)ret=(long long)ret*a%M;
			a=(long long)a*a%M;
			b/=2;
		}
		return ret;
	}
	int ABS(int a){return max(a,-a);}
	long long ABS(long long a){return max(a,-a);}
	double ABS(double a){return max(a,-a);}
	int sig(double r) { return (r < -EPS) ? -1 : (r > +EPS) ? +1 : 0; }
}
// ここから編集しろ
int x[210000];
int y[210000];
int UF[210000];
vector<int>g[210000];
int FIND(int a){
	if(UF[a]<0)return a;return UF[a]=FIND(UF[a]);
}
void UNION(int a,int b){
	a=FIND(a);b=FIND(b);if(a==b)return;UF[a]+=UF[b];UF[b]=a;
}
int h[210000];
int v[210000];
vector<int>cy;
vector<int>cur;
int dfs(int a,int b){
	v[a]=-1;
	
	cur.push_back(a);
	int ret=-2;
	
	for(int i=0;i<g[a].size();i++){
		if(g[a][i]==b)continue;
		if(v[g[a][i]]==-2){
			int val=dfs(g[a][i],a);
			if(val!=-2){
				cy.push_back(a);
				if(a==val)val=-2;
			}
			if(val!=-2) ret=val;
		}else if(v[g[a][i]]!=a){
			cy.push_back(a);
			ret=g[a][i];
		}
	}
	v[a]=ret;
	// printf("%d %d: %d\n",a,b,ret);fflush(stdout);
	return ret;
}
vector<int>g2[210000];
int sz[210000];
int ou[210000];
long long calc(int a){
	long long ret=1;
	sz[a]=0;
	for(int i=0;i<g2[a].size();i++){
		ret=ret*calc(g2[a][i])%mod;
		ret=ret*Comb(sz[a]+sz[g2[a][i]],sz[a])%mod;
		sz[a]+=sz[g2[a][i]];
	}
	sz[a]++;
	// printf("%d: %lld\n",a,ret);
	return ret;
}
int tl[210000];
int n;
void dfs2(int a,int b){
	for(int i=0;i<g[a].size();i++){
		if(b==g[a][i])continue;
		// printf("%d %d\n",a,g[a][i]);
		tl[g[a][i]]=a;
		if(g[a][i]==cy[0])continue;
		dfs2(g[a][i],a);
	}
}
int main(){
	int a;scanf("%d",&a);
	n=a;
	for(int i=0;i<a*2;i++)UF[i]=-1;
	for(int i=0;i<a*2;i++)v[i]=-2;
	init_C(410000);
	for(int i=0;i<2*a;i++){
		scanf("%d%d",x+i,y+i);
		x[i]--;y[i]--;
		g[x[i]].push_back(y[i]+a);
		g[y[i]+a].push_back(x[i]);
		UNION(x[i],y[i]+a);
	}
	for(int i=0;i<a*2;i++){
		h[FIND(x[i])]++;
	}
	long long ret=1;
	long long ks=1;
	int tn=0;
	for(int i=0;i<a*2;i++){
		if(UF[i]<0&&h[i]!=-UF[i]){
			printf("0\n");return 0;
		}
		if(UF[i]<0){
			cur.clear();
			ks=ks*Comb(tn+h[i],h[i])%mod;
			tn+=h[i];
			cy.clear();
			dfs(i,-1);
			// for(int j=0;j<cy.size();j++)printf("%d ",cy[j]);printf("\n");
			long long tmp=0;
			for(int k=0;k<2;k++){

				if(k==0)dfs2(cy[0],cy[1]);
				else dfs2(cy[0],cy[cy.size()-1]);
				for(int j=0;j<cur.size();j++){
					for(int l=0;l<g[cur[j]].size();l++){
						int to=g[cur[j]][l];
						if(to<tl[cur[j]]){
							g2[cur[j]].push_back(to);
							// printf("%d %d\n",cur[j],to);
							ou[to]++;
						}
					}
				}
				for(int j=0;j<cur.size();j++){
					if(!ou[cur[j]]){
						g2[a*2].push_back(cur[j]);
					}
				}
				tmp+=calc(a*2);

				for(int j=0;j<cur.size();j++){
					g2[cur[j]].clear();
					ou[cur[j]]=0;
				}
				g2[a*2].clear();
			}
			tmp%=mod;
			// printf("%d: %lld\n",i,tmp);
			ret=ret*tmp%mod;
		}
	}

	ret=ret*ks%mod;
	printf("%lld\n",ret);
}

Measured memory_bytes: 45304