login
The sequence 2, floor(a), floor(a^2), floor(a^3), ..., with a = 1+sqrt(5).
1

%I #15 Sep 08 2022 08:45:02

%S 2,3,10,33,109,354,1148,3716,12026,38918,125943,407562,1318899,

%T 4268047,13811692,44695576,144637922,468058148,1514667986,4901568568,

%U 15861809082,51329892437,166107021206,537533612162,1739495309150,5629125066952,18216231370504

%N The sequence 2, floor(a), floor(a^2), floor(a^3), ..., with a = 1+sqrt(5).

%H G. Harman, <a href="http://at.yorku.ca/cgi-bin/amca/cadx-39">One hundred years of normal numbers</a>

%F a(n) = floor( (1+sqrt(5))^n ) for n>0, a(0)=2. - _Artur Jasinski_, Nov 22 2006

%t Join[{2}, Table[Floor[(Sqrt[5] + 1)^n], {n, 1, 30}]] (* _Artur Jasinski_, Nov 22 2006 *)

%o (PARI) a(n) = if (n==0, 2, floor((1+sqrt(5))^n)); \\ _Michel Marcus_, Feb 19 2016

%o (Magma) [2] cat [Floor((1+Sqrt(5))^n): n in [1..30]]; // _Vincenzo Librandi_, Feb 20 2016

%Y Cf. A014217.

%K nonn

%O 0,1

%A _N. J. A. Sloane_, Sep 12 2000, May 14 2007