N=int(input()) X=list(map(int, input().split())) Y=list(map(int, input().split())) Z=[] for i in range(N): Z.append((X[i],1,i)) Z.append((Y[i],2,i)) Z=sorted(Z) x,y=0,0 for _ ,a,i in Z[:N]: if a==1: x+=1 else: y+=1 D=[1] mod=998244353 for i in range(1,10**5+10): D.append(D[-1]*i) D[-1]%=mod ans=D[x]*D[y]%mod print(ans)