結果

問題 No.3332 Consecutive Power Sum (Small)
コンテスト
ユーザー wasd314
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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)
0