#ifdef MAIN bool __multi__ = 0; namespace XK { const ll mod = 998244353; struct mm { ll x; mm(ll x_ = 0) : x(x_ % mod) { if(x < 0) x += mod; } friend mm operator+(mm a, mm b) { return a.x + b.x; } friend mm operator-(mm a, mm b) { return a.x - b.x; } friend mm operator*(mm a, mm b) { return a.x * b.x; } friend mm operator/(mm a, mm b) { return a * b.inv(); } friend mm& operator+=(mm& a, mm b) { return a = a.x + b.x; } friend mm& operator-=(mm& a, mm b) { return a = a.x - b.x; } friend mm& operator*=(mm& a, mm b) { return a = a.x * b.x; } friend mm& operator/=(mm& a, mm b) { return a = a * b.inv(); } mm inv() const { return pow(mod - 2); } mm pow(ll b) const { mm a = *this, c = 1; while(b) { if(b & 1) c *= a; a *= a; b >>= 1; } return c; } }; struct Comb { private: vector fact; vector invf; public: Comb(int _n) { fact.resize(_n + 1); invf.resize(_n + 1); fact[0] = 1; Rep(i, 1, _n + 1) fact[i] = fact[i - 1] * i; invf[_n] = fact[_n].inv(); for(int i = _n; i > 0; i --) invf[i - 1] = invf[i] * i; } mm C(int n, int k) { if(k < 0 || k > n) return 0; return fact[n] * invf[k] * invf[n - k]; } }; void solve() { int N, M, K; cin >> N >> M >> K; Comb comb(N * N + 1); mm ans = 0; rep(x, N + 1) rep(y, N + 1) { if(x * y < M) continue; if(2 * N - x - y < K) continue; mm cur = comb.C(N, x) * comb.C(N, y) * comb.C(2 * N - x - y, K) * comb.C(x * y, M); if(K + x + y & 1) ans -= cur; else ans += cur; } cout << ans.x; } }; #else #include "cassert" #include "cmath" #include "cstdint" #include "cstdio" #include "cstdlib" #include "cstring" #include "algorithm" #include "bitset" #include "chrono" #include "complex" #include "deque" #include "functional" #include "iostream" #include "limits" #include "map" #include "numeric" #include "queue" #include "random" #include "set" #include "sstream" #include "string" #include "unordered_map" #include "unordered_set" #include "utility" #include "vector" #include "array" using namespace std; #define int long long using ll = long long; using ull = unsigned long long; const ll INF = 1ll << 60; const ll LINF = 0x1fffffffffffffff; const ll MINF = 0x7fffffffffff; template bool chmax(A& l, const B& r){ return r > l ? l = r, 1 : 0; } template bool chmin(A& l, const B& r){ return r < l ? l = r, 1 : 0; } #define sz(x) ssize(x) #define rep(i, a) for(ll i = 0; i < (a); i ++) #define Rep(i, a, b) for(ll i = (a); i < (b); i ++) #define rrep(i, a, b) for(ll i = (b); i --> (a); ) #define all(x) begin(x), end(x) #define fst first #define snd second #define pb push_back template using V = vector; template using AR = array; #define MAIN #include __FILE__ signed main() {ios::sync_with_stdio(0);cin.tie(0);fixed(cout).precision(12);int t = 1;if(__multi__) cin >> t;while(t --) XK::solve();} #endif