結果
問題 | No.391 CODING WAR |
ユーザー |
|
提出日時 | 2020-06-20 14:09:42 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
WA
|
実行時間 | - |
コード長 | 850 bytes |
コンパイル時間 | 155 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 67,832 KB |
最終ジャッジ日時 | 2024-07-03 17:15:37 |
合計ジャッジ時間 | 19,931 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 7 WA * 9 |
ソースコード
import numpy as npN = 10 ** 6Nsq = 10 ** 3mod = 10 ** 9 + 7fac = np.arange(N, dtype=np.int64).reshape(Nsq, Nsq)fac[0, 0] = 1for i in range(1, Nsq):fac[:, i] *= fac[:, i - 1]fac[:, i] %= modfor i in range(1, Nsq):fac[i] *= fac[i - 1, -1]fac[i] %= modfac = 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] %= modfor i in range(1, Nsq):finv[i] *= finv[i - 1, -1]finv[i] %= modfinv = finv.ravel()[::-1]def nCr(a, b):comb = fac[a] * finv[b] % mod * finv[a - b] % modreturn combn, m = map(int, input().split())ans = pow(m, n, mod)for i in range(m, 1, -1):ans -= (-1)**(i + 1) * nCr(m, i - 1) * pow(i - 1, n, mod)ans %= modprint(ans)