|
%I #5 May 11 2014 04:46:32
%S 0,2,-6,42,-468,7080,-133128,2938824,-73169568,1997384832,
%T -58814501760,1868053207680,-65311214042880,2585560450337280,
%U -115344597684718080,5424254194395456000,-244310147229735014400,10256126830544041574400
%N E.g.f. satisfies: A(x) = (1+x)/(1+3*x) * A(x*(1+x)^2).
%F E.g.f. A=A(x) satisfies:
%F . (1+x)^2 = 1 + A + A*Dx(A)/2! + A*Dx(A*Dx(A))/3! + A*Dx(A*Dx(A*Dx(A)))/4! +...
%F . (1+x)^2*(1+x*(1+x)^2)^2 = 1 + 2*A + 2^2*A*Dx(A)/2! + 2^3*A*Dx(A*Dx(A))/3! + 2^4*A*Dx(A*Dx(A*Dx(A)))/4! +...
%F . G001764(-x)^2 = 1 - A + A*Dx(A)/2! - A*Dx(A*Dx(A))/3! + A*Dx(A*Dx(A*Dx(A)))/4! -+...; G001764(x) = g.f. of A001764;
%F where Dx(F) = d/dx(x*F).
%F INVERSION FORMULA:
%F More generally, if A(x) = A(G(x)) * G(x)/(x*G'(x)) with G(0)=0, G'(0)=1,
%F then G(x) can be obtained from A=A(x) by the series:
%F . G(x)/x = 1 + A + A*Dx(A)/2! + A*Dx(A*Dx(A))/3! + A*Dx(A*Dx(A*Dx(A)))/4! +... where Dx(F) = d/dx(x*F).
%F ITERATION FORMULA:
%F Let G_n(x) denote the n-th iteration of G(x) = x*(1+x)^2, and A=A(x), then:
%F . A(x) = A(G_n(x)) * G_n(x)/(x*G_n'(x)) for all n;
%F . G_n(x)/x = 1 + n*A + n^2*A*Dx(A)/2! + n^3*A*Dx(A*Dx(A))/3! + n^4*A*Dx(A*Dx(A*Dx(A)))/4! +... where Dx(F) = d/dx(x*F).
%F ...
%F MATRIX LOG OF RIORDAN ARRAY (G(x)/x, G(x)) where G(x) = x*(1+x)^2:
%F . k*A(x) = e.g.f. of column k of the matrix log of triangle A116088 for k>=0.
%e E.g.f.: A(x) = 2*x - 6*x^2/2! + 42*x^3/3! - 468*x^4/4! + 7080*x^5/5! - 133128*x^6/6! + 2938824*x^7/7! - 73169568*x^8/8! + 1997384832*x^9/9! - 58814501760*x^10/10! + 1868053207680*x^11/11! - 65311214042880*x^12/12! +...
%e ...
%e A(x*(1+x)^2) = 2*x + 2*x^2/2! - 18*x^3/3! + 108*x^4/4! - 480*x^5/5! - 2808*x^6/6! + 162792*x^7/7! - 3940128*x^8/8! + 57267648*x^9/9! + 534366720*x^10/10! - 78703384320*x^11/11! + 2883142045440*x^12/12! +...
%e ...
%e where A(x*(1+x)^2) = (1+3*x)/(1+x) * A(x).
%e ...
%e Related expansions begin:
%e . A = 2*x - 6*x^2/2! + 42*x^3/3! - 468*x^4/4! + 7080*x^5/5! +...
%e . A*Dx(A)/2! = 8*x^2/2! - 90*x^3/3! + 1332*x^4/4! - 25200*x^5/5! +...
%e . A*Dx(A*Dx(A))/3! = 48*x^3/3! - 1248*x^4/4! + 32760*x^5/5! -+...
%e . A*Dx(A*Dx(A*Dx(A)))/4! = 384*x^4/4! - 18480*x^5/5! + 770400*x^6/6! -+...
%e . A*Dx(A*Dx(A*Dx(A*Dx(A))))/5! = 3840*x^5/5! - 300672*x^6/6! +-...
%e ...
%e Sums of which generate the square of the g.f. of A001764:
%e . G001764(-x)^2 = 1 - A + A*Dx(A)/2! - A*Dx(A*Dx(A))/3! + A*Dx(A*Dx(A*Dx(A)))/4! -+...
%e . G001764(-x)^2 = 1 - 2*x + 7*x^2 - 30*x^3 + 143*x^4 - 728*x^5 +...+ A006013(n)*(-x)^n +...
%e ...
%e The Riordan array ((1+x)^2, x*(1+x)^2) (cf. A116088) begins:
%e 1;
%e 2, 1;
%e 1, 4, 1;
%e 0, 6, 6, 1;
%e 0, 4, 15, 8, 1;
%e 0, 1, 20, 28, 10, 1;
%e 0, 0, 15, 56, 45, 12, 1; ...
%e The matrix log of Riordan array ((1+x)^2, x*(1+x)^2) begins:
%e 0;
%e 2, 0;
%e -6/2!, 4, 0;
%e 42/3!, -12/2!, 6, 0;
%e -468/4!, 84/3!, -18/2!, 8, 0;
%e 7080/5!, -936/4!, 126/3!, -24/2!, 10, 0;
%e -133128/6!, 14160/5!, -1404/4!, 168/3!, -30/2!, 12, 0; ...
%e where the g.f. of the leftmost column equals the e.g.f. of this sequence.
%o (PARI) /* E.g.f. satisfies: A(x) = (1+x)/(1+3*x)*A(x*(1+x)^2): */
%o {a(n)=local(A=2*x, B); for(m=2, n, B=(1+x)/(1+3*x+O(x^(n+3)))*subst(A,x,x*(1+x)^2+O(x^(n+3))); A=A-polcoeff(B, m+1)*x^m/(m-1)/2); n!*polcoeff(A, n)}
%o (PARI) /* (1+x)^2 = 1 + A + A*Dx(A)/2! + A*Dx(A*Dx(A))/3! +...: */
%o {a(n)=local(A=0+sum(m=1,n-1,a(m)*x^m/m!),D=1,R=0);R=-((1+x)^2+x*O(x^n))+1+sum(m=1,n,(D=A*deriv(x*D+x*O(x^n)))/m!);-n!*polcoeff(R,n)}
%o (PARI) /* First column of the matrix log of triangle A116088: */
%o {a(n)=local(M=matrix(n+1, n+1, r, c, if(r>=c, polcoeff(((1+x)^2+x*O(x^n))^c,r-c))), LOG, ID=M^0); LOG=sum(m=1, n+1, -(ID-M)^m/m); n!*LOG[n+1, 1]}
%Y Cf. A179331, variants: A179320, A179420.
%K sign
%O 0,2
%A _Paul D. Hanna_, Jul 21 2010
| 