結果
| 問題 |
No.1973 Divisor Sequence
|
| コンテスト | |
| ユーザー |
rlangevin
|
| 提出日時 | 2023-04-23 11:33:10 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 513 ms / 2,000 ms |
| コード長 | 859 bytes |
| コンパイル時間 | 884 ms |
| コンパイル使用メモリ | 82,176 KB |
| 実行使用メモリ | 76,416 KB |
| 最終ジャッジ日時 | 2024-11-07 19:14:00 |
| 合計ジャッジ時間 | 5,888 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 22 |
ソースコード
def factorization(n):
arr = []
temp = n
for i in range(2, int(-(-n**0.5//1))+1):
if temp%i==0:
cnt=0
while temp%i==0:
cnt+=1
temp //= i
arr.append([i, cnt])
if temp!=1:
arr.append([temp, 1])
return arr
def solve(K):
pre = [0] * (K + 1)
pre[0] = 1
Acp = [1] * (K + 2)
Acp[0] = 0
for i in range(N):
dp = [0] * (K + 1)
Ac = [0] * (K + 2)
for j in range(K + 1):
dp[j] += Acp[K + 1 - j]
Ac[j + 1] = Ac[j] + dp[j]
dp[j] %= mod
Ac[j + 1] %= mod
dp, pre = pre, dp
Ac, Acp = Acp, Ac
return sum(pre) % mod
N, M = map(int, input().split())
mod = 10 ** 9 + 7
ans = 1
for p, m in factorization(M):
ans *= solve(m)
ans %= mod
print(ans)
rlangevin