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 15 results. Next

A006480 De Bruijn's S(3,n): (3n)!/(n!)^3.

Original entry on oeis.org

1, 6, 90, 1680, 34650, 756756, 17153136, 399072960, 9465511770, 227873431500, 5550996791340, 136526995463040, 3384731762521200, 84478098072866400, 2120572665910728000, 53494979785374631680, 1355345464406015082330, 34469858696831179429500, 879619727485803060256500, 22514366432046593564460000
Offset: 0

Views

Author

Keywords

Comments

Number of paths of length 3n in an n X n X n grid from (0,0,0) to (n,n,n), using steps (0,0,1), (0,1,0), and (1,0,0).
Appears in Ramanujan's theory of elliptic functions of signature 3.
S(s,n) = Sum_{k=0..2n} (-1)^(k+n) * binomial(2n, k)^s. The formula S(3,n) = (3n)!/(n!)^3 is due to Dixon (according to W. N. Bailey 1935). - Charles R Greathouse IV, Dec 28 2011
a(n) is the number of ballot results that end in a 3-way tie when 3n voters each cast two votes for two out of three candidates vying for 2 slots on a county board; in such a tie, each of the three candidates receives 2n votes. Note there are C(3n,2n) ways to choose the voters who cast a vote for the youngest candidate. The n voters who did note vote for the youngest candidate voted for the two older candidates. Then there are C(2n,n) ways to choose the other n voters who voted for both the youngest and the second youngest candidate. The remaining voters vote for the oldest candidate. Hence there are C(3n,2n)*C(2n,n)=(3n)!/(n!)^3 ballot results. - Dennis P. Walsh, May 02 2013
a(n) is the constant term of (X+Y+1/(X*Y))^(3*n). - Mark van Hoeij, May 07 2013
For n > 2 a(n) is divisible by (n+2)*(n+1)^2, a(n) = (n+1)^2*(n+2)*A161581(n). - Alexander Adamchuk, Dec 27 2013
a(n) is the number of permutations of the multiset {1^n, 2^n, 3^n}, the number of ternary words of length 3*n with n of each letters. - Joerg Arndt, Feb 28 2016
Diagonal of the rational function 1/(1 - x - y - z). - Gheorghe Coserea, Jul 06 2016
At least two families of elliptic curves, x = 2*H1 = (p^2+q^2)*(1-q) and x = 2*H2 = p^2+q^2-3*p^2*q+q^3 (0Bradley Klee, Feb 25 2018
The ordinary generating function also determines periods along a family of tetrahedral-symmetric sphere curves ("du troisième ordre"). Compare links to Goursat "Étude des surfaces..." and "Proof Certificate". - Bradley Klee, Sep 28 2018

Examples

			G.f.: 1 + 6*x + 90*x^2 + 1680*x^3 + 34650*x^4 + 756756*x^5 + 17153136*x^6 + ...
		

References

  • L. A. Aizenberg and A. P. Yuzhakov, "Integral representations and residues in multidimensional complex analysis", American Mathematical Society, 1983, p. 194.
  • Louis Comtet, Advanced Combinatorics, Reidel, 1974, p. 174.
  • N. G. de Bruijn, Asymptotic Methods in Analysis, North-Holland Publishing Co., 1958. See chapters 4 and 6.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Row 3 of A187783.
Related to diagonal of rational functions: A268545-A268555. Elliptic Integrals: A002894, A113424, A000897. Factors: A005809, A000984. Integrals: A007004, A024486. Sphere Curves: A318245, A318495.

Programs

  • GAP
    List([0..20],n->Factorial(3*n)/Factorial(n)^3); # Muniru A Asiru, Mar 31 2018
    
  • Magma
    [Factorial(3*n)/(Factorial(n))^3: n in [0..20] ]; // Vincenzo Librandi, Aug 20 2011
    
  • Maple
    seq((3*n)!/(n!)^3, n=0..16); # Zerinvary Lajos, Jun 28 2007
  • Mathematica
    Sum [ (-1)^(k+n) Binomial[ 2n, k ]^3, {k, 0, 2n} ]
    a[ n_] := If[ n < 0, 0, (-1)^n HypergeometricPFQ[ {-2 n, -2 n, -2 n}, {1, 1}, 1]]; (* Michael Somos, Oct 22 2014 *)
    Table[Multinomial[n, n, n], {n, 0, 100}] (* Emanuele Munarini, Oct 25 2016 *)
    CoefficientList[Series[Hypergeometric2F1[1/3,2/3,1,27*x],{x,0,5}],x] (* Bradley Klee, Feb 28 2018 *)
    Table[(3n)!/(n!)^3,{n,0,20}] (* Harvey P. Dale, Mar 09 2025 *)
  • Maxima
    makelist(multinomial_coeff(n,n,n),n,0,24); /* Emanuele Munarini, Oct 25 2016 */
    
  • PARI
    {a(n) = if( n<0, 0, (3*n)! / n!^3)}; /* Michael Somos, Dec 03 2002 */
    
  • PARI
    {a(n) = my(A, m); if( n<1, n==0, m=1; A = 1 + O(x); while( m<=n, m*=3; A = subst( (1 + 2*x) * subst(A, x, (x/3)^3), x, serreverse(x * (1 + x + x^2) / (1 + 2*x)^3 / 3 + O(x^m)))); polcoeff(A, n))}; /* Michael Somos, Dec 03 2002 */
    
  • Python
    from math import factorial
    def A006480(n): return factorial(3*n)//factorial(n)**3 # Chai Wah Wu, Oct 04 2022

Formula

Using Stirling's formula in A000142 it is easy to get the asymptotic expression a(n) ~ 1/2 * sqrt(3) * 27^n / (Pi*n) - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001
From Karol A. Penson, Nov 21 2001: (Start)
O.g.f.: hypergeom([1/3, 2/3], [1], 27*x).
E.g.f.: hypergeom([1/3, 2/3], [1, 1], 27*x).
Integral representation as n-th moment of a positive function on [0, 27]:
a(n) = int( x^n*(-1/24*(3*sqrt(3)*hypergeom([2/3, 2/3], [4/3], 1/27*x)* Gamma(2/3)^6*x^(1/3) - 8*hypergeom([1/3, 1/3], [2/3], 1/27*x)*Pi^3)/Pi^3 /x^(2/3)/Gamma(2/3)^3), x=0..27). This representation is unique. (End)
a(n) = Sum_{k=-n..n} (-1)^k*binomial(2*n, n+k)^3. - Benoit Cloitre, Mar 02 2005
a(n) = C(2n,n)*C(3n,n) = A104684(2n,n). - Paul Barry, Mar 14 2006
G.f. satisfies: A(x^3) = A( x*(1+3*x+9*x^2)/(1+6*x)^3 )/(1+6*x). - Paul D. Hanna, Oct 29 2010
D-finite with recurrence: n^2*a(n) - 3*(3*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Dec 04 2012
a(n) = (n+1)^2*(n+2)*A161581(n) for n>2. - Alexander Adamchuk, Dec 27 2013
0 = a(n)^2*(472392*a(n+1)^2 - 83106*a(n+1)*a(n+2) + 3600*a(n+2)^2) + a(n)*a(n+1)*(-8748*a(n+1)^2 + 1953*a(n+1)*a(n+2) - 120*a(n+2)^2) + a(n+1)^2*(36*a(n+1)^2 - 12*a(n+1)*a(n+2) + a(n+2)^2) for all n in Z. - Michael Somos, Oct 22 2014
0 = x*(27*x-1)*y'' + (54*x-1)*y' + 6*y, where y is g.f. - Gheorghe Coserea, Jul 06 2016
From Peter Bala, Jul 15 2016: (Start)
a(n) = 3*binomial(2*n - 1,n)*binomial(3*n - 1,n) = 3*[x^n] 1/(1 - x)^n * [x^n] 1/(1 - x)^(2*n) for n >= 1.
a(n) = binomial(2*n,n)*binomial(3*n,n) = ([x^n](1 + x)^(2*n)) *([x^n](1 + x)^(3*n)) = [x^n](F(x)^(6*n)), where F(x) = 1 + x + 2*x^2 + 14*x^3 + 127*x^4 + 1364*x^5 + 16219*x^6 + ... appears to have integer coefficients. Cf. A002894.
This sequence occurs as the right-hand side of several binomial sums:
Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n,k)^3 = a(n) (Dixon's identity).
Sum_{k = 0..n} binomial(n,k)*binomial(2*n,n - k)*binomial(3*n + k,k) = a(n) (Gould, Vol. 4, 6.86)
Sum_{k = 0..n} (-1)^(n+k)*binomial(n,k)*binomial(2*n + k,n)*binomial(3*n + k,n) = a(n).
Sum_{k = 0..n} binomial(n,k)*binomial(2*n + k,k)*binomial(3*n,n - k) = a(n).
Sum_{k = 0..n} (-1)^(k)*binomial(n,k)*binomial(3*n - k,n)*binomial(4*n - k,n) = a(n).
Sum_{k = 0..2*n} (-1)^(n+k)*binomial(2*n + k,2*n - k)*binomial(2*k,k)*binomial(4*n - k,2*n) = a(n) (see Gould, Vol.5, 9.23).
Sum_{k = 0..2*n} (-1)^k*binomial(3*n,k)*binomial(3*n - k,n)^3 = a(n) (Sprugnoli, Section 2.9, Table 10, p. 123). (End)
From Bradley Klee, Feb 28 2018: (Start)
a(n) = A005809(n)*A000984(n).
G.f.: F(x) = 1/(2*Pi) Integral_{z=0..2*Pi} 2F1(1/3,2/3; 1/2; 27*x*sin^2(z)) dz.
With G(x) = x*2F1(1/3,2/3; 2; 27*x): F(x) = d/dx G(x). (Cf. A007004) (End)
F(x)*G(1/27-x) + F(1/27-x)*G(x) = 1/(4*Pi*sqrt(3)). - Bradley Klee, Sep 29 2018
Sum_{n>=0} 1/a(n) = A091683. - Amiram Eldar, Nov 15 2020
From Peter Bala, Sep 20 2021: (Start)
a(n) = Sum_{k = n..2*n} binomial(2*n,k)^2 * binomial(k,n). Cf. A001459.
a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for any prime p >= 5 and any positive integers n and k (write a(n) as C(3*n,2*n)*C(2*n,n) and apply Mestrovic, equation 39, p. 12). (End)
a(n) = 6*A060542(n). - R. J. Mathar, Jun 21 2023
Occurs on the right-hand side of the binomial sum identities Sum_{k = -n..n} (-1)^k * (n + x - k) * binomial(2*n, n+k)^3 = (x + n)*a(n) and Sum_{k = -n..n} (-1)^k * (n + x - k)^3 * binomial(2*n, n+k)^3 = x*(x + n)*(x + 2*n)*a(n) (x arbitrary). Compare with Dixon's identity: Sum_{k = -n..n} (-1)^k * binomial(2*n, n+k)^3 = a(n). - Peter Bala, Jul 31 2023
From Peter Bala, Aug 14 2023: (Start)
a(n) = (-1)^n * [x^(2*n)] ( (1 - x)^(4*n) * Legendre_P(2*n, (1 + x)/(1 - x)) ).
Row 1 of A364509. (End)
From Peter Bala, Oct 10 2024: (Start)
The following hold for n >= 1:
a(n) = Sum_{k = 0.. 2*n} (-1)^(n+k) * k/n * binomial(2*n, k)^3 = 3/2 * Sum_{k = 0.. 2*n} (-1)^(n+k) * (k/n)^2 * binomial(2*n, k)^3.
a(n) = 3/2 * Sum_{0..2*n-1} (-1)^(n+k) * k/n * binomial(2*n, k)^2*binomial(2*n-1, k).
a(n) = 3 * Sum_{0..2*n-1} (-1)^(n+k) * k/n * binomial(2*n, k)*binomial(2*n-1, k)^2. (End)
a(n) = Sum_{k = 0..n} (-1)^(n+k) * binomial(n, k) * A108625(2*n, k) (verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). - Peter Bala, Oct 15 2024

Extensions

a(14)-a(16) from Eric W. Weisstein
Terms a(17) and beyond from T. D. Noe, Jun 29 2008

A000897 a(n) = (4*n)! / ((2*n)!*n!^2).

Original entry on oeis.org

1, 12, 420, 18480, 900900, 46558512, 2498640144, 137680171200, 7735904619300, 441233078286000, 25467973278667920, 1484298740174927040, 87202550985276963600, 5157850293780050462400, 306839461354466267304000, 18344908596179023234548480
Offset: 0

Views

Author

Keywords

Comments

Appears in Ramanujan's theory of elliptic functions of signature 4.
H. A. Verrill proves that a(n) = Sum_{p + q + r = 3n} w^(p-q) * {(3n)!/(p! q! r!)}^2, with p, q, r >= 0 and w = primitive 3rd root of unity.
The family of elliptic curves "x=2*H1=p^2+q^2-(1/4)*q^4, 0sqrt(-1)*q" to H1 produces "x=2*H2=p^2-q^2-(1/4)*q^4, 0Bradley Klee, Feb 25 2018
Even-order terms in the diagonal of rational function 1/(1 - (x^2 + y^2 + z)). - Gheorghe Coserea, Aug 09 2018

Examples

			G.f.: 1 + 12*x + 420*x^2 + 18480*x^3 + 900900*x^4 + 46558512*x^5 + 2498640144*x^6 + ...
		

References

  • E. R. Hansen, A Table of Series and Products, Prentice-Hall, Englewood Cliffs, NJ, 1975, p. 96.

Crossrefs

Cf. A002897, A008977, A186420, A188662. Elliptic Integrals: A002894, A113424, A006480. Factors: A005809, A005810, A000984, A001448.

Programs

  • GAP
    a:=n->Sum([0..3*n],k->(-1)^k*Binomial(3*n,k)*Binomial(6*n-k,3*n)*
    Binomial(2*k,k));;
    A000897:=List([0..14],n->a(n)); # Muniru A Asiru, Feb 11 2018
  • Maple
    seq((4*n)!/(n!)^4/binomial(2*n,n), n=0..14); # Zerinvary Lajos, Jun 28 2007
  • Mathematica
    Table[(4n)!/((2n)! n!^2), {n, 0, 30}] (* Stefan Steinerberger, Apr 14 2006 *)
    a[ n_] := Binomial[ 4 n, 2 n] Binomial[ 2 n, n]; (* Michael Somos, Mar 24 2013 *)
    a[ n_] := SeriesCoefficient[ Hypergeometric2F1[ 1/4, 3/4, 1, 64 x], {x, 0, n}]; (* Michael Somos, Mar 24 2013 *)
    a[ n_] := If[ n < 0, 0, With[{m = 4 n}, (-1)^n m! SeriesCoefficient[ BesselI[ 0, 2 x] BesselJ[ 0, 2 x], {x, 0, m}]]]; (* Michael Somos, Aug 12 2014 *)
    a[ n_] := 64^n Pochhammer[1/4, n] Pochhammer[3/4, n] / n!^2; (* Michael Somos, Aug 12 2014 *)
  • PARI
    {a(n) = if( n<0, 0, (4*n)! / ((2*n)! * n!^2))}; /* Michael Somos, Oct 31 2005 */
    

Formula

E.g.f.: Sum_{k>=0} (-1)^k * a(k) * x^(4*k) / (4*k)! = BesselI(0, 2x) * BesselJ(0, 2x).
G.f.: F(1/4, 3/4; 1; 64*x). - Michael Somos, Oct 31 2005
a(n) = A008977(n)/A000984(n) - Zerinvary Lajos, Jun 28 2007
Sum_{k>=0} a(k) * x^(3k)/(3k)!^2 = f(x)*f(x*w)*f(x/w) where f(x) = BesselI(0, 2*sqrt(x)) and w = primitive 3rd root of unity. - Michael Somos, Jul 25 2007
In general, for (BesselI(b, 2x))*(BesselJ(b, 2x))=((x^(2*b))/((GAMMA(b+1))^2)*(1-(x^4)/(Q(0)+(x^4))); Q(k)=(k+1)*(k+b+1)*(2*k+b+1)*(2*k+b+2)-(x^4)+(x^4)*(k+1)*(k+b+1)*(2*k+b+1)*(2*k+b+2)/Q(k+1)) ; (continued fraction). - Sergei N. Gladkovskii, Nov 24 2011
D-finite with recurrence 0 = a(n)*4*(4*n + 1)*(4*n + 3) - a(n+1)*(n + 1)^2 for all n in Z. - Michael Somos, Aug 12 2014
0 = a(n)*(-4026531840*a(n+2) +2005401600*a(n+3) -103896576*a(n+4) +1251948*a(n+5)) + a(n+1)*(+41418752*a(n+2) -30435328*a(n+3) +1863228*a(n+4) -24604*a(n+5)) + a(n+2)*(-16896*a(n+2) +75608*a(n+3) -6740*a(n+4) +105*a(n+5)) for all n in Z. - Michael Somos, Aug 12 2014
From Peter Bala, Jul 12 2016: (Start)
a(n) = binomial(3*n,n)*binomial(4*n,n) = A005809(n)*A005810(n) = ( [x^n](1 + x)^(3*n) ) * ( [x^n](1 + x)^(4*n) ) = [x^n](F(x)^(12*n)), where F(x) = 1 + x + 6*x^2 + 105*x^3 + 2448*x^4 + 67043*x^5 + 2028307*x^6 + ... appears to have integer coefficients. Cf. A002894, A002897, A006480, A008977, A186420 and A188662. (End)
a(n) ~ 2^(6*n-1/2)/(Pi*n). - Ilya Gutkovskiy, Jul 12 2016
G.f.: 2*EllipticK(sqrt((sqrt(1-64*x)-1)/(2*sqrt(1-64*x))))/(Pi*(1-64*x)^(1/4)) where EllipticK is the complete elliptic integral of the first kind (in Maple's notation). - Robert Israel, Jul 12 2016
a(n) = Sum_{k = 0..3*n} (-1)^k*C(3*n,k)*C(6*n-k,3*n)*C(2*k,k). - Peter Bala, Feb 10 2018
From Bradley Klee, Feb 27 2018: (Start)
a(n) = A000984(n)*A001448(n).
G.f.: (1/(sqrt(2)*Pi))*Integral_{q=-oo..oo} 1/sqrt(q^2+(1/4)*q^4+(1-64*x)) dq.
G.f.: (1/(2*Pi))*Integral_{phi=0..2*Pi} 1/sqrt(1-64*x*sin^4(phi)) dphi. (End)
From Peter Bala, Mar 20 2022: (Start)
Right-hand side of the following identities valid for n >= 1:
Sum_{k = 0..2*n} 2*n*(2*n+k-1)!/(k!*n!^2) = (4*n)!/((2*n)!*n!^2);
(3/2)*Sum_{k = 0..n} 2*n*(3*n+k-1)!/(k!*n!*(2*n)!) = (4*n)!/((2*n)!*n!^2).
Cf. A001451. (End)
a(n) = (4^n/n!^2)*Product_{k = 0..2*n-1} (2*k + 1). - Peter Bala, Feb 26 2023
a(n) = Sum_{k = 0..n} (-1)^(n+k) * binomial(n, k) * A108625(3*n, k) (verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). - Peter Bala, Oct 15 2024

A001451 a(n) = (5*n)!/((3*n)!*n!*n!).

Original entry on oeis.org

1, 20, 1260, 100100, 8817900, 823727520, 79919739900, 7962100660800, 808906548235500, 83426304143982800, 8707404737345073760, 917663774856743842200, 97491279924241456098300, 10427604345391237790688000, 1121786259855036145008408000
Offset: 0

Views

Author

Keywords

Comments

From Peter Bala, Jul 15 2016: (Start)
This sequence occurs as the right-hand side of the binomial sum identity Sum_{k = 0..n} (-1)^(k)*binomial(n,k)*binomial(2*n + k,n)*binomial(3*n - k,n) = (-1)^m*a(m) for n = 2*m. The sum vanishes for n odd. Cf. A273628 and A273629.
Note the similar results:
Sum_{k = 0..n} (-1)^k*binomial(n,k)* binomial(2*n + k,n)*binomial(3*n + k,n) = (-1)^n*(3*n)!/n!^3 = (-1)^n*A006480(n);
Sum_{k = 0..n} (-1)^k*binomial(n,k)*binomial(2*n - k,n)*binomial(3*n - k,n) = binomial(2*n,n)^2 = A002894(n);
Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(2*n + k,n)*binomial(3*n + k,n) = Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(2*n - k,n)*binomial(3*n - k,n) = binomial(2*n,n) = A000984(n);
Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(2*n + k,n)*binomial(3*n - k,n) = Sum_{k = 0..2*n} (-1)^k*binomial(2*n,k)*binomial(2*n - k,n)*binomial(3*n + k,n) = (-1)^n*binomial(2*n,n) = (-1)^n*A000984(n). (End)
Choose three noncollinear step vectors to satisfy the zero sum, 3*v_1 + v_2 + v_3 = 0. Then a(n) is the number of loop plane walks of length 5*n which depart from and return to the origin. Equivalently, a(n) counts distinct permutations of a (5*n)-digit integer with digits 1,2,3 of multiplicity 3*n,n,n respectively. - Bradley Klee, Aug 12 2018

Examples

			G.f. = 1 + 20*x + 1260*x^2 + 100100*x^3 + 8817900*x^4 + 823727520*x^5 + ... - _Michael Somos_, Aug 12 2018
		

Crossrefs

Programs

  • GAP
    List([0..15],n->Factorial(5*n)/(Factorial(3*n)*Fact0rial(n)*Factorial(n))); # Muniru A Asiru, Aug 12 2018
  • Magma
    [Factorial(5*n)/(Factorial(3*n)*Factorial(n)*Factorial(n)): n in [0..30]]; // Vincenzo Librandi, May 22 2011
    
  • Maple
    f := n->(5*n)!/((3*n)!*n!*n!);
  • Mathematica
    Table[(5*n)!/((3*n)!*n!*n!), {n, 0, 20}] (* Vincenzo Librandi, Sep 04 2012 *)

Formula

a(n) = binomial(4*n,n)*binomial(5*n,n) = ( [x^n](1 + x)^(4*n) ) * ( [x^n](1 + x)^(5*n) ) = [x^n](F(x)^(20*n)), where F(x) = 1 + x + 12*x^2 + 390*x^3 + 16984*x^4 + 867042*x^5 + 48848541*x^6 + ... appears to have integer coefficients. For similar results see A000897, A002894, A002897, A006480, A008977, A186420 and A188662. - Peter Bala, Jul 14 2016
a(n) ~ 3^(-3*n-1/2)*5^(5*n+1/2)/(2*Pi*n). - Ilya Gutkovskiy, Jul 13 2016
G.f.: G(x) = 4F3(1/5,2/5,3/5,4/5;1/3,2/3,1;(5^5/3^3)*x). Let G^(n)(x) = d^n/dx^n G(x), and c = {120, 15000*x-6, 45000*x^2-114*x, 25000*x^3-135*x^2, 3125*x^4-27*x^3}, then Sum_{n=0..4} c_n*G^(n)(x) = 0. - Bradley Klee, Aug 12 2018
From Peter Bala, Mar 20 2022: (Start)
Right-hand side of the following identities valid for n >= 1:
Sum_{k = 0..3*n} 2*n*(2*n+k-1)!/(k!*n!^2) = (5*n)!/((3*n)!*n!^2);
Sum_{k = 0..n} 4*n*(4*n+k-1)!/(k!*n!*(3*n)!) = (5*n)!/((3*n)!*n!^2). Cf. A000897 and A113424. (End)
a(n) = Sum_{k = 0..n} (-1)^(n+k) * binomial(n, k) * A108625(4*n, k) (verified using the MulZeil procedure in Doron Zeilberger's MultiZeilberger package). - Peter Bala, Oct 12 2024

A318417 Scaled g.f. T(u) = Sum_{n>=0} a(n)*(3*u/48)^n satisfies 3*(2*u-1)*T + d/du(4*u*(2*u-1)*(u-1)*T') = 0, and a(0)=1; sequence gives a(n).

Original entry on oeis.org

1, 12, 228, 5040, 121380, 3093552, 82047504, 2240162496, 62508328740, 1773580002480, 50988042273168, 1481392081181376, 43413834762798864, 1281498837550545600, 38059165854011995200, 1136249610240102992640, 34076899109906247654180, 1026061759878805529676720
Offset: 0

Views

Author

Bradley Klee, Aug 26 2018

Keywords

Comments

The defining differential equation determines the period function T(u) along a square-symmetric Hamiltonian family of algebraic plane curves, u=2*H=x^2+y^2-(1/2)*(x^4+y^4). The separatrix with u=1/2 is a product of two congruent ellipses. Transformation u->1-u leaves T(u) invariant.

Examples

			Singular Points of u=x^2+y^2-(1/2)*(x^4+y^4).
   u             (x,y)              type
=============================================
   0             (0,0)            circular
  1/2      (0,+/-1) (+/-1,0)     hyperbolic
   1       +/-(1,1) +/-(1,-1)     circular
G.f. = 1 + 12*x + 228*x^2 + 5040*x^3 + 121380*x^4 + 3093552*x^5 + 82047504*x^6 + ... - _Michael Somos_, Aug 27 2018
		

Crossrefs

Factors: A000984, A098410. Quartic Periods: A113424, A002894, A318245.

Programs

  • Mathematica
    RecurrenceTable[{n^2*a[n]-12*(2*n-1)^2*a[n-1] + 128*(2*n-3)*(2*n-1)*a[n-2] == 0, a[0] == 1, a[1] == 12}, a, {n,0,1000}]
    PeriodIntegral[n_]:= CoefficientList[Series[1/Sqrt[1-2*x*((z+y)^4+(z-y)^4)], {x,0,n}], {x, y, z}][[#+1,2*#+1,2*#+1]]&/@Range[0,n];
    PeriodIntegral[10]
    HadamardProduct[n_]:=Times@@Map[CoefficientList[Normal[Series[#,{x,0,n}]],x,n+1]&, {1/Sqrt[1-12*x+32*x^2], 1/Sqrt[1-4*x]}];
    HadamardProduct[10]
    a[ n_] := Binomial[2 n, n] SeriesCoefficient[ (1 - 12 x + 32 x^2)^(-1/2), {x, 0, n}]; (* Michael Somos, Aug 27 2018 *)
  • PARI
    {a(n) = if( n<0, 0, binomial(2*n, n) * polcoeff( (1 - 12*x + 32*x^2 + x * O(x^n))^(-1/2), n))}; /* Michael Somos, Aug 27 2018 */

Formula

Define t(u,z) = 1/sqrt(1-2*u*(cos(z)^4+sin(z)^4)), then certify that:
3*(2*u-1)*t + d/du(4*u*(2*u-1)*(u-1)*dt/du) = d/dz((1/4)*sin(4*z)*t^3).
G.f.: G(x) = T(48*x/3), T(u) = 1/(2*Pi)*Integral_{z=0..2*Pi} t(u)*dz.
n^2*a(n) - 12*(2*n-1)^2*a(n-1) + 128*(2*n-3)*(2*n-1)*a(n-2) = 0.
a(n) = A000984(n)*A098410(n).
0 = +a(n)*(+a(n+1)*(+134217728*a(n+2) -25165824*a(n+3) +819200*a(n+4)) +a(n+2)*(-6291456*a(n+2) +1998848*a(n+3) -78336*a(n+4)) +a(n+3)*(-29184*a(n+3) +1472*a(n+4))) +a(n+1)*(+a(n+1)*(-6291456*a(n+2) +1179648*a(n+3) -38400*a(n+4)) +a(n+2)*(+327680*a(n+2) -110592*a(n+3) +4448*a(n+4)) +a(n+3)*(+1728*a(n+3) -90*a(n+4))) +a(n+2)*(+a(n+2)*(-1536*a(n+2) +800*a(n+3) -36*a(n+4)) +a(n+3)*(-18*a(n+3) +a(n+4))) for all n in Z. - Michael Somos, Aug 27 2018
G.f.: hypergeom([1/2, 1/2],[1],16*x/(32*x-1))/sqrt(1-32*x). - Mark van Hoeij, Dec 13 2024

A178819 Pascal's prism (3-dimensional array) read by folded antidiagonal cross-sections: (h+i; h, i-j, j), h >= 0, i >= 0, 0 <= j <= i.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 1, 2, 1, 1, 3, 3, 1, 3, 6, 3, 3, 3, 1, 1, 4, 4, 6, 12, 6, 4, 12, 12, 4, 1, 4, 6, 4, 1, 1, 5, 10, 10, 5, 1, 5, 20, 30, 20, 5, 10, 30, 30, 10, 10, 20, 10, 5, 5, 1, 1, 6, 6, 15, 30, 15, 20, 60, 60, 20, 15, 60, 90, 60, 15, 6, 30, 60, 60, 30, 6, 1, 6, 15, 20, 15, 6, 1
Offset: 0

Views

Author

Harlan J. Brothers, Jun 16 2010

Keywords

Comments

P_h = level h of Pascal's prism where P_1 = Pascal's triangle (A007318) and P_2 = denominators of Leibniz harmonic triangle (A003506). A sequence of length k through P is defined by P for n = {1, 2, 3, ..., k}.

Examples

			Prism begins (levels 1-4):
1
1 1
1 2 1
1 3 3 1
1
2 2
3 6 3
4 12 12 4
1
3 3
6 12 6
10 30 30 10
1
4 4
10 20 10
20 60 60 20
		

Crossrefs

Level 1 = A007318.
Level 2 = A003506.
Level 3 = A094305.
Level 4 = A178820.
Level 5 = A178821.
Level 6 = A178822.
Sums of shallow diagonals for each level correspond to rows of square A037027.
Contains A109649 and A046816.
P = A000984.
P = A006480.
P = A000897.
P<3n-2, 3n-2, n> = A113424.

Programs

  • Mathematica
    end = 5; Column/@Table[Multinomial[h, i-j, j], {h, 0, end}, {i, 0, end}, {j, 0, i}]

Formula

a_(h, i, j) = (h+i-2; h-1, i-j, j-1), h >= 1, i >= 1, 1 <= j <= i.
Recurrence:
For P_h, element a is given by: a_(1, 1) = 1; a_(i, j) = ((i+h-2)/(i-1)) (a_(i-1, j) + a_(i-1, j-1)).

Extensions

Keyword tabf by Michel Marcus, Oct 22 2017

A295870 a(n) = binomial(3n,n)*CQC(n), where CQC(n) = A005721(n) = A005190(2n) is a central quadrinomial coefficient.

Original entry on oeis.org

1, 12, 660, 48720, 4005540, 349260912, 31626298704, 2940502593600, 278788387440420, 26831860080682800, 2613367831568654160, 257012469788428710720, 25479526081439438845200, 2543092744417831625342400, 255292245777771431285140800, 25755871314484468746363582720
Offset: 0

Views

Author

Bradley Klee, Feb 23 2018

Keywords

Comments

Compare with EllipticK A002894 and the Ramanujan period-energy functions A113424, A006480, A000897. The series expansion "T(x) = 2*Pi*Sum_{n>=0} a_n*x^n" determines the real period T of elliptic curves in the family "x=p^2+q^2-4*(q^2-p^2)*q, 0 < x < 1/108". This sequence serves as a counterexample to the naive idea that elliptic integrals will always evaluate to a hypergeometric function such as 2F1(a,b;c;x).
A300058 is the complex period-energy function, after scaling energy and time dimensions such that all a(n) are integers and a(0)=1. The Picard-Fuchs equation is "(12-288*x+9216*x^2)*T(x) + (-1+232*x-8160*x^2+82944*x^3)*T'(x) + (-x+164*x^2-6432*x^3+41472*x^4)*T''(x)".
Although the sequence is not generated by a hypergeometric function, it can be formulated in terms of Hypergeometric numbers, specifically the binomial coefficients. Then Zeilberger's algorithm outputs a second order recurrence with polynomial coefficients.
The contour plot is nice to look at, with reflection symmetry, three critical points, and two separatrices dividing the phase plane into eight distinct regions.
Hyperbolic Critical points are located at (q,p) locations (1/6,0) and (-1/4,sqrt(5)/4) and (-1/4,-sqrt(5)/4). Is it possible to use chord-and-tangent addition rules to produce an exponentially-convergent Diophantine approximation to sqrt(5) that moves along the upper separatrix x=1/8?
Does there exist a period-preserving transformation that takes any one of the curves with 0 < x < 1/108 into a particular Weierstrass curve from the L-function and Modular Forms Database?

References

  • D. Husemöller, Elliptic Curves, 2nd ed., New York: Springer, 2004.
  • J. H. Silverman, The Arithmetic of Elliptic Curves, 2nd ed., New York: Springer, 2009.

Crossrefs

Factors: A005190, A005809, A005721.
Complex Period: A300058.

Programs

  • Mathematica
    b[NN_]:=Total/@Table[((-1)^k)*Binomial[3*n,n]*Binomial[2*n,k]*Binomial[5*n-4*k-1,3*n-4*k],{n,0,NN},{k,0,Floor[3*n/4]}];
    c1=8*(-30+201*n-319*n^2+145*n^3);c2=-8640*(n-5/3)*(n-4/3)*(n-1/5);c3=10*(n-6/5)*n^2;a[0]=1;a[1]=12;a[n0_]:=ReplaceAll[(c1/c3)*a[n0-1]+(c2/c3)*a[n0-2],{n->n0}];
    ({#,SameQ[a/@Range[0, 15],#]}&@b[15])[[1]]

Formula

a(n) = A005809(n)*A005721(n).
a(n) = Sum_{k=0..floor(3n/4)} ((-1)^k)*binomial(3*n,n)*binomial(2 *n, k)*binomial(5*n - 4*k - 1, 3*n - 4*k).
c1 = 8 *(-30 + 201*n - 319*n^2 + 145*n^3); c2 = -8640*(n - 5/3)*(n - 4/3)*(n - 1/5); c3 = 10*(n - 6/5)*n^2; a(0)=1; a(1)=12; a(n) = (c1/c3)*a(n-1) + (c2/c3)*a(n-2).

A300058 a(n) = binomial(3*n,n)/(2*Pi)*Integral_{x=0..2*Pi} (12*cos^2(x)*sin(x) + 20*sin^3(x))^(2*n) dx.

Original entry on oeis.org

1, 492, 707220, 1298204880, 2654173160100, 5765723073622512, 13021894087331233104, 30217387890886676251200, 71532102917478013611243300, 171944976047709681477985038000, 418347201888204996027087975427920
Offset: 0

Views

Author

Bradley Klee, Feb 23 2018

Keywords

Comments

Compare with A295870. The series expansion "T(x)=2*Pi*sqrt(3/5)*Sum_{n>=0} a_n*(x/25)^n" determines the period T of anharmonic oscillation along a contour of the Hamiltonian energy surface "x=2H=(5/3)*p^2+q^2+4*(p^2+q^2)*q,0
The period-energy function T(x) satisfies the Picard-Fuchs equation "(2460+28512*x+2239488*x^2)*T(x)-(125-24840*x-1423008*x^2-20155392*x^3)*T'(x)+(-125*x+1620*x^2+1189728*x^3+10077696 x^4)*T''(x)", also the P.F.Eq. of A295870 under transformation x->x'=1/108-x.
A300057 has a similar definition to A005721, with a couple of extra integers appearing in the integrand. This makes a nice analogy between real and complex periods A295870, A300058. Second-order recurrences with polynomial coefficients define both sequences.

Crossrefs

Cf. A002894, A113424, A006480, A000897. Factors: A005809, A300057. Real Period: A295870.

Programs

  • Maple
    a := n -> 36^n*(3*n)!/n!^3*hypergeom([-2*n, n+1/2], [n+1], -2/3):
    seq(simplify(a(n)), n=0..10); # Peter Luschny, Apr 19 2018
  • Mathematica
    c1=12*(-230+2259*n-3933*n^2+1863*n^3);c2=5248800*(n-5/3)*(n-4/3)*(n-1/9);c3=9*n^2*(n-10/9);a[0]=1;a[1]=492;a[n0_]:=ReplaceAll[(c1/c3)*a[n0-1]+(c2/c3)*a[n0-2],n->n0];
    b[NN_]:=Expand[Total[Flatten[#]]&/@Table[Binomial[3*n,n]*Binomial[2*n,k2]*Binomial[2*n,k1]*Binomial[2*n,3*n-k1-k2]*((4+Sqrt[15])^(2*n-k1))*((4-Sqrt[15])^(2*n-k2)),{n,0,NN},{k1,0,2*n},{k2,0,2*n}]]; ({#,SameQ[#,a/@Range[0,10]]}&@b[10])[[1]]
    Table[Binomial[3*n, n] * SeriesCoefficient[(1 + 9*z + 9*z^2 + z^3)^(2*n), {z, 0, 3*n}], {n, 0, 15}] (* Vaclav Kotesovec, Apr 18 2018 *)

Formula

a(n) = A005809(n)*A300057(n).
a(n) = Sum_{k1=0..2n} Sum_{k2=0..2n} binomial(3*n,n)*binomial(2*n,k2)*binomial(2*n,k1)*binomial(2*n,3*n-k1-k2)*((4+sqrt(15))^(2*n-k1))*((4-sqrt(15))^(2*n-k2))
a(0) = 1; a(1) = 492; a(n):=(c1/c3)*a(n-1)+(c2/c3)*a(n-2); with
c1 = 12*(-230+2259*n-3933*n^2+1863*n^3);
c2 = 5248800*(n-5/3)*(n-4/3)*(n-1/9);
c3 = 9*n^2*(n-10/9);
a(n) ~ 2^(2*n - 1) * 3^(3*n - 1/2) * 5^(2*n + 1/2) / (Pi*n). - Vaclav Kotesovec, Apr 18 2018

A308837 The nome q=exp(T_C/T_R)=Sum_{n>=0} a(n)*(x/432)^n follows from the series solutions of 5*T-d/dx(36*(1-x)*x*dT/dx)=0.

Original entry on oeis.org

0, 1, 312, 107604, 39073568, 14645965026, 5609733423408, 2182717163349896, 859521859502348352, 341679883727799750159, 136868519056531319862408, 55173969942211048781835468, 22360181278518828446785034976, 9103073677708423854325869548662
Offset: 0

Author

Bradley Klee, Jun 27 2019

Keywords

Comments

Also appears in Ramanujan's theory of elliptic functions, signature 6 (cf. A113424). Almkvist et al. give a real and complex Ansatz for the second-order, ordinary differential equation: T_R = 1 + x*{Z[[x]]}, T_C = T_R*log(x) + x*{Z[[x]]}.

References

  • B.C. Berndt, "Ramanujan's Notebooks Part II", Springer, 2012, pages 80-82.

Crossrefs

Programs

  • Mathematica
    G[nMax_] := Dot[RecurrenceTable[{Dot[{(6*n - 11)^2 (6*n - 7)^2 (18*n - 5), -36 (n - 1) (385 - 2426*n + 4968*n^2 - 4248*n^3 + 1296*n^4), 1296 (n - 1) n^3 (18*n - 23)},
    a[n - #] & /@ Reverse[Range[0, 2]]] == 0, a[0] == 0, a[1] == 13/18}, a, {n, 0, nMax}], x^Range[0, nMax]];
    qSer[nMax_] := Expand[Times[x,Normal[ Series[Exp[Divide[G[nMax], Hypergeometric2F1[1/6, 5/6, 1, x]]], {x, 0, nMax}]]]];
    CoefficientList[(1/k)*qSer[12] /. {x -> k*x}, x] /. {k -> 432}

A368545 a(n) = (6*n)!*(9*n)!/((2*n)!*(3*n)!*((5*n)!)^2).

Original entry on oeis.org

1, 1512, 13477464, 156037527360, 2018989675062360, 27791880541579666512, 397974308267103144857280, 5857301482132559510090140920, 87960819866039245485380114158680, 1341476995616395775866469804150505600, 20709814202640945368698086671062023292464, 322897098547446926110838713561401344838974976
Offset: 0

Author

Karol A. Penson, Dec 29 2023

Keywords

Crossrefs

Programs

  • Maple
    seq((6*n)!*(9*n)!/((2*n)!*(3*n)!*((5*n)!)^2),n=0..12);

Formula

G.f.: hypergeometric10F9([1/9, 1/6, 2/9, 1/3, 4/9, 5/9, 2/3, 7/9, 5/6, 8/9], [1/5, 1/5, 2/5, 2/5, 3/5, 3/5, 4/5, 4/5, 1], (167365651248*z)/9765625).
E.g.f.: hypergeometric10F10([1/9, 1/6, 2/9, 1/3, 4/9, 5/9, 2/3, 7/9, 5/6, 8/9], [1/5, 1/5, 2/5, 2/5, 3/5, 3/5, 4/5, 4/5, 1,1], (167365651248*z)/9765625).
a(n) = Integral_{x=0..167365651248/9765625} x^n*W(x) dx, n>=0, where W(x) = (1953125*MeijerG([[], [-4/5, -4/5, -3/5, -3/5, -2/5, -2/5, -1/5, -1/5, 0, 0]], [[-1/9, -1/6, -2/9, -1/3, -4/9, -5/9, -2/3, -7/9, -5/6, -8/9], []], (9765625*x)/167365651248))/(111577100832*Pi). MeijerG is the Meijer G - function. This integral representation as the n-th power moment of the W(x) function in the interval [0, 167365651248/9765625] is verified by direct integration. However, we are unable at present to determine the geometric shape of W(x).

A347304 a(n) = n!/(floor(n/2)!*floor(n/3)!*floor(n/6)!).

Original entry on oeis.org

1, 1, 2, 6, 12, 60, 60, 420, 840, 2520, 5040, 55440, 13860, 180180, 360360, 1081080, 2162160, 36756720, 4084080, 77597520, 155195040, 465585120, 931170240, 21416915520, 1338557220, 33463930500, 66927861000, 200783583000, 401567166000, 11645447814000, 465817912560
Offset: 0

Author

N. J. A. Sloane, Aug 28 2021

Keywords

Crossrefs

a(6*n) = A113424(n).

Programs

  • PARI
    a(n) = n!/((n\2)!*(n\3)!*(n\6)!); \\ Seiichi Manyama, Aug 28 2021
    
  • Python
    from math import factorial
    def A347304(n): return factorial(n)//factorial(n//2)//factorial(n//3)//factorial(n//6) # Chai Wah Wu, Aug 28 2021

Extensions

a(0)=1 prepended by Seiichi Manyama, Aug 28 2021
Showing 1-10 of 15 results. Next