|
%I #6 Jul 30 2023 02:14:27
%S 0,1,3,6,13,19,28,37,48,64,77,94,110,131,152,172,199,226,253,281,308,
%T 343,377,412,447,488,528,567,612,654,703,750,796,847,902,957,1013,
%U 1068,1129,1187,1252,1313,1378,1446,1511,1582,1650,1725,1800,1877,1955,2034
%N Number of unit square lattice cells inside quadrant of origin centered circle of diameter 2n+1.
%C Number of unit square lattice cells inside quadrant of origin centered circle of radius n+1/2.
%H G. C. Greubel, <a href="/A136484/b136484.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) = Sum_{k=1..n} floor(sqrt((n+1/2)^2 - k^2)).
%F a(n) = (1/2) * A136515(n).
%F a(n) = (1/4) * A136486(n).
%F a(n) = A136483(2*n+1).
%F Lim_{n -> oo} a(n)/(n^2) -> Pi/16 (A019683).
%e a(2) = 3 because a circle of radius 2+1/2 in the first quadrant encloses (2,1), (1,1), (1,2).
%t Table[Sum[Floor[Sqrt[(n+1/2)^2 - k^2]], {k,n}], {n,0,100}]
%o (Magma)
%o A136484:= func< n | n eq 0 select 0 else (&+[Floor(Sqrt((n+1/2)^2-j^2)): j in [1..n]]) >;
%o [A136484(n): n in [0..100]]; // _G. C. Greubel_, Jul 29 2023
%o (SageMath)
%o def A136484(n): return sum(floor(sqrt((n+1/2)^2-k^2)) for k in range(1, n+1))
%o [A136484(n) for n in range(101)] # _G. C. Greubel_, Jul 29 2023
%Y Cf. A019683, A136483, A136486, A136515.
%K easy,nonn
%O 0,3
%A Glenn C. Foster (gfoster(AT)uiuc.edu), Jan 02 2008
|
