OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..345
FORMULA
a(n) ~ 5^(-1/4) * ((1+sqrt(5))/2)^(3*n-1/2) * n^n / exp(2*n/(1+sqrt(5))). - Vaclav Kotesovec, Aug 07 2014
a(n) = Sum_{k = 1..n} A060281(n,k) n^k. - David Einstein, Oct 31 2016
a(n) = n! * [x^n] 1/(1 + LambertW(-x))^n. - Ilya Gutkovskiy, Oct 03 2017
MAPLE
b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
add(b(n-j, i-1)*binomial(n, j)*j^j, j=0..n)))
end:
a:= n-> b(n$2):
seq(a(n), n=0..20); # Alois P. Heinz, Jul 17 2014
MATHEMATICA
f4[n_] := Sum[n^k Sum[Binomial[n - 1, j]*n^(n - 1 - j)*StirlingS1[j + 1, k] *(-1)^(j + k + 1), {j, 0, n - 1}], {k, 1, n}] (* David Einstein, Oct 31 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Chad Brewbaker, Mar 26 2014
EXTENSIONS
a(6)-a(7) from Giovanni Resta, Mar 28 2014
a(8)-a(17) from Alois P. Heinz, Jul 17 2014
STATUS
approved