login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A222468 Numerator sequence of the n-th convergent of the continued fraction 1/(1+2/(2+2/(3+2/(4+... 4
1, 2, 8, 36, 196, 1248, 9128, 75520, 697936, 7130400, 79830272, 972224064, 12798573376, 181124475392, 2742464277632, 44241677392896, 757593444234496, 13725165351006720, 262293328557596672, 5273316901853946880 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
The corresponding denominator sequence is A222467(n).
a(n) = Phat(n,2) with the numerator polynomials Phat of A221913. All the given formulas follow from there and the comments given under A084950. The limit of the continued fraction (0 + K_{k>=1} (2/k))/2 = 1/(1+2/(2+2/(3+2/(4+... is (1/2)*sqrt(2)*BesselI(1,2*sqrt(2))/BesselI(0,2*sqrt(2)) = 0.5631786198117... See A222466 for more decimals.
For a combinatorial interpretation in terms of labeled Morse codes see a comment on A221913. Here each dash has label x=2, and the dots have label j if they are at position j. Labels are multiplied and all codes on [2,...,n+1] are summed.
LINKS
FORMULA
Recurrence: a(n) = n*a(n-1) + 2*a(n-2), a(-1) = 1/2, a(0) = 0, n >= 1.
As a sum: a(n) = Sum_{m=0..floor(n/2)} a(n-m,m)*2^m, n >= 1, with a(n,m) = binomial(n-1-m,m)*(n-m)!/(m+1)! = |A066667(n,m)| (Laguerre coefficients, parameter alpha = 1).
Explicit form: a(n) = 2*(w/2)^n*(BesselI(1,w)*BesselK(n+1,w) + BesselK(1,w)*BesselI(n+1,w)*(-1)^(n+1)), n >= 1, where w := -2*sqrt(2).
E.g.f.: Pi*(BesselJ(1, 2*i*sqrt(2)*sqrt(1-z))* BesselY(1, 2*i*sqrt(2)) - BesselY(1, (2*i)*sqrt(2)*sqrt(1-z))*BesselJ(1, 2*i*sqrt(2)))/sqrt(1-z) with Bessel functions and the imaginary unit i = sqrt(-1). Phat(0,x) = 0.
Asymptotics: lim_{n -> infinity} a(n)/n! = BesselI(1,2*sqrt(2)) /sqrt(2) = 2.3948330992734...
EXAMPLE
a(4) = 4*a(3) + 2*a(2) = 4*8 + 2*2 = 36.
Continued fraction convergent: 1/(1+2/(2+2/(3+2/4))) = 9/16 = 36/64 = a(4)/A222467(4).
Morse code a(5) = 196 from the sum of all 5 labeled codes on [2,3,4,5], one with no dash, three with one dash and one with two dashes: 5!/1 + (4*5 + 2*5 + 2*3)*2 +2^2 = 196.
MATHEMATICA
RecurrenceTable[{a[1] == 1, a[2] == 2, a[n] == n*a[n - 1] + 2 a[n - 2]}, a[n], {n, 20}] (* Ray Chandler, Jul 30 2015 *)
PROG
(PARI) m=30; v=concat([1, 2], vector(m-2)); for(n=3, m, v[n]=n*v[n-1] +2*v[n-2]); v \\ G. C. Greubel, May 17 2018
(Magma) I:=[1, 2]; [n le 2 select I[n] else n*Self(n-1) + 2*Self(n-2): n in [1..30]]; // G. C. Greubel, May 17 2018
CROSSREFS
Cf. A084950, A221913, A222467, A001040(n+1) (x=1), A058797 (x=-1).
Sequence in context: A129044 A052582 A367920 * A020021 A213076 A233744
KEYWORD
nonn,easy
AUTHOR
Gary Detlefs and Wolfdieter Lang, Mar 21 2013
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 17:25 EDT 2024. Contains 371254 sequences. (Running on oeis4.)