%I #8 Jul 19 2013 05:07:17
%S 1,1,3,9,34,139,643,3198,17186,98438,598551,3842893,25957607,
%T 183858842,1361853786,10523285935,84648613474,707461189513,
%U 6132612218266,55050390426042,510994289786018,4898133459331104,48424584171850411
%N G.f. satisfies: A(x) = x*exp( Sum_{n>=1} A(x^n/(1-x^n))/n ).
%C Compare the g.f. to the g.f. P(x) of the partition numbers (A000041):
%C * P(x) = exp( Sum_{n>=1} (x^n/(1-x^n))/n ),
%C and to the g.f. R(x) of rooted trees with n nodes (A000081):
%C * R(x) = x*exp( Sum_{n>=1} R(x^n)/n ).
%C Consider also the trivial identity:
%C * B(x) = x*exp( Sum_{n>=1} B(x^n/(1+x^n))/n ) where B(x) = x/(1-x).
%e G.f.: A(x) = x + x^2 + 3*x^3 + 9*x^4 + 34*x^5 + 139*x^6 + 643*x^7 +...
%e where the g.f. A(x) satisfies:
%e log(A(x)/x) = A(x/(1-x)) + A(x^2/(1-x^2))/2 + A(x^3/(1-x^3))/3 + A(x^4/(1-x^4))/4 +...
%e Explicitly, the logarithmic series begins:
%e log(A(x)/x) = x + 5*x^2/2 + 19*x^3/3 + 93*x^4/4 + 466*x^5/5 + 2633*x^6/6 + 15534*x^7/7 + 97645*x^8/8 +...
%o (PARI) {a(n)=local(A=x+x^2);for(i=1,n,A=x*exp(sum(m=1,n,subst(A,x,x^m/(1-x^m+x*O(x^n)))/m)));polcoeff(A,n)}
%K nonn
%O 1,3
%A _Paul D. Hanna_, Jun 01 2011