#include using namespace std; typedef long long int ll; typedef pair P; typedef vector VI; typedef vector VVI; #define REP(i,n) for(int i=0;i<(n);i++) #define ALL(v) v.begin(),v.end() template bool chmax(T &x, const T &y){return (x bool chmin(T &x, const T &y){return (x>y)?(x=y,true):false;}; constexpr ll MOD=998244353; constexpr ll INF=2e18; long long inv(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(MOD-a[MOD%i]*(MOD/i)%MOD); } return a[n]; } long long fact(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(a.back()*i%MOD); } return a[n]; } long long finv(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(a.back()*inv(i)%MOD); } return a[n]; } long long com(int n, int k){ if(n> n; VI x(n), y(n); REP(i,n) cin >> x[i]; REP(i,n) cin >> y[i]; sort(ALL(x)); sort(ALL(y),greater()); y.push_back(-1); int a=0, b=0; for(int i=n-1;i>=0;i--){ if(x[i]>y[i]) a++; if(y[i]!=y[i+1]&&x[i]<=y[i]) break; b++; } cout << com(b,a)*fact(a)%MOD*fact(n-a)%MOD << endl; return 0; }