login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A341962 G.f. B(x) satisfies: B(x) = (1 - x^2*B(x)^2) / (1 - 2*x*B(x))^2. 4
1, 4, 27, 224, 2070, 20444, 211239, 2255200, 24680862, 275408456, 3121711758, 35842765872, 415999482912, 4872611960268, 57524376899871, 683778868508352, 8176889078761590, 98303183244150968, 1187413394405415498, 14403838768083054208, 175394125615368312900 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

LINKS

Table of n, a(n) for n=0..20.

FORMULA

G.f.: B(x) = (1/x)*Series_Reversion( x*(1 - 2*x)^2 / (1 - x^2) ).

G.f. B = B(x) and related functions A = A(x),  C = C(x),  D = D(x), satisfy:

(1.a) A = 1/((1 - 2*x*B)*(1 - 3*x*C)).

(1.b) B = 1/((1 - x*A)*(1 - 3*x*C)).

(1.c) C = 1/((1 - x*A)*(1 - 2*x*B)).

(1.d) D = 1/((1 - x*A)*(1 - 2*x*B)*(1 - 3*x*C)).

(1.e) D = sqrt(A*B*C).

(2.a) A = (1 + 2*x*D)*(1 + 3*x*D).

(2.b) B = (1 + x*D)*(1 + 3*x*D).

(2.c) C = (1 + x*D)*(1 + 2*x*D).

(2.d) D = (sqrt(24*A + 1) - 5)/(12*x) = (sqrt(12*B + 4) - 4)/(6*x) = (sqrt(8*C + 1) - 3)/(4*x).

(3.a) A = B/(1 - x*B) = C/(1 - 2*x*C) = D/(1 + x*D).

(3.b) B = C/(1 - x*C) = A/(1 + x*A) = D/(1 + 2*x*D).

(3.c) C = A/(1 + 2*x*A) = B/(1 + x*B) = D/(1 + 3*x*D).

(3.d) D = A/(1 - x*A) = B/(1 - 2*x*B) = C/(1 - 3*x*C).

(3.e) 1 = (1 + x*A)*(1 - x*B) = (1 + 2*x*A)*(1 - 2*x*C) = (1 + x*B)*(1 - x*C).

(3.f) 1 = (1 - x*A)*(1 + x*D) = (1 - 2*x*B)*(1 + 2*x*D) = (1 - 3*x*C)*(1 + 3*x*D).

(4.a) A = (1 + x*A)*(1 + 2*x*A)/(1 - x*A)^2.

(4.b) B = (1 - x^2*B^2)/(1 - 2*x*B)^2.

(4.c) C = (1 - x*C)*(1 - 2*x*C)/(1 - 3*x*C)^2.

(4.d) D = (1 + x*D)*(1 + 2*x*D)*(1 + 3*x*D).

(5.a) A = (1/x)*Series_Reversion( x*(1 - x)^2 / ((1 + x)*(1 + 2*x)) ).

(5.b) B = (1/x)*Series_Reversion( x*(1 - 2*x)^2 / (1 - x^2) ).

(5.c) C = (1/x)*Series_Reversion( x*(1 - 3*x)^2 / ((1 - x)*(1 - 2*x)) ).

(5.d) D = (1/x)*Series_Reversion( x / ((1 + x)*(1 + 2*x)*(1 + 3*x)) ).

a(n) ~ sqrt((2*s - 5*r*s^2 + 2*r^2*s^3) / (Pi*(22 - 8*r*s))) / (n^(3/2)*r^(n + 1/2)), where r = 0.07627811703169412709742160523783922642030319519275992338... and s = 2.2734876474240885308453076415547165781643109176277594587178... are positive real roots of the system of equations (1 - r^2*s^2)/(1 - 2*r*s)^2 = s, -1 + 8*r^3*s^3 - 2*r^2*s*(1 + 6*s) + r*(4 + 6*s) = 0. - Vaclav Kotesovec, Mar 02 2021

EXAMPLE

G.f. B(x) = 1 + 4*x + 27*x^2 + 224*x^3 + 2070*x^4 + 20444*x^5 + 211239*x^6 + 2255200*x^7 + 24680862*x^8 + 275408456*x^9 + 3121711758*x^10 + ...

such that B(x) = 1/((1 - x*A(x))*(1 - 3*x*C(x))) where

A(x) = 1 + 5*x + 36*x^2 + 307*x^3 + 2880*x^4 + 28714*x^5 + 298620*x^6 + 3203183*x^7 + 35181792*x^8 + 393697030*x^9 + 4472679816*x^10 + ...

C(x) = 1 + 3*x + 20*x^2 + 165*x^3 + 1520*x^4 + 14982*x^5 + 154588*x^6 + 1648713*x^7 + 18029456*x^8 + 201063402*x^9 + 2277890472*x^10 + ...

RELATED SERIES.

D(x) = sqrt(A(x)*B(x)*C(x)) = 1 + 6*x + 47*x^2 + 420*x^3 + 4058*x^4 + 41286*x^5 + 435739*x^6 + 4726644*x^7 + 52373294*x^8 + 590247900*x^9 + 6744908118*x^10 + ...

D(x)^2 = A(x)*B(x)*C(x) = 1 + 12*x + 130*x^2 + 1404*x^3 + 15365*x^4 + 170748*x^5 + 1924762*x^6 + 21971760*x^7 + 253573386*x^8 + 2954377800*x^9 + ...

B(x)*C(x) = D(x) + x*D(x)^2 = 1 + 7*x + 59*x^2 + 550*x^3 + 5462*x^4 + 56651*x^5 + 606487*x^6 + 6651406*x^7 + 74345054*x^8 + ...

A(x)*C(x) = D(x) + 2*x*D(x)^2 = 1 + 8*x + 71*x^2 + 680*x^3 + 6866*x^4 + 72016*x^5 + 777235*x^6 + 8576168*x^7 + 96316814*x^8 + ...

A(x)*B(x) = D(x) + 3*x*D(x)^2 = 1 + 9*x + 83*x^2 + 810*x^3 + 8270*x^4 + 87381*x^5 + 947983*x^6 + 10500930*x^7 + 118288574*x^8 + ...

MATHEMATICA

CoefficientList[1/x * InverseSeries[Series[x*(1 - 2*x)^2/((1 - x^2)), {x, 0, 20}], x], x] (* Vaclav Kotesovec, Mar 02 2021 *)

PROG

(PARI) {b(n) = my(A=1, B=1, C=1); for(i=1, n,

A = 1/((1-2*x*B)*(1-3*x*C) +x*O(x^n));

B = 1/((1-1*x*A)*(1-3*x*C) +x*O(x^n));

C = 1/((1-1*x*A)*(1-2*x*B) +x*O(x^n)); );

polcoeff(B, n)}

for(n=0, 30, print1(b(n), ", "))

(PARI) /* By Series Reversion: */

{b(n) = my(B = 1/x*serreverse( x*(1 - 2*x)^2/((1 - x^2) +x*O(x^n)))); polcoeff(B, n)}

for(n=0, 30, print1(b(n), ", "))

CROSSREFS

Cf. A341961 (A(x)), A341963 (C(x)), A071878 (D(x)).

Sequence in context: A319518 A304045 A317103 * A276029 A160883 A328978

Adjacent sequences:  A341959 A341960 A341961 * A341963 A341964 A341965

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Feb 27 2021

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified November 27 20:49 EST 2021. Contains 349395 sequences. (Running on oeis4.)