結果
問題 | No.1581 Multiple Sequence |
ユーザー |
![]() |
提出日時 | 2021-07-02 23:13:30 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 398 ms / 2,000 ms |
コード長 | 744 bytes |
コンパイル時間 | 152 ms |
コンパイル使用メモリ | 82,116 KB |
実行使用メモリ | 86,012 KB |
最終ジャッジ日時 | 2024-06-29 13:01:35 |
合計ジャッジ時間 | 7,264 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 21 |
ソースコード
def primeFactor(N): i, n, ret, d, sq = 2, N, {}, 2, 99 while i <= sq: k = 0 while n % i == 0: n, k, ret[i] = n//i, k+1, k+1 if k > 0 or i == 97: sq = int(n**(1/2)+0.5) if i < 4: i = i * 2 - 1 else: i, d = i+d, d^6 if n > 1: ret[n] = 1 return ret def divisors(N): pf = primeFactor(N) ret = [1] for p in pf: ret_prev = ret ret = [] for i in range(pf[p]+1): for r in ret_prev: ret.append(r * (p ** i)) return sorted(ret) N = int(input()) X = [0] * (10 ** 5 + 1) X[0] = X[1] = 1 for i in range(2, N + 1): a = 0 for d in divisors(i): b = i // d - 1 a += X[b] X[i] = a print(X[N] % (10 ** 9 + 7))