OFFSET
1,3
LINKS
Zhuorui He, Table of n, a(n) for n = 1..10000
PROG
(Python)
def A004886_upto(n):
a=set()
for i in range(n):
if 2*(i**9)>n: break
for j in range(i, n):
if i**9+j**9<=n: a.add(i**9+j**9)
else: break
return sorted(a) # Zhuorui He, Jun 30 2025
(PARI) lista(nn) = setbinop((x, y)->x^9+y^9, [0..nn]); \\ Michel Marcus, Jul 02 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
