結果
問題 |
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])