|
%I #15 Oct 15 2017 12:47:45
%S 1,1,5,31,265,2621,31621,426595,6574961,111673945,2092318021,
%T 42552808871,937495160185,22150499622421,559765402811525,
%U 15039597200385451,428293292251548001,12875707199330296625,407547173842501629061
%N E.g.f.: exp(f(x)-1), where f(x) = o.g.f. for partitions (A000041), Product_{k>=1} 1/(1-x^k).
%H Seiichi Manyama, <a href="/A058892/b058892.txt">Table of n, a(n) for n = 0..422</a>
%F a(0) = 1 and a(n) = (n-1)! * Sum_{k=1..n} k*A000041(k)*a(n-k)/(n-k)! for n > 0. - _Seiichi Manyama_, Oct 15 2017
%t nmax = 30; CoefficientList[Series[1/E*Exp[Product[1/(1 - x^k), {k, 1, nmax}]], {x, 0, nmax}], x] * Range[0, nmax]!(* _Vaclav Kotesovec_, Aug 19 2015 *)
%o (PARI)
%o N=66; q='q+O('q^N);
%o f=exp( 1/prod(n=1,N, 1-q^n ) - 1 );
%o egf=serlaplace(f);
%o Vec(egf)
%o /* _Joerg Arndt_, Oct 06 2012 */
%Y Cf. A000041, A038048.
%K nonn
%O 0,3
%A _N. J. A. Sloane_, Jan 08 2001
|
