OFFSET
0,3
LINKS
Eric Weisstein's World of Mathematics, Bell Polynomial.
FORMULA
a(n) = Sum_{k=0..n} 5^(n-k) * |Stirling1(n,k)| * A004211(k) = Sum_{k=0..n} 2^k * 5^(n-k) * |Stirling1(n,k)| * Bell_k(1/2), where Bell_n(x) is n-th Bell polynomial.
a(n) = (1/exp(1/2)) * (-5)^n * n! * Sum_{k>=0} binomial(-2*k/5,n)/(2^k * k!).
From Vaclav Kotesovec, Feb 04 2026: (Start)
a(n) = 35*(n-3)*a(n-1) - 105*(5*n^2 - 35*n + 63)*a(n-2) + 875*(n-4)*(5*n^2 - 40*n + 84)*a(n-3) - 7*(3125*n^4 - 56250*n^3 + 383125*n^2 - 1170000*n + 1351287)*a(n-4) + (65625*n^5 - 1640625*n^4 + 16471875*n^3 - 83015625*n^2 + 210003885*n - 213300674)*a(n-5) - 175*(n-6)*(n-5)*(625*n^4 - 13750*n^3 + 113750*n^2 - 419375*n + 581361)*a(n-6) + 125*(n-7)*(n-6)*(n-5)*(5*n - 33)*(5*n - 31)*(5*n - 29)*(5*n - 27)*a(n-7).
a(n) ~ 5^(n + 1/7) * n^(n - 5/14) / (sqrt(7) * exp(n + 1/2 - 7*5^(-5/7)*n^(2/7)/2)) * (1 + 1/(2*5^(10/7)*n^(3/7))). (End)
MATHEMATICA
CoefficientList[Series[Exp[ (1/(1-5*x)^(2/5) - 1)/2 ], {x, 0, 18}], x]Range[0, 18]! (* Stefano Spezia, Mar 31 2025 *)
PROG
(PARI) my(N=20, x='x+O('x^N)); Vec(serlaplace(exp((1/(1-5*x)^(2/5)-1)/2)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jan 18 2025
STATUS
approved
