%I #24 Jul 14 2021 17:07:48
%S 1,4,1,32,1,196,1,1152,1,6724,1,39200,1,228484,1,1331712,1,7761796,1,
%T 45239072,1,263672644,1,1536796800,1,8957108164,1,52205852192,1,
%U 304278004996,1,1773462177792,1,10336495061764,1,60245508192800,1,351136554095044,1
%N Continued fraction expansion for exp( Sum_{n>=1} 1/(n*A003499(n)) ), where A003499(n) = (3+sqrt(8))^n + (3-sqrt(8))^n.
%H Colin Barker, <a href="/A174501/b174501.txt">Table of n, a(n) for n = 1..1000</a>
%H Peter Bala, <a href="/A174500/a174500_2.pdf">Some simple continued fraction expansions for an infinite product, Part 1</a>
%H <a href="/index/Rec#order_06">Index entries for linear recurrences with constant coefficients</a>, signature (0,7,0,-7,0,1).
%F a(2n-1) = 1, a(2n) = A003499(n) - 2, for n>=1 [conjecture].
%F The above conjectures are correct. See the Bala link for details. - _Peter Bala_, Jan 08 2013
%F a(n) = 7*a(n-2)-7*a(n-4)+a(n-6). G.f.: -x*(x^4+4*x^3-6*x^2+4*x+1) / ((x-1)*(x+1)*(x^2-2*x-1)*(x^2+2*x-1)). - _Colin Barker_, Jan 20 2013
%F a(n) = (((-1-sqrt(2))^n+(1-sqrt(2))^n+(sqrt(2)-1)^n+(1+sqrt(2))^n-4))/2 for n even. - _Colin Barker_, May 11 2016
%e Let L = Sum_{n>=1} 1/(n*A003499(n)) or, more explicitly,
%e L = 1/6 + 1/(2*34) + 1/(3*198) + 1/(4*1154) + 1/(5*6726) +...
%e so that L = 0.1833074113563494600094468694966574405706183998044...
%e then exp(L) = 1.2011836088120841844713993433258934531421726294252...
%e equals the continued fraction given by this sequence:
%e exp(L) = [1;4,1,32,1,196,1,1152,1,6724,1,39200,1,...]; i.e.,
%e exp(L) = 1 + 1/(4 + 1/(1 + 1/(32 + 1/(1 + 1/(196 + 1/(1 +...)))))).
%e Compare these partial quotients to A003499(n), n=1,2,3,...:
%e [6,34,198,1154,6726,39202,228486,1331714,7761798,45239074,...].
%t LinearRecurrence[{0,7,0,-7,0,1},{1,4,1,32,1,196},50] (* _Harvey P. Dale_, Jul 14 2021 *)
%o (PARI) {a(n)=local(L=sum(m=1,2*n+1000,1./(m*round((3+sqrt(8))^m+(3-sqrt(8))^m))));contfrac(exp(L))[n]}
%o (PARI) Vec(-x*(x^4+4*x^3-6*x^2+4*x+1)/((x-1)*(x+1)*(x^2-2*x-1)*(x^2+2*x-1)) + O(x^50)) \\ _Colin Barker_, May 11 2016
%Y Cf. A003499, A174500, A174502.
%K cofr,nonn,easy
%O 1,2
%A _Paul D. Hanna_, Mar 20 2010
|