editing
approved
editing
approved
a(n) ~ sqrt(2*Pi) * n^(n + 1/2) / ((1 + LambertW(1)) * exp(n) * (1/LambertW(1) - 1)^(n+1)). - Vaclav Kotesovec, Nov 11 2023
approved
editing
proposed
approved
editing
proposed
editing
proposed
(PARI) my(N=20, 30, x='x+O('x^N)); Vec(serlaplace(1/(1-(1+x)*log(1+x))))
1, 1, 3, 11, 60, 384, 3062, 27838, 293416, 3447768, 45277392, 651587760, 10254900048, 174557518992, 3203361670896, 62938642659504, 1319693558377728, 29390794198726656, 693223221342879360, 17256288944072200320, 452215395177034040064, 12442278072964931318016
1, 1, 3, 11, 60, 384, 3062, 27838, 293416, 3447768, 45277392, 651587760, 10254900048, 174557518992, 3203361670896, 62938642659504, 1319693558377728, 29390794198726656, 693223221342879360, 17256288944072200320, 452215395177034040064, 12442278072964931318016
a(n) = Sum_{k=0..n} |Stirling1(n,k)| * A006153(k).
a(n) = Sum_{k=0..n} |Stirling1(n,k)| * A006153(k).