login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(n) = Lucas(2^n - 2).
3

%I #9 Jan 03 2014 00:52:33

%S 2,3,18,843,1860498,9062201101803,215002084978043708894524818,

%T 121020968315000050139390193037122554865361969834971243,

%U 38343921554607207587938114587587818441864732465057252794474861753545122655196096751375348482086938743684498

%N a(n) = Lucas(2^n - 2).

%C Let phi := 1/2*(1 + sqrt(5)) denote the golden ratio. This sequence, apart from the initial term, gives an Engel expansion of 1 to the base phi^2 (see A230601 for a definition of this term). The associated Engel series expansion of 1 to the base phi^2 begins 1 = phi^2/3 + phi^4/(3*18) + phi^6/(3*18*843) + phi^8/(3*18*843*1860498) + .... This result can be extended in two ways. Firstly, the sequence Lucas(2*k*(2^n - 1)) for k = 1,2,3,... is an Engel expansion of 1 to the base phi^(2*k). Secondly, for n = 1,2,3,... the sequence [a(n),a(n+1),a(n+2),...] is an Engel expansion of phi^(4 - 2^n) to the base phi^2. Some example are given below.

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/Engel_expansion">Engel Expansion</a>

%F a(n) = A000032(2^n-2) = phi^(2^n-2) + (1/phi)^(2^n-2), where phi := 1/2*(1 + sqrt(5)).

%F Recurrence equation: a(1) = 2, a(2) = 3 and a(n) = floor(phi^2*a(n-1)^2) - 5 for n >= 3.

%e Engel series expansion of phi^(4 - 2^n) to the base phi^2 for n = 1 to 5.

%e n = 1:

%e phi^2 = phi^2/2 + phi^4/(2*3) + phi^6/(2*3*18) + phi^8/(2*3*18*843) + ...

%e n = 2:

%e 1 = phi^2/3 + phi^4/(3*18) + phi^6/(3*18*843) + phi^8/(3*18*843*1860498) + ...

%e n = 3:

%e 1/phi^4 = phi^2/18 + phi^4/(18*843) + phi^6/(18*843*1860498) + ...

%e n = 4:

%e 1/phi^12 = phi^2/843 + phi^4/(843*1860498) + phi^6/(843*1860498*9062201101803) + ...

%e n = 5:

%e 1/phi^28 = phi^2/1860498 + phi^4/(1860498*9062201101803) + ...

%t Table[LucasL[2^n - 2], {n, 1, 10}]

%Y Cf. A000032, A001622, A192223, A230600, A230601.

%K nonn,easy

%O 1,1

%A _Peter Bala_, Oct 28 2013