結果
| 問題 |
No.3332 Consecutive Power Sum (Small)
|
| コンテスト | |
| ユーザー |
norioc
|
| 提出日時 | 2025-11-02 22:56:59 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,564 bytes |
| コンパイル時間 | 403 ms |
| コンパイル使用メモリ | 82,596 KB |
| 実行使用メモリ | 210,140 KB |
| 最終ジャッジ日時 | 2025-11-02 22:57:14 |
| 合計ジャッジ時間 | 13,258 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 2 |
| other | AC * 12 WA * 40 |
ソースコード
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 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:
res.append((1, a, a+m-1))
res.append((1, n, n))
return 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)
norioc