#include #include #include #include #include #include #include #include #include #include #include #include #include #include #define fi first #define se second #define rep(i,n) for(ll i=0;i<(n);i++) #define rrep(i,n) for(ll i=(n)-1;i>=0;i--) #define orep(i,n) for(ll i=1;i<=(n);i++) #define nfor(i,s,n) for(ll i=(s);i<(n);i++) #define dfor(i,s,n) for(ll i=(s)-1;i>=n;i--) #define INF 9e18//9223372036854775807 #define all(a) a.begin(),a.end() #define rall(a) a.rbegin(),a.rend() #define chmax(x,y) x = max(x,y) #define chmin(x,y) x = min(x,y) #define pb push_back #define pob pop_back #define vc vector #define YES cout << "Yes" << endl; #define NO cout << "No" << endl; #define YN {cout << "Yes" << endl;}else{cout << "No" << endl;} #define dame cout << -1 << endl; #define vc_unique(v) v.erase(unique(v.begin(), v.end()), v.end()) #define vc_remove(v,x) v.erase(remove(v.begin(), v.end(),x), v.end()) #define vc_rotate(v) rotate(v.begin(), v.begin()+1, v.end()) #define pop_cnt(s) ll(popcount(uint64_t(s))) #define next_p(v) next_permutation(v.begin(),v.end()) #ifndef ONLINE_JUDGE #define _GLIBCXX_DEBUG #endif using namespace std; using ll = long long; using ld = long double; using pll = pair; using vvvvvl = vector > > > >; using vvvvl = vector > > >; using vvvl = vector > >; using vvl = vector >; using vl = vector; using Graph = vector >; template using pq = priority_queue,less >; template using pq_g = priority_queue,greater >; ll dx[4] = {0,1,0,-1};//vl dx = {1,1,0,-1,-1,-1,0,1}; ll dy[4] = {1,0,-1,0};//vl dy = {0,1,1,1,0,-1,-1,-1}; bool out_grid(ll i, ll j, ll h, ll w){//trueならcontinueする return(!(0 <= i && i < h && 0 <= j && j 0){ if((1ll << cnt) & y){ y ^= (1ll << cnt); ans *= a; ans %= MOD; } a = a*a; a %= MOD; cnt++; } return ans; } ll npow(ll x,ll y){ ll a = x; ll cnt = 0; ll ans = 1; while(y > 0){ if((1ll << cnt) & y){ y ^= (1ll << cnt); ans *= a; } a = a*a; cnt++; } return ans; } void print(ld x){printf("%.20Lf\n", x);} ////////////////////////////////////////////////////////////////// int main(){ ll N; cin >> N; vl X(N); vl Y(N); rep(i,N)cin >> X[i]; rep(i,N)cin >> Y[i]; vl Z; rep(i,N)Z.pb(X[i]); rep(i,N)Z.pb(Y[i]); sort(all(X)); sort(all(Y)); sort(all(Z)); ll cntX = 0; ll cntY = 0; rep(i,N){ if(X[cntX] == Z[i])cntX++; else cntY++; } ll ans = 1; while(cntX > 0){ ans *= cntX; ans %= MOD; cntX--; } while(cntY > 0){ ans *= cntY; ans %= MOD; cntY--; } cout << ans << endl; }