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-10 of 22 results. Next

A007297 Number of connected graphs on n labeled nodes on a circle with straight-line edges that don't cross.

Original entry on oeis.org

1, 1, 4, 23, 156, 1162, 9192, 75819, 644908, 5616182, 49826712, 448771622, 4092553752, 37714212564, 350658882768, 3285490743987, 30989950019532, 294031964658430, 2804331954047160, 26870823304476690, 258548658860327880
Offset: 1

Views

Author

Keywords

Comments

Apart from the initial 1, reversion of g.f. for A162395 (squares with signs): see A263843.

Examples

			G.f. = x*(1 + x + 4*x^2 + 23*x^3 + 156*x^4 + 1162*x^5 + 9192*x^6 + 75819*x^7 + ...).
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A162395, A000290. 4th row of A107111. Row sums of A089434.
See A263843 for a variant.
Cf. A000108 (non-crossing set partitions), A001006, A001187, A054726 (non-crossing graphs), A054921, A099947, A194560, A293510, A323818, A324167, A324169, A324173.

Programs

  • Maple
    A007297:=proc(n) if n = 1 then 1 else add(binomial(3*n - 3, n + j)*binomial(j - 1, j - n + 1), j = n - 1 .. 2*n - 3)/(n - 1); fi; end;
  • Mathematica
    CoefficientList[ InverseSeries[ Series[(x-x^2)/(1+x)^3, {x, 0, 20}], x], x] // Rest (* From Jean-François Alcover, May 19 2011, after PARI prog. *)
    Table[Binomial[3n, 2n+1] Hypergeometric2F1[1-n, n, 2n+2, -1]/n, {n, 1, 20}] (* Vladimir Reshetnikov, Oct 25 2015 *)
  • PARI
    a(n)=if(n<0,0,polcoeff(serreverse((x-x^2)/(1+x)^3+O(x^(n+2))),n+1)) /* Ralf Stephan */

Formula

Apart from initial term, g.f. is the series reversion of (x-x^2)/(1+x)^3 (A162395). See A263843. - Vladimir Kruchinin, Feb 08 2013
G.f.: (g-z)/z, where g=-1/3+(2/3)*sqrt(1+9z)*sin((1/3)*arcsin((2+27z+54z^2)/2/(1+9*z)^(3/2))). - Emeric Deutsch, Dec 02 2002
a(n) = (1/n)*Sum_{k=0..n} binomial(3n, n-k-1)*binomial(n+k-1, k). - Paul Barry, May 11 2005
a(n) = 4^(n-1)*(Gamma(3*n/2-1)/Gamma(n/2+1)/Gamma(n) -Gamma((3*n-1)/2)/ Gamma( (n+1)/2)/Gamma(n+1)). - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = 4^n * binomial(3*n/2, n/2) / (9*n-6) - 4^(n-1) * binomial(3*(n-1)/2, (n-1)/2 ) / n. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
D-finite with recurrence: n*(n-1)*(3*n-4)*a(n) +36*(n-1)*a(n-1) -12*(3*n-8)*(3*n-1)*(3*n-7)*a(n-2)=0. - Mark van Hoeij, Aug 27 2005, adapted to offset Feb 21 2020 by R. J. Mathar
a(n) = (1/n)*Sum_{k=0..n} C(3n, k)*C(2n-k-2, n-1). - Paul Barry, Sep 27 2005
a(n) ~ (2-sqrt(3)) * 6^n * 3^(n/2) / (sqrt(2*Pi) * n^(3/2)). - Vaclav Kotesovec, Mar 17 2014
a(n) = binomial(3*n,2*n+1)*hypergeom([1-n,n], [2*n+2], -1)/n. - Vladimir Reshetnikov, Oct 25 2015
a(n) = 2*A078531(n) - A085614(n+1). - Vladimir Reshetnikov, Apr 24 2016

Extensions

Better description from Philippe Flajolet, Apr 20 2000
More terms from James Sellers, Aug 21 2000
Definition revised and initial a(1)=1 added by N. J. A. Sloane, Nov 05 2015 at the suggestion of Axel Boldt. Some of the formulas may now need to be adjusted slightly.

A048990 Catalan numbers with even index (A000108(2*n), n >= 0): a(n) = binomial(4*n, 2*n)/(2*n+1).

Original entry on oeis.org

1, 2, 14, 132, 1430, 16796, 208012, 2674440, 35357670, 477638700, 6564120420, 91482563640, 1289904147324, 18367353072152, 263747951750360, 3814986502092304, 55534064877048198, 812944042149730764, 11959798385860453492, 176733862787006701400
Offset: 0

Views

Author

Keywords

Comments

With interpolated zeros, this is C(n)*(1+(-1)^n)/2 with g.f. given by 2/(sqrt(1+4x) + sqrt(1-4x)). - Paul Barry, Sep 09 2004
Self-convolution of a(n)/4^n gives Catalan numbers (A000108). - Vladimir Reshetnikov, Oct 10 2016
a(n) is the number of grand Dyck paths from (0,0) to (4n,0) that avoid vertices (2k,0) for all odd k > 0. - Alexander Burstein, May 11 2021
a(n) is the number of lattice paths from (0,0) to (2n,2n) with steps (1,0) and (0,1) that avoid the points (1,1), (3,3), (5,5), ..., (2n-1,2n-1). This is Example 2.5 of the Shapiro reference. - Lucas A. Brown, Jul 24 2025

Examples

			sqrt(2*x^-1*(1-sqrt(1-x))) = 1 + (1/8)*x + (7/128)*x^2 + (33/1024)*x^3 + ...
		

Crossrefs

Programs

  • Mathematica
    a[n_] := CatalanNumber[2n]; Array[a, 18, 0] (* Or *)
    CoefficientList[ Series[ Sqrt[2]/Sqrt[1 + Sqrt[1 - 16 x]], {x, 0, 17}], x] (* Robert G. Wilson v *)
    CatalanNumber[Range[0,40,2]] (* Harvey P. Dale, Mar 19 2015 *)
  • MuPAD
    combinat::dyckWords::count(2*n) $ n = 0..28 // Zerinvary Lajos, Apr 14 2007
    
  • PARI
    /* G.f.: A(x) = exp( x*A(x)^4 + Integral(A(x)^4 dx) ): */
    {a(n)=local(A=1+x); for(i=1, n, A=exp(x*A^4 + intformal(A^4 +x*O(x^n)))); polcoeff(A, n)} \\ Paul D. Hanna, Nov 09 2013
    for(n=0, 30, print1(a(n), ", "))
    
  • Sage
    A048990 = lambda n: hypergeometric([1-2*n,-2*n],[2],1)
    [Integer(A048990(n).n()) for n in range(20)] # Peter Luschny, Sep 22 2014

Formula

a(n) = 2 * A065097(n) - A000007(n).
G.f.: A(x) = sqrt((1/8)*x^(-1)*(1-sqrt(1-16*x))).
G.f.: 2F1( (1/4, 3/4); (3/2))(16*x). - Olivier Gérard Feb 17 2011
D-finite with recurrence n*(2*n+1)*a(n) - 2*(4*n-1)*(4*n-3)*a(n-1) = 0. - R. J. Mathar, Nov 30 2012
E.g.f: 2F2(1/4, 3/4; 1, 3/2; 16*x). - Vladimir Reshetnikov, Apr 24 2013
G.f. A(x) satisfies: A(x) = exp( x*A(x)^4 + Integral(A(x)^4 dx) ). - Paul D. Hanna, Nov 09 2013
G.f. A(x) satisfies: A(x) = sqrt(1 + 4*x*A(x)^4). - Paul D. Hanna, Nov 09 2013
a(n) = hypergeom([1-2*n,-2*n],[2],1). - Peter Luschny, Sep 22 2014
a(n) ~ 2^(4*n-3/2)/(sqrt(Pi)*n^(3/2)). - Ilya Gutkovskiy, Oct 10 2016
From Peter Bala, Feb 27 2020: (Start)
a(n) = (4^n)*binomial(2*n + 1/2, n)/(4*n + 1).
O.g.f.: A(x) = sqrt(c(4*x)), where c(x) = (1 - sqrt(1 - 4*x))/(2*x) is the o.g.f. of the Catalan numbers. Cf. A228411. (End)
Sum_{n>=0} 1/a(n) = A276483. - Amiram Eldar, Nov 18 2020
Sum_{n>=0} a(n)/4^n = sqrt(2). - Amiram Eldar, Mar 16 2022
From Peter Bala, Feb 22 2023: (Start)
a(n) = (1/2^(2*n-1)) * Product_{1 <= i <= j <= 2*n-1} (i + j + 2)/(i + j - 1) for n >= 1.
a(n) = Product_{1 <= i <= j <= 2*n-1} (3*i + j + 2)/(3*i + j - 1). Cf. A024492. (End)
a(n) = Sum_{k = 0..2*n-1} (-1)^k * 4^(2*n-k-1)*binomial(2*n-1, k)*Catalan(k+1). - Peter Bala, Apr 29 2024
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^6). - Seiichi Manyama, Jun 20 2025

A091527 a(n) = ((3*n)!/n!^2)*(Gamma(1+n/2)/Gamma(1+3n/2)).

Original entry on oeis.org

1, 4, 30, 256, 2310, 21504, 204204, 1966080, 19122246, 187432960, 1848483780, 18320719872, 182327718300, 1820797698048, 18236779032600, 183120225632256, 1842826521244230, 18581317012684800, 187679234340049620, 1898554215471513600, 19232182592635611060
Offset: 0

Views

Author

Michael Somos, Jan 18 2004

Keywords

Comments

Sequence terms are given by [x^n] ( (1 + x)^(k+2)/(1 - x)^k )^n for k = 1. See the crossreferences for related sequences obtained from other values of k. - Peter Bala, Sep 29 2015
Let a > b be nonnegative integers. Then the ratio of factorials ((2*a + 1)*n)!*((b + 1/2)*n)!/(((a + 1/2)*n)!*((2*b + 1)*n)!*((a - b)*n)!) is an integer for n >= 0. This is the case a = 1, b = 0. - Peter Bala, Aug 28 2016

References

  • R. P. Stanley, Enumerative Combinatorics Volume 2, Cambridge Univ. Press, 1999, Theorem 6.33, p. 197.

Crossrefs

Cf. A061162(n) = a(2n), A007297, A000984 (k = 0), A001448 (k = 2), A262732 (k = 3), A211419 (k = 4), A262733 (k = 5), A211421 (k = 6), A276098, A276099.

Programs

  • Maple
    a := n -> 4^n * `if`(n<2, 1, (2*(n+1)*binomial((3*n-1)/2, n + 1))/(n-1)):
    seq(a(n), n=0..18); # Peter Luschny, Feb 03 2020
  • Mathematica
    Table[((3 n)!/n!^2) Gamma[1 + n/2]/Gamma[1 + 3 n/2], {n, 0, 18}] (* Michael De Vlieger, Oct 02 2015 *)
    Table[4^n Sum[Binomial[k - 1 + (n - 1)/2, k], {k, 0, n}], {n, 0, 18}] (* Michael De Vlieger, Aug 28 2016 *)
  • Maxima
    B(x):=(-1/3+(2/3)*sqrt(1+9*x)*sin((1/3)*asin((2+27*x+54*x^2)/2/(1+9*x)^(3/2))))/x-1;
    taylor(x*diff(B(x),x)/B(x),x,0,10); /* Vladimir Kruchinin, Oct 02 2015 */
    
  • PARI
    a(n)=4^n*sum(i=0,n,binomial(i-1+(n-1)/2,i))
    
  • PARI
    vector(30, n, sum(k=0, n, binomial(3*n-3, k)*binomial(2*n-k-3, n-k-1))) \\ Altug Alkan, Oct 04 2015
    
  • Python
    from math import factorial
    from sympy import factorial2
    def A091527(n): return int((factorial(3*n)*factorial2(n)<Chai Wah Wu, Aug 10 2023

Formula

D-finite with recurrence n*(n - 1)*a(n) = 12*(3*n - 1)*(3*n - 5)*a(n-2).
From Peter Bala, Sep 29 2015: (Start)
a(n) = Sum_{i = 0..n} binomial(3*n,i) * binomial(2*n-i-1,n-i).
a(n) = [x^n] ( (1 + x)^3/(1 - x) )^n.
exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 4*x + 23*x^2 + 156*x^3 + 1162*x^4 + 9192*x^5 + ... is the o.g.f. for A007297 (but with an offset of 0). (End)
a(n) = (n+1)*A078531(n). [Barry, JIS (2011)]
G.f.: x*B'(x)/B(x), where x*B(x)+1 is g.f. of A007297. - Vladimir Kruchinin, Oct 02 2015
From Peter Bala, Aug 22 2016: (Start)
a(n) = Sum_{k = 0..floor(n/2)} binomial(4*n,n-2*k)*binomial(n+k-1,k).
O.g.f.: A(x) = Hypergeom([5/6, 1/6], [1/2], 108*x^2) + 4*x*Hypergeom([4/3, 2/3], [3/2], 108*x^2).
The o.g.f. is the diagonal of the bivariate rational function 1/(1 - t*(1 + x)^3/(1 - x)) and hence is algebraic by Stanley 1999, Theorem 6.33, p. 197. (End)
a(n) ~ 2^n*3^(3*n/2)/sqrt(2*Pi*n). - Ilya Gutkovskiy, Aug 22 2016
a(n) = 4^n*2*(n+1)*binomial((3*n-1)/2, n+1)/(n-1) for n >= 2. - Peter Luschny, Feb 03 2020
From Peter Bala, Mar 04 2022: (Start)
The o.g.f. A(x) satisfies the algebraic equation (1 - 108*x^2)*A(x)^3 - A(x) = 8*x. Cf. A244039.
The Gauss congruences a(n*p^k) == a(n*p^(k-1)) (mod p^k) hold for all primes p and positive integers n and k.
Conjecture: the stronger supercongruences a(n*p^k) == a(n*p^(k-1)) (mod p^(3*k)) hold for primes p >= 5 and positive integers n and k. (End)
From Seiichi Manyama, Aug 09 2025: (Start)
a(n) = [x^n] 1/((1-x)^(n+1) * (1-2*x)^n).
a(n) = Sum_{k=0..n} 2^k * (-1)^(n-k) * binomial(3*n,k) * binomial(2*n-k,n-k).
a(n) = Sum_{k=0..n} 2^k * binomial(n+k-1,k) * binomial(2*n-k,n-k).
a(n) = 4^n * binomial((3*n-1)/2,n).
a(n) = [x^n] 1/(1-4*x)^((n+1)/2).
a(n) = [x^n] (1+4*x)^((3*n-1)/2). (End)

A078532 Coefficients of power series that satisfies A(x)^3 - 9*x*A(x)^4 = 1, A(0)=1.

Original entry on oeis.org

1, 3, 27, 315, 4158, 59049, 880308, 13586859, 215233605, 3479417370, 57168561996, 951892141473, 16026585711660, 272383068872700, 4666865660812044, 80521573261807755, 1397858693681272230, 24398716826612190447, 427921056863230599900, 7537621933880388620010
Offset: 0

Views

Author

Paul D. Hanna, Nov 28 2002

Keywords

Comments

If A(x) = Sum_{k>=1} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2) (conjecture).
If A(x) = Sum_{k>=1} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(k)=n^(2k)*binomial(k/n+1/n+k-1,k)/(k+1) and, consequently, a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2). - Emeric Deutsch, Dec 10 2002
A generalization of the Catalan sequence (A000108) since for n = 1 the equation A(x)^n -(n^2)*x*A(x)^(n+1) = 1 reduces to A(x)=1+xA(x)^2. - Emeric Deutsch, Dec 10 2002
Radius of convergence of g.f. A(x) is r = 1/(3*4^(4/3)) where A(r) = 4^(1/3). - Paul D. Hanna, Jul 24 2012
Self-convolution cube yields A214668.

Examples

			A(x)^3 - 9x*A(x)^4 = 1 since A(x)^3 = 1 +9x +108x^2 +1458x^3 +21060x^4 +... and A(x)^4 = 1 +12x +162x^2 +2340x^3 +... also a(2)=3^3, a(5)=3^10.
		

Crossrefs

Programs

  • Mathematica
    Table[3^(2n) Binomial[(4n-2)/3,n]/(n+1),{n,0,20}] (* Harvey P. Dale, Nov 03 2011 *)
  • PARI
    for(n=0,25, print1(9^n * binomial((4*n-2)/3, n)/(n+1), ", ")) \\ G. C. Greubel, Jan 26 2017

Formula

a(n) = 3^(2n)*binomial(4n/3-2/3, n)/(n+1). - Emeric Deutsch, Dec 10 2002
Sequence with offset 1 is expansion of reversion of g.f. x*(1-9*x)^(1/3), which equals x times the g.f. of A004990.
a(n) ~ 2^(8*n/3-5/6) * 3^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Dec 03 2014
D-finite with recurrence n*(n-1)*(n+1)*a(n) -216*(4*n-5)*(2*n-1)*(4*n-11)*a(n-3)=0. - R. J. Mathar, Mar 24 2023
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^5). - Seiichi Manyama, Jun 20 2025

Extensions

More terms from Harvey P. Dale, Nov 03 2011

A078534 Coefficients of power series that satisfies A(x)^5 - 25*x*A(x)^6 = 1, A(0)=1.

Original entry on oeis.org

1, 5, 100, 2625, 78125, 2502500, 84150000, 2929265625, 104646953125, 3814697265625, 141323284375000, 5305403695312500, 201382633183593750, 7715985752343750000, 298023223876953125000, 11591412585295166015625, 453601640704152832031250
Offset: 0

Views

Author

Paul D. Hanna, Nov 28 2002

Keywords

Comments

If A(x) = Sum_{k>=1} a(k)*x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(n-1) = n^(2*n-3) and a(2*n-1) = n^(4*n-2) (conjecture).
From Emeric Deutsch, Dec 10 2002: (Start)
If A(x) = Sum_{k>=1} a(k)*x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(k) = n^(2*k)*binomial(k/n+1/n+k-1,k)/(k+1) and, consequently, a(n-1) = n^(2*n-3) and a(2*n-1) = n^(4*n-2).
A generalization of the Catalan sequence (A000108) since for n = 1 the equation A(x)^n - (n^2)*x*A(x)^(n+1) = 1 reduces to A(x) = 1 + x*A(x)^2. (End)

Examples

			A(x)^5 - 25*x*A(x)^6 = 1 since A(x)^5 = 1 + 25*x + 750*x^2 + 24375*x^3 + 831250*x^4 + ... and A(x)^6 = 1 + 30*x + 975*x^2 + 33250*x^3 + ... also a(4) = 5^7, a(9) = 5^18 = 3814697265625.
		

Crossrefs

Programs

  • Mathematica
    Table[5^(2n) Binomial[(6n-4)/5,n]/(n+1),{n,0,25}]  (* Harvey P. Dale, Mar 27 2011 *)
  • PARI
    for(n=0,50, print1(5^(2*n)*binomial((6*n-4)/5, n)/(n+1), ", ")) \\ G. C. Greubel, Jan 30 2017

Formula

a(n) = 5^(2*n)*binomial(6*n/5 - 4/5, n)/(n+1). - Emeric Deutsch, Dec 10 2002
a(n) ~ sqrt(3) * 6^(6*n/5 - 4/5) * 5^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Dec 03 2014
From Seiichi Manyama, Jun 21 2025: (Start)
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^7).
G.f.: ( (1/x) * Series_Reversion(x/(1+25*x)^(6/5)) )^(1/6). (End)

Extensions

More terms from Harvey P. Dale, Mar 27 2011

A182122 Expansion of exp( arcsinh( 2*x ) ).

Original entry on oeis.org

1, 2, 2, 0, -2, 0, 4, 0, -10, 0, 28, 0, -84, 0, 264, 0, -858, 0, 2860, 0, -9724, 0, 33592, 0, -117572, 0, 416024, 0, -1485800, 0, 5348880, 0, -19389690, 0, 70715340, 0, -259289580, 0, 955277400, 0, -3534526380, 0, 13128240840, 0, -48932534040, 0, 182965127280
Offset: 0

Views

Author

Michael Somos, Apr 13 2012

Keywords

Examples

			G.f. = 1 + 2*x + 2*x^2 - 2*x^4 + 4*x^6 - 10*x^8 + 28*x^10 - 84*x^12 + ...
		

Crossrefs

Programs

  • Magma
    m:=50; R:=PowerSeriesRing(Rationals(), m); Coefficients(R!(Exp(Argsinh(2*x)))); // G. C. Greubel, Aug 12 2018
  • Maple
    s := proc(n) option remember; `if`(n<2, n+1, -4*(n-2)*s(n-2)/(n+1)) end: A127846 := n -> `if`(n<2,n+1,s(n-1)); seq(A127846(n), n=0..47); # Peter Luschny, Sep 23 2014
  • Mathematica
    CoefficientList[Series[Exp[ArcSinh[2x]],{x,0,50}],x] (* Harvey P. Dale, Aug 18 2012 *)
    Table[2 HypergeometricPFQ[{-n+1,2-n},{2},-1],{n,0,46}] (* Peter Luschny, Sep 23 2014 *)
  • PARI
    {a(n) = if( n<2, (n>=0) + (n>0), n = n-2; if( n%2, 0, (-1)^(n/2) * 4 * binomial( n, n/2) / (n + 2)))};
    
  • PARI
    {a(n) = if( n<0, 0, polcoeff( sqrt( 1 + 4*x^2 + x*O(x^n) ) + 2*x, n ) )};
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = 1 + O(x); for( k=1, n, A = sqrt( 1 + 4*x * A)); polcoeff( A, n))};
    
  • Sage
    def A182122(n):
        if n < 2: return n+1
        if n % 2 == 1: return 0
        return (-1)^(n/2-1)*binomial(n,n/2)/(n-1)
    [A182122(n) for n in range(47)] # Peter Luschny, Sep 23 2014
    

Formula

G.f.: 2*x + sqrt( 1 + 4*x^2 ) = 1 / (1 - 2*x / (1 + x / (1 - x / (1 + x / ... )))).
The g.f. A(x) satisfies: A(x) = sqrt(1 + 4*x * A(x)).
a(n) = (-1)^n * A104624(n). Convolution inverse of A104624.
Conjecture : n*(n+1)*a(n) + (n+2)*(n-1)*a(n-1) +4*(n+1)*(n-3)*a(n-2) +4*(n+2)*(n-4)*a(n-3) = 0.- R. J. Mathar, Jul 24 2012
a(n) = 2*hypergeom([-n+1,2-n],[2],-1). - Peter Luschny, Sep 23 2014
0 = a(n)*(+16*a(n+2) + 10*a(n+4)) + a(n+2)*(-2*a(n+2) + a(n+4)) if n>=0. - Michael Somos, Jan 10 2017
a(n+4) = 2 * a(n+2) * (a(n+2) - 8*a(n)) / (a(n+2) + 10*a(n)) if n>=0 is even. - Michael Somos, Jan 10 2017
G.f. A(x) satisfies A(x) = 1/A(-x). - Seiichi Manyama, Jun 20 2025

A214377 G.f. A(x) satisfies A(x) = 1 + 4*x*A(x)^(3/2).

Original entry on oeis.org

1, 4, 24, 168, 1280, 10296, 86016, 739024, 6488064, 57946200, 524812288, 4808643120, 44493176832, 415146189360, 3901709352960, 36902658748320, 350980432461824, 3354743017001880, 32207616155320320, 310446853795570800, 3003167577200394240, 29146910264615460240
Offset: 0

Views

Author

Paul D. Hanna, Jul 14 2012

Keywords

Comments

Radius of convergence of g.f. A(x) is r = sqrt(3)/18 where A(r) = 3.

Examples

			G.f.: A(x) = 1 + 4*x + 24*x^2 + 168*x^3 + 1280*x^4 + 10296*x^5 + 86016*x^6 + ... where A(x) = 1 + 4*x*A(x)^(3/2).
Radius of convergence: r = 1/(2*3^(3/2)) = 0.09622504486...
Related expansions:
A(x)^(3/2) = 1 + 6*x + 42*x^2 + 320*x^3 + 2574*x^4 + 21504*x^5 + 184756*x^6 + ...
A(x)^(1/2) = 1 + 2*x + 10*x^2 + 64*x^3 + 462*x^4 + 3584*x^5 + 29172*x^6 + ... + A078531(n)*x^n + ...
		

References

  • Bruce C. Berndt, Ramanujan's Notebooks Part I, Springer Verlag, 1985, p. 305.

Crossrefs

Programs

  • Mathematica
    Table[4^n*Binomial[3*n/2, n]*2/(n+2), {n, 0, 20}] (* Vaclav Kotesovec, Oct 20 2015 *)
  • PARI
    {a(n)=4^n*binomial(3/2*n,n)/(n/2+1)}
    for(n=0,30,print1(a(n),", "))

Formula

a(n) = 2^(2*n+1) * binomial(3*n/2, n) / (n+2).
Self-convolution of A078531.
A(-x) = 1/x * series reversion( x*(2*x + sqrt(1 + 4*x^2))^2 ) follows from the Lagrange inversion formula and equation 1.13, p. 305 in Berndt. Cf. A098616. - Peter Bala, Oct 19 2015
a(n) ~ 2^(n + 1/2) * 3^(3*n/2 + 1/2) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Oct 20 2015
G.f.: 4*x*(sin(asin(216*x^2-1)/3)/(6*x)+1/(12*x))^3+1. - Vladimir Kruchinin, Sep 30 2022
From Paul D. Hanna, Feb 03 2023: (Start)
G.f. A(x) satisfies: A(x) = 1/(1 - 4*x*A(x)^(1/2)).
G.f. A(x) satisfies: A(x) = B(x*A(x)) where B(x) = A(x/B(x)) = 1 + 8*x^2 + 4*x*sqrt(1 + 4*x^2) is the g.f. of A135863. (End)
From Karol A. Penson, Mar 23 2024: (Start)
G.f. = 1/(48*z^2) - 2F1([-2/3, -1/3], [-1/2], 108*z^2)/(48*z^2) + 4*z*2F1([5/6, 7/6],[5/2],108*z^2); a(n) = Integral_{x=0..sqrt(108)} x^n*W(x), with W(x) = ((72*(g1(x) - g2(x)) + x^2*(-g1(x) + g2(x)) + 4*sqrt(-3*x^2 + 324)*(g1(x) + g2(x)))*3^(1/6))/(96*Pi*(x^2)^(5/6)),
where g1(x) = (18 - sqrt(324 - 3*x^2))^(2/3) and
g2(x) = (18 + sqrt(324 - 3*x^2))^(2/3).
This integral representation is unique as W(x) is the solution of the Hausdorff power moment problem on x = (0, sqrt(108)). Using only the definition of a(n), W(x) can be proven to be positive. W(x) is singular at x = 0, with singularity x^(-1/3), and for x > 0 is monotonically decreasing to zero at x = sqrt(108). For x -> sqrt(108), W'(x) tends to -infinity. (End)
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^2). - Seiichi Manyama, Jun 17 2025
D-finite with recurrence -(n+2)*(n-1)*a(n) +12*(3*n-2)*(3*n-4)*a(n-2)=0. - R. J. Mathar, Jul 30 2025

A078533 Coefficients of power series that satisfies A(x)^4 - 16x*A(x)^5 = 1, A(0)=1.

Original entry on oeis.org

1, 4, 56, 1024, 21216, 473088, 11075328, 268435456, 6677665280, 169514369024, 4373549027328, 114349209288704, 3023068543631360, 80675644291153920, 2170389180446539776, 58798996734949195776, 1602737048880933109760, 43924199383151211970560
Offset: 0

Views

Author

Paul D. Hanna, Nov 28 2002

Keywords

Comments

If A(x) = Sum_{k>=1} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2) (conjecture).
If A(x) = Sum_{k>=1} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(k) = n^(2k)*binomial(k/n + 1/n + k - 1, k)/(k+1) and, consequently, a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2). - Emeric Deutsch, Dec 10 2002
A generalization of the Catalan sequence (A000108) since for n = 1 the equation A(x)^n - (n^2)*x*A(x)^(n+1) = 1 reduces to A(x) = 1 + xA(x)^2. - Emeric Deutsch, Dec 10 2002

Examples

			A(x)^4 - 16x*A(x)^5 = 1 since A(x)^4 = 1 + 16x + 320x^2 + 7040x^3 + 163840x^4 + ... and A(x)^5 = 1 + 20x + 440x^2 + 10240x^3 + ... also a(3) = 4^5, a(7) = 4^14 = 268435456.
		

Crossrefs

Programs

  • Mathematica
    Table[4^(2*n)*Binomial[5*n/4-3/4, n]/(n+1),{n,0,20}] (* Vaclav Kotesovec, Dec 03 2014 *)
  • PARI
    for(n=0,50, print1(2^(4*n)*binomial((5*n-3)/4,n)/(n+1), ", ")) \\ G. C. Greubel, Jan 30 2017

Formula

a(n) = 4^(2n)*binomial(5n/4 - 3/4, n)/(n+1). - Emeric Deutsch, Dec 10 2002
a(n) ~ 5^(5*n/4 - 1/4) * 2^(2*n - 1/2) / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Dec 03 2014
From Seiichi Manyama, Jun 21 2025: (Start)
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^6).
G.f.: ( (1/x) * Series_Reversion(x/(1+16*x)^(5/4)) )^(1/5). (End)

A078535 Coefficients of power series that satisfies A(x)^6 - 36x*A(x)^7 = 1, A(0)=1.

Original entry on oeis.org

1, 6, 162, 5760, 232254, 10077696, 458960580, 21634449408, 1046465787510, 51644846702592, 2590092194793948, 131621703842267136, 6762649550214036780, 350714987252652441600, 18334388441036020419720, 965148007553698721955840, 51116742846877582931249574
Offset: 0

Views

Author

Paul D. Hanna, Nov 28 2002

Keywords

Comments

If A(x)=sum_{k=1..inf} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2) (conjecture).
If A(x)=sum_{k=1..inf} a(k)x^k satisfies A(x)^n - (n^2)*x*A(x)^(n+1) = 1, then a(k)=n^(2k)*binomial(k/n+1/n+k-1,k)/(k+1) and, consequently, a(n-1) = n^(2n-3) and a(2n-1) = n^(4n-2). - Emeric Deutsch, Dec 10 2002
A generalization of the Catalan sequence (A000108) since for n = 1 the equation A(x)^n -(n^2)*x*A(x)^(n+1) = 1 reduces to A(x)=1+xA(x)^2. - Emeric Deutsch, Dec 10 2002

Examples

			A(x)^6 - 36x*A(x)^7 = 1 since A(x)^6 = 1 +36x +1512x^2 +68040x^3 +3193344x^4 +... and A(x)^7 = 1 +42x +1890x^2 +88704x^3 +... also a(5)=6^9, a(11)=6^22 = 131621703842267136.
		

Crossrefs

Programs

  • Mathematica
    Table[6^(2*n)*Binomial[7*n/6-5/6, n]/(n+1),{n,0,20}] (* Vaclav Kotesovec, Dec 03 2014 *)
  • PARI
    a(n) = {6^(2*n)*binomial((7*n-5)/6, n)/(n+1)} \\ Andrew Howroyd, Nov 05 2019

Formula

a(n) = 6^(2n)*binomial(7n/6-5/6, n)/(n+1). - Emeric Deutsch, Dec 10 2002
a(n) ~ 7^(7*n/6-1/3) * 6^n / (sqrt(2*Pi) * n^(3/2)). - Vaclav Kotesovec, Dec 03 2014
From Seiichi Manyama, Jun 21 2025: (Start)
G.f. A(x) satisfies A(x) = 1/A(-x*A(x)^8).
G.f.: ( (1/x) * Series_Reversion(x/(1+36*x)^(7/6)) )^(1/7). (End)

Extensions

Terms a(13) and beyond from Andrew Howroyd, Nov 05 2019

A138020 G.f. satisfies A(x) = sqrt( (1 + 2*x*A(x)) / (1 - 2*x*A(x)) ).

Original entry on oeis.org

1, 2, 6, 24, 110, 544, 2828, 15232, 84246, 475648, 2730068, 15882240, 93438540, 554967040, 3323125528, 20039827456, 121597985254, 741871845376, 4548193111428, 28004975116288, 173113004348580, 1073893324357632, 6683288759506856, 41715337804120064
Offset: 0

Views

Author

Paul D. Hanna, Feb 28 2008

Keywords

Crossrefs

Programs

  • Maple
    A138020 := proc(n)
        option remember ;
        if n < 5 then
            op(n+1,[1,2,6,24,110]) ;
        else
            4*(-55*n^3 +231*n^2 -263*n +51)*procname(n-2) -16*(n-3)*(n-4)*(5*n-1)*procname(n-4) ;
            -%/n/(n+1)/(5*n-11)
        end if;
    end proc:
    seq(A138020(n),n=0..30) ; # R. J. Mathar, Sep 27 2024
  • Mathematica
    CoefficientList[y/.AsymptoticSolve[y^2-1-2x(y+y^3) ==0,y->1,{x,0,23}][[1]],x] (* Alexander Burstein, Nov 26 2021 *)
  • PARI
    a(n)=polcoeff((1/x)*serreverse(x*sqrt((1-2*x)/(1+2*x+x^2*O(x^n)))),n)

Formula

a(n) ~ 2^(n - 1/2) * phi^((5*n + 3)/2) / (sqrt(Pi) * 5^(1/4) * n^(3/2)), where phi = A001622 = (1+sqrt(5))/2 is the golden ratio. - Vaclav Kotesovec, Oct 04 2020
From Alexander Burstein, Nov 26 2021: (Start)
G.f.: A(x) = 1 + 2*x*A(x)*(1 + A(x)^2)/(1 + A(x)).
G.f.: A(-x*A(x)^2) = 1/A(x). (End)
D-finite with recurrence +n*(n+1)*(5*n-11) *a(n) +4*(-55*n^3 +231*n^2 -263*n +51)*a(n-2) -16*(n-3)*(n-4)*(5*n-1)*a(n-4)=0. - R. J. Mathar, Mar 25 2024
From Seiichi Manyama, Dec 22 2024: (Start)
a(n) = (2^n/(n+1)) * Sum_{k=0..n} binomial(n/2+k-1/2,k) * binomial(n/2+1/2,n-k).
a(n) = 2^n * Sum_{k=0..n} binomial(n,k) * binomial(n/2+k+1/2,n)/(n+2*k+1). (End)
Showing 1-10 of 22 results. Next