login
A Beatty sequence: a(n) = floor(n * (Pi-1)/(Pi-2)).
2

%I #28 Sep 04 2020 05:22:49

%S 1,3,5,7,9,11,13,15,16,18,20,22,24,26,28,30,31,33,35,37,39,41,43,45,

%T 46,48,50,52,54,56,58,60,61,63,65,67,69,71,73,75,76,78,80,82,84,86,88,

%U 90,91,93,95,97,99,101,103,105,106,108,110,112

%N A Beatty sequence: a(n) = floor(n * (Pi-1)/(Pi-2)).

%H Harry J. Smith, <a href="/A063460/b063460.txt">Table of n, a(n) for n = 1..1000</a>

%H J. Tamura, <a href="http://matwbn.icm.edu.pl/ksiazki/aa/aa71/aa7141.pdf">A class of transcendental numbers having explicit 0-adic and Jacobi-Perron expansions of arbitrary dimension</a>, Acta Arithmetica, LXXI.4 (1995), 301-329 (see Example p. 327).

%p Digits:=100:t:=evalf(1+sech(Pi/6)):A:=n->(t*n):seq(floor((t*n)),n=1..60); # _Zerinvary Lajos_, Mar 27 2009

%t Table[Floor[n (Pi - 1)/(Pi - 2)], {n, 50}] (* _Wesley Ivan Hurt_, Jul 04 2014 *)

%o (PARI) { default(realprecision, 50); f=(Pi - 1)/(Pi - 2); for (n=1, 1000, write("b063460.txt", n, " ", floor(n*f)) ) } \\ _Harry J. Smith_, Aug 22 2009

%o (Maxima) makelist(floor(n*(%pi - 1)/(%pi - 2)), n, 1, 65); /* _Martin Ettl_, Oct 17 2012 */

%K nonn

%O 1,2

%A _Robert G. Wilson v_, Jul 25 2001

%E Incorrect comment removed by _Georg Fischer_, Sep 04 2020