#include #ifndef MODULO_HPP #define MODULO_HPP #include #include template class modint { int val; public: constexpr modint() noexcept : val{0} {} constexpr modint(long long x) noexcept : val((x %= mod) < 0 ? mod + x : x) {} constexpr long long value() const noexcept { return val; } constexpr modint &operator+=(const modint &other) noexcept { return (val += other.val) < mod ? 0 : val -= mod, *this; } constexpr modint &operator++() noexcept { return ++val, *this; } constexpr modint operator++(int) noexcept { modint t = *this; return ++val, t; } constexpr modint &operator-=(const modint &other) noexcept { return (val += mod - other.val) < mod ? 0 : val -= mod, *this; } constexpr modint &operator--() noexcept { return --val, *this; } constexpr modint operator--(int) noexcept { modint t = *this; return --val, t; } constexpr modint &operator*=(const modint &other) noexcept { return val = (long long)val * other.val % mod, *this; } constexpr modint &operator/=(const modint &other) noexcept { return *this *= inverse(other); } constexpr modint operator-() const noexcept { return modint(-val); } constexpr modint operator+(const modint &other) const noexcept { return modint(*this) += other; } constexpr modint operator-(const modint &other) const noexcept { return modint(*this) -= other; } constexpr modint operator*(const modint &other) const noexcept { return modint(*this) *= other; } constexpr modint operator/(const modint &other) const noexcept { return modint(*this) /= other; } constexpr bool operator==(const modint &other) const noexcept { return val == other.val; } constexpr bool operator!=(const modint &other) const noexcept { return val != other.val; } constexpr bool operator!() const noexcept { return !val; } template friend constexpr modint operator+(T x, modint y) noexcept { return modint(x) + y; } template friend constexpr modint operator-(T x, modint y) noexcept { return modint(x) - y; } template friend constexpr modint operator*(T x, modint y) noexcept { return modint(x) * y; } template friend constexpr modint operator/(T x, modint y) noexcept { return modint(x) / y; } friend constexpr modint inverse(const modint &other) noexcept { 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) noexcept { 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) noexcept { return s << other.val; } friend std::istream &operator>>(std::istream &s, modint &other) noexcept { long long val; other = modint{(s >> val, val)}; return s; } }; // class modint template <> class modint<2> { bool val; public: constexpr modint(bool x = false) noexcept : val{x} {} constexpr modint(int x) noexcept : val(x & 1) {} constexpr modint(long long x) noexcept : val(x & 1) {} constexpr bool value() const noexcept { return val; } constexpr modint &operator+=(const modint &other) noexcept { return val ^= other.val, *this; } constexpr modint &operator++() noexcept { return val = !val, *this; } constexpr modint operator++(int) noexcept { modint t = *this; return val = !val, t; } constexpr modint &operator-=(const modint &other) noexcept { return val ^= other.val, *this; } constexpr modint &operator--() noexcept { return val = !val, *this; } constexpr modint operator--(int) { modint t = *this; return val = !val, t; } constexpr modint &operator*=(const modint &other) noexcept { return val &= other.val, *this; } constexpr modint &operator/=(const modint &other) noexcept { return *this; } constexpr modint operator-() const noexcept { return *this; } constexpr modint operator+(const modint &other) const noexcept { return val != other.val; } constexpr modint operator-(const modint &other) const noexcept { return val != other.val; } constexpr modint operator*(const modint &other) const noexcept { return val && other.val; } constexpr modint operator/(const modint &other) const noexcept { return *this; } constexpr bool operator==(const modint &other) const noexcept { return val == other.val; } constexpr bool operator!=(const modint &other) const noexcept { return val != other.val; } constexpr bool operator!() const noexcept { return !val; } operator bool() const noexcept { return val; } template friend constexpr modint operator+(T x, modint y) noexcept { return modint(x) + y; } template friend constexpr modint operator-(T x, modint y) noexcept { return modint(x) - y; } template friend constexpr modint operator*(T x, modint y) noexcept { return modint(x) * y; } template friend constexpr modint operator/(T x, modint y) noexcept { return modint(x) / y; } friend constexpr modint inverse(const modint &other) noexcept { return other; } friend constexpr modint pow(const modint &other, long long exp) noexcept { return other; } friend std::ostream &operator<<(std::ostream &os, const modint &other) noexcept { return os << other.val; } friend std::istream &operator>>(std::istream &is, modint &other) noexcept { long long val; other.val = (is >> val, val & 1); return is; } }; // class modint specialization #endif class binom_runtime { const int mod; int n; std::vector _fact, _inv, _invfact; public: binom_runtime(const int _mod) noexcept : mod(_mod), n{2}, _fact{1, 1}, _inv{0, 1}, _invfact{1, 1} {} void build(const int m) noexcept { 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(const int x) noexcept { return build(x), x < 0 ? 0 : _fact[x]; } long long inv(const int x) noexcept { return build(x), x < 0 ? 0 : _inv[x]; } long long invfact(const int x) noexcept { return build(x), x < 0 ? 0 : _invfact[x]; } long long binom(const int x, int y) noexcept { return fact(x) * invfact(y) % mod * invfact(x - y) % mod; } }; using mint = modint<1000000007>; using namespace std; main() { binom_runtime gen(1000000007); gen.build(2200010); int X,Y,Z; while(cin>>X>>Y>>Z) { mint ans,s{1}; for(int k=X+Y+Z,f=1; k>=0; --k,f^=1) { mint 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"; } }