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

%I #20 Apr 15 2023 10:17:45

%S 1,2,18,224,3230,50688,840420,14483456,256856886,4656988160,

%T 85929839996,1608379269120,30463651429484,582796191989760,

%U 11245047027447240,218581150665277440,4276257634911525670,84135742205488791552,1663738200769421021580,33047906167191995678720

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

%C Radius of convergence of g.f. A(x) is r = (3/5)^(5/2) / 6 where A(r) = sqrt(5/3).

%H Gi-Sang Cheon, S.-T. Jin, and L. W. Shapiro, <a href="https://doi.org/10.1016/j.laa.2015.03.015">A combinatorial equivalence relation for formal power series</a>, Linear Algebra and its Applications, Volume 491, 15 February 2016, Pages 123-137.

%F a(n) = 4^n * binomial((5*n - 1)/2, n) / (3*n + 1).

%F G.f. A(x) satisfies: A(x) = sqrt(1 + 4*x^2*A(x)^8) + 2*x*A(x)^4.

%F Self convolution yields A214553.

%F G.f. A(x) = 1/x * series reversion of x*sqrt(1 - 4*x*C(4*x)), where C(x) = (1 - sqrt(1 - 4*x))/(2*x) is the generating function of the Catalan numbers A000108. See A024491. - _Peter Bala_, Mar 27 2023

%e G.f.: A(x) = = 1 + 2*x + 18*x^2 + 224*x^3 + 3230*x^4 + 50688*x^5 +...

%e where A(x)^2 = 1 + 4*x*A(x)^5:

%e A(x)^2 = 1 + 4*x + 40*x^2 + 520*x^3 + 7680*x^4 + 122360*x^5 +...

%e A(x)^5 = 1 + 10*x + 130*x^2 + 1920*x^3 + 30590*x^4 + 512512*x^5 +...

%e Related series:

%e A(x)^4 = 1 + 8*x + 96*x^2 + 1360*x^3 + 21120*x^4 + 347760*x^5 +...

%e A(x)^8 = 1 + 16*x + 256*x^2 + 4256*x^3 + 73216*x^4 + 1294560*x^5 +...

%e where A(x) = sqrt(1 + 4*x^2*A(x)^8) + 2*x*A(x)^4.

%o (PARI) /* From A(x)^2 = 1 + 4*x*A(x)^5 : */

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

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

%o (PARI) {a(n) = 4^n * binomial((5*n - 1)/2, n) / (3*n + 1)}

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

%o (PARI) /* From A(x) = sqrt(1 + 4*x^2*A(x)^8) + 2*x*A(x)^4 : */

%o {a(n) = local(A=1+x);for(i=1,n,A = sqrt(1 + 4*x^2*A^8 +x*O(x^n)) + 2*x*A^4);polcoeff(A,n)}

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

%Y Cf. A000108, A024491, A214553, A245113.

%K nonn,easy

%O 0,2

%A _Paul D. Hanna_, Jul 31 2014