|
%I #4 Dec 31 2012 17:26:03
%S 1,1,4,33,440,8380,211824,6771422,264621216,12329282160,671882721600,
%T 42181858413552,3013915118776704,242626985772839616,
%U 21821596448977248000,2176989083049432207600,239420370429753669425664
%N E.g.f. A(x) satisfies: x*A(x) equals column 0 in the matrix log of the Riordan array (A(x), x*A(x)).
%H Paul D. Hanna, <a href="/A179421/b179421.txt">Table of n, a(n) for n = 0..150</a>
%F a(n) = A179420(n+1)/(n+1).
%F E.g.f. A(x) satisfies: A(x) = x*A'(x)/(A(x*A(x)) - 1).
%F Let L_{n}(x) denote the n-th iteration of L=L(x)=x*A(x) then:
%F . L(L(x)) = x*L'(x),
%F . L_{n+1}(x) = L(x)*x*[d/dx L_n(x)]/L_n(x),
%F . L(x)/x = 1 + L + L*Dx(L)/2! + L*Dx(L*Dx(L))/3! + L*Dx(L*Dx(L*Dx(L)))/4! +...
%F . L_{n}(x)/x = 1 + n*L + n^2*L*Dx(L)/2! + n^3*L*Dx(L*Dx(L))/3! + n^4*L*Dx(L*Dx(L*Dx(L)))/4! +...
%F where Dx(F) = d/dx(x*F).
%e E.g.f.: A(x) = 1 + x + 4*x^2/2! + 33*x^3/3! + 440*x^4/4! +...
%e x*A(x) = x + 2*x^2/2! + 12*x^3/3! + 132*x^4/4! + 2200*x^5/5! +...
%e ...
%e The RIORDAN ARRAY (A(x), x*A(x)) begins:
%e . 1;
%e . 1, 1;
%e . 4/2!, 2, 1;
%e . 33/3!, 10/2!, 3, 1;
%e . 440/4!, 90/3!, 18/2!, 4, 1;
%e . 8380/5!, 1240/4!, 177/3!, 28/2!, 5, 1;
%e . 211824/6!, 23800/5!, 2544/4!, 300/3!, 40/2!, 6, 1;
%e . 6771422/7!, 598788/6!, 49680/5!, 4520/4!, 465/3!, 54/2!, 7, 1; ...
%e where the e.g.f. of column k = A(x)^(k+1) for k>=0.
%e ...
%e The MATRIX LOG of the above Riordan array (A(x), x*A(x)) begins:
%e . 0;
%e . 1, 0;
%e . 2/2!, 2, 0;
%e . 12/3!, 4/2!, 3, 0;
%e . 132/4!, 24/3!, 6/2!, 4, 0;
%e . 2200/5!, 264/4!, 36/3!, 8/2!, 5, 0;
%e . 50280/6!, 4400/5!, 396/4!, 48/3!, 10/2!, 6, 0;
%e . 1482768/7!, 100560/6!, 6600/5!, 528/4!, 60/3!, 12/2!, 7, 0; ...
%e where the e.g.f. of column k = (k+1)*x*A(x) for k>=0.
%o (PARI) {a(n)=local(A=1+x+sum(m=2,n-1,a(m)*x^m/m!)+x*O(x^n));if(n<2,n!*polcoeff(A,n),n!*polcoeff(subst(x*A,x,x*A)/x,n)/(n-1))}
%Y Cf. A179420, A179422, A179423, A179424.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Jul 15 2010
|
