結果

問題 No.895 MESE
ユーザー H3PO4
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0