OFFSET
1,2
LINKS
MAPLE
N:= 1000: # to get all terms <= N
V:= Vector(N):
for x from 1 to floor(sqrt(N/3)) do
for y from x+1 to floor(sqrt((N-x^2)/2)) do
zs:= [$(y+1).. floor(sqrt(N-x^2-y^2))];
V[map(z -> x^2 + y^2 + z^2, zs)]:= 1;
od
od:
select(i -> V[i] = 0, [$1..N]); # Robert Israel, Dec 31 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved