OFFSET
1,1
LINKS
MAPLE
N:= 10^6; # to get all terms <= N
S:= Vector(N):
for a from 1 to floor(sqrt(N/3)) do
for b from a+1 to floor(sqrt((N-a^2)/2)) do
c:= [$(b+1) .. floor(sqrt(N-a^2-b^2))]:
v:= map(t -> a^2 + b^2 + t^2, c):
S[v]:= map(`+`, S[v], 1)
od od: select(t -> S[t]=6, [$1..N]); # Robert Israel, Jan 03 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved