結果
| 問題 | No.1973 Divisor Sequence |
| コンテスト | |
| ユーザー |
学ぶマン
|
| 提出日時 | 2025-12-20 13:56:05 |
| 言語 | PyPy3 (7.3.17) |
| 結果 |
AC
|
| 実行時間 | 1,207 ms / 2,000 ms |
| コード長 | 846 bytes |
| 記録 | |
| コンパイル時間 | 362 ms |
| コンパイル使用メモリ | 82,896 KB |
| 実行使用メモリ | 364,920 KB |
| 最終ジャッジ日時 | 2025-12-20 13:56:15 |
| 合計ジャッジ時間 | 9,368 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 22 |
ソースコード
MOD = 10**9+7
from collections import Counter
from itertools import accumulate
N, M = map(int, input().split())
def prime_factorize(n): # O(√n)
a = []
while n % 2 == 0:
a.append(2)
n //= 2
f = 3
while f * f <= n:
if n % f == 0:
a.append(f)
n //= f
else:
f += 2
if n != 1:
a.append(n)
return a
cnt = Counter(prime_factorize(M))
ans = 1
for k, v in cnt.items():
dp = [[0] * (v + 1) for _ in range(N + 1)]
dp[0][0] = 1
# 累積和+もらう
for i in range(1, N + 1):
p = dp[i - 1]
r = [p[0]]
for j in range(v):
r.append((r[-1] + p[j + 1])% MOD)
for j in range(v + 1):
dp[i][j] = r[v - j]
sub = sum(dp[-1]) % MOD
#print(dp)
ans *= sub
ans %= MOD
print(ans)
学ぶマン