OFFSET
1,1
COMMENTS
{a(n)} UNION A376291 = positive numbers k such that k^k is not a cube and can be expressed as (x^3 + y^3)/2 for nonnegative integers x, y.
All terms are even.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..897
PROG
(Python)
from itertools import count, islice
from sympy import factorint
def A376315_gen(startvalue=2): # generator of terms >= startvalue
for k in count(max(startvalue+(startvalue&1), 2), 2):
f = {p:k*e for p, e in factorint(k).items()}
f[2] += 1
if not any(v%3 for v in f.values()):
yield k
CROSSREFS
KEYWORD
nonn
AUTHOR
Chai Wah Wu, Sep 20 2024
STATUS
approved
