nmax = 21; CoefficientList[Series[Product[(1 + Log[1 + x]^k), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
nmax = 21; CoefficientList[Series[Exp[Sum[(-1)^(k + 1) Log[1 + x]^k/(k (1 - Log[1 + x]^k)), {k, 1, nmax}]], {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[StirlingS1[n, k] PartitionsQ[k] k!, {k, 0, n}], {n, 0, 21}]