#include #ifndef MODULO_HPP #define MODULO_HPP #include namespace modulo { template class modint { int val; public: constexpr modint() : val{0} {} constexpr modint(long long x) : val((x %= mod) < 0 ? mod + x : x) {} constexpr long long get() const { return val; } constexpr modint &operator+=(const modint &other) { return (val += other.val) < mod ? 0 : val -= mod, *this; } constexpr modint &operator++() { return ++val, *this; } constexpr modint operator++(int) { modint t = *this; return ++val, t; } constexpr modint &operator-=(const modint &other) { return (val += mod - other.val) < mod ? 0 : val -= mod, *this; } constexpr modint &operator--() { return --val, *this; } constexpr modint operator--(int) { modint t = *this; return --val, t; } constexpr modint &operator*=(const modint &other) { return val = (long long)val * other.val % mod, *this; } constexpr modint &operator/=(const modint &other) { return *this *= inverse(other); } constexpr modint operator-() const { return modint(-val); } constexpr modint operator+(const modint &other) const { return modint(*this) += other; } constexpr modint operator-(const modint &other) const { return modint(*this) -= other; } constexpr modint operator*(const modint &other) const { return modint(*this) *= other; } constexpr modint operator/(const modint &other) const { return modint(*this) /= other; } constexpr bool operator==(const modint &other) const { return val == other.val; } constexpr bool operator!=(const modint &other) const { return val != other.val; } constexpr bool operator!() const { return !val; } friend constexpr modint inverse(const modint &other) { assert(other != 0); int a{mod}, b{other.val}, u{}, v{1}, t{}; while(b) t = a / b, a ^= b ^= (a -= t * b) ^= b, u ^= v ^= (u -= t * v) ^= v; return modint{u}; } friend constexpr modint pow(modint other, long long e) { if(e < 0) e = e % (mod - 1) + mod - 1; modint res{1}; while(e) { if(e & 1) res *= other; other *= other, e >>= 1; } return res; } friend std::ostream &operator<<(std::ostream &s, const modint &other) { return s << other.val; } friend std::istream &operator>>(std::istream &s, modint &other) { long long val; other = modint{(s >> val, val)}; return s; } }; // class modint template <> class modint<2> { bool val; public: constexpr modint(bool x = false) : val{x} {} constexpr modint(int x) : val(x & 1) {} constexpr modint(long long x) : val(x & 1) {} constexpr bool get() const { return val; } constexpr modint &operator+=(const modint &other) { return val ^= other.val, *this; } constexpr modint &operator++() { return val = !val, *this; } constexpr modint operator++(int) { modint t = *this; return val = !val, t; } constexpr modint &operator-=(const modint &other) { return val ^= other.val, *this; } constexpr modint &operator--() { return val = !val, *this; } constexpr modint operator--(int) { modint t = *this; return val = !val, t; } constexpr modint &operator*=(const modint &other) { return val &= other.val, *this; } constexpr modint &operator/=(const modint &other) { return *this; } constexpr modint operator-() const { return *this; } constexpr modint operator+(const modint &other) const { return val != other.val; } constexpr modint operator-(const modint &other) const { return val != other.val; } constexpr modint operator*(const modint &other) const { return val && other.val; } constexpr modint operator/(const modint &other) const { return *this; } constexpr bool operator==(const modint &other) const { return val == other.val; } constexpr bool operator!=(const modint &other) const { return val != other.val; } constexpr bool operator!() const { return !val; } operator bool() const { return val; } friend constexpr modint inverse(const modint &other) { return other; } friend constexpr modint pow(const modint &other, long long exp) { return other; } friend std::ostream &operator<<(std::ostream &os, const modint &other) { return os << other.val; } friend std::istream &operator>>(std::istream &is, modint &other) { long long val; other.val = (is >> val, val & 1); return is; } }; // class modint specialization } // namespace modulo #endif class binom_runtime { const int mod; int n; std::vector _fact, _inv, _invfact; public: binom_runtime(int _mod) : mod(_mod), n{2}, _fact{1, 1}, _inv{0, 1}, _invfact{1, 1} {} void build(int m) { if(m < n) return; int _n = n; while(m >= n) n <<= 1; _fact.resize(n), _inv.resize(n), _invfact.resize(n); for(int i = _n; i < n; ++i) { _fact[i] = (long long)_fact[i - 1] * i % mod; _inv[i] = mod - (long long)mod / i * _inv[mod % i] % mod; _invfact[i] = (long long)_invfact[i - 1] * _inv[i] % mod; } } long long fact(int x) { return build(x), x < 0 ? 0 : _fact[x]; } long long inv(int x) { return build(x), x < 0 ? 0 : _inv[x]; } long long invfact(int x) { return build(x), x < 0 ? 0 : _invfact[x]; } long long binom(int x, int y) { return fact(x) * invfact(y) % mod * invfact(x - y) % mod; } }; using modint = modulo::modint<1000000007>; using namespace std; main() { binom_runtime gen(1000000007); gen.build(2200010); int X,Y,Z; while(cin>>X>>Y>>Z) { modint ans,s{1}; for(int k=X+Y+Z,f=1; k>=0; --k,f^=1) { modint tmp=s; if(X) { tmp*=gen.binom(X+k-1,k-1); } if(Y) { tmp*=gen.binom(Y+k-1,k-1); } if(Z) { tmp*=gen.binom(Z+k-1,k-1); } ans+=tmp; s*=2; if(f) { s-=gen.binom(X+Y+Z+1,k); } else { s+=gen.binom(X+Y+Z+1,k); } } cout << ans << "\n"; } }