OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..30
MAPLE
b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
b(n, i-1, k)+b(n-i, min(n-i, i), k)*((i-1)!*binomial(n, i))^k))
end:
a:= n-> b(n$3):
seq(a(n), n=0..12); # Alois P. Heinz, Jul 27 2023
MATHEMATICA
Table[(n!)^n SeriesCoefficient[Product[1/(1 - x^k/k^n), {k, 1, n}], {x, 0, n}], {n, 0, 10}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Jul 16 2020
STATUS
approved