%I #30 Nov 22 2022 10:46:47
%S 1,1,7,7,19,19,37,43,61,73,91,109,127,151,187,199,241,253,301,313,367,
%T 397,439,475,517,571,613,661,721,757,823,859,931,979,1045,1111,1165,
%U 1237,1303,1381,1459,1519,1615,1663,1765,1813,1921,1993,2083,2173,2263
%N Circle numbers (version 4): a(n)= number of points (i+j/2,j*sqrt(3)/2), i,j integers (triangular grid) contained in a circle of diameter n, centered at (0,0).
%C In other words, number of points in a hexagonal lattice covered by a circle of diameter n if the center of the circle is chosen at a grid point. - _Hugo Pfoertner_, Jan 07 2007
%C Same as above but "number of disks (r = 1)" instead of "number of points". See illustration in links. - _Kival Ngaokrajang_, Apr 06 2014
%H H. v. Eitzen, <a href="/A053416/b053416.txt">Table of n, a(n) for n = 0..1000</a>
%H Kival Ngaokrajang, <a href="/A053416/a053416_1.pdf">Illustration of initial terms</a>
%H <a href="/index/Aa#A2">Index entries for sequences related to A2 = hexagonal = triangular lattice</a>
%F a(n)/(n/2)^2->Pi*2/sqrt(3).
%F a(n) >= A053458(n). - _R. J. Mathar_, Nov 22 2022
%F a(2*n) = A308685(n). - _R. J. Mathar_, Nov 22 2022
%p A053416 := proc(d)
%p local a,j,imin,imax ;
%p a := 0 ;
%p for j from -floor(d/sqrt(3)) do
%p if j^2*3 > d^2 and j> 0 then
%p break ;
%p end if;
%p imin := ceil((-j-sqrt(d^2-3*j^2))/2) ;
%p imax := floor((-j+sqrt(d^2-3*j^2))/2) ;
%p a := a+imax-imin+1 ;
%p end do:
%p a ;
%p end proc:
%p seq(A053416(d),d=0..30) ; # _R. J. Mathar_, Nov 22 2022
%t a[n_] := Sum[Boole[4*(i^2 + i*j + j^2) <= n^2], {i, -n, n}, {j, -n, n}];
%t Table[a[n], {n, 0, 100}] (* _Jean-François Alcover_, Jun 06 2013, updated Apr 08 2022 to correct a discrepancy wrt b-file noticed by Georg Fischer *)
%Y Cf. A003215, A125849, A125850, A125851, A125852.
%Y Cf. A053411, A053414, A053415, A053417, A053458 (open disk), A308685 (bisection).
%K easy,nonn
%O 0,3
%A Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Jan 10 2000
%E Edited by _N. J. A. Sloane_, Jul 03 2008 at the suggestion of _R. J. Mathar_
|