OFFSET
1,3
COMMENTS
No more terms up to 10^6.
MATHEMATICA
n = 10^5;
t = Union@Flatten@Table[x^2 + y^2 + z^2 + w^5, {x, (n/3)^(1/2)}, {y, x, ((n-x^2)/2)^(1/2)}, {z, y, (n-x^2-y^2)^(1/2)}, {w, 1, (n-x^2-y^2-z^2)^(1/5)}];
Complement[Range[0, n], t]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
XU Pingya, Jun 17 2020
STATUS
approved