結果
問題 | No.823 Many Shifts Easy |
ユーザー |
![]() |
提出日時 | 2021-01-28 23:42:34 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 93 ms / 2,000 ms |
コード長 | 1,437 bytes |
コンパイル時間 | 392 ms |
コンパイル使用メモリ | 82,308 KB |
実行使用メモリ | 75,768 KB |
最終ジャッジ日時 | 2024-06-26 07:24:49 |
合計ジャッジ時間 | 1,856 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
class Combination:"""SIZEが10**6程度以下の二項係数を何回も呼び出したいときに使う使い方:comb = Combination(SIZE, MOD)comb(10, 3) => 120"""def __init__(self, N, MOD=10 ** 9 + 7):self.MOD = MODself.__make_factorial_list(N)def __call__(self, n, k):if k < 0 or k > n:return 0res = self.fact[n] * self.inv[k] % self.MODres = res * self.inv[n - k] % self.MODreturn resdef nPk(self, n, k):if k < 0 or k > n:return 0return self.fact[n] * self.inv[n - k] % self.MODdef nHk(self, n, k):if k == 0:return 1return self.__call__(n + k - 1, k)def __make_factorial_list(self, N):self.fact = [1] * (N + 1)self.inv = [1] * (N + 1)MOD = self.MODfor i in range(1, N + 1):self.fact[i] = (self.fact[i - 1] * i) % MODself.inv[N] = pow(self.fact[N], MOD - 2, MOD)for i in range(N, 0, -1):self.inv[i - 1] = (self.inv[i] * i) % MODreturnmod = 10 ** 9 + 7comb = Combination(10 ** 6, mod)inv2 = pow(2, mod - 2, mod)N, K = map(int, input().split())def sum(x): return x * (x + 1) * inv2 % modans = comb.nPk(N-1, K) * sum(N) % modif K >= 2:tmp = sum(N - 1) * comb(N - 2, K - 2) * comb.fact[K] * inv2 % modans = (ans + tmp) % modprint(ans)