I,R=input,range
J=lambda:map(int,I().split())
d=N=int(I())
X=sorted(J())
Y=sorted(J())
for i in R(N-1,-1,-1):
	if X[i]>Y[N-i-1]:d=i
f=[1]
P=998244353
for i in R(N):f+=[f[-1]*(i+1)%P]
print(f[d]*f[N-d]%P)