OFFSET
1,1
COMMENTS
First differs from A003338 at term 64: A003338(64) = 1393 is also a term of A003337, so not a term here. - Michael S. Branicky, Apr 19 2021
FORMULA
Equals A003338 - A344188 - A344187 - A000583, where "-" denotes "set difference". - Sean A. Irvine, May 15 2021
PROG
(Python)
limit = 1153
from functools import lru_cache
qd = [k**4 for k in range(1, int(limit**.25)+2) if k**4 + 3 <= limit]
qds = set(qd)
@lru_cache(maxsize=None)
def findsums(n, m):
if m == 1: return {(n, )} if n in qds else set()
return set(tuple(sorted(t+(q, ))) for q in qds for t in findsums(n-q, m-1))
A003338s = set(n for n in range(4, limit+1) if len(findsums(n, 4)) >= 1)
A003337s = set(n for n in range(3, limit+1) if len(findsums(n, 3)) >= 1)
A003336s = set(n for n in range(2, limit+1) if len(findsums(n, 2)) >= 1)
print(sorted(A003338s - A003337s - A003336s - qds)) # Michael S. Branicky, Apr 19 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Arlin Anderson (starship1(AT)gmail.com)
STATUS
approved