OFFSET
0,4
MATHEMATICA
nmax = 20; CoefficientList[Series[Product[1/(1+x^k)^(1/k), {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Oct 29 2017 *)
PROG
(PARI) N=66; x='x+O('x^N); Vec(serlaplace(prod(k=1, N, 1/(1+x^k)^(1/k))))
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Oct 29 2017
STATUS
approved