結果
問題 | No.802 だいたい等差数列 |
ユーザー | ptotq |
提出日時 | 2020-04-18 14:09:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 158 ms / 2,000 ms |
コード長 | 604 bytes |
コンパイル時間 | 311 ms |
コンパイル使用メモリ | 82,916 KB |
実行使用メモリ | 111,596 KB |
最終ジャッジ日時 | 2024-04-14 20:48:54 |
合計ジャッジ時間 | 4,814 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 91 ms
92,108 KB |
testcase_01 | AC | 80 ms
91,796 KB |
testcase_02 | AC | 82 ms
92,512 KB |
testcase_03 | AC | 80 ms
89,748 KB |
testcase_04 | AC | 83 ms
91,956 KB |
testcase_05 | AC | 79 ms
90,308 KB |
testcase_06 | AC | 82 ms
90,792 KB |
testcase_07 | AC | 83 ms
91,704 KB |
testcase_08 | AC | 80 ms
91,732 KB |
testcase_09 | AC | 79 ms
90,668 KB |
testcase_10 | AC | 84 ms
92,316 KB |
testcase_11 | AC | 83 ms
91,776 KB |
testcase_12 | AC | 82 ms
92,820 KB |
testcase_13 | AC | 82 ms
91,260 KB |
testcase_14 | AC | 108 ms
107,168 KB |
testcase_15 | AC | 119 ms
111,596 KB |
testcase_16 | AC | 98 ms
98,724 KB |
testcase_17 | AC | 79 ms
90,388 KB |
testcase_18 | AC | 79 ms
90,484 KB |
testcase_19 | AC | 101 ms
105,492 KB |
testcase_20 | AC | 157 ms
107,072 KB |
testcase_21 | AC | 157 ms
106,696 KB |
testcase_22 | AC | 82 ms
91,948 KB |
testcase_23 | AC | 119 ms
107,048 KB |
testcase_24 | AC | 80 ms
90,408 KB |
testcase_25 | AC | 80 ms
91,188 KB |
testcase_26 | AC | 81 ms
92,292 KB |
testcase_27 | AC | 83 ms
91,128 KB |
testcase_28 | AC | 117 ms
107,296 KB |
testcase_29 | AC | 158 ms
106,944 KB |
testcase_30 | AC | 79 ms
90,472 KB |
testcase_31 | AC | 79 ms
90,364 KB |
testcase_32 | AC | 80 ms
89,624 KB |
testcase_33 | AC | 87 ms
94,164 KB |
ソースコード
import sys readline = sys.stdin.readline N, M, d1, d2 = map(int, input().split()) mod = 10**9 + 7 max_n = 2 * 10**6 fac, inv = [1]*(max_n+1), [0]*(max_n+1) for i in range(2, max_n+1): fac[i] = fac[i-1] * i % mod inv[-1] = pow(fac[-1], mod-2, mod) for i in range(max_n, 0, -1): inv[i-1] = inv[i] * i % mod d = d2 - d1 x = M - 1 - d1 * (N-1) if x < 0: print(0) exit() ans = 0 for i in range(N): sign = -1 if i & 1 else 1 if x+N-(d+1)*i >= N: ans += fac[x+N-(d+1)*i] * inv[N] * inv[x-(d+1)*i]\ * fac[N-1] * inv[i] * inv[N-1-i] * sign % mod print(ans % mod)