OFFSET
0,2
COMMENTS
3rd binomial transform of A001850.
Diagonal of the rational function 1 / (1 - x^2 - 6*x*y - 2*y^2).
FORMULA
E.g.f.: exp(6*x) * BesselI(0,2*sqrt(2)*x).
Coefficient of x^n in (1 + 6*x + 2*x^2)^n.
n * a(n) = 6 * (2*n-1) * a(n-1) - 28 * (n-1) * a(n-2) with a(0) = 1, a(1) = 6.
a(n) = Sum_{k=0..floor(n/2)} binomial(n,2*k) * binomial(2*k,k) * 2^(n-k) * 3^(n-2*k).
MATHEMATICA
nmax = 22; CoefficientList[Series[1/Sqrt[1 - 12 x + 28 x^2], {x, 0, nmax}], x]
nmax = 22; CoefficientList[Series[Exp[6 x] BesselI[0, 2 Sqrt[2] x], {x, 0, nmax}], x] Range[0, nmax]!
a[0] = 1; a[1] = 6; a[n_] := a[n] = (6 (2 n - 1) a[n - 1] - 28 (n - 1) a[n - 2])/n; Table[a[n], {n, 0, 22}]
Table[Sum[Binomial[n, 2 k] Binomial[2 k, k] 2^(n - k) 3^(n - 2 k), {k, 0, Floor[n/2]}], {n, 0, 22}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Mar 13 2026
STATUS
approved
