結果

問題 No.391 CODING WAR
ユーザー Ryushi-techRyushi-tech
提出日時 2020-06-20 14:12:00
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 1,130 ms / 2,000 ms
コード長 867 bytes
コンパイル時間 81 ms
コンパイル使用メモリ 12,672 KB
実行使用メモリ 67,848 KB
最終ジャッジ日時 2024-07-03 17:16:34
合計ジャッジ時間 18,863 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 568 ms
67,052 KB
testcase_01 AC 579 ms
66,916 KB
testcase_02 AC 595 ms
67,496 KB
testcase_03 AC 565 ms
67,172 KB
testcase_04 AC 567 ms
67,568 KB
testcase_05 AC 609 ms
67,500 KB
testcase_06 AC 570 ms
67,620 KB
testcase_07 AC 550 ms
67,108 KB
testcase_08 AC 557 ms
67,552 KB
testcase_09 AC 1,130 ms
67,192 KB
testcase_10 AC 1,041 ms
67,152 KB
testcase_11 AC 915 ms
67,236 KB
testcase_12 AC 576 ms
67,092 KB
testcase_13 AC 1,102 ms
67,848 KB
testcase_14 AC 1,011 ms
67,552 KB
testcase_15 AC 1,099 ms
67,400 KB
testcase_16 AC 886 ms
67,172 KB
testcase_17 AC 934 ms
67,404 KB
testcase_18 AC 817 ms
67,156 KB
testcase_19 AC 808 ms
67,236 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import numpy as np

N = 10 ** 6
Nsq = 10 ** 3
mod = 10 ** 9 + 7

fac = np.arange(N, dtype=np.int64).reshape(Nsq, Nsq)
fac[0, 0] = 1
for i in range(1, Nsq):
    fac[:, i] *= fac[:, i - 1]
    fac[:, i] %= mod
for i in range(1, Nsq):
    fac[i] *= fac[i - 1, -1]
    fac[i] %= mod
fac = fac.ravel()

finv = np.arange(1, N + 1, dtype=np.int64)[::-1].reshape(Nsq, Nsq)
finv[0, 0] = pow(int(fac[N - 1]), mod - 2, mod)
for i in range(1, Nsq):
    finv[:, i] *= finv[:, i - 1]
    finv[:, i] %= mod
for i in range(1, Nsq):
    finv[i] *= finv[i - 1, -1]
    finv[i] %= mod
finv = finv.ravel()[::-1]


def nCr(a, b):
    comb = fac[a] * finv[b] % mod * finv[a - b] % mod
    return comb


n, m = map(int, input().split())

ans = pow(m, n, mod)
sign = m & 1

for i in range(m, 1, -1):
    ans -= (-1)**(i + sign) * nCr(m, i - 1) * pow(i - 1, n, mod)
    ans %= mod
print(ans)
0