cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A213591 G.f. A(x) satisfies A( x - A(x)^2 ) = x.

Original entry on oeis.org

1, 1, 4, 24, 178, 1512, 14152, 142705, 1528212, 17211564, 202460400, 2474708496, 31310415376, 408815254832, 5495451727376, 75907303147652, 1075685334980240, 15618612118252960, 232102241507321384, 3526880759915999016, 54755450619399484512, 867928449982022915984
Offset: 1

Views

Author

Paul D. Hanna, Jun 14 2012

Keywords

Comments

Unsigned version of A139702.
Self-convolution is A276370.
Row sums of triangle A277295.

Examples

			G.f.: A(x) = x + x^2 + 4*x^3 + 24*x^4 + 178*x^5 + 1512*x^6 + 14152*x^7 +...
where A(x) = x + A(A(x))^2:
A(A(x)) = x + 2*x^2 + 10*x^3 + 69*x^4 + 568*x^5 + 5250*x^6 + 52792*x^7 +...
A(A(x))^2 = x^2 + 4*x^3 + 24*x^4 + 178*x^5 + 1512*x^6 + 14152*x^7 +...
The g.f. satisfies the series:
A(x) = x + A(x)^2 + d/dx A(x)^4/2! + d^2/dx^2 A(x)^6/3! + d^3/dx^3 A(x)^8/4! +...
Logarithmic series:
log(A(x)/x) = A(x)^2/x + [d/dx A(x)^4/x]/2! + [d^2/dx^2 A(x)^6/x]/3! + [d^3/dx^3 A(x)^8/x]/4! +...
Also, A(x) = x*G(A(x)^2/x) where G(x) = x/A(x/G(x)^2) is the g.f. of A212411:
G(x) = 1 + x + 2*x^2 + 7*x^3 + 36*x^4 + 235*x^5 + 1792*x^6 + 15261*x^7 +...
Also, A(x)^2 = x*F(A(x)) where F(x) is the g.f. of A213628:
F(x) = x + x^2 + 3*x^3 + 14*x^4 + 85*x^5 + 616*x^6 + 5072*x^7 + 46013*x^8 +...

Links

Crossrefs

Cf. A220379, A139702, A212411, A138740, A213628, A213639.
Cf. A088714, A088717, A091713, A120971, A140094, A140095.
Cf. A143426, A087949, A143435, A182969.
Cf. A275765, A276360, A276361, A276362, A276363, A276370.
Cf. A277295 (triangle).

Programs

  • Mathematica
    terms = 22; A[] = 0; Do[A[x] = x + A[A[x]]^2 + O[x]^(terms+1) // Normal, terms+1]; CoefficientList[A[x], x] // Rest (* Jean-François Alcover, Jan 09 2018 *)
  • PARI
    {a(n)=local(A=x); if(n<1, 0, for(i=1, n, A=serreverse(x - A^2+x*O(x^n))); polcoeff(A, n))}
  • PARI
    {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x+x^2+x*O(x^n)); for(i=1, n, A=x+sum(m=1, n, Dx(m-1, A^(2*m))/m!)+x*O(x^n)); polcoeff(A, n)}
  • PARI
    {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x+x^2+x*O(x^n)); for(i=1, n, A=x*exp(sum(m=1, n, Dx(m-1, A^(2*m)/x)/m!)+x*O(x^n))); polcoeff(A, n)}
for(n=1,21,print1(a(n),", "))
  • PARI
    b(n, k) = if(k==0, 0^n, k*sum(j=0, n, binomial(n+j+k, j)/(n+j+k)*b(n-j, 2*j)));
a(n) = b(n-1, 1); \\ Seiichi Manyama, Jun 05 2025

Formula

G.f. satisfies:
(1) A(x) = x + A(A(x))^2.
(2) A(x) = x + Sum_{n>=1} d^(n-1)/dx^(n-1) A(x)^(2*n) / n!.
(3) A(x) = x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) A(x)^(2*n)/x / n! ).
(4) A(x) = x*G(A(x)^2/x) where G(x) = 1 + x*G(1-1/G(x))^2 is the g.f. of A212411.
(5) A(x)^2 = x*F(A(x)) where F(x) = 1 - (x^2/F(x))/F(x^2/F(x)) is the g.f. of A213628.
(6) x = A(A( x-x^2 - A(x)^2 )). - Paul D. Hanna, Jul 01 2012
(7) A(x) is the unique solution to variable A in the infinite system of simultaneous equations starting with:
A = x + B^2;
B = A + C^2;
C = B + D^2;
D = C + E^2; ...
where B = A(A(x)), C = A(A(A(x))), D = A(A(A(A(x)))), etc.
...
a(n) = Sum_{k=0..n-1} A277295(n,k).
From Seiichi Manyama, Jun 05 2025: (Start)
Let b(n,k) = [x^n] (A(x)/x)^k.
b(n,0) = 0^n; b(n,k) = k * Sum_{j=0..n} binomial(n+j+k,j)/(n+j+k) * b(n-j,2*j).
a(n) = b(n-1,1). (End)

A139702 G.f. satisfies: x = A( x + A(x)^2 ).

Original entry on oeis.org

1, -1, 4, -24, 178, -1512, 14152, -142705, 1528212, -17211564, 202460400, -2474708496, 31310415376, -408815254832, 5495451727376, -75907303147652, 1075685334980240, -15618612118252960, 232102241507321384, -3526880759915999016
Offset: 1

Views

Author

Paul D. Hanna, Apr 30 2008, May 20 2008

Keywords

Comments

Signed version of A213591.

Examples

			G.f.: A(x) = x - x^2 + 4*x^3 - 24*x^4 + 178*x^5 - 1512*x^6 +-...
A(x)^2 = x^2 - 2*x^3 + 9*x^4 - 56*x^5 + 420*x^6 - 3572*x^7 +-...
where A(x + A(x)^2) = x.
Let G(x) = Series_Reversion( A(x) ) = x + A(x)^2, then:
G(x) = x + x^2 - 2*x^3 + 9*x^4 - 56*x^5 + 420*x^6 -+... and
G(G(x)) = x + 2*x^2 - 2*x^3 + 9*x^4 - 56*x^5 + 420*x^6 -+...
so that G(x) = G(G(x)) - x^2 = g.f. of A138740.
Logarithmic series:
log(A(x)/x) = -A(x)^2/x + [d/dx A(x)^4/x]/2! - [d^2/dx^2 A(x)^6/x]/3! + [d^3/dx^3 A(x)^8/x]/4! -+...

Links

Crossrefs

Cf. A138740, A212411, A213591.
Cf. A088714, A088717, A091713, A120971, A140094, A140095.
Cf. A143426, A087949, A143435, A182969.

Programs

  • Mathematica
    nmax = 20; sol = {a[1] -> 1}; nmin = Length[sol]+1;
Do[A[x_] = Sum[a[k] x^k, {k, 0, n}] /. sol; eq = CoefficientList[x - A[x + A[x]^2] + O[x]^(n+1), x][[nmin;;]] == 0 /. sol; sol = sol ~Join~ Solve[eq][[1]], {n, nmin, nmax}];
a /@ Range[nmax] /. sol (* Jean-François Alcover, Nov 06 2019 *)
  • PARI
    {a(n)=local(A=x); if(n<1, 0, for(i=1,n, A=serreverse(x + (A+x*O(x^n))^2)); polcoeff(A, n))}
  • PARI
    /* n-th Derivative: */
{Dx(n,F)=local(D=F);for(i=1,n,D=deriv(D));D}
/* G.f.: [Paul D. Hanna, Dec 18 2010] */
{a(n)=local(A=x-x^2+x*O(x^n));for(i=1,n,
A=x*exp(sum(m=0,n,(-1)^(m+1)*Dx(m,A^(2*m+2)/x)/(m+1)!)+x*O(x^n)));polcoeff(A,n)}

Formula

Let G(x) = Series_Reversion( A(x) ) = x + A(x)^2, then G(x) = G(G(x)) - x^2 = g.f. of A138740.
G.f. satisfies: A(x) = x*G(-A(x)^2/x) where G(x) = 1 + x*G(1-1/G(x))^2 is the g.f. of A212411.
G.f.: A(x)/x is the unique solution to variable A in the infinite system of simultaneous equations starting with:
A = 1 - x*B^2;
B = A - x*C^2;
C = B - x*D^2;
D = C - x*E^2;
E = D - x*F^2; ...
G.f. satisfies: A(x) = x*exp( Sum_{n>=0} (-1)^(n+1)*[d^n/dx^n A(x)^(2n+2)/x]/(n+1)! ). [Paul D. Hanna, Dec 18 2010]

A140094 G.f. satisfies: A(x) = x/(1 - A(A(A(x)))).

Original entry on oeis.org

1, 1, 4, 25, 199, 1855, 19387, 221407, 2717782, 35455981, 487672243, 7029980797, 105732907498, 1653377947393, 26805765569863, 449568735630517, 7785116448484318, 138980739891821269, 2554369130466577138
Offset: 1

Views

Author

Paul D. Hanna, May 08 2008, May 20 2008

Keywords

Examples

			G.f.: A(x) = x + x^2 + 4*x^3 + 25*x^4 + 199*x^5 + 1855*x^6 + 19387*x^7 +...
Iterations A_{n+1}(x) = A( A_{n}(x) ) are related as follows.
A_2(x) = 1 - Series_Reversion( A(x) ) / x;
A_3(x) = 1 - x / A(x);
A_4(x) = 1 - A(x) / A_2(x);
A_5(x) = 1 - A_2(x) / A_3(x);
A_6(x) = 1 - A_3(x) / A_4(x); ...
where the iterations of A(x) begin:
A_2(x) = x + 2*x^2 + 10*x^3 + 71*x^4 + 616*x^5 + 6119*x^6 + 67210*x^7 +...;
A_3(x) = x + 3*x^2 + 18*x^3 + 144*x^4 + 1365*x^5 + 14544*x^6 +...;
A_4(x) = x + 4*x^2 + 28*x^3 + 250*x^4 + 2584*x^5 + 29584*x^6 +...;
A_5(x) = x + 5*x^2 + 40*x^3 + 395*x^4 + 4435*x^5 + 54515*x^6 +...;
A_6(x) = x + 6*x^2 + 54*x^3 + 585*x^4 + 7104*x^5 + 93555*x^6 +...;
...
Iterations are also related by continued fractions:
A(x) = x/(1 - A_2(x)/(1 - A_4(x)/(1 - A_6(x)/(1 -...)))) ;
A_2(x) = A(x)/(1 - A_3(x)/(1 - A_5(x)/(1 - A_7(x)/(1 -...)))).

Links

Crossrefs

Cf. A140095, A088714.
Cf. A088717, A091713, A120971, A139702.

Programs

  • PARI
    {a(n)=local(A); if(n<0, 0, n++; A=x+O(x^2); for(i=2, n, A=x/(1-subst(A, x, subst(A, x, A)))); polcoeff(A, n))}

Formula

G.f. A(x) satisfies:
(1) A(x) = Series_Reversion(x - x*A(A(x))).
(2) A(x) = x + Sum_{n=1} d^(n-1)/dx^(n-1) x^n * A(A(x))^n / n!.
(3) A(x) = x*exp( Sum_{n=1} d^(n-1)/dx^(n-1) x^n * A(A(x))^n/x / n! ).
Define A_{n} such that A_{n+1}(x) = A( A_{n}(x) ) with A_0(x) = x,
then A_{n}(x) = A_{n-1}/[1 - A_{n+2}(x)] ;
thus A_{n}(x) = 1 - A_{n-3}(x) / A_{n-2}(x).
G.f. A(x)/x is the unique solution to variable A in the infinite system of simultaneous equations starting with:
A = 1 + x*A*C;
B = A + x*B*D;
C = B + x*C*E;
D = C + x*D*F;
E = D + x*E*G; ...

Extensions

Name, formulas, and examples revised by Paul D. Hanna, Feb 03 2013
Showing 1-3 of 3 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.