OFFSET
0,3
FORMULA
a(n) ~ sqrt(2*Pi) * 3^n * n^(n - 1/6) / (Gamma(1/3) * (1/r + sqrt(1 - r^2))^(1/3) * exp(n) * r^(n + 1/3)), where r = A069814. - Vaclav Kotesovec, Jun 24 2025
PROG
(PARI) a185951(n, k) = binomial(n, k)/2^k*sum(j=0, k, (2*j-k)^(n-k)*binomial(k, j));
a007559(n) = prod(k=0, n-1, 3*k+1);
a(n) = sum(k=0, n, a007559(k)*3^(n-k)*a185951(n, k));
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jun 24 2025
STATUS
approved
