from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline

N = int(input())
X = list(map(int,input().split()))
Y = list(map(int,input().split()))
mod = 998244353

Z = [(x,1) for x in X] + [(y,0) for y in Y]
Z.sort()
n = sum(flg for x,flg in Z[:N])
m = N - n

ans = 1
for i in range(1,n+1):
    ans = (ans*i)%mod

for i in range(1,m+1):
    ans = (ans*i)%mod
print(ans)