OFFSET
0,3
MATHEMATICA
r[n_] := Reduce[n == x^2 + y^2 + 7z^2, {x, y, z}, Integers]; Reap[For[n = 0, n <= 100, n++, If[r[n] =!= False, Print[n]; Sow[n]]]][[2, 1]] (* Jean-François Alcover, Feb 08 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved