結果
問題 | No.895 MESE |
ユーザー |
|
提出日時 | 2021-11-07 08:56:53 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 169 ms / 2,000 ms |
コード長 | 819 bytes |
コンパイル時間 | 376 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 89,084 KB |
最終ジャッジ日時 | 2024-11-08 19:24:16 |
合計ジャッジ時間 | 4,721 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
a, b, c = map(int, input().split()) N = a + b + c MOD = 10 ** 9 + 7 MAX = 3 * 10 ** 5 + 10 # https://tjkendev.github.io/procon-library/python/math/factorial.html から拝借しています。 fact = [1] * (MAX + 1) rfact = [1] * (MAX + 1) r = 1 for i in range(1, MAX + 1): fact[i] = r = r * i % MOD rfact[MAX] = r = pow(fact[MAX], MOD - 2, MOD) for i in range(MAX, 0, -1): rfact[i - 1] = r = r * i % MOD def comb(n, k, error=0): if not 0 <= k <= n: return error return fact[n] * rfact[k] * rfact[n - k] % MOD ans = 0 p2 = 1 for i in range(N - 2): ct = (comb(a + b - 1, b) - comb(i + 1 - c, b)) % MOD val_sum = comb(i, c - 1) * p2 % MOD if i >= 1: val_sum += comb(i - 1, c - 2) * (p2 - 1) % MOD p2 *= 2 p2 %= MOD ans += ct * val_sum ans %= MOD print(ans)