#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define repr(i, a, b) for (int i = (int)(a); i < (int)(b); i++) #define rep(i, n) repr(i, 0, n) #define INF 2e9 #define MOD 1000000007 //#define MOD 998244353 #define LINF (long long)4e18 #define jck 3.141592 const double EPS = 1e-10; using ll = long long; using Pi = pair; using Pl = pair; ll fac[2000000]; ll finv[2000000]; ll inv[2000000]; void COMinit(){ fac[0] = fac[1] = 1; finv[0] = finv[1] = 1; inv[1] = 1; repr(i,2,2000000){ fac[i] = fac[i-1] * i % MOD; inv[i] = MOD - inv[MOD%i] * (MOD/i) % MOD; finv[i] = finv[i-1] * inv[i] % MOD; } } ll COM(int n, int k){ if(n < k) return 0; if(n < 0 || k < 0) return 0; return fac[n] * (finv[k] * finv[n-k] % MOD) % MOD; } ll modpow(ll a, ll n, ll m){ if(n == 0) return 1; ll half = modpow(a,n/2,m); ll res = half * half % m; if(n & 1) res = res * a % m; return res; } int main(){ COMinit(); int X,Y,Z; cin >> X >> Y >> Z; vector now(X+Y+Z); vector ans(X+Y+Z); ll res = 0; rep(i,X+Y+Z){ if(i >= 1){ repr(j,1,i+1){ now[i] += COM(i+1,j)*ans[j-1]; now[i] %= MOD; } } ans[i] += COM(X+i,i)*COM(Y+i,i)%MOD*COM(Z+i,i)%MOD-(i == 0?0:now[i]); ans[i] %= MOD; if(ans[i] < 0) ans[i] += MOD; res += ans[i]; res %= MOD; } cout << res << endl; }