%I #11 Dec 14 2019 12:04:31
%S 1,4,18,112,810,7144,73458,850672,11069370,161190904,2575237698,
%T 44571447232,836188737930,16970931765064,368985732635538,
%U 8524290269083792,208874053200038490,5428866923032585624,149250273758730282978,4318265042184721248352
%N Expansion of e.g.f. Sum_{k>=1} (exp(x) - 1)^k / (k * (1 - (exp(x) - 1)^k)).
%H Vaclav Kotesovec, <a href="/A330353/b330353.txt">Table of n, a(n) for n = 1..420</a>
%F E.g.f.: -Sum_{k>=1} log(1 - (exp(x) - 1)^k).
%F E.g.f.: A(x) = log(B(x)), where B(x) = e.g.f. of A167137.
%F G.f.: Sum_{k>=1} (k - 1)! * sigma(k) * x^k / Product_{j=1..k} (1 - j*x), where sigma = A000203.
%F exp(Sum_{n>=1} a(n) * log(1 + x)^n / n!) = g.f. of the partition numbers (A000041).
%F a(n) = Sum_{k=1..n} Stirling2(n,k) * (k - 1)! * sigma(k).
%F a(n) ~ n! * Pi^2 / (12 * (log(2))^(n+1)). - _Vaclav Kotesovec_, Dec 14 2019
%t nmax = 20; CoefficientList[Series[Sum[(Exp[x] - 1)^k/(k (1 - (Exp[x] - 1)^k)), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]! // Rest
%t Table[Sum[StirlingS2[n, k] (k - 1)! DivisorSigma[1, k], {k, 1, n}], {n, 1, 20}]
%Y Cf. A000041, A000203, A000629, A002745, A008277, A038048, A167137, A308555, A330351, A330352, A330354.
%K nonn
%O 1,2
%A _Ilya Gutkovskiy_, Dec 11 2019
|