結果
| 問題 |
No.1581 Multiple Sequence
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2022-01-11 22:44:37 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 418 ms / 2,000 ms |
| コード長 | 841 bytes |
| コンパイル時間 | 372 ms |
| コンパイル使用メモリ | 81,900 KB |
| 実行使用メモリ | 79,124 KB |
| 最終ジャッジ日時 | 2024-11-14 12:16:27 |
| 合計ジャッジ時間 | 7,651 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 21 |
ソースコード
from collections import defaultdict
import itertools
def Sieve_of_Eratosthenes(maxA):
lst = [-1]*(maxA+1)
lst[0] = 0
lst[1] = 1
for i in range(2, maxA+1):
if lst[i] == -1:
for g in range(i, maxA+1, i):
lst[g] = i
return lst
def factorization(n):
d = defaultdict(int)
now = n
while now != 1:
r = lst[now]
now //= r
d[r] += 1
divs = []
for l in itertools.product(*[range(v+1) for v in d.values()]):
div = 1
for idx, k in enumerate(d):
div *= k**l[idx]
divs.append(div)
return divs
m = int(input())
mod = 10**9+7
lst = Sieve_of_Eratosthenes(m+1)
dp = [0]*(m+1)
dp[0] = 1
for i in range(1, m+1):
f = factorization(i)
for j in f:
dp[i] += dp[i//j-1]
dp[i] %= mod
print(dp[m])