#include #include using namespace std; #define int long long const int mod = 998244353; const int max_fact = 9000001; vector fact(max_fact); vector fact_inve(max_fact,1); void precompute() { fact[0] = 1; for (int i = 1; i < max_fact; ++i) fact[i] = (fact[i - 1] * i) % mod; fact_inve[max_fact - 1] = 921885471; for (int i = max_fact - 2; i > 0; --i) fact_inve[i] = (fact_inve[i + 1] * (i + 1)) % mod; } long long Comb(int a, int b) { if (a < b || b < 0) return 0; return (fact[a] * fact_inve[b] % mod * fact_inve[a - b]) % mod; } signed main() { int N, M, K; cin >> N >> M >> K; precompute(); long long ans = 0; vector C(N + 1); vector CC(2 * N + 1); for (int i = 0; i <= N; ++i) C[i] = Comb(N, i); for (int i = 0; i <= 2 * N; ++i) CC[i] = Comb(i, K); for (int a = 0; a < N; ++a) { for (int b = 0; b < N; ++b) { if ((a + b) % 2 == K % 2) ans = (ans + (C[a] * C[b] % mod * Comb((N - a) * (N - b), M) % mod * CC[a + b]) % mod) % mod; else ans = (ans - (C[a] * C[b] % mod * Comb((N - a) * (N - b), M) % mod * CC[a + b]) % mod + mod) % mod; } } cout << ans << endl; return 0; }