OFFSET
3,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 3..1000
FORMULA
a(n) = 2*A025179(n-1).
From G. C. Greubel, Mar 01 2017: (Start)
a(n) = 2*Sum_{k=0..floor(n/2)} binomial(n, 2*k)*binomial(2*k+1, k+1), for n>=1.
O.g.f.: ((1-x)^2-(1-x+2*x^2)*sqrt(1-2*x-3*x^2)) / sqrt(1-2*x-3*x^2) [corrected by Charles R Greathouse IV, Mar 05 2017]
E.g.f.: 2*exp(x)*(BesselI(0, 2*x) + BesselI(2, 2*x)). (End)
MATHEMATICA
Rest[Differences[CoefficientList[Series[1/Sqrt[(1 + x) (1 - 3 x)], {x, 0, 30}], x], 2]] (* Harvey P. Dale, May 11 2013 *)
Table[2 Sum[Binomial[n, 2 k] Binomial[2 k + 1, k + 1], {k, 0, Floor[n/2]}], {n, 1, 25}] (* G. C. Greubel, Mar 01 2017 *)
Rest[Rest[CoefficientList[Series[((1 - x)^2 - (1 - x) Sqrt[1 - 2 x - 3 x^2])/(x Sqrt[1 - 2 x - 3 x^2]), {x, 0, 15}], x]]] (* G. C. Greubel, Mar 02 2017 *)
PROG
(PARI) x='x +O('x^50); Vec(((1-x)^2-(1-x +2*x^2)*sqrt(1-2*x-3*x^2)) /(x*sqrt(1 - 2*x -3*x^2))) \\ G. C. Greubel, Mar 01 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved