OFFSET
1,2
LINKS
MATHEMATICA
nn = 100; lim = nn^2; t = Table[0, {lim}]; Do[k = x^2 + y^2 + z^2; If[0 < k <= lim, t[[k]]++], {x, 0, nn}, {y, x, nn}, {z, y, nn}]; u = Union[t]; c = Complement[Range[u[[-1]]], u]; If[c == {}, mx = u[[-1]], mx = c[[1]] - 1]; Join[{0}, Flatten[Table[Position[t, _?(# >= n &), 1, 1], {n, 2, mx}]]] (* T. D. Noe, Jun 20 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved