MOD = 998244353
N,M,K = map(int,input().split())
P = {}
for _ in range(K):
    x,y = map(int,input().split())
    if x not in P:
        P[x] = set()
    if y not in P:
        P[y] = set()
    P[y].add(x)
V = set(list(P.keys()))
listV = list(V)
L = len(V)
J = 0
for v in P:
    if v<=M:
        J += 1
dpf = [0]*(1<<L)
dpf[0] = 1
for i in range(1<<L):
    S = set()
    for k in range(L):
        if i>>k & 1:
            S.add(listV[k])
    Sc = V-S
    flag = True
    for v in S:
        if not(P[v] <= S):
            flag = False
            break
    if not flag:continue
    nextV = set()
    for v in Sc:
        if P[v] <= S:
            nextV.add(v)
    for v in nextV:
        indv = listV.index(v)
        dpf[i|1<<indv] = (dpf[i|1<<indv]+dpf[i])%MOD
dpg = [0]*(1<<L)
dpg[0] = 1
for i in range(1<<L):
    S = set()
    for k in range(L):
        if i>>k & 1:
            S.add(listV[k])
    Sc = V-S
    flag = True
    for v in S:
        if not(P[v] <= S):
            flag = False
            break
    if not flag:continue
    nextV = set()
    for v in Sc:
        if i==0 and v<=M:continue
        else:
            if P[v] <= S:
                nextV.add(v)
    for v in nextV:
        indv = listV.index(v)
        dpg[i|1<<indv] = (dpg[i|1<<indv]+dpg[i])%MOD
f = dpf[(1<<L)-1]
g = dpg[(1<<L)-1]
A = [1]*(N+1)
for i in range(2,N+1):
    A[i] = (A[i-1]*i)%MOD
B = [1]*(N+1)
B[N] = pow(A[N],MOD-2,MOD)
for i in range(N-1,1,-1):
    B[i] = (B[i+1]*(i+1))%MOD
def comb(n,k):
    if k<0 or k>n:return 0
    return A[n]*B[k]*B[n-k]%MOD
ans = (comb(N-1,L)*f*(N-L-(M-J))*A[N-L-1])%MOD+(comb(N-1,L-1)*g*A[N-L])%MOD
ans = ans%MOD
print(ans)