login
G.f. satisfies: A(x) = 1 + x*A(x)^4*A(x*A(x)^4).
3

%I #5 Aug 13 2012 23:21:27

%S 1,1,5,39,387,4516,59522,864981,13641761,230931591,4162594198,

%T 79412758356,1595971048438,33662065586127,742852904575909,

%U 17107642358401929,410242691108675077,10223836926997402818,264339609703677876178,7079673476173815446742

%N G.f. satisfies: A(x) = 1 + x*A(x)^4*A(x*A(x)^4).

%F a(n) = coefficient of x^n in (1+x*A(x))^(4*n+1)/(4*n+1) where A(x) = Sum_{n=0} a(n)*x^n.

%F Recurrence:

%F Let A(x)^m = Sum_{n>=0} a(n,m)*x^n with a(0,m)=1, then

%F a(n,m) = Sum_{k=0..n} m*binomial(4*n+m,k)/(4*n+m) * a(n-k,k).

%F G.f. A(x) = F(x,1) where F(x,n) satisfies: F(x,n) = F(x,n-1)*(1 + x*F(x,n)^3*F(x,n+1)) for n>0 with F(x,0)=1.

%e G.f.: A(x) = 1 + x + 5*x^2 + 39*x^3 + 387*x^4 + 4516*x^5 + 59522*x^6 +...

%e G.f. A(x) is the unique solution to variable A in the infinite system of simultaneous equations:

%e A = 1 + x*A^3*B;

%e B = A*(1 + x*B^3*C);

%e C = B*(1 + x*C^3*D);

%e D = C*(1 + x*D^3*E);

%e E = D*(1 + x*E^3*F); ...

%e where B(x) = A(x)*A(x*A(x)^4), C(x) = A(x)*B(x*A(x)^4), D(x) = A(x)*C(x*A(x)^4), ...

%e Expansions of a few of the functions described above begin:

%e B(x) = 1 + 2*x + 15*x^2 + 158*x^3 + 1991*x^4 + 28283*x^5 + 439880*x^6 +...

%e C(x) = 1 + 3*x + 30*x^2 + 402*x^3 + 6257*x^4 + 107405*x^5 + 1983027*x^6 +...

%e D(x) = 1 + 4*x + 50*x^2 + 816*x^3 + 15215*x^4 + 308986*x^5 + 6679230*x^6 +...

%e E(x) = 1 + 5*x + 75*x^2 + 1445*x^3 + 31480*x^4 + 741075*x^5 + 18447355*x^6 +...

%e ALTERNATE GENERATING METHOD.

%e Suppose functions A=A(x), B=B(x), C=C(x), etc., satisfy:

%e A = 1 + x*A^4*B,

%e B = 1 + x*(A*B)^4*C,

%e C = 1 + x*(A*B*C)^4*D,

%e D = 1 + x*(A*B*C*D)^4*E, etc.,

%e then B(x) = A(x*A(x)^4), C(x) = B(x*A(x)^4), D(x) = C(x*A(x)^4), etc.,

%e where A(x) = 1 + x*A(x)^4*A(x*A(x)^4) is the g.f. of this sequence.

%e Expansions of a few of the functions described above begin:

%e B(x) = 1 + x + 9*x^2 + 105*x^3 + 1415*x^4 + 21089*x^5 + 340100*x^6 +...

%e C(x) = 1 + x + 13*x^2 + 203*x^3 + 3507*x^4 + 65018*x^5 + 1274266*x^6 +...

%e D(x) = 1 + x + 17*x^2 + 333*x^3 + 7047*x^4 + 157359*x^5 + 3665076*x^6 +...

%e E(x) = 1 + x + 21*x^2 + 495*x^3 + 12419*x^4 + 325312*x^5 + 8813506*x^6 +...

%o (PARI) {a(n)=local(A=1+x); for(i=0, n, A=1+x*A^4*subst(A, x, x*A^4+x*O(x^n))); polcoeff(A, n)}

%o (PARI) /* a(n) = [x^n] (1+x*A(x))^(4*n+1)/(4*n+1): */

%o {a(n)=local(A=1+x); for(i=0, n, A=sum(m=0,n,polcoeff((1+x*A+x*O(x^m))^(4*m+1)/(4*m+1),m)*x^m)+x*O(x^n));polcoeff(A,n)}

%o (PARI) {a(n, m=1)=if(n==0, 1, if(m==0, 0^n, sum(k=0, n, m*binomial(4*n+m, k)/(4*n+m)*a(n-k, k))))}

%o for(n=0,41,print1(a(n),", "))

%Y Cf. A088717, A215505, A215507.

%K nonn

%O 0,3

%A _Paul D. Hanna_, Aug 13 2012