結果
| 問題 | No.895 MESE |
| ユーザー |
 jell
|
| 提出日時 | 2019-10-01 23:07:57 |
| 言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
| 結果 |
CE
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 21,883 bytes |
| コンパイル時間 | 18,982 ms |
| コンパイル使用メモリ | 434,932 KB |
| 最終ジャッジ日時 | 2024-11-14 21:43:37 |
| 合計ジャッジ時間 | 20,010 ms |
|
ジャッジサーバーID （参考情報） |
judge1 / judge5 |
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
ただし、clay言語の場合は開発者のデバッグのため、公開されます。
コンパイルメッセージ
main.cpp:635:25: in 'constexpr' expansion of 'math::fact_impl()'
main.cpp:635:25: error: 'constexpr' evaluation operation count exceeds limit of 33554432 (use '-fconstexpr-ops-limit=' to increase the limit)
635 | constexpr fact_impl _fact_impl;
| ^~~~~~~~~~
ソースコード
#ifdef stderr_path
#define LOCAL
#endif
#ifdef LOCAL
#define _GLIBCXX_DEBUG
#else
#pragma GCC optimize("Ofast")
#endif
#include <algorithm>
#include <bitset>
#include <cassert>
#include <chrono>
#include <complex>
#include <cstring>
#include <deque>
#include <functional>
#include <iomanip>
#include <iostream>
#include <map>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <unordered_map>
#include <unordered_set>
// #define NDEBUG
#define debug_stream std::cerr
#define iostream_untie true
#define __precision__ 10
#define all(v) std::begin(v), std::end(v)
#define rall(v) std::rbegin(v), std::rend(v)
#define __odd(n) ((n)&1)
#define __even(n) (not __odd(n))
#define __popcount(n) __builtin_popcountll(n)
#define __clz32(n) __builtin_clz(n)
#define __clz64(n) __builtin_clzll(n)
#define __ctz32(n) __builtin_ctz(n)
#define __ctz64(n) __builtin_ctzll(n)
using i64 = int_fast64_t;
using pii = std::pair<int, int>;
using pll = std::pair<int_fast64_t, int_fast64_t>;
template <class T>
using heap = std::priority_queue<T>;
template <class T>
using minheap = std::priority_queue<T, std::vector<T>, std::greater<T>>;
template <class T>
constexpr T inf = std::numeric_limits<T>::max() / T(2) - T(1123456);
namespace execution
{
std::chrono::system_clock::time_point start_time, end_time;
void print_elapsed_time()
{
end_time = std::chrono::system_clock::now();
std::cerr << "\n----- Exec time : ";
std::cerr << std::chrono::duration_cast<std::chrono::milliseconds>(
end_time - start_time)
.count();
std::cerr << " ms -----\n\n";
}
struct setupper
{
setupper()
{
if(iostream_untie)
{
std::ios::sync_with_stdio(false);
std::cin.tie(nullptr);
}
std::cout << std::fixed << std::setprecision(__precision__);
#ifdef stderr_path
if(freopen(stderr_path, "a", stderr))
{
std::cerr << std::fixed << std::setprecision(__precision__);
}
#endif
#ifdef stdout_path
if(not freopen(stdout_path, "w", stdout))
{
freopen("CON", "w", stdout);
std::cerr << "Failed to open the stdout file\n\n";
}
std::cout << "";
#endif
#ifdef stdin_path
if(not freopen(stdin_path, "r", stdin))
{
freopen("CON", "r", stdin);
std::cerr << "Failed to open the stdin file\n\n";
}
#endif
#ifdef LOCAL
std::cerr << "----- stderr at LOCAL -----\n\n";
atexit(print_elapsed_time);
start_time = std::chrono::system_clock::now();
#else
fclose(stderr);
#endif
}
} __setupper;
} // namespace execution
class myclock_t
{
std::chrono::system_clock::time_point built_pt, last_pt;
int built_ln, last_ln;
std::string built_func, last_func;
bool is_built;
public:
explicit myclock_t() : is_built(false)
{
}
void build(int crt_ln, const std::string &crt_func)
{
is_built = true;
last_pt = built_pt = std::chrono::system_clock::now();
last_ln = built_ln = crt_ln, last_func = built_func = crt_func;
}
void set(int crt_ln, const std::string &crt_func)
{
if(is_built)
{
last_pt = std::chrono::system_clock::now();
last_ln = crt_ln, last_func = crt_func;
}
else
{
debug_stream << "[ " << crt_ln << " : " << crt_func << " ] "
<< "myclock_t::set failed (yet to be built!)\n";
}
}
void get(int crt_ln, const std::string &crt_func)
{
if(is_built)
{
std::chrono::system_clock::time_point crt_pt(
std::chrono::system_clock::now());
int64_t diff =
std::chrono::duration_cast<std::chrono::milliseconds>(crt_pt -
last_pt)
.count();
debug_stream << diff << " ms elapsed from"
<< " [ " << last_ln << " : " << last_func << " ]";
if(last_ln == built_ln) debug_stream << " (when built)";
debug_stream << " to"
<< " [ " << crt_ln << " : " << crt_func << " ]"
<< "\n";
last_pt = built_pt, last_ln = built_ln, last_func = built_func;
}
else
{
debug_stream << "[ " << crt_ln << " : " << crt_func << " ] "
<< "myclock_t::get failed (yet to be built!)\n";
}
}
};
#ifdef LOCAL
myclock_t __myclock;
#define build_clock() __myclock.build(__LINE__, __func__)
#define set_clock() __myclock.set(__LINE__, __func__)
#define get_clock() __myclock.get(__LINE__, __func__)
#else
#define build_clock() ((void)0)
#define set_clock() ((void)0)
#define get_clock() ((void)0)
#endif
namespace std
{
template <class RAitr>
void rsort(RAitr __first, RAitr __last)
{
sort(__first, __last, greater<>());
}
template <class T>
size_t hash_combine(size_t seed, T const &key)
{
return seed ^ (hash<T>()(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2));
}
template <class T, class U>
struct hash<pair<T, U>>
{
size_t operator()(pair<T, U> const &pr) const
{
return hash_combine(hash_combine(0, pr.first), pr.second);
}
};
template <class tuple_t, size_t index = tuple_size<tuple_t>::value - 1>
struct tuple_hash_calc
{
static size_t apply(size_t seed, tuple_t const &t)
{
return hash_combine(
tuple_hash_calc<tuple_t, index - 1>::apply(seed, t),
get<index>(t));
}
};
template <class tuple_t>
struct tuple_hash_calc<tuple_t, 0>
{
static size_t apply(size_t seed, tuple_t const &t)
{
return hash_combine(seed, get<0>(t));
}
};
template <class... T>
struct hash<tuple<T...>>
{
size_t operator()(tuple<T...> const &t) const
{
return tuple_hash_calc<tuple<T...>>::apply(0, t);
}
};
template <class T, class U>
istream &operator>>(std::istream &s, pair<T, U> &p)
{
return s >> p.first >> p.second;
}
template <class T, class U>
ostream &operator<<(std::ostream &s, const pair<T, U> p)
{
return s << p.first << " " << p.second;
}
template <class T>
istream &operator>>(istream &s, vector<T> &v)
{
for(T &e : v)
{
s >> e;
}
return s;
}
template <class T>
ostream &operator<<(ostream &s, const vector<T> &v)
{
bool is_front = true;
for(const T &e : v)
{
if(not is_front)
{
s << ' ';
}
else
{
is_front = false;
}
s << e;
}
return s;
}
template <class tuple_t, size_t index>
struct tupleos
{
static ostream &apply(ostream &s, const tuple_t &t)
{
tupleos<tuple_t, index - 1>::apply(s, t);
return s << " " << get<index>(t);
}
};
template <class tuple_t>
struct tupleos<tuple_t, 0>
{
static ostream &apply(ostream &s, const tuple_t &t)
{
return s << get<0>(t);
}
};
template <class... T>
ostream &operator<<(ostream &s, const tuple<T...> &t)
{
return tupleos<tuple<T...>, tuple_size<tuple<T...>>::value - 1>::apply(
s, t);
}
template <>
ostream &operator<<(ostream &s, const tuple<> &t)
{
return s;
}
string revstr(string str)
{
reverse(str.begin(), str.end());
return str;
}
} // namespace std
#ifdef LOCAL
#define dump(...) \
debug_stream << "[ " << __LINE__ << " : " << __FUNCTION__ << " ]\n", \
dump_func(#__VA_ARGS__, __VA_ARGS__)
template <class T>
void dump_func(const char *ptr, const T &x)
{
debug_stream << '\t';
for(char c = *ptr; c != '\0'; c = *++ptr)
{
if(c != ' ') debug_stream << c;
}
debug_stream << " : " << x << '\n';
}
template <class T, class... rest_t>
void dump_func(const char *ptr, const T &x, rest_t... rest)
{
debug_stream << '\t';
for(char c = *ptr; c != ','; c = *++ptr)
{
if(c != ' ') debug_stream << c;
}
debug_stream << " : " << x << ",\n";
dump_func(++ptr, rest...);
}
#else
#define dump(...) ((void)0)
#endif
template <class P>
void read_range(P __first, P __second)
{
for(P i = __first; i != __second; ++i)
std::cin >> *i;
}
template <class P>
void write_range(P __first, P __second)
{
for(P i = __first; i != __second;
std::cout << (++i == __second ? '\n' : ' '))
{
std::cout << *i;
}
}
// substitute y for x.
template <class T>
void subst(T &x, const T &y)
{
x = y;
}
// substitue y for x iff x > y.
template <class T>
bool chmin(T &x, const T &y)
{
return x > y ? x = y, true : false;
}
// substitue y for x iff x < y.
template <class T>
bool chmax(T &x, const T &y)
{
return x < y ? x = y, true : false;
}
template <class T>
constexpr T minf(const T &x, const T &y)
{
return std::min(x, y);
}
template <class T>
constexpr T maxf(const T &x, const T &y)
{
return std::max(x, y);
}
// binary search.
template <class int_t, class F>
int_t bin(int_t ok, int_t ng, const F &f)
{
while(std::abs(ok - ng) > 1)
{
int_t mid = (ok + ng) / 2;
(f(mid) ? ok : ng) = mid;
}
return ok;
}
// be careful that val is type-sensitive.
template <class T, class A, size_t N>
void init(A (&array)[N], const T &val)
{
std::fill((T *)array, (T *)(array + N), val);
}
void reset()
{
}
template <class A, class... rest_t>
void reset(A &array, rest_t... rest)
{
memset(array, 0, sizeof(array));
reset(rest...);
}
// a integer uniformly and randomly chosen from the interval [l, r).
template <typename int_t>
int_t rand_int(int_t l, int_t r)
{
static std::random_device seed_gen;
static std::mt19937 engine(seed_gen());
std::uniform_int_distribution<int_t> unid(l, r - 1);
return unid(engine);
}
// a real number uniformly and randomly chosen from the interval [l, r).
template <typename real_t>
real_t rand_real(real_t l, real_t r)
{
static std::random_device seed_gen;
static std::mt19937 engine(seed_gen());
std::uniform_real_distribution<real_t> unid(l, r);
return unid(engine);
}
/* The main code follows. */
namespace math
{
template <int mod>
struct modint
{
int rep;
constexpr modint() : rep(0)
{}
constexpr modint(int_fast64_t y)
: rep(y >= 0 ? y % mod : (mod - (-y) % mod) % mod)
{}
constexpr modint &operator+=(const modint &p)
{
return (rep += p.rep) < mod ? 0 : rep -= mod, *this;
}
constexpr modint &operator++()
{
return ++rep, *this;
}
constexpr modint operator++(int)
{
modint t = *this;
return ++rep, t;
}
constexpr modint &operator-=(const modint &p)
{
return (rep += mod - p.rep) < mod ? 0 : rep -= mod, *this;
}
constexpr modint &operator--()
{
return --rep, *this;
}
constexpr modint operator--(int)
{
modint t = *this;
return --rep, t;
}
constexpr modint &operator*=(const modint &p)
{
return rep = (int_fast64_t)rep * p.rep % mod, *this;
}
constexpr modint &operator/=(const modint &p)
{
return *this *= inverse(p);
}
// constexpr modint &operator%=(int m) { return rep %= m, *this; }
constexpr modint operator-() const
{
return modint(-rep);
}
constexpr modint operator+(const modint &p) const
{
return modint(*this) += p;
}
constexpr modint operator-(const modint &p) const
{
return modint(*this) -= p;
}
constexpr modint operator*(const modint &p) const
{
return modint(*this) *= p;
}
constexpr modint operator/(const modint &p) const
{
return modint(*this) /= p;
}
// constexpr modint operator%(int m) const { return modint(*this) %= m;
// }
constexpr bool operator==(const modint &p) const
{
return rep == p.rep;
}
constexpr bool operator!=(const modint &p) const
{
return rep != p.rep;
}
constexpr bool operator!() const
{
return !rep;
}
// constexpr bool operator>(const modint &p) const { return rep > p.rep; }
// constexpr bool operator<(const modint &p) const { return rep < p.rep; }
// constexpr bool operator>=(const modint &p) const { return rep >= p.rep; }
// constexpr bool operator<=(const modint &p) const { return rep <= p.rep; }
constexpr friend modint<mod> inverse(const modint<mod> &p)
{
int a = p.rep, b = mod, u = 1, v = 0;
while(b > 0)
{
int t = a / b;
a -= t * b;
a ^= b ^= a ^= b;
u -= t * v;
u ^= v ^= u ^= v;
}
return modint(u);
}
constexpr friend modint pow(modint p, int_fast64_t e)
{
if(e < 0) e = (e % (mod - 1) + mod - 1) % (mod - 1);
modint ret = 1;
while(e)
{
if(e & 1) ret *= p;
p *= p;
e >>= 1;
}
return ret;
}
friend std::ostream &operator<<(std::ostream &s, const modint &p)
{
return s << p.rep;
}
friend std::istream &operator>>(std::istream &s, modint &p)
{
int_fast64_t rep;
p = modint((s >> rep, rep));
return s;
}
};
} // namespace math
// require modint
namespace math
{
constexpr int mod = 1000000007;
constexpr size_t N = 2.5e5, ext_N = 2.2e6;
class fact_impl
{
int _fact[N], _invfact[N], _inv[N];
static int _ext_fact[], _ext_invfact[], _ext_inv[];
static size_t fact_itr, invfact_itr, inv_itr;
public:
constexpr fact_impl() : _fact(), _invfact(), _inv()
{
_fact[0] = 1;
for(size_t i = 1; i < N; ++i)
_fact[i] = (int_fast64_t)_fact[i - 1] * i % mod;
_inv[1] = 1;
for(size_t i = 2; i < N && i < mod; ++i)
_inv[i] = mod - (int_fast64_t)_inv[mod % i] * (mod / i) % mod;
_invfact[0] = 1;
for(size_t i = 1; i < N && i < mod; ++i)
_invfact[i] = (int_fast64_t)_invfact[i - 1] * _inv[i] % mod;
}
int ext_fact(int x) const
{
assert(x < (int)ext_N);
if(x < 0) return 0;
if(not fact_itr)
{
_ext_fact[N - 1] = _fact[N - 1];
fact_itr = N;
}
for(size_t &i = fact_itr; (int)i <= x; ++i)
{
_ext_fact[i] = (int_fast64_t)_ext_fact[i - 1] * i % mod;
}
return x < (int)N ? _fact[x] : _ext_fact[x];
}
int ext_inv(int x) const
{
assert(x < (int)ext_N and x % mod);
for(size_t &i = inv_itr, r; (int)i <= x; ++i)
{
r = mod % i,
_ext_inv[i] = mod - (r < N ? _inv[r] : _ext_inv[r]) *
int_fast64_t(mod / i) % mod;
}
return x < (int)N ? _inv[x] : _ext_inv[x];
}
int ext_invfact(int x) const
{
assert(x < (int)ext_N and x < mod);
if(not invfact_itr)
{
_ext_invfact[N - 1] = _invfact[N - 1];
invfact_itr = N;
}
for(size_t &i = inv_itr, r; (int)i <= x; ++i)
{
r = mod % i,
_ext_inv[i] = mod - (r < N ? _inv[r] : _ext_inv[r]) *
int_fast64_t(mod / i) % mod;
}
for(size_t &i = invfact_itr; (int)i <= x; ++i)
{
_ext_invfact[i] = i < N
? _invfact[i]
: _ext_invfact[i - 1] * _ext_inv[i] % mod;
}
return x < (int)N ? _invfact[x] : _ext_invfact[x];
}
};
constexpr fact_impl _fact_impl;
int fact_impl::_ext_fact[ext_N], fact_impl::_ext_inv[ext_N],
fact_impl::_ext_invfact[ext_N];
size_t fact_impl::fact_itr, fact_impl::inv_itr = N, fact_impl::invfact_itr;
template <>
struct modint<mod>
{
int rep;
constexpr modint() : rep(0)
{}
constexpr modint(int_fast64_t y)
: rep(y >= 0 ? y % mod : (mod - (-y) % mod) % mod)
{}
constexpr modint &operator+=(const modint &p)
{
return (rep += p.rep) < mod ? 0 : rep -= mod, *this;
}
constexpr modint &operator++()
{
return ++rep, *this;
}
constexpr modint operator++(int)
{
modint t = *this;
return ++rep, t;
}
constexpr modint &operator-=(const modint &p)
{
return (rep += mod - p.rep) < mod ? 0 : rep -= mod, *this;
}
constexpr modint &operator--()
{
return --rep, *this;
}
constexpr modint operator--(int)
{
modint t = *this;
return --rep, t;
}
constexpr modint &operator*=(const modint &p)
{
return rep = (int_fast64_t)rep * p.rep % mod, *this;
}
modint &operator/=(const modint &p)
{
return *this *= inverse(p);
}
// constexpr modint &operator%=(int m) { return rep %= m, *this; }
constexpr modint operator-() const
{
return modint(-rep);
}
constexpr modint operator+(const modint &p) const
{
return modint(*this) += p;
}
constexpr modint operator-(const modint &p) const
{
return modint(*this) -= p;
}
constexpr modint operator*(const modint &p) const
{
return modint(*this) *= p;
}
modint operator/(const modint &p) const
{
return modint(*this) /= p;
}
// constexpr modint operator%(int m) const
// {
// return modint(*this) %= m;
// }
constexpr bool operator==(const modint &p) const
{
return rep == p.rep;
}
constexpr bool operator!=(const modint &p) const
{
return rep != p.rep;
}
constexpr bool operator!() const
{
return !rep;
}
// constexpr bool operator>(const modint &p) const { return rep > p.rep; }
// constexpr bool operator<(const modint &p) const { return rep < p.rep; }
// constexpr bool operator>=(const modint &p) const { return rep >= p.rep; }
// constexpr bool operator<=(const modint &p) const { return rep <= p.rep; }
friend modint inverse(const modint &p)
{
if(p.rep <= (int)ext_N) return _fact_impl.ext_inv(p.rep);
int a = p.rep, b = mod, u = 1, v = 0;
while(b > 0)
{
int t = a / b;
a -= t * b;
a ^= b ^= a ^= b;
u -= t * v;
u ^= v ^= u ^= v;
}
return modint(u);
}
constexpr friend modint pow(modint p, int_fast64_t e)
{
if(e < 0) e = (e % (mod - 1) + mod - 1) % (mod - 1);
modint ret = 1;
while(e)
{
if(e & 1) ret *= p;
p *= p;
e >>= 1;
}
return ret;
}
friend std::ostream &operator<<(std::ostream &s, const modint &p)
{
return s << p.rep;
}
friend std::istream &operator>>(std::istream &s, modint &p)
{
int_fast64_t rep;
p = modint((s >> rep, rep));
return s;
}
};
modint<mod> fact(int x)
{
return modint<mod>(_fact_impl.ext_fact(x));
}
modint<mod> invfact(int x)
{
return modint<mod>(_fact_impl.ext_invfact(x));
}
modint<mod> binom(int x, int y)
{
if(y < 0 or x < y) return 0;
return fact(x) * invfact(y) * invfact(x - y);
}
modint<mod> perm(int x, int y)
{
return binom(x, y) * fact(y);
}
} // namespace math
using namespace std;
using namespace math;
signed main()
{
void __solve();
void __precalc();
unsigned int t = 1;
// cin >> t;
// __precalc();
#ifdef LOCAL
t = 3;
#endif
while(t--)
{
__solve();
}
}
using mint=modint<mod>;
void __solve()
{
int a, b, c;
cin >> a >> b >> c;
const int n=a+b+c;
mint ans=0;
mint po=pow(mint(2),c);
for(int i = c; i <= n-2; ++i,po*=2)
{
ans+=(po-1)*binom(i-1,c-1)*binom(i-c,b-1);
}
std::cout << ans << "\n";
}