OFFSET
0,4
COMMENTS
For a guide to related sequences, see A211422.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
MAPLE
b:= proc(n) option remember; 1+floor(sqrt(n)) end:
a:= proc(n) local c, x, y, w;
c:= 0;
for x to n do
for y from x to n do
w:= b(x^2+y^2);
if w>n then break fi;
c:= c+ (n-w+1)*`if`(x=y, 1, 2)
od
od: c
end:
seq (a(n), n=0..50); # Alois P. Heinz, May 10 2012
MATHEMATICA
t = Compile[{{n, _Integer}}, Module[{s = 0},
(Do[If[w^2 > x^2 + y^2, s = s + 1],
{w, 1, #}, {x, 1, #}, {y, 1, #}] &[n]; s)]];
Map[t[#] &, Range[0, 50]] (* A211637 *)
(* Peter J. C. Moses, Apr 13 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Apr 18 2012
STATUS
approved