結果
問題 | No.1973 Divisor Sequence |
ユーザー | qib |
提出日時 | 2023-01-13 00:39:59 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,663 ms / 2,000 ms |
コード長 | 932 bytes |
コンパイル時間 | 389 ms |
コンパイル使用メモリ | 86,876 KB |
実行使用メモリ | 350,772 KB |
最終ジャッジ日時 | 2023-08-25 04:30:08 |
合計ジャッジ時間 | 13,306 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 75 ms
71,220 KB |
testcase_01 | AC | 75 ms
71,480 KB |
testcase_02 | AC | 789 ms
244,944 KB |
testcase_03 | AC | 133 ms
84,560 KB |
testcase_04 | AC | 500 ms
174,972 KB |
testcase_05 | AC | 174 ms
97,120 KB |
testcase_06 | AC | 518 ms
181,564 KB |
testcase_07 | AC | 178 ms
93,284 KB |
testcase_08 | AC | 328 ms
130,028 KB |
testcase_09 | AC | 442 ms
162,932 KB |
testcase_10 | AC | 535 ms
187,268 KB |
testcase_11 | AC | 170 ms
90,216 KB |
testcase_12 | AC | 455 ms
160,308 KB |
testcase_13 | AC | 199 ms
99,380 KB |
testcase_14 | AC | 366 ms
138,504 KB |
testcase_15 | AC | 673 ms
223,320 KB |
testcase_16 | AC | 684 ms
220,260 KB |
testcase_17 | AC | 483 ms
169,448 KB |
testcase_18 | AC | 143 ms
81,568 KB |
testcase_19 | AC | 598 ms
198,048 KB |
testcase_20 | AC | 188 ms
95,888 KB |
testcase_21 | AC | 582 ms
205,784 KB |
testcase_22 | AC | 571 ms
196,240 KB |
testcase_23 | AC | 752 ms
284,600 KB |
testcase_24 | AC | 1,663 ms
350,772 KB |
ソースコード
import math MOD = 10 ** 9 + 7 n, m = map(int, input().split()) ans = 1 cur = m for x in range(2, int(math.sqrt(m)) + 1): if cur % x != 0: continue ei = 0 while cur % x == 0: ei += 1 cur //= x dp = [[0 for _ in range(ei + 1)] for _ in range(n)] for j in range(ei + 1): dp[0][j] = 1 for i in range(1, n): dp[i][ei] = dp[i - 1][0] for j in range(ei - 1, -1, -1): dp[i][j] = (dp[i][j + 1] + dp[i - 1][ei - j]) % MOD cnt = 0 for j in range(ei + 1): cnt = (cnt + dp[-1][j]) % MOD ans = (ans * cnt) % MOD if cur > 1: ei = 1 dp = [[0 for _ in range(ei + 1)] for _ in range(n)] for j in range(ei + 1): dp[0][j] = 1 for i in range(1, n): dp[i][ei] = dp[i - 1][0] for j in range(ei - 1, -1, -1): dp[i][j] = (dp[i][j + 1] + dp[i - 1][ei - j]) % MOD cnt = 0 for j in range(ei + 1): cnt = (cnt + dp[-1][j]) % MOD ans = (ans * cnt) % MOD print(ans)