|
%I M0914 #41 Mar 19 2024 11:58:32
%S 1,1,2,3,14,129,25298,420984147,269425140741515486,
%T 47749585090209528873482531562977121,
%U 3466137915373323052799848584927709551269254572949111609037058632767202
%N Denominators of convergents to Cahen's constant: a(n+2) = a(n)^2*a(n+1) + a(n).
%C Shifted square roots of partial quotients in continued fraction expansion of Cahen's constant: a(n) = sqrt(A006280(n+2)). - _Jonathan Sondow_, Aug 20 2014
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Amiram Eldar, <a href="/A006279/b006279.txt">Table of n, a(n) for n = 0..13</a>
%H J. L. Davison and Jeffrey O. Shallit, <a href="https://doi.org/10.1007/BF01332350">Continued Fractions for Some Alternating Series</a>, Monatshefte für Mathematik, Vol. 111 (1991), pp. 119-126; <a href="http://gdz.sub.uni-goettingen.de/dms/load/img/?PPN=PPN362162050_0111&DMDID=DMDLOG_0013&LOGID=LOG_0013&PHYSID=PHYS_0126">alternative link</a>.
%p A006279 := proc(n) option remember; if n <= 1 then 1 else A006279(n-2)^2*A006279(n-1)+A006279(n-2) fi end:
%p seq(A006279(n), n=0..10);
%t a[n_] := a[n] = If[n < 2, 1, a[n-2]^2*a[n-1] + a[n-2]];
%t Table[a[n], {n, 0, 9}] (* _Jean-François Alcover_, Sep 23 2022 *)
%o (Python)
%o from itertools import islice
%o def A006279_gen(): # generator of terms
%o a, b = 1, 1
%o yield a
%o while True:
%o yield b
%o a, b = b, a*(a*b+1)
%o A006279_list = list(islice(A006279_gen(),10)) # _Chai Wah Wu_, Mar 19 2024
%Y Cf. A118227, A006280, A006281.
%K nonn,easy,frac
%O 0,3
%A _N. J. A. Sloane_.
%E Definition clarified by _Jonathan Sondow_, Aug 20 2014
| 