OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..150
FORMULA
E.g.f.: 1/(1- T(x) + T(x)^2/2) where T(x) is e.g.f. for A000169.
a(n) = Sum_{m=1..n} A105599(n,m)*m!.
a(n) ~ 4*n^(n-2). - Vaclav Kotesovec, Aug 16 2013
a(0) = 1; a(n) = Sum_{k=1..n} binomial(n,k) * k^(k-2) * a(n-k). - Ilya Gutkovskiy, Jan 26 2020
MAPLE
T:= -LambertW(-x):
egf:= 1/(1-T+T^2/2):
a:= n-> n! * coeff(series(egf, x, n+1), x, n):
seq(a(n), n=0..30); # Alois P. Heinz, Nov 04 2012
MATHEMATICA
nn=20; t=Sum[n^(n-1)x^n/n!, {n, 1, nn}]; Range[0, nn]!CoefficientList[ Series[1/(1-t+t^2/2), {x, 0, nn}], x]
PROG
(PARI) A218688_vec(n, A=List(1))={until(#A>n, listput(A, sum(k=1, #A, binomial(#A, k)*k^(k-2)*A[#A-k+1]))); Vec(A)} \\ M. F. Hasler, Jan 26 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Geoffrey Critzer, Nov 04 2012
STATUS
approved