結果

問題 No.665 Bernoulli Bernoulli
ユーザー H3PO4H3PO4
提出日時 2021-02-17 10:05:34
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 721 bytes
コンパイル時間 164 ms
コンパイル使用メモリ 82,156 KB
実行使用メモリ 83,068 KB
最終ジャッジ日時 2024-09-13 20:28:01
合計ジャッジ時間 3,955 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 63 ms
66,688 KB
testcase_01 AC 60 ms
61,184 KB
testcase_02 TLE -
testcase_03 -- -
testcase_04 -- -
testcase_05 -- -
testcase_06 -- -
testcase_07 -- -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

MOD = 10 ** 9 + 7
MAX = 10 ** 4 + 10

modinv_table = [pow(i, MOD - 2, MOD) for i in range(MAX)]

fac, inv = [1] * MAX, [1] * MAX
for i in range(1, MAX):
    fac[i] = fac[i - 1] * i % MOD
    inv[i] = pow(fac[i], MOD - 2, MOD)


def comb(n, k, error=0):
    if not 0 <= k <= n:
        return error
    return fac[n] * inv[k] * inv[n - k] % MOD


N, K = map(int, input().split())

B = [1] * (K + 1)
for i in range(1, K + 1):
    tmp = 0
    for k in range(i):
        tmp += comb(i + 1, k) * B[k]
        tmp %= MOD
    B[i] = -tmp * modinv_table[i + 1] % MOD

ans = 0
P = N + 1
for j in range(K + 1):
    ans += comb(K + 1, K - j) * B[K - j] * P
    P = P * (N + 1) % MOD
ans *= modinv_table[K + 1]
ans %= MOD
print(ans)
0