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)