#include #include #include #include #include #include #include using namespace std; #define int long long #define endl "\n" constexpr long long INF = (long long)1e18; constexpr long long MOD = 1'000'000'007; struct fast_io { fast_io(){ std::cin.tie(nullptr); std::ios::sync_with_stdio(false); }; } fio; class binomial_coefficients { long long MAX_VAL; public: vector fac, mmi; binomial_coefficients(){ } binomial_coefficients(long long num){ init(num); } ~binomial_coefficients(){ } void init(long long num){ MAX_VAL = num+1; fac.resize(MAX_VAL); mmi.resize(MAX_VAL); factorial_mod(); modular_multiplicatibe_inverse(); } void factorial_mod(){ fac[0] = 1; for(long long i = 1; i < MAX_VAL; fac[i] %= MOD, i++) fac[i] = fac[i - 1] * (i % MOD); } long long power(long long x, long long n){ long long ans = 1; for(;n;n >>= 1, x *= x, ans %= MOD, x %= MOD) if(n&1)ans*=x; return ans % MOD; } void exgcd(long long a, long long b, long long &x, long long &y){ if(b == 0){ x = 1; y = 0; return ; } exgcd(b, a % b, y, x); y -= a / b * x; } void modular_multiplicatibe_inverse(){ long long x, y; exgcd(fac[MAX_VAL - 1], MOD, x, y); mmi[MAX_VAL-1] = x; // mmi[MAX_VAL-1] = power(fac[MAX_VAL-1], MOD-2); for(long long i = MAX_VAL - 2; i >= 0; mmi[i]%=MOD, i--) mmi[i] = mmi[i + 1] * ((i + 1) % MOD); } long long combination(long long n, long long r){ return n < r ? 0 :fac[n] * (mmi[r] * mmi[n-r] % MOD) % MOD; } }; signed main(){ cout<>N>>M>>K; BC.init((N+M) * 2); for(int i = 2; i <= N + M; i += 2){ int A = 0, B = 0; int a = (i) / 2, b = (i) / 2; if(N + M - i < K) continue; if(N - a >= 0) A = BC.combination(N - a + (a - 1), a - 1) ; if(M - b >= 0) B = BC.combination(M - b + (b - 1), b - 1); ans += A * B % MOD * BC.power(i/2, MOD-2) % MOD; ans %= MOD; } ans *= BC.fac[N]; ans %= MOD; ans *= BC.fac[M]; ans %= MOD; cout<