#include #include #include using namespace std; using i32 = int32_t; using i64 = int64_t; #define rep(i,n) for(int i=0; i<(n); i++) const i32 MOD = 998244353; struct imod{ i32 val; imod(i32 v=0){ val = v; } imod operator+(const imod& r) const { imod res; res.val = (val + r.val) % MOD; return res; } imod operator*(const imod& r) const { imod res; res.val = (i64) val * r.val % MOD; return res; } }; vector F, I, iF; void initComb(int Z){ F.assign(Z+1,1); for(int i=1; i<=Z; i++) F[i] = (i64)F[i-1] * i % MOD; I.assign(Z+1,1); for(int i=2; i<=Z; i++) I[i] = MOD - (i64)MOD / i * I[MOD % i] % MOD; iF.assign(Z+1,1); for(int i=1; i<=Z; i++) iF[i] = (i64)iF[i-1] * I[i] % MOD; } i64 Comb(int n, int r){ return (i64)F[n] * iF[r] % MOD * iF[n-r] % MOD; } int N,M,K; int main(){ cin >> N >> M >> K; initComb(N+M+K); imod ans = 0; // o-o o-o o o o o if(N >= K-2 && M >= K && K >= 4 && N >= K-2 && M >= K) ans = ans + imod(Comb(N,K-2)) * Comb(M,K) * Comb(K-2,2) * F[K] * I[2] * I[2]; // o=o o=o o o o o if(N >= K && M >= K-2 && K >= 4 && N >= K && M >= K-2) ans = ans + imod(Comb(N,K)) * Comb(M,K-2) * Comb(K-2,2) * F[K] * I[2] * I[2]; cerr << ans.val << endl; // o-o-o o o o o o if(N >= K-2 && M >= K && K >= 3 && N >= K-2 && M >= K) ans = ans + imod(Comb(N,K-2)) * Comb(M,K) * Comb(K-2,1) * F[K] * I[2] * I[3]; // o=o=o o o o o o if(N >= K && M >= K-2 && K >= 3 && N >= K && M >= K-2) ans = ans + imod(Comb(N,K)) * Comb(M,K-2) * Comb(K-2,1) * F[K] * I[2] * I[3]; cerr << ans.val << endl; // o-o-o=o o o o o if(N >= K-2 && M >= K-1 && K >= 4 && N >= K-2 && M >= K-1) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-1) * Comb(K-2,1) * Comb(K-1,1) * (K-3) * F[K-2] * I[2]; // o=o=o-o o o o o if(N >= K-1 && M >= K-2 && K >= 4 && N >= K-1 && M >= K-2) ans = ans + imod(Comb(N,K-1)) * Comb(M,K-2) * Comb(K-2,1) * Comb(K-1,1) * (K-3) * F[K-2] * I[2]; cerr << ans.val << endl; // o-o=o o o o o // | // o=o if(N >= K-2 && M >= K-2 && K >= 5 && N >= K-2 && M >= K-2) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-2) * Comb(K-2,1) * Comb(K-2,2) * (K-3) * (K-4) * F[K-4]; // o=o=o o o o o // | | // o o if(N >= K-2 && M >= K-2 && K >= 5 && N >= K-2 && M >= K-2) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-2) * Comb(K-2,1) * Comb(K-2,2) * (K-3) * (K-4) * F[K-4]; cerr << ans.val << endl; // o-o o=o o o o o if(N >= K-1 && M >= K-1 && K >= 4 && N >= K-1 && M >= K-1) ans = ans + imod(Comb(N,K-1)) * Comb(M,K-1) * Comb(K-1,1) * Comb(K-1,1) * Comb(K-2,2) * F[K-2] * I[2]; cerr << ans.val << endl; // o-o=o-o o o o o if(N >= K-2 && M >= K-1 && K >= 4 && N >= K-2 && M >= K-1) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-1) * Comb(K-2,2) * Comb(K-1,1) * F[K-2]; // o=o-o=o o o o o if(N >= K-1 && M >= K-2 && K >= 4 && N >= K-1 && M >= K-2) ans = ans + imod(Comb(N,K-1)) * Comb(M,K-2) * Comb(K-2,2) * Comb(K-1,1) * F[K-2]; cerr << ans.val << endl; // o=o o o o o // | | // o=o if(N >= K-2 && M >= K-2 && K >= 4 && N >= K-2 && M >= K-2) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-2) * Comb(K-2,2) * Comb(K-2,2) * F[K-4]; cerr << ans.val << endl; // o-o=o-o=o o o o if(N >= K-2 && M >= K-2 && K >= 5 && N >= K-2 && M >= K-2) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-2) * Comb(K-2,2) * Comb(K-2,2) * 4 * (K-4) * F[K-4]; cerr << ans.val << endl; // o-o=o o-o=o o o if(N >= K-2 && M >= K-2 && K >= 6 && N >= K-2 && M >= K-2) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-2) * Comb(K-2,2) * Comb(K-2,2) * 2 * (K-4) * (K-5) * F[K-4]; cerr << ans.val << endl; // o-o o-o=o o o o if(N >= K-2 && M >= K-1 && K >= 5 && N >= K-2 && M >= K-1) ans = ans + imod(Comb(N,K-2)) * Comb(M,K-1) * Comb(K-2,2) * Comb(K-1,1) * 2 * (K-4) * F[K-2] * I[2]; // o=o o-o=o o o o if(N >= K-1 && M >= K-2 && K >= 5 && N >= K-1 && M >= K-2) ans = ans + imod(Comb(N,K-1)) * Comb(M,K-2) * Comb(K-2,2) * Comb(K-1,1) * 2 * (K-4) * F[K-2] * I[2]; cout << ans.val << "\n"; return 0; } struct ios_do_not_sync{ ios_do_not_sync(){ ios::sync_with_stdio(false); cin.tie(nullptr); } } ios_do_not_sync_instance;