OFFSET
0,3
LINKS
Robert Israel, Table of n, a(n) for n = 0..445
FORMULA
E.g.f.: Product_{k>0} exp(x^(3*k-1)).
a(n) ~ exp(2*sqrt(3*n)/3 - n - 1/6) * n^(n-1/4) / (sqrt(2) * 3^(1/4)). - Vaclav Kotesovec, Oct 10 2017
D-finite with recurrence: n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*a(n) - (n + 5)*(n + 4)*(n + 3)*(n + 2)*a(n + 1) - 2*(n + 5)*(n + 4)*(n + 3)*a(n + 3) - 2*a(n + 4)*(n + 5) + a(n + 6) = 0. - Robert Israel, Feb 22 2026
MAPLE
f:= gfun:-rectoproc({n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*a(n) - (n + 5)*(n + 4)*(n + 3)*(n + 2)*a(n + 1) - 2*(n + 5)*(n + 4)*(n + 3)*a(n + 3) - 2*a(n + 4)*(n + 5) + a(n + 6), a(0)=1, a(1)=0, a(2)=2, a(3)=0, a(4)=12, a(5)=120}, a(n), remember):
map(f, [$0..30]); # Robert Israel, Feb 22 2026
PROG
(PARI) N=66; x='x+O('x^N); Vec(serlaplace(exp(x^2/(1-x^3))))
(PARI) N=66; x='x+O('x^N); Vec(serlaplace(prod(k=1, N, exp(x^(3*k-1)))))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Oct 10 2017
STATUS
approved
