A059551 Beatty sequence for Gamma(1/3).
2, 5, 8, 10, 13, 16, 18, 21, 24, 26, 29, 32, 34, 37, 40, 42, 45, 48, 50, 53, 56, 58, 61, 64, 66, 69, 72, 75, 77, 80, 83, 85, 88, 91, 93, 96, 99, 101, 104, 107, 109, 112, 115, 117, 120, 123, 125, 128, 131, 133, 136, 139, 141, 144, 147, 150, 152, 155, 158, 160, 163
Offset: 1
Links
- Harry J. Smith, Table of n, a(n) for n = 1..2000
- Aviezri S. Fraenkel, Jonathan Levitt, Michael Shimshoni, Characterization of the set of values f(n)=[n alpha], n=1,2,..., Discrete Math. 2 (1972), no.4, 335-345.
- Eric Weisstein's World of Mathematics, Beatty Sequence
- Index entries for sequences related to Beatty sequences
Programs
-
Magma
[Floor(n*Gamma(1/3)): n in [1..80]]; // Vincenzo Librandi, Jan 07 2015
-
Mathematica
Table[Floor[n Gamma[1/3]], {n, 70}] (* Vincenzo Librandi, Jan 07 2015 *)
-
PARI
{ default(realprecision, 100); b=gamma(1/3); for (n = 1, 2000, write("b059551.txt", n, " ", floor(n*b)); ) } \\ Harry J. Smith, Jun 28 2009
Formula
a(n) = floor(n*Gamma(1/3)). - Michel Marcus, Jan 04 2015