login
E.g.f.: Product_{k>0} (1+x^k)^(-1/k).
6

%I #13 Oct 30 2017 10:00:02

%S 1,-1,1,-5,23,-119,619,-4759,48145,-476657,4249961,-48286061,

%T 691331431,-9132207655,117900772963,-2025161870159,37607411624609,

%U -628236985455329,10768798391659345,-215626810984559317,4751529623277906871,-105427459848063440471

%N E.g.f.: Product_{k>0} (1+x^k)^(-1/k).

%t 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 *)

%o (PARI) N=66; x='x+O('x^N); Vec(serlaplace(prod(k=1, N, 1/(1+x^k)^(1/k))))

%Y Cf. A028342, A028343, A168243.

%K sign

%O 0,4

%A _Seiichi Manyama_, Oct 29 2017