結果
問題 | No.802 だいたい等差数列 |
ユーザー |
![]() |
提出日時 | 2021-04-28 00:54:56 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 104 ms / 2,000 ms |
コード長 | 1,469 bytes |
コンパイル時間 | 308 ms |
コンパイル使用メモリ | 82,472 KB |
実行使用メモリ | 85,256 KB |
最終ジャッジ日時 | 2024-07-06 20:12:21 |
合計ジャッジ時間 | 3,943 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 30 |
ソースコード
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) % MODreturnN, M, D1, D2 = map(int, input().split())mod = 10**9+7comb = Combination(N+M, mod)ans = 0for i in range(N):power = D1*(N-1) + i*(D2-D1+1)need = M-1-powerif need < 0:continuecnt = comb(N-1, i) * comb(need + N, N) % modif i % 2 == 0:ans += cntelse:ans += mod - cntans %= modprint(ans)