OFFSET
1,2
COMMENTS
A(x) is E.g.f. for a(n)/n. - Vaclav Kotesovec, Jan 24 2014
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..360
Vladimir Kruchinin, D. V. Kruchinin, Composita and their properties, arXiv:1103.2582 [math.CO], 2013.
FORMULA
a(n) = (1/n)*Sum_{m=1..(n-1)} m!*binomial(n,m)*Sum_{i=2*m..(n+m-1)} binomial(n-1,i-m)*Stirling1(i-m,m)*(-1)^(i-2*m)*(n-m)^(n+m-i-1))) + n^(n-2), n > 0.
a(n) ~ (1-s)*sqrt(s/(1-(2-r)*s+s^2)) * n^(n-2) / (exp(n) * r^n), where s = 0.6276084459806374335... is the root of the equation (exp(s) + 1/(1-s))*s = exp(s) + log(1/(1-s)), and r = 1/(exp(s) + 1/(1-s)) = 0.2193717986293565123... - Vaclav Kotesovec, Jan 24 2014
MATHEMATICA
Table[(Rest[CoefficientList[InverseSeries[Series[x/(E^x + Log[(1 - x)^(-1)]), {x, 0, 20}], x], x] * Range[0, 20]!])[[n]]/n, {n, 1, 20}] (* Vaclav Kotesovec, Jan 24 2014 *)
Table[1/n*Sum[m!*Binomial[n, m]*Sum[Binomial[n-1, i-m]*StirlingS1[i-m, m]*(-1)^(i-2*m)*(n-m)^(n+m-i-1), {i, 2*m, n+m-1}], {m, 1, n-1}]+n^(n-2), {n, 1, 20}] (* Vaclav Kotesovec after Vladimir Kruchinin, Jan 24 2014 *)
PROG
(Maxima)
a(n):=1/n*sum(m!*binomial(n, m)*sum(binomial(n-1, i-m)*stirling1(i-m, m)*(-1)^(i-2*m)*(n-m)^(n+m-i-1), i, 2*m, n+m-1), m, 1, n-1)+n^(n-2);
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, Feb 21 2011
EXTENSIONS
Offset corrected by Vaclav Kotesovec, Jan 24 2014
STATUS
approved