A138020 G.f. satisfies A(x) = sqrt( (1 + 2x*A(x)) / (1 - 2x*A(x)) ).

1, 2, 6, 24, 110, 544, 2828, 15232, 84246, 475648, 2730068, 15882240, 93438540, 554967040, 3323125528, 20039827456, 121597985254, 741871845376, 4548193111428, 28004975116288, 173113004348580, 1073893324357632

Offset

0,2

Links

Vaclav Kotesovec, Table of n, a(n) for n = 0..1000

Formula

a(n) ~ 2^(n - 1/2) * phi^((5*n + 3)/2) / (sqrt(Pi) * 5^(1/4) * n^(3/2)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Oct 04 2020

From Alexander Burstein, Nov 26 2021: (Start)

G.f.: A(x) = 1 + 2*x*A(x)*(1 + A(x)^2)/(1 + A(x)).

G.f.: A(-x*A(x)^2) = 1/A(x). (End)

D-finite with recurrence +n*(n+1)*(5*n-11) *a(n) +4*(-55*n^3 +231*n^2 -263*n +51)*a(n-2) -16*(n-3)*(n-4)*(5*n-1)*a(n-4)=0. - R. J. Mathar, Mar 25 2024

Mathematica

CoefficientList[y/.AsymptoticSolve[y^2-1-2x(y+y^3) ==0, y->1, {x, 0, 21}][[1]], x] - Alexander Burstein, Nov 26 2021

Prog

(PARI) a(n)=polcoeff((1/x)*serreverse(x*sqrt((1-2*x)/(1+2*x+x^2*O(x^n)))), n)

Crossrefs

Sequence in context: A141254 A366706 A216879 * A046646 A342284 A174195

Adjacent sequences: A138017 A138018 A138019 * A138021 A138022 A138023

Keyword

nonn

Author

Paul D. Hanna, Feb 28 2008

Status

approved