STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
(PARI) a(n) = sum(k=1, n, (-1)^(n-k)*stirling(n, k, 1)*(k-1)!*numdiv(k)); \\ Michel Marcus, Dec 16 2019
approved
editing
editing
approved
Vaclav Kotesovec, <a href="/A330493/b330493.txt">Table of n, a(n) for n = 1..400</a>
approved
editing
editing
approved
E.g.f.: -Sum_{k>=1} log(1 - log(1/(1 - x))^k) / k.
approved
editing
editing
approved
nmax = 20; Rest[CoefficientList[Series[-Sum[Log[1 - Log[1/(1 - x)]^k]/k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]!]
approved
editing