OFFSET
0,2
FORMULA
G.f.: exp(Sum_{k>=1} ( Sum_{d|k} d * (6 - x^d)^(k/d) ) * x^k / (2*k)).
G.f.: A(x) = Product_{k>=1} B(x^k), where B(x) = g.f. of A001850 (central Delannoy numbers).
a(n) ~ sqrt(2) * (1 + sqrt(2))^(2*n - 1/2) / (c * sqrt(Pi*n)), where c = QPochhammer[1/(1 + sqrt(2))^2] = 0.799142925985081767883272500537236047... - Vaclav Kotesovec, Nov 21 2019
MATHEMATICA
nmax = 24; CoefficientList[Series[Product[1/(1 - 6 x^k + x^(2 k))^(1/2), {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 24; CoefficientList[Series[Exp[Sum[Sum[d (6 - x^d)^(k/d), {d, Divisors[k]}] x^k/(2 k), {k, 1, nmax}]], {x, 0, nmax}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Nov 21 2019
STATUS
approved