|
%I #5 Mar 30 2012 18:37:27
%S 1,1,6,48,977,13378,680739,16253410,638406407,21983092893,
%T 992808221997,42510501798854,20888455486970142,1827747975355320051,
%U 96671149737735542775,4002066332008089249033,193076352317855172664077,13301772601630930169228482
%N G.f. satisfies: A(x) = exp( Sum_{n>=} (x^n/n)*[Sum_{d|n} d*A(x)^d]^n ).
%e G.f.: A(x) = 1 + x + 6*x^2 + 48*x^3 + 977*x^4 + 13378*x^5 +...
%e The g.f. A = A(x) satisfies:
%e log(A) = x*A + x^2*(A + 2*A^2)^2/2 + x^3*(A + 3*A^3)^3/3 + x^4*(A + 2*A^2 + 4*A^4)^4/4 + x^5*(A + 5*A^5)^5/5 + x^6*(A + 2*A^2 + 3*A^3 + 6*A^6)^6/6 +...
%o (PARI) {a(n)=local(A=1+x);for(i=1,n,A=exp(sum(m=1,n,(x^m/m)*sumdiv(m,d,d*A^d)^m+x*O(x^n))));polcoeff(A,n)}
%K nonn
%O 0,3
%A _Paul D. Hanna_, Jul 09 2011
|
