%I #14 Jul 06 2020 11:25:57
%S 4,8,12,16,20,25,29,33,37,41,46,50,54,58,62,67,71,75,79,83,87,92,96,
%T 100,104,108,113,117,121,125,129,134,138,142,146,150,154,159,163,167,
%U 171,175,180,184,188,192,196,201,205,209,213,217,222,226,230,234,238,242
%N Beatty sequence for 4*Pi/3 = 4.1887902... .
%C a(n) = A109238(n) for n <= 20;
%C complement of A164087;
%C a(n) = A164088(A164087(n)) and A164088(a(n)) = A164087(n);
%C a(A000578(n)) = A066645(n).
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/BeattySequence.html">Beatty Sequence</a>
%H <a href="/index/Be#Beatty">Index entries for sequences related to Beatty sequences</a>
%F a(n) = floor(4*n*Pi/3).
%e a(3^3) = a(27) = 113 = (integer part of volume of sphere with radius=3) = A066645(3).
%t With[{c=4 \[Pi]/3},Floor[c #]&/@Range[70]] (* _Harvey P. Dale_, Mar 19 2011 *)
%o (Maxima)
%o fprec:100$
%o A164086:4*%pi/3$
%o ev(bfloat(A164086)); /* _Martin Ettl_, Nov 03 2012 */
%Y Cf. A022844, A019699.
%K nonn
%O 1,1
%A _Reinhard Zumkeller_, Aug 11 2009