M=998244353
n,m,k=map(int,input().split())
F=[c:=1]+[c:=c*-~i%M for i in range(n*n+3)]
iF=([c:=pow(F[-1],-1,M)]+[c:=c*i%M for i in range(n*n+3,0,-1)])[::-1]
C=lambda n,k:n>=k and F[n]*iF[n-k]*iF[k]%M
nC=[C(n,i)for i in range(n+1)]
Ck=[C(i,k)for i in range(2*n+2)]
r=0
for i in range(1,n+1):
    for j in range(1,n+1):
        d=C(i*j,m)*nC[i]*nC[j]*Ck[n+n-i-j]%M
        if i+j+k&1:r-=d
        else:r+=d
print(r%M)