OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..385
FORMULA
E.g.f.: exp(T(x)/(1-T(x))) where T(x) is the e.g.f. for A000169.
a(n) ~ n^(n-1/3) * exp(3/2*n^(1/3) - 2/3) / sqrt(3). - Vaclav Kotesovec, Sep 24 2013
a(n) = Sum_{k=0..n} n^(n-k)*binomial(n-1,k-1)*A000262(k). - Fabian Pereyra, Jul 12 2024
The above formula can be written with the Abel polynomials: a(n) = Sum_{k=0..n} (-1)^(n - k) * A137452(n, k) * A000262(k). - Peter Luschny, Jul 13 2024
MAPLE
with(combinat):
b:= proc(n, i) option remember; `if`(n=0, 1,
`if`(i<1, 0, add(i^(i*j)*b(n-i*j, i-1)*
multinomial(n, n-i*j, i$j)/j!, j=0..n/i)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..20); # Alois P. Heinz, Mar 29 2016
MATHEMATICA
nn = 20; t = Sum[n^(n - 1) x^n/n!, {n, 1, nn}] ;
Range[0, nn]! CoefficientList[Series[Exp[t/(1 - t)], {x, 0, nn}], x]
CROSSREFS
KEYWORD
nonn
AUTHOR
Geoffrey Critzer, Dec 19 2011
STATUS
approved