%I #13 Oct 28 2023 09:23:54
%S 0,0,1,5,13,28,50,80,121,175,244,327,425,544,683,845,1028,1236,1470,
%T 1733,2027,2349,2706,3096,3520,3985,4489,5034,5619,6247,6922,7641,
%U 8411,9230,10102,11030,12007,13043,14133,15288,16504,17778,19117
%N Number of ordered triples (w,x,y) with all terms in {1,...,n} and w^2>=x^2+y^2.
%C For a guide to related sequences, see A211422.
%H Seiichi Manyama, <a href="/A211636/b211636.txt">Table of n, a(n) for n = 0..10000</a>
%t t = Compile[{{n, _Integer}}, Module[{s = 0},
%t (Do[If[w^2 >= x^2 + y^2, s = s + 1],
%t {w, 1, #}, {x, 1, #}, {y, 1, #}] &[n]; s)]];
%t Map[t[#] &, Range[0, 50]] (* A211636 *)
%t (* _Peter J. C. Moses_, Apr 13 2012 *)
%o (PARI) a(n) = sum(j=0, n, sum(k=1, j-1, sqrtint(j^2-k^2))); \\ _Seiichi Manyama_, Oct 28 2023
%Y Cf. A211422, A211651.
%K nonn
%O 0,4
%A _Clark Kimberling_, Apr 18 2012
|