結果
問題 | No.1973 Divisor Sequence |
ユーザー | ああいい |
提出日時 | 2022-06-12 09:43:07 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 483 ms / 2,000 ms |
コード長 | 1,214 bytes |
コンパイル時間 | 269 ms |
コンパイル使用メモリ | 81,500 KB |
実行使用メモリ | 84,172 KB |
最終ジャッジ日時 | 2023-10-24 02:33:53 |
合計ジャッジ時間 | 6,368 ms |
ジャッジサーバーID (参考情報) |
judge11 / judge13 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 87 ms
75,992 KB |
testcase_01 | AC | 84 ms
75,992 KB |
testcase_02 | AC | 207 ms
84,016 KB |
testcase_03 | AC | 109 ms
83,352 KB |
testcase_04 | AC | 170 ms
84,072 KB |
testcase_05 | AC | 116 ms
83,808 KB |
testcase_06 | AC | 159 ms
83,468 KB |
testcase_07 | AC | 123 ms
84,168 KB |
testcase_08 | AC | 147 ms
84,028 KB |
testcase_09 | AC | 166 ms
84,028 KB |
testcase_10 | AC | 200 ms
83,484 KB |
testcase_11 | AC | 116 ms
83,336 KB |
testcase_12 | AC | 163 ms
84,028 KB |
testcase_13 | AC | 133 ms
84,160 KB |
testcase_14 | AC | 170 ms
83,540 KB |
testcase_15 | AC | 180 ms
83,440 KB |
testcase_16 | AC | 183 ms
83,464 KB |
testcase_17 | AC | 198 ms
84,064 KB |
testcase_18 | AC | 122 ms
83,584 KB |
testcase_19 | AC | 218 ms
84,100 KB |
testcase_20 | AC | 123 ms
84,172 KB |
testcase_21 | AC | 193 ms
84,028 KB |
testcase_22 | AC | 183 ms
84,132 KB |
testcase_23 | AC | 383 ms
83,280 KB |
testcase_24 | AC | 483 ms
83,488 KB |
ソースコード
N,M = map(int,input().split()) C = 10 ** 6 + 2 P = 10 ** 9 + 7 dat = [0] * C ans = 1 for i in range(2,C): if dat[i] == 0: for j in range(2 * i,C,i): dat[j] = 1 if M % i == 0: M //= i count = 1 while M % i == 0: count += 1 M //= i dp = [1] * (count + 1) for _ in range(N-1): nx = [0] * (count + 1) S = [0] * (count + 1) S[0] = dp[0] for k in range(1,count + 1): S[k] = dp[k] + S[k-1] S[k] %= P for k in range(count + 1): nx[k] = S[count - k] dp = nx ans = ans * sum(dp) % P from functools import lru_cache @lru_cache def calc(count): dp = [1] * (count + 1) for _ in range(N - 1): nx = [0] * (count + 1) S = [0] * (count + 1) S[0] = dp[0] for k in range(1,count + 1): S[k] = dp[k] + S[k-1] S[k] %= P for k in range(count + 1): nx[k] = S[count - k] dp = nx return sum(dp) % P if M != 1: ans = ans * calc(1) % P print(ans)