OFFSET
1,1
COMMENTS
These are the ordered numbers for which A223730 is not zero. The multiplicity for the number a(n) is A223730(a(n)).
According to the Halter-Koch reference the present sequence lists the ordered positive integers satisfying i) n not 0, 4, or 7 (mod 8) (see p.10, formula for r_3(n) attributed to A. Schinzel) and ii) n not from the set {1,2,5,10,13,25,37,58,85,130} with possibly one more positive integer member of this set which has to be >= 5*10^10 (if it exists at all). (Korollar 1. (b), p. 13). For this set see also A051952.
The first members with multiplicity 1 (precisely one representation) are 3, 6, 9, 11, 14, 17, 18, 19, 21, 22, 26, 27, 29, 30, 34, 35, 42, 43, 45, 46, 49, 50, 53, 61, 65, 67 ... A223732.
The first members with multiplicity 2 are 33, 38, 41, 51, 54, 57, 59, 62, 69, 74, 77, 81, 83, 90, 94, 98, 99, ... A223733.
The first members with multiplicity 3 are 66, 86, 89, 101, 110, 114, 131, 149, 153, 166, 171, 173, ... A223734.
For the complement see A223735.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
F. Halter-Koch, Darstellung natürlicher Zahlen als Summe von Quadraten, Acta Arith. 42 (1982) 11-20.
FORMULA
The sequence a(n) is obtained from the ordered set
{m positive integer | m = a^2 + b^2 + c^2 , a,b,c integer, 0 < a <= b <= c, gcd(a,b,c) = 1} with entries appearing only once.
Conjectured g.f.: (x^77 +2*x^76 -2*x^75 +x^74 -x^73 -x^72 +2*x^50 -x^49 +2*x^47 -2*x^46 -x^45 +x^34 +2*x^33 -2*x^32 +x^31 -x^30 -x^29 +2*x^22 -x^21 +2*x^19 -2*x^18 -x^17 +3*x^15 -2*x^14 +x^13 -x^12 -x^10 +2*x^9 +2*x^7 +2*x^6 -3*x^4 -2*x^3 -3*x^2 -3*x -3)*x / (-x^6 +x^5 +x -1). - Alois P. Heinz, Apr 06 2013
EXAMPLE
a(12) = 27 because 27 is the 12th number for which A223730 is nonzero. Because A223730(27) = 1 there is only one primitive sum of three nonzero squares which is 27 denoted by [1,1,5]:
1^2 + 1^2 + 5^2 = 27.
a(28) = 54 has two primitive representations in question, namely [1,2,7] and [2,5,5]. A223730(54) = 2. The representation [3,3,6] is not primitive because gcd(3,3,6) = 3 not 1.
a(34) = 66 has three representations in question, namely [1,1,8], [1,4,7] and [4,5,5].
MATHEMATICA
threeSquaresQ[n_] := Select[ PowersRepresentations[n, 3, 2], Times @@ #1 != 0 && GCD @@ #1 == 1 & ] != {}; Select[Range[120], threeSquaresQ] (* Jean-François Alcover, Jun 21 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Wolfdieter Lang, Apr 05 2013
STATUS
approved