# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/

Search: id:a144683
Showing 1-1 of 1

%I A144683 #5 Mar 30 2012 18:37:14
%S A144683 1,1,7,154,6625,446776,42088609,5171653432,792466370161,
%T A144683 146737621112464,32079442602647521,8134165156479090136,
%U A144683 2358873194743497162889,773523517692799444058632,284154419348084944647780289
%N A144683 E.g.f. satisfies: A(x/A(x)^3) = exp(x).
%F A144683 E.g.f. satisfies: A(x) = exp( x*A(log A(x))^3 ).
%F A144683 E.g.f. satisfies: a(n+1) = [x^n/n!] exp(x)*A(x)^(3*n+3).
%F A144683 E.g.f. A(x) = G(3x)^(1/3) where G(x/G(x)) = exp(x) and G(x) is the g.f. of A144681.
%e A144683 E.g.f. A(x) = 1 + x + 7*x^2/2! + 154*x^3/3! + 6625*x^4/4! +...
%e A144683 ILLUSTRATE FORMULA a(n+1) = [x^n/n!] exp(x)*A(x)^(3*n+3) as follows.
%e A144683 Form a table of coefficients of x^k/k! in exp(x)*A(x)^(3*n) for n>=1, k>=0:
%e A144683 exp(x)*A(x)^3: [(1), 4, 34, 685, 27256, 1747159, 159049000, ...];
%e A144683 exp(x)*A(x)^6: [1, (7), 85, 1909, 75193, 4654417, 410053357, ...];
%e A144683 exp(x)*A(x)^9: [1, 10,(154), 3835, 153208, 9284725, 795356632, ...];
%e A144683 exp(x)*A(x)^12:[1, 13, 241, (6625), 272641, 16385713, 1373165425, ...];
%e A144683 exp(x)*A(x)^15:[1, 16, 346, 10441,(446776), 26918851, 2221660936, ...];
%e A144683 exp(x)*A(x)^18:[1, 19, 469, 15445, 690841, (42088609), 3443635405, ...];
%e A144683 exp(x)*A(x)^21:[1, 22, 610, 21799, 1022008, 63371617,(5171653432), ...]; ...
%e A144683 then the terms along the main diagonal form this sequence shift left.
%o A144683 (PARI) {a(n)=local(A=1+x+x*O(x^n));for(n=0,n,A=exp(serreverse(x/A^3)));n!*polcoeff(A,n)}
%o A144683 (PARI) {a(n)=local(A=1+x+sum(k=2, n-1, a(k)*x^k/k!)+x*O(x^n)); if(n==0, 1, (n-1)!*polcoeff(exp(x+x*O(x^n))*A^(3*n), n-1))}
%Y A144683 Cf. A144681, A144682, A144684.
%K A144683 nonn
%O A144683 0,3
%A A144683 _Paul D. Hanna_, Sep 19 2008

# Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE