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”).

A119810
Partial quotients of the continued fraction of the constant defined by binary sums involving Beatty sequences: c = Sum_{n>=1} 1/2^A049472(n) = Sum_{n>=1} A001951(n)/2^n.
1
2, 3, 10, 132, 131104, 2199023259648, 633825300114114700748888473600, 883423532389192164791648750371459257913741948437810659652423818057613312
OFFSET
1,1
COMMENTS
Convergents A119811: [2/1,7/3,72/31,9511/4095,1246930216/536870911,...], where the denominators of the convergents are equal to [2^A000129(n-1)-1] and A000129 is the Pell numbers. The number of digits in these partial quotients are (beginning at n=1): [1,1,2,3,6,13,30,72,174,420,1013,2445,5901,14246,34391,83027,...].
LINKS
W. W. Adams and J. L. Davison, A remarkable class of continued fractions, Proc. Amer. Math. Soc. 65 (1977), 194-198.
P. G. Anderson, T. C. Brown, P. J.-S. Shiue, A simple proof of a remarkable continued fraction identity, Proc. Amer. Math. Soc. 123 (1995), 2005-2009.
J. L. Davison, A series and its associated continued fraction, Proc. Amer. Math. Soc., 63 (1977), 29-32.
FORMULA
a(n) = 2^A001333(n-1) + 2^A000129(n-2) for n>1, with a(1)=2.
EXAMPLE
c = 2.32258852258806773012144068278798408011950250800432925665718...
The partial quotients start:
a(1) = 2^1; a(2) = 2^1 + 2^0; a(3) = 2^3 + 2^1;
a(4) = 2^7 + 2^2; a(5) = 2^17 + 2^5; a(6) = 2^41 + 2^12;
and continue as a(n) = 2^A001333(n-1) + 2^A000129(n-2) where
A001333(n) = ( (1+sqrt(2))^n + (1-sqrt(2))^n )/2;
A000129(n) = ( (1+sqrt(2))^n - (1-sqrt(2))^n )/(2*sqrt(2)).
MATHEMATICA
(* b = A001333 *) b[0] = 1; b[1] = 1; b[n_] := b[n] = 2 b[n-1] + b[n-2]; a[1] = 2; a[n_] := 2^b[n-1] + 2^Fibonacci[n-2, 2]; Array[a, 10] (* Jean-François Alcover, May 04 2017 *)
PROG
(PARI) {a(n)=if(n==1, 2, 2^round(((1+sqrt(2))^(n-1)+(1-sqrt(2))^(n-1))/2) +2^round(((1+sqrt(2))^(n-2)-(1-sqrt(2))^(n-2))/(2*sqrt(2))))}
CROSSREFS
Cf. A119809 (decimal expansion), A119811 (convergents); A119812 (dual constant).
Sequence in context: A155148 A076927 A064647 * A055708 A162647 A368063
KEYWORD
cofr,nonn
AUTHOR
Paul D. Hanna, May 26 2006
STATUS
approved