OFFSET
0,3
COMMENTS
Generally, for e.g.f. 1/(1-sinh(p*x))^(1/p) we have a(n) ~ n! * p^n / (Gamma(1/p) * 2^(1/(2*p)) * n^(1-1/p) * (arcsinh(1))^(n+1/p)).
FORMULA
a(n) ~ n! * 3^n / (Gamma(1/3) * 2^(1/6) * n^(2/3) * (log(1+sqrt(2)))^(n+1/3)).
MATHEMATICA
CoefficientList[Series[1/(1-Sinh[3*x])^(1/3), {x, 0, 20}], x] * Range[0, 20]!
PROG
(PARI) a136630(n, k) = 1/(2^k*k!)*sum(j=0, k, (-1)^(k-j)*(2*j-k)^n*binomial(k, j));
a007559(n) = prod(k=0, n-1, 3*k+1);
a(n) = sum(k=0, n, a007559(k)*3^(n-k)*a136630(n, k)); \\ Seiichi Manyama, Jun 24 2025
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vaclav Kotesovec, Jan 03 2014
STATUS
approved
