OFFSET
1,1
FORMULA
{n: A025442(n) >= 2}.
MAPLE
isA024804 := proc(n) local ws; ws := 0 ; for x from 1 do if x^2 > n then break; end if; for y from x+1 do if x^2+y^2>n then break; end if; z := n-x^2-y^2 ; if issqr(z) then z := sqrt(z) ; if z > y then ws := ws+1 ; end if; end if ; end do: end do: ws >= 2 ; end proc:
for n from 1 to 230 do if isA024804(n) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Dec 07 2010
MATHEMATICA
nn=230; t=Table[0, {nn}]; Do[n=i^2+j^2+k^2; If[n<=nn, t[[n]]++], {i, Sqrt[nn]}, {j, i+1, Sqrt[nn]}, {k, j+1, Sqrt[nn]}]; Flatten[Position[t, _?(#>1&)]]
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved