結果
| 問題 | No.1973 Divisor Sequence |
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2026-09-16 00:35:47 |
| 言語 | PyPy3 (7.3.23 + ACL) |
| 結果 |
AC
不安定
|
| 実行時間 | 1,112 ms / 2,000 ms |
| + 158µs | |
| コード長 | 938 bytes |
| 記録 | |
| コンパイル時間 | 71 ms |
| コンパイル使用メモリ | 80,384 KB |
| 実行使用メモリ | 83,200 KB |
| 最終ジャッジ日時 | 2026-09-16 00:35:52 |
| 合計ジャッジ時間 | 4,871 ms |
|
ジャッジサーバーID (参考情報) |
judge2_1 / judge1_0 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 22 |
ソースコード
## https://yukicoder.me/problems/no/1973
import math
MOD = 10 ** 9 + 7
def solve(N, e):
dp = [0] * (e + 1)
dp[0] = 1
for _ in range(N):
new_dp = [0] * (e + 1)
for j in range(e + 1):
for k in range(e + 1 - j):
new_dp[k] += dp[j]
new_dp[k] %= MOD
dp = new_dp
ans = 0
for i in range(e + 1):
ans += dp[i]
ans %= MOD
return ans
def main():
N, M = map(int, input().split())
# 素因数分解
sqrt_m = int(math.sqrt(M))
p_map = {}
for p in range(2 ,sqrt_m + 1):
if M % p == 0:
p_map[p] = 0
while M % p == 0:
p_map[p] += 1
M //= p
if M > 1:
p_map[M] = 1
answer = 1
for value in p_map.values():
ans = solve(N, value)
answer *= ans
answer %= MOD
print(answer)
if __name__ == "__main__":
main()