login
A341227
Expansion of (-1 + Product_{k>=1} 1 / (1 - x^k))^8.
11
1, 16, 136, 824, 4004, 16608, 61076, 204200, 631714, 1831752, 5027312, 13159104, 33049090, 80030808, 187613348, 427201176, 947520103, 2051989360, 4347996772, 9030416704, 18412343832, 36905322248, 72807201940, 141525042736, 271321432489, 513454659312
OFFSET
8,2
LINKS
MAPLE
b:= proc(n, k) option remember; `if`(k<2, `if`(n=0, 1-k, combinat[
numbpart](n)), (q-> add(b(j, q)*b(n-j, k-q), j=0..n))(iquo(k, 2)))
end:
a:= n-> b(n, 8):
seq(a(n), n=8..33); # Alois P. Heinz, Feb 07 2021
MATHEMATICA
nmax = 33; CoefficientList[Series[(-1 + Product[1/(1 - x^k), {k, 1, nmax}])^8, {x, 0, nmax}], x] // Drop[#, 8] &
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Feb 07 2021
STATUS
approved