OFFSET
0,3
MATHEMATICA
Table[SeriesCoefficient[Product[(1 + x^k)^Binomial[k + n - 1, n - 1], {k, 1, n}], {x, 0, n}], {n, 0, 21}]
A[n_, k_] := A[n, k] = If[n == 0, 1, (1/n) Sum[Sum[(-1)^(j/d + 1) d Binomial[d + k - 1, k - 1], {d, Divisors[j]}] A[n - j, k], {j, 1, n}]]; a[n_] := A[n, n]; Table[a[n], {n, 0, 21}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 09 2021
STATUS
approved