def f(a:list,x:int): a=set(a) while x in a:x+=1 return x def g(a:list,m:int): ret=1;mod=998244353 for i in range(1,m+1): ret*=f(a,i) ret%=mod return ret import itertools def naive(N,M,K): ret=0 mod=998244353 for v in itertools.product(range(1,M+1),repeat=N): ret+=g(list(v),K) ret%=mod return ret def fast(N,M,K): DP=[[[[0,0]for _ in range(89)]for _ in range(89)]for _ in range(89)] bi=[[0]*89 for _ in range(89)] bi[0][0]=1 mod=998244353 for i in range(85): for j in range(85): bi[i+1][j]+=bi[i][j] bi[i+1][j]%=mod bi[i+1][j+1]+=bi[i][j] bi[i+1][j+1]%=mod fit=lambda l,r,x:0 if x0:continue tmp=(DP[i][j][k][0]+DP[i][j][k][1])*pow(i,fit(i-k,i,K),mod) for l in range(i+1,K+1): tmp*=l tmp%=mod DP[i+1][0][0][1]+=tmp DP[i+1][0][0][1]%=mod return DP[M+2][0][0][1] def fast2(N,M,K): DP=[[[0]*89 for _ in range(89)]for _ in range(89)] bi=[[0]*89 for _ in range(89)] bi[0][0]=1 mod=998244353 for i in range(85): for j in range(85): bi[i+1][j]+=bi[i][j] bi[i+1][j]%=mod bi[i+1][j+1]+=bi[i][j] bi[i+1][j+1]%=mod fit=lambda l,r,x:0 if x0:continue tmp=(DP[i][j][k])*pow(i,fit(i-k,i,K),mod) for l in range(i+1,K+1): tmp*=l tmp%=mod DP[i+1][0][0]+=tmp DP[i+1][0][0]%=mod return DP[M+2][0][0] def main(): N,M,K=map(int,input().split()) print(fast2(N,M,K)) main()