login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A169580
Squares of the form x^2+y^2+z^2 with x,y,z positive integers.
10
9, 36, 49, 81, 121, 144, 169, 196, 225, 289, 324, 361, 441, 484, 529, 576, 625, 676, 729, 784, 841, 900, 961, 1089, 1156, 1225, 1296, 1369, 1444, 1521, 1681, 1764, 1849, 1936, 2025, 2116, 2209, 2304, 2401, 2500, 2601, 2704, 2809, 2916, 3025, 3136, 3249, 3364
OFFSET
1,1
COMMENTS
Integer solutions of a^2 = b^2 + c^2 + d^2, i.e., Pythagorean Quadruples. - Jon Perry, Oct 06 2012
Also null (or light-like, or isotropic) vectors in Minkowski 4-space. - Jon Perry, Oct 06 2012
REFERENCES
T. Nagell, Introduction to Number Theory, Wiley, 1951, p. 194.
EXAMPLE
9 = 1 + 4 + 4,
36 = 16 + 16 + 4,
49 = 36 + 9 + 4,
81 = 49 + 16 + 16,
so these are in the sequence.
16 cannot be written as the sum of 3 squares if zero is not allowed, therefore 16 is not in the sequence.
Also we can see that 49-36-9-4=0, so (7,6,3,2) is a null vector in the signatures (+,-,-,-) and (-,+,+,+). - Jon Perry, Oct 06 2012
MAPLE
M:= 10000: # to get all terms <= M
sort(convert(select(issqr, {seq(seq(seq(x^2 + y^2 + z^2,
z=y..floor(sqrt(M-x^2-y^2))), y=x..floor(sqrt((M-x^2)/2))),
x=1..floor(sqrt(M/3)))}), list)); # Robert Israel, Jan 28 2016
MATHEMATICA
Select[Range[60]^2, Resolve@ Exists[{x, y, z}, Reduce[# == x^2 + y^2 + z^2, {x, y, z}, Integers], And[x > 0, y > 0, z > 0]] &] (* Michael De Vlieger, Jan 27 2016 *)
CROSSREFS
For the square roots see A005767. Cf. A000378, A000419.
Cf. A217554.
Sequence in context: A192610 A319958 A329808 * A258844 A294952 A068810
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Mar 02 2010
STATUS
approved