結果
| 問題 |
No.3332 Consecutive Power Sum (Small)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-10-24 21:14:26 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 151 ms / 2,025 ms |
| コード長 | 1,085 bytes |
| コンパイル時間 | 408 ms |
| コンパイル使用メモリ | 82,224 KB |
| 実行使用メモリ | 76,920 KB |
| 最終ジャッジ日時 | 2025-11-02 21:15:38 |
| 合計ジャッジ時間 | 6,932 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 52 |
ソースコード
def solve_1(n: int):
"""
e = 1 の解を列挙する
試し割り法により Θ(n^{1/2}) 時間
"""
ans = []
n2 = n * 2
for w in range(1, n2):
if w * w > n2:
break
wl, r = divmod(n2, w)
if r or w >= wl or (w + wl) % 2 == 0:
continue
l = (wl - w + 1) // 2
ans.append((1, l, l + w - 1))
ans.reverse()
return ans
def solve_e(n: int, e: int):
"""
e = e の解を列挙する
尺取り法を用いて Θ(n^{1/e}) 時間
"""
# floor(n^{1/e}) + 1
m = next(i for i in range(n + 2) if i**e > n)
assert (m - 1)**e <= n < m**e
ans = []
# S(e, l, r-1) == s になるようにする
r = 1
s = 0
for l in range(1, m):
while s + r**e <= n:
s += r**e
r += 1
if s == n:
ans.append((e, l, r - 1))
s -= l**e
return ans
n = int(input())
assert 2 <= n <= 10**12
ans = solve_1(n)
for e in range(2, n.bit_length()):
ans.extend(solve_e(n, e))
print(len(ans))
for t in ans:
print(*t)