OFFSET
1,3
MAPLE
A:= proc(n) option remember; series(x*`if`(n=1, 1, mul(mul(
1+a(j)*x^(j*k), k=1..(n-1)/j), j=1..n-1)), x, n+1)
end:
a:= n-> coeff(A(n), x, n):
seq(a(n), n=1..35); # Alois P. Heinz, May 14 2019
MATHEMATICA
a[n_] := a[n] = SeriesCoefficient[x Product[Product[(1 + a[j] x^(j k)), {k, 1, n - 1}], {j, 1, n - 1}], {x, 0, n}]; Table[a[n], {n, 1, 32}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, May 14 2019
STATUS
approved