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-5 of 5 results.

A000172 The Franel number a(n) = Sum_{k = 0..n} binomial(n,k)^3.

Original entry on oeis.org

1, 2, 10, 56, 346, 2252, 15184, 104960, 739162, 5280932, 38165260, 278415920, 2046924400, 15148345760, 112738423360, 843126957056, 6332299624282, 47737325577620, 361077477684436, 2739270870994736, 20836827035351596, 158883473753259752, 1214171997616258240
Offset: 0

Views

Author

Keywords

Comments

Cusick gives a general method of deriving recurrences for the r-th order Franel numbers (this is the sequence of third-order Franel numbers), with floor((r+3)/2) terms.
This is the Taylor expansion of a special point on a curve described by Beauville. - Matthijs Coster, Apr 28 2004
An identity of V. Strehl states that a(n) = Sum_{k = 0..n} C(n,k)^2 * binomial(2*k,n). Zhi-Wei Sun conjectured that for every n = 2,3,... the polynomial f_n(x) = Sum_{k = 0..n} binomial(n,k)^2 * binomial(2*k,n) * x^(n-k) is irreducible over the field of rational numbers. - Zhi-Wei Sun, Mar 21 2013
Conjecture: a(n) == 2 (mod n^3) iff n is prime. - Gary Detlefs, Mar 22 2013
a(p) == 2 (mod p^3) for any prime p since p | C(p,k) for all k = 1,...,p-1. - Zhi-Wei Sun, Aug 14 2013
a(n) is the maximal number of totally mixed Nash equilibria in games of 3 players, each with n+1 pure options. - Raimundas Vidunas, Jan 22 2014
This is one of the Apéry-like sequences - see Cross-references. - Hugo Pfoertner, Aug 06 2017
Diagonal of rational functions 1/(1 - x*y - y*z - x*z - 2*x*y*z), 1/(1 - x - y - z + 4*x*y*z), 1/(1 + y + z + x*y + y*z + x*z + 2*x*y*z), 1/(1 + x + y + z + 2*(x*y + y*z + x*z) + 4*x*y*z). - Gheorghe Coserea, Jul 04 2018
a(n) is the constant term in the expansion of ((1 + x) * (1 + y) + (1 + 1/x) * (1 + 1/y))^n. - Seiichi Manyama, Oct 27 2019
Diagonal of rational function 1 / ((1-x)*(1-y)*(1-z) - x*y*z). - Seiichi Manyama, Jul 11 2020
Named after the Swiss mathematician Jérôme Franel (1859-1939). - Amiram Eldar, Jun 15 2021
It appears that a(n) is equal to the coefficient of (x*y*z)^n in the expansion of (1 + x + y - z)^n * (1 + x - y + z)^n * (1 - x + y + z)^n. Cf. A036917. - Peter Bala, Sep 20 2021

Examples

			O.g.f.: A(x) = 1 + 2*x + 10*x^2 + 56*x^3 + 346*x^4 + 2252*x^5 + ...
O.g.f.: A(x) = 1/(1-2*x) + 3!*x^2/(1-2*x)^4 + (6!/2!^3)*x^4/(1-2*x)^7 + (9!/3!^3)*x^6/(1-2*x)^10 + (12!/4!^3)*x^8/(1-2*x)^13 + ... - _Paul D. Hanna_, Oct 30 2010
Let g.f. A(x) = Sum_{n >= 0} a(n)*x^n/n!^3, then
A(x) = 1 + 2*x + 10*x^2/2!^3 + 56*x^3/3!^3 + 346*x^4/4!^3 + ... where
A(x) = [1 + x + x^2/2!^3 + x^3/3!^3 + x^4/4!^3 + ...]^2. - _Paul D. Hanna_
		

References

  • Matthijs Coster, Over 6 families van krommen [On 6 families of curves], Master's Thesis (unpublished), Aug 26 1983.
  • Jérôme Franel, On a question of Laisant, Intermédiaire des Mathématiciens, vol 1 1894 pp 45-47
  • H. W. Gould, Combinatorial Identities, Morgantown, 1972, (X.14), p. 56.
  • Murray Klamkin, ed., Problems in Applied Mathematics: Selections from SIAM Review, SIAM, 1990; see pp. 148-149.
  • John Riordan, An Introduction to Combinatorial Analysis, Wiley, 1958, p. 193.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A002893, A052144, A005260, A096191, A033581, A189791. Second row of array A094424.
The Apéry-like numbers [or Apéry-like sequences, Apery-like numbers, Apery-like sequences] include A000172, A000984, A002893, A002895, A005258, A005259, A005260, A006077, A036917, A063007, A081085, A093388, A125143 (apart from signs), A143003, A143007, A143413, A143414, A143415, A143583, A183204, A214262, A219692,A226535, A227216, A227454, A229111 (apart from signs), A260667, A260832, A262177, A264541, A264542, A279619, A290575, A290576. (The term "Apery-like" is not well-defined.)
For primes that do not divide the terms of the sequences A000172, A005258, A002893, A081085, A006077, A093388, A125143, A229111, A002895, A290575, A290576, A005259 see A260793, A291275-A291284 and A133370 respectively.
Sum_{k = 0..n} C(n,k)^m for m = 1..12: A000079, A000984, A000172, A005260, A005261, A069865, A182421, A182422, A182446, A182447, A342294, A342295.
Column k=3 of A372307.

Programs

  • Haskell
    a000172 = sum . map a000578 . a007318_row
    -- Reinhard Zumkeller, Jan 06 2013
    
  • Maple
    A000172 := proc(n)
        add(binomial(n,k)^3,k=0..n) ;
    end proc:
    seq(A000172(n),n=0..10) ; # R. J. Mathar, Jul 26 2014
    A000172_list := proc(len) series(hypergeom([], [1, 1], x)^2, x, len);
    seq((n!)^3*coeff(%, x, n), n=0..len-1) end:
    A000172_list(21); # Peter Luschny, May 31 2017
  • Mathematica
    Table[Sum[Binomial[n,k]^3,{k,0,n}],{n,0,30}] (* Harvey P. Dale, Aug 24 2011 *)
    Table[ HypergeometricPFQ[{-n, -n, -n}, {1, 1}, -1], {n, 0, 20}]  (* Jean-François Alcover, Jul 16 2012, after symbolic sum *)
    a[n_] := Sum[ Binomial[2k, n]*Binomial[2k, k]*Binomial[2(n-k), n-k], {k, 0, n}]/2^n; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Mar 20 2013, after Zhi-Wei Sun *)
    a[ n_] := SeriesCoefficient[ Hypergeometric2F1[ 1/3, 2/3, 1, 27 x^2 / (1 - 2 x)^3] / (1 - 2 x), {x, 0, n}]; (* Michael Somos, Jul 16 2014 *)
  • PARI
    {a(n)=polcoeff(sum(m=0,n,(3*m)!/m!^3*x^(2*m)/(1-2*x+x*O(x^n))^(3*m+1)),n)} \\ Paul D. Hanna, Oct 30 2010
    
  • PARI
    {a(n)=n!^3*polcoeff(sum(m=0,n,x^m/m!^3+x*O(x^n))^2,n)} \\ Paul D. Hanna, Jan 19 2011
    
  • PARI
    A000172(n)={sum(k=0,(n-1)\2,binomial(n,k)^3)*2+if(!bittest(n,0),binomial(n,n\2)^3)} \\ M. F. Hasler, Sep 21 2015
    
  • Sage
    def A000172():
        x, y, n = 1, 2, 1
        while True:
            yield x
            n += 1
            x, y = y, (8*(n-1)^2*x + (7*n^2-7*n + 2)*y) // n^2
    a = A000172()
    [next(a) for i in range(21)]   # Peter Luschny, Oct 12 2013

Formula

A002893(n) = Sum_{m = 0..n} binomial(n, m)*a(m) [Barrucand].
Sum_{k = 0..n} C(n, k)^3 = (-1)^n*Integral_{x = 0..infinity} L_k(x)^3 exp(-x) dx. - from Askey's book, p. 43.
D-finite with recurrence (n + 1)^2*a(n+1) = (7*n^2 + 7*n + 2)*a(n) + 8*n^2*a(n-1) [Franel]. - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jan 31 2001
a(n) ~ 2*3^(-1/2)*Pi^-1*n^-1*2^(3*n). - Joe Keane (jgk(AT)jgk.org), Jun 21 2002
O.g.f.: A(x) = Sum_{n >= 0} (3*n)!/n!^3 * x^(2*n)/(1 - 2*x)^(3*n+1). - Paul D. Hanna, Oct 30 2010
G.f.: hypergeom([1/3, 2/3], [1], 27 x^2 / (1 - 2x)^3) / (1 - 2x). - Michael Somos, Dec 17 2010
G.f.: Sum_{n >= 0} a(n)*x^n/n!^3 = [ Sum_{n >= 0} x^n/n!^3 ]^2. - Paul D. Hanna, Jan 19 2011
G.f.: A(x) = 1/(1-2*x)*(1+6*(x^2)/(G(0)-6*x^2)),
with G(k) = 3*(x^2)*(3*k+1)*(3*k+2) + ((1-2*x)^3)*((k+1)^2) - 3*(x^2)*((1-2*x)^3)*((k+1)^2)*(3*k+4)*(3*k+5)/G(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Dec 03 2011
In 2011 Zhi-Wei Sun found the formula Sum_{k = 0..n} C(2*k,n)*C(2*k,k)*C(2*(n-k),n-k) = (2^n)*a(n) and proved it via the Zeilberger algorithm. - Zhi-Wei Sun, Mar 20 2013
0 = a(n)*(a(n+1)*(-2048*a(n+2) - 3392*a(n+3) + 768*a(n+4)) + a(n+2)*(-1280*a(n+2) - 2912*a(n+3) + 744*a(n+4)) + a(n+3)*(+288*a(n+3) - 96*a(n+4))) + a(n+1)*(a(n+1)*(-704*a(n+2) - 1232*a(n+3) + 288*a(n+4)) + a(n+2)*(-560*a(n+2) - 1372*a(n+3) + 364*a(n+4)) + a(n+3)*(+154*a(n+3) - 53*a(n+4))) + a(n+2)*(a(n+2)*(+24*a(n+2) + 70*a(n+3) - 20*a(n+4)) + a(n+3)*(-11*a(n+3) + 4*a(n+4))) for all n in Z. - Michael Somos, Jul 16 2014
For r a nonnegative integer, Sum_{k = r..n} C(k,r)^3*C(n,k)^3 = C(n,r)^3*a(n-r), where we take a(n) = 0 for n < 0. - Peter Bala, Jul 27 2016
a(n) = (n!)^3 * [x^n] hypergeom([], [1, 1], x)^2. - Peter Luschny, May 31 2017
From Gheorghe Coserea, Jul 04 2018: (Start)
a(n) = Sum_{k=0..floor(n/2)} (n+k)!/(k!^3*(n-2*k)!) * 2^(n-2*k).
G.f. y=A(x) satisfies: 0 = x*(x + 1)*(8*x - 1)*y'' + (24*x^2 + 14*x - 1)*y' + 2*(4*x + 1)*y. (End)
a(n) = [x^n] (1 - x^2)^n*P(n,(1 + x)/(1 - x)), where P(n,x) denotes the n-th Legendre polynomial. See Gould, p. 56. - Peter Bala, Mar 24 2022
a(n) = (2^n/(4*Pi^2)) * Integral_{x,y=0..2*Pi} (1+cos(x)+cos(y)+cos(x+y))^n dx dy = (8^n/(Pi^2)) * Integral_{x,y=0..Pi} (cos(x)*cos(y)*cos(x+y))^n dx dy (Pla, 1995). - Amiram Eldar, Jul 16 2022
a(n) = Sum_{k = 0..n} m^(n-k)*binomial(n,k)*binomial(n+2*k,n)*binomial(2*k,k) at m = -4. Cf. A081798 (m = 1), A006480 (m = 0), A124435 (m = -1), A318109 (m = -2) and A318108 (m = -3). - Peter Bala, Mar 16 2023
From Bradley Klee, Jun 05 2023: (Start)
The g.f. T(x) obeys a period-annihilating ODE:
0=2*(1 + 4*x)*T(x) + (-1 + 14*x + 24*x^2)*T'(x) + x*(1 + x)*(-1 + 8*x)*T''(x).
The periods ODE can be derived from the following Weierstrass data:
g2 = (4/243)*(1 - 8*x + 240*x^2 - 464*x^3 + 16*x^4);
g3 = -(8/19683)*(1 - 12*x - 480*x^2 + 3080*x^3 - 12072*x^4 + 4128*x^5 +
64*x^6);
which determine an elliptic surface with four singular fibers. (End)
From Peter Bala, Oct 31 2024: (Start)
For n >= 1, a(n) = 2 * Sum_{k = 0..n-1} binomial(n, k)^2 * binomial(n-1, k). Cf. A361716.
For n >= 1, a(n) = 2 * hypergeom([-n, -n, -n + 1], [1, 1], -1). (End)

A081656 a(n) = 2*6^n - 2^n.

Original entry on oeis.org

1, 10, 68, 424, 2576, 15520, 93248, 559744, 3358976, 20154880, 120931328, 725592064, 4353560576, 26121379840, 156728311808, 940369936384, 5642219749376, 33853318758400, 203119913074688, 1218719479496704, 7312316879077376, 43873901278658560, 263243407680339968
Offset: 0

Views

Author

Paul Barry, Mar 26 2003

Keywords

Comments

Binomial transform of A081655.
Inverse binomial transform of A081657.

Crossrefs

Third column of array A094424.

Programs

  • Magma
    [2*6^n-2^n: n in [0..30]]; // Vincenzo Librandi, Aug 10 2013
    
  • Mathematica
    CoefficientList[Series[(1 + 2 x) / ((1 - 6 x) (1 - 2 x)), {x, 0, 40}], x] (* Vincenzo Librandi, Aug 10 2013 *)
    Table[2*6^n-2^n,{n,0,20}] (* or *) LinearRecurrence[{8,-12},{1,10},20] (* Harvey P. Dale, Mar 17 2019 *)
  • PARI
    a(n)=2*6^n-2^n \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: (1+2*x)/((1-6*x)(1-2*x)).
E.g.f.: 2*exp(6*x) - exp(2*x).
a(0)=1, a(2)=10; for n>2, a(n) = 8*a(n-1) - 12*a(n-2). - Vincenzo Librandi, Aug 10 2013

A000658 Strehl's sequence "C_n^(3)".

Original entry on oeis.org

1, 4, 68, 1732, 51076, 1657904, 57793316, 2117525792, 80483121028, 3147565679824, 125937573689968, 5133632426499152, 212530848994367524, 8914634034287235856, 378138515326996979168, 16196097181014298854032
Offset: 0

Views

Author

Keywords

References

  • Volker Strehl, Binomial identities - combinatorial and algorithmic aspects. Trends in discrete mathematics. Discrete Math. 136 (1994), no. 1-3, 309-346.

Crossrefs

Third row of array A094424.
Cf. A007318.

Programs

  • Haskell
    a000658 n = sum $ map c3 [0..n] where
       c3 k = (a007318' n k)^2 * (a007318' (2*k) k)^2 * a007318' (2*k) (n-k)
    -- Reinhard Zumkeller, May 20 2015
  • Maple
    A000658:=n->add(binomial(n,k)^2*binomial(2*k,k)^2*binomial(2*k,n-k), k=0..n): seq(A000658(n), n=0..15); # Wesley Ivan Hurt, Sep 19 2014
  • Mathematica
    Table[Sum[Binomial[n,k]^2 Binomial[2k,k]^2 Binomial[2k,n-k],{k,0,n}], {n,0,25}] (* Harvey P. Dale, Oct 19 2011 *)
  • PARI
    a(n)=sum(k=1,n, binomial(n, k)^2 * binomial(2k, k)^2 * binomial(2k, n-k)) \\ Charles R Greathouse IV, Sep 19 2014
    

Formula

Sum binomial(n, k)^2 * binomial(2k, k)^2 * binomial(2k, n-k); k=0..n.
a(n) ~ 7^(2*n+5/2) / (20 * sqrt(15) * Pi^2 * n^2). - Vaclav Kotesovec, Mar 09 2014

A092868 Coefficients c[r,n] in Schmidt's problem Sum[Binomial[n,k]^r Binomial[n+k,k]^r,{k,0,n}] == Sum[Binomial[n,k]Binomial[n+k,k]c[r,k],{k,0,n}] for r=4.

Original entry on oeis.org

1, 8, 424, 48896, 6672232, 1022309408, 176808084544, 33055112886272, 6507475475389288, 1336577286762538496, 284198765977135568224, 62135041920796512325952, 13901968841738902540019776
Offset: 0

Views

Author

Eric W. Weisstein, Mar 08 2004

Keywords

Crossrefs

Fourth row of array A094424.

Programs

  • Mathematica
    c[4, n_] := Sum[Binomial[2j, j]^3Binomial[n, j]Binomial[k+j, k-j]Binomial[j, n-k]Binomial[k, j]Binomial[2j, k-j], {k, 0, n}, {j, 0, n}]

Formula

a(n) ~ sqrt(3) * 2^(5*n+6) * 3^(2*n+3) / (17^(5/2) * Pi^3 * n^3). - Vaclav Kotesovec, Mar 09 2014

A379610 a(n) = Sum_{j=0..n} binomial(2*j,j)^4 * binomial(n,j)^2 * Sum_{k=0..n} binomial(k+j,k-j)^2 * binomial(2*j,n-k) * binomial(2*j,k-j).

Original entry on oeis.org

1, 16, 2576, 1383568, 873960976, 615833930816, 526152430612496, 502263183380000576, 507670642053018634768, 542989589742114444434176, 609997002301177503142835776, 710665335270156912049219940096, 853208134042451055649274133396496, 1051685299255545900736463773134099328
Offset: 0

Views

Author

Seiichi Manyama, Dec 27 2024

Keywords

Crossrefs

Fifth row of array A094424.
Cf. A092815.

Programs

  • PARI
    a(n) = sum(j=0, n, binomial(2*j, j)^4*binomial(n, j)^2*sum(k=0, n, binomial(k+j, k-j)^2*binomial(2*j, n-k)*binomial(2*j, k-j)));

Formula

A092815(n) = Sum_{k=0..n} binomial(n,k) * binomial(n+k,k) * a(k).
Showing 1-5 of 5 results.