OFFSET
3,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 3..1000
FORMULA
G.f.: x^3*C(x)^7/(1 - x/Sqrt(1-4*x)) = x^3*(1-2*x*C(x))*C(x)^9/(1-x*C(x)^3), where C(x) is the g.f. of A000108. - G. C. Greubel, Jul 17 2019
a(n) ~ (2 + sqrt(5))^(n+3) * (3 - sqrt(5))^7 / (128*sqrt(5)). - Vaclav Kotesovec, Jul 18 2019
MATHEMATICA
Drop[CoefficientList[Series[Sqrt[1-4*x]*(1-Sqrt[1-4*x])^9/(64*x^4*(8*x^2 -(1-Sqrt[1-4*x])^3)), {x, 0, 40}], x], 3] (* G. C. Greubel, Jul 17 2019 *)
PROG
(PARI) my(x='x+O('x^40)); Vec(sqrt(1-4*x)*(1-sqrt(1-4*x))^9/(64*x^4*(8*x^2 -(1 - sqrt(1-4*x))^3 ))) \\ G. C. Greubel, Jul 17 2019
(Magma) R<x>:=PowerSeriesRing(Rationals(), 30); Coefficients(R!( Sqrt(1-4*x)*(1-Sqrt(1-4*x))^9/(64*x^4*(8*x^2 -(1-Sqrt(1-4*x))^3 )) )); // G. C. Greubel, Jul 17 2019
(Sage) a=(sqrt(1-4*x)*(1-sqrt(1-4*x))^9/(64*x^4*(8*x^2 -(1-sqrt(1-4*x))^3 ))).series(x, 45).coefficients(x, sparse=False); a[3:40] # G. C. Greubel, Jul 17 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved