|
%I #20 Sep 13 2018 03:04:34
%S 1,1,4,18,114,810,7140,68880,766920,9304680,125086080,1814015280,
%T 28588356720,481128888240,8678237087520,166041500264640,
%U 3371031116893440,72153115744469760,1627441316510929920,38500269726897538560,954533425718494702080
%N E.g.f.: 1/Product_{k>0} (1 - x^k/k)^k.
%H Seiichi Manyama, <a href="/A294469/b294469.txt">Table of n, a(n) for n = 0..443</a>
%F E.g.f.: exp(Sum_{k>=1} Sum_{j>=1} x^(j*k)/(k*j^(k-1))). - _Ilya Gutkovskiy_, Sep 12 2018
%t nmax = 20; CoefficientList[Series[Product[1/(1-x^k/k)^k, {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* _Vaclav Kotesovec_, Nov 01 2017 *)
%o (PARI) N=66; x='x+O('x^N); Vec(serlaplace(1/prod(k=1, N, (1-x^k/k)^k)))
%Y Cf. A007841, A181541, A294470, A294471.
%K nonn
%O 0,3
%A _Seiichi Manyama_, Oct 31 2017
|
