結果

問題 No.3332 Consecutive Power Sum (Small)
コンテスト
ユーザー norioc
提出日時 2025-11-02 22:58:44
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 542 ms / 2,025 ms
コード長 1,590 bytes
コンパイル時間 317 ms
コンパイル使用メモリ 82,684 KB
実行使用メモリ 209,808 KB
最終ジャッジ日時 2025-11-02 22:58:58
合計ジャッジ時間 13,130 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 52
権限があれば一括ダウンロードができます

ソースコード

diff # 

from itertools import accumulate


def divisors(n: int) -> list[int]:
    """n の約数を求める"""
    s = set()
    p = 1
    while p * p <= n:
        if n % p == 0:
            s.add(p)
            s.add(n // p)
        p += 1
    return sorted(s)


def accum(a: list):
    acc = list(accumulate(a))
    return lambda l, r: acc[r] - (acc[l-1] if l > 0 else 0)


def solve(n: int, k: int):
    assert k > 1

    a = []
    i = 0
    while True:
        x = pow(i, k)
        if x > n: break
        a.append(x)
        i += 1

    acc = accum(a)
    ans = []
    for i in range(1, len(a)):
        if a[i] > n: break
        lo = i
        hi = len(a)-1
        while lo <= hi:
            m = (lo + hi) // 2
            x = acc(i, m)
            if x == n:
                ans.append((k, i, m))
                break
            elif x > n:
                hi = m - 1
            else:
                lo = m + 1

    return sorted(ans)


def solve_k0(n: int):
    divs = divisors(n)
    res = []
    for d in divs:
        if d == 1: continue

        m = d
        a = N//d - (d-1)//2
        s = m * (2*a + m - 1) // 2
        if s == n and a > 0:
            res.append((1, a, a+m-1))

        m = 2*N//d
        a = (d - 2*N//d + 1) // 2
        s = m * (2*a + m - 1) // 2
        if s == n and a > 0:
            res.append((1, a, a+m-1))

    res.append((1, n, n))
    return sorted(res)


N = int(input())

ans = solve_k0(N)
for i in range(2, 41):
    if pow(2, i) > N: break
    res = solve(N, i)
    ans.extend(res)

print(len(ans))
for k, l, r in ans:
    print(k, l, r)
0