%I #23 Jul 05 2024 10:29:47
%S 3,7,11,15,19,22,26,30,34,38,42,45,49,53,57,61,65,68,72,76,80,84,87,
%T 91,95,99,103,107,110,114,118,122,126,130,133,137,141,145,149,152,156,
%U 160,164,168,172,175,179,183,187,191,195,198,202,206,210,214,217,221
%N Beatty sequence for Gamma(2/3)/(Gamma(2/3)-1).
%H Harry J. Smith, <a href="/A059554/b059554.txt">Table of n, a(n) for n = 1..2000</a>
%H Aviezri S. Fraenkel, Jonathan Levitt and Michael Shimshoni, <a href="http://dx.doi.org/10.1016/0012-365X(72)90012-X">Characterization of the set of values f(n)=[n alpha], n=1,2,...</a>, Discrete Math. 2 (1972), no.4, 335-345.
%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(n*gamma(2/3)/(gamma(2/3)-1)). - _Michel Marcus_, Jan 05 2015
%t Floor[Range[100]*(1 + 1/(Gamma[2/3] - 1))] (* _Paolo Xausa_, Jul 05 2024 *)
%o (PARI) { default(realprecision, 100); b=gamma(2/3)/(gamma(2/3) - 1); for (n = 1, 2000, write("b059554.txt", n, " ", floor(n*b)); ) } \\ _Harry J. Smith_, Jun 28 2009
%Y Beatty complement is A059553.
%Y Cf. A073006.
%K nonn,easy
%O 1,1
%A _Mitch Harris_, Jan 22 2001