OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..356
FORMULA
a(n) = Sum_{k=0..n} (2*k+1)^(k-1) * Stirling1(n,k).
E.g.f. A(x) = Sum_{k>=0} (2*k+1)^(k-1) * (log(1+x))^k / k!.
From Vaclav Kotesovec, Nov 20 2021: (Start)
E.g.f.: sqrt(-LambertW(-2*log(1 + x)) / (2*log(1 + x))).
a(n) ~ n^(n-1) / (sqrt(2) * (exp(exp(-1)/2) - 1)^(n - 1/2) * exp(n + exp(-1)/4 - 1)). (End)
MATHEMATICA
nmax = 20; CoefficientList[Series[Sqrt[-LambertW[-2*Log[1 + x]]/(2*Log[1 + x])], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Nov 20 2021 *)
PROG
(PARI) a(n) = sum(k=0, n, (2*k+1)^(k-1)*stirling(n, k, 1));
(PARI) N=20; x='x+O('x^N); Vec(serlaplace(sum(k=0, N, (2*k+1)^(k-1)*log(1+x)^k/k!)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Nov 20 2021
STATUS
approved