OFFSET
0,4
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..10000
MAPLE
b:= proc(n, i) option remember; `if`(i*(i+1)/2<n, 0,
`if`(n=0, 1, b(n, i-1)+(i-1)*b(n-i, min(n-i, i-1))))
end:
a:= n-> b(n$2):
seq(a(n), n=0..42); # Alois P. Heinz, Aug 15 2019
MATHEMATICA
nmax = 50; CoefficientList[Series[Product[1+(k-1)*x^k, {k, 1, nmax}], {x, 0, nmax}], x]
nmax = 50; poly = ConstantArray[0, nmax+1]; poly[[1]] = 1; poly[[2]] = 0; Do[Do[poly[[j+1]] += (k-1)*poly[[j-k+1]], {j, nmax, k, -1}]; , {k, 2, nmax}]; poly
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Jan 08 2016
STATUS
approved