A211638 Number of ordered triples (w, x, y) with all terms in {1, ..., n} and w^2 + x^2 + y^2 < n.

0, 0, 0, 0, 1, 1, 1, 4, 4, 4, 7, 7, 10, 11, 11, 17, 17, 17, 20, 23, 26, 26, 32, 35, 35, 38, 38, 44, 48, 48, 54, 60, 60, 60, 66, 69, 75, 78, 78, 87, 87, 87, 96, 102, 105, 108, 114, 120, 120, 121, 127, 133, 139, 139, 145, 157, 157, 163, 169, 169, 178, 178, 184

Offset

0,8

Comments

For a guide to related sequences, see A211422.

Links

David A. Corneth, Table of n, a(n) for n = 0..10000

Formula

a(n) + A063691(n) = A211639(n). - R. J. Mathar, Jun 16 2023

a(n) = A211639(n-1). - R. J. Mathar, Jun 16 2023

Mathematica

t = Compile[{{n, _Integer}}, Module[{s = 0},
    (Do[If[w^2 + x^2 + y^2 < n, s = s + 1],
        {w, 1, #}, {x, 1, #}, {y, 1, #}] &[n]; s)]];
Map[t[#] &, Range[0, 80]]     (* A211638 *)
(* Peter J. C. Moses, Apr 13 2012 *)

Prog

(PARI) first(n) = {n = max(n, 2); n-=2; my(res = vector(n), v = vector(n)); forvec(x = vector(3, i, [1, sqrtint(n)]), c = sum(i = 1, 3, x[i]^2); if(c <= n, v[c]++)); for(i = 2, #v, v[i]+=v[i-1]); concat([0, 0], v)} \\ David A. Corneth, Jun 16 2023

Crossrefs

Cf. A211422.

Sequence in context: A023958 A137751 A244481 * A211639 A010658 A128939

Adjacent sequences: A211635 A211636 A211637 * A211639 A211640 A211641

Keyword

nonn

Author

Clark Kimberling, Apr 18 2012

Status

approved