A109971 -id:A109971 - cp's OEIS frontend

A005809 a(n) = binomial(3n,n).

1, 3, 15, 84, 495, 3003, 18564, 116280, 735471, 4686825, 30045015, 193536720, 1251677700, 8122425444, 52860229080, 344867425584, 2254848913647, 14771069086725, 96926348578605, 636983969321700, 4191844505805495, 27619435402363035

Offset: 0

N. J. A. Sloane

Number of paths in Z X Z starting at (0,0) and ending at (3n,0) using steps in {(1,1),(1,-2)}.
Number of even trees with 2n edges and one distinguished vertex. Even trees are rooted plane trees where every vertex (including root) has even degree.
Hankel transform is 3^n*A051255(n), where A051255 is the Hankel transform of C(3n,n)/(2n+1). - Paul Barry, Jan 21 2007
a(n) is the number of stack polyominoes inscribed in an (n+1) X (n+1) box. Equivalently, a(n) is the number of unimodal compositions with n+1 parts in which the maximum value of the parts is n+1. For instance, for n = 2, we have the following compositions: (3,3,3), (2,3,3), (1,3,3), (3,3,1), (3,3,2), (2,2,3), (1,2,3), (2,3,1), (1,1,3), (1,3,1), (3,1,1), (2,3,2), (1,3,2), (3,2,1), (3,2,2). - Emanuele Munarini, Apr 07 2011
Conjecture: a(n)==3 (mod n^3) iff n is an odd prime. - Gary Detlefs, Mar 23 2013. The congruence a(p) = binomial(3*p,p) = 3 (mod p^3) for odd prime p is a known generalization of Wolstenholme's theorem. See Mestrovic, Section 6, equation 35. - Peter Bala, Dec 28 2014
In general, C(k*n,n) = C(k*n-1,n-1)*C((k*n)^2,2)/(3*n*C(k*n+1,3)), n>0. - Gary Detlefs, Jan 02 2014

			G.f. = 1 + 3*x + 15*x^2 + 84*x^3 + 495*x^4 + 3003*x^5 + 18564*x^6 + ... - _Michael Somos_, Jan 30 2019

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 828.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

G. C. Greubel, Table of n, a(n) for n = 0..1000[Terms 0 to 100 computed by T. D. Noe; terms 101 to 1000 by G. C. Greubel, Jan 14 2017]
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Paul Barry, On the Central Coefficients of Riordan Matrices, Journal of Integer Sequences, 16 (2013), Article 13.5.1.
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
Naiomi Tuere Cameron, Random walks, trees and extensions of Riordan group techniques, Dissertation, Howard University, 2002.
Tom C. Copeland, Discriminating Deltas, Depressed Equations, and Generalized Catalan Numbers
Maciej Dziemianczuk, On Directed Lattice Paths With Additional Vertical Steps, arXiv:1410.5747 [math.CO], 2014.
Milan Janjic, Two Enumerative Functions.
Vaclav Kotesovec, Non-attacking chess pieces, 6ed, 2013, p. 436.
Yaping Liu, On the Recursiveness of Pascal Sequences, Global J. of Pure and Appl. Math. (2022) Vol. 18, No. 1, 71-80.
Mathematics Stack Exchange, Ordinary generating function for binomial(3n,n), Nov 2013.
Romeo Meštrović, Wolstenholme's theorem: Its Generalizations and Extensions in the last hundred and fifty years (1862-2011), arXiv:1111.3057 [math.NT], 2011.
Wojciech Młotkowski and Karol A. Penson, Probability distributions with binomial moments, Infinite Dimensional Analysis, Quantum Probability and Related Topics, Vol. 17, No. 2 (2014), 1450014; arXiv preprint, arXiv:1309.0595 [math.PR], 2013.
Khodabakhsh Hessami Pilehrood and Tatiana Hessami Pilehrood, Jacobi Polynomials and Congruences Involving Some Higher-Order Catalan Numbers and Binomial Coefficients, Journal of Integer Sequences, Vol. 18 (2015), Article 15.11.7.
Yash Puri and Thomas Ward, Arithmetic and growth of periodic orbits, Journal of Integer Sequences, Vol. 4 (2001), Article 01.2.1.

binomial(k*n,n): A000984 (k = 2), A005810 (k = 4), A001449 (k = 5), A004355 (k = 6), A004368 (k = 7), A004381 (k = 8), A169958 - A169961 (k = 9 thru 12).

Cf. A001764, A007318, A110608, A229705, A254157, A268196.

a005809 n = a007318 (3*n) n  -- Reinhard Zumkeller, May 06 2012

[ Binomial(3*n,n): n in [0..150] ]; // Vincenzo Librandi, Apr 21 2011

A005809:=n->binomial(3*n,n); seq(A005809(n), n=0..40); # Wesley Ivan Hurt, Mar 21 2014

R[ z_ ] := ((2-18*z + 27*z^2 + 3^(3/2)*z^(3/2)*(27*z-4)^(1/2))/2)^(1/3); f[ z_ ] := ( (R[ z ])^3 + (1-3*z)*(R[ z ])^2 + (1-6*z)*R[ z ] )/( (R[ z ])^4 + (1-6*z)*(R[ z ])^2 + (6*z-1)^2 )
Table[Binomial[3*n,n],{n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Mar 03 2011 *)

makelist(binomial(3*n,n),n,0,100); /* Emanuele Munarini, Apr 07 2011 */

B(x):=(2/sqrt(3*x))*sin((1/3)*asin(sqrt(27*x/4)))-1;
taylor(x*diff(B(x),x)/B(x),x,0,10); /* Vladimir Kruchinin, Oct 02 2015 */

a(n)=binomial(3*n,n) \\ Charles R Greathouse IV, Nov 20 2012

[binomial(3*n,n) for n in range(0, 22)] # Zerinvary Lajos, Dec 16 2009

The g.f. R[ z_ ] below (in the Mathematica field) was found by Kurt Persson (kurt(AT)math.chalmers.se) and communicated by Einar Steingrimsson (einar(AT)math.chalmers.se).
Using Stirling's formula in A000142, it is easy to get the asymptotic expression a(n) ~ (1/2) * (27/4)^n / sqrt(Pi*n / 3). - Dan Fux (dan.fux(AT)OpenGaia.com or danfux(AT)OpenGaia.com), Apr 07 2001
a(n) = Sum_{k=0..n} C(n, k)*C(2n, k). - Paul Barry, May 15 2003
G.f.: 1/(1-3zg^2), where g=g(z) is given by g=1+zg^3, g(0)=1, i.e., (in Maple notation) g := 2*sin(arcsin(3*sqrt(3*z)/2)/3)/sqrt(3*z). - Emeric Deutsch, May 22 2003
G.f.: x*B'(x)/B(x), where B(x)+1 is the g.f. for A001764. -  Vladimir Kruchinin, Oct 02 2015
a(n) ~ (1/2)*3^(1/2)*Pi^(-1/2)*n^(-1/2)*2^(-2*n)*3^(3*n)*(1 - 7/72*n^-1 + 49/10368*n^-2 + 6425/2239488*n^-3 - ...). - Joe Keane (jgk(AT)jgk.org), Nov 07 2003
a(n) = A006480(n)/A000984(n). - Lior Manor, May 04 2004
a(n) = Sum_{i_1=0..n, i_2=0..n} binomial(n, i_1)*binomial(n, i_2)*binomial(n, i_1+i_2). - Benoit Cloitre, Oct 14 2004
a(n) = Sum_{k=0..n} A109971(k)*3^k; a(0)=1, a(n) = Sum_{k=0..n} 3^k*C(3n-k,n-k)2k/(3n-k), n>0. - Paul Barry, Jan 21 2007
a(n) = A085478(2n,n). - Philippe Deléham, Sep 17 2009
E.g.f.: 2F2(1/3,2/3;1/2,1;27*x/4), where F(a1,a2;b1,b2;z) is a hypergeometric series. - Emanuele Munarini, Apr 12 2011
a(n) = Sum_{k=0..n} binomial(2*n+k-1,k). - Arkadiusz Wesolowski, Apr 02 2012
G.f.: cos((1/3)*asin(sqrt(27x/4)))/sqrt(1-27x/4). - Tom Copeland, May 24 2012
G.f.: A(x) = 1 + 6*x/(G(0)-6*x) where G(k) = (2*k+2)*(2*k+1) + 3*x*(3*k+1)*(3*k+2) - 6*x*(k+1)*(2*k+1)*(3*k+4)*(3*k+5)/G(k+1); (continued fraction). - Sergei N. Gladkovskii, Jun 30 2012
D-finite with recurrence: 2*n*(2*n-1)*a(n) - 3*(3*n-1)*(3*n-2)*a(n-1) = 0. - R. J. Mathar, Feb 05 2013
a(n) = (2n+1)*A001764(n). - Johannes W. Meijer, Aug 22 2013
a(n) = C(3*n-1,n-1)*C(9*n^2,2)/(3*n*C(3*n+1,3)), n>0. - Gary Detlefs, Jan 02 2014
a(n) = [x^n] 1/(1 - x)^(2*n+1). - Ilya Gutkovskiy, Oct 03 2017
a(n) = hypergeom([-2*n, -n], [1], 1). - Peter Luschny, Mar 19 2018
a(n) = Sum_{k=0..n} binomial(n, k) * binomial(2*n, n-k) = row sums of A110608. - Michael Somos, Jan 30 2019
0 = a(n)*(-3188646*a(n+2) +7322076*a(n+3) -2805111*a(n+4) +273585*a(n+5)) +a(n+1)*(+413343*a(n+2) -1252017*a(n+3) +538344*a(n+4) -55940*a(n+5)) +a(n+2)*(-4131*a(n+2) +38733*a(n+3) -21628*a(n+4) +2528*a(n+5)) for all n in Z. - Michael Somos, Jan 30 2019
Sum_{n>=1} 1/a(n) = A229705. - Amiram Eldar, Nov 14 2020
From Peter Bala, Feb 20 2022: (Start)
The o.g.f. A(x) satisfies the differential equation (4*x - 27*x^2)*A''(x) + (2 - 54*x)*A'(x) - 6*A(x) = 0, with A(0) = 1 and A'(0) = 3.
Algebraic equation: (1 - A(x))*(1 + 2*A(x))^2 +  27*x*A(x)^3 = 0.
Sum_{n >= 1} a(n)*( x*(2*x + 3)^2/(27*(1 + x)^3) )^n = x. (End)
From Vaclav Kotesovec, May 13 2022: (Start)
Sum_{n>=0} a(n) / 3^(2*n) = 2*cos(Pi/9).
Sum_{n>=0} a(n) / (27/2)^n = (1 + sqrt(3))/2.
Sum_{n>=0} a(n) / 3^(3*n) = 2*cos(Pi/18) / sqrt(3).
In general, for k > 27/4, Sum_{n>=0} a(n)/k^n = 2*cos(arccos(1 - 27/(2*k))/6) / sqrt(4 - 27/k). (End)
G.f.: hypergeom([1/3, 2/3], [1/2], 27*z/4), the Gauss hypergeometric function 2F1. - Karol A. Penson, Dec 12 2023
a(n) = 1/4^n * Sum_{k = n..3*n} binomial(k, n)*binomial(3*n, k). - Peter Bala, Jun 29 2025

A047098 a(n) = 2binomial(3n, n) - Sum_{k=0..n} binomial(3*n, k).

Original entry on oeis.org

1, 2, 8, 38, 196, 1062, 5948, 34120, 199316, 1181126, 7080928, 42860534, 261542752, 1607076200, 9934255472, 61732449648, 385393229460, 2415935640198, 15200964233864, 95962904716402, 607640599286276, 3858198001960438, 24559243585545644, 156692889782067712

Offset: 0

Clark Kimberling, Aug 15 1998

nonn

T(2n,n), array T as in A047089. [Corrected Dec 08 2006]
Let B_3^+ denote the semigroup with presentation . Let D=aba be the 'fundamental word'. Then this sequence is also equal to the number of words in B_3^+ equal in B_3^+ to D^n, n >= 0. - Stephen P. Humphries, Jan 20 2004
In the language of Riordan arrays, row sums of (1/(1+x), x/(1+x)^3)^-1, where (1/(1+x), x/(1+x)^3) has general term (-1)^(n-k)*binomial(n+2k, 3k). - Paul Barry, May 09 2005
Hankel transform is 2^n*A051255(n) where A051255 is the Hankel transform of C(3n,n)/(2n+1). - Paul Barry, Jan 21 2007

Michael De Vlieger, Table of n, a(n) for n = 0..1211
Paul Barry, Centered polygon numbers, heptagons and nonagons, and the Robbins numbers, arXiv:2104.01644 [math.CO], 2021.
Christopher R. Cornwell and Stephen P. Humphries, Counting fundamental paths in certain Garside semigroups, Journal of Knot Theory and Its Ramifications, Vol. 17 (2008), No. 02, pp. 191-211.

Column k=2 of A213028.

Cf. A047089, A047099, A107026, A321957.

A047098 := n -> 2*binomial(3*n, n)-add(binomial(3*n, k), k=0..n);

Table[2Binomial[3n,n]-Sum[Binomial[3n,k],{k,0,n}],{n,0,35}] (* Harvey P. Dale, Jul 27 2011 *)

a(n)=if(n<0,0,polcoeff((((1+10*x-2*x^2)+(1-4*x)*sqrt(1-4*x+x*O(x^n)))/2)^n,n))

a(n)=if(n<0,0, 2*binomial(3*n,n)-sum(k=0,n,binomial(3*n,k)))

G.f. A(x)=y satisfies (8x-1)y^3-y^2+y+1=0. - Michael Somos, Jan 28 2004
Coefficient of x^n in ((1+10x-2x^2+(1-4x)^(3/2))/2)^n. - Michael Somos, Sep 25 2003
a(n) = Sum_{k = 0..n} A109971(k)*2^k; a(0) = 1, a(n) = Sum_{k = 0..n} 2^k*C(3n-k,n-k)*2*k/(3*n-k), n > 0. - Paul Barry, Jan 21 2007
Conjecture: 2*n*(2*n-1)*a(n) +(-71*n^2+112*n-48)*a(n-1) +3*(131*n^2-391*n+296)*a(n-2) -72*(3*n-7)*(3*n-8)*a(n-3)=0. - R. J. Mathar, Nov 30 2012
a(n) = A321957(n) + 2*binomial(3*n, n) - 8^n. - Peter Luschny, Nov 22 2018
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. - Peter Bala, Mar 05 2022

Clark Kimberling, Dec 08 2006, changed "T(3n,2n)" to "T(2n,n)" in the comment line, but observes that some of the other comments seem to apply to the sequence T(3n,2n) rather than to the sequence T(2n,n).

Edited by N. J. A. Sloane, Dec 21 2006, replacing the old definition in terms of A047089 by an explicit formula supplied by Benoit Cloitre, Oct 25 2003.

A116088 Riordan array (1, x*(1+x)^2).

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 1, 4, 1, 0, 0, 6, 6, 1, 0, 0, 4, 15, 8, 1, 0, 0, 1, 20, 28, 10, 1, 0, 0, 0, 15, 56, 45, 12, 1, 0, 0, 0, 6, 70, 120, 66, 14, 1, 0, 0, 0, 1, 56, 210, 220, 91, 16, 1, 0, 0, 0, 0, 28, 252, 495, 364, 120, 18, 1

Offset: 0

Paul Barry, Feb 04 2006

			Triangle begins as:
  1;
  0, 1;
  0, 2, 1;
  0, 1, 4,  1;
  0, 0, 6,  6,  1;
  0, 0, 4, 15,  8,  1;
  0, 0, 1, 20, 28, 10,  1;
  0, 0, 0, 15, 56, 45, 12, 1;

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (Rows 0 <= n <= 150).
Milan Janjić, Pascal Matrices and Restricted Words, J. Int. Seq., Vol. 21 (2018), Article 18.5.2.

Row sums are A002478. Diagonal sums are A094686. Inverse is (-1)^(n-k) * A109971(n,k). Unsigned version of A109970.

Flat(List([0..10], n->List([0..n], k-> Binomial(2*k, n-k) ))); # G. C. Greubel, May 09 2019

[[Binomial(2*k, n-k): k in [0..n]]: n in [0..10]]; // G. C. Greubel, May 09 2019

Flatten[Table[Binomial[2k,n-k],{n,0,20},{k,0,n}]] (* Harvey P. Dale, Oct 22 2012 *)

{T(n,k) = binomial(2*k, n-k)}; \\ G. C. Greubel, May 09 2019

[[binomial(2*k, n-k) for k in (0..n)] for n in (0..10)] # G. C. Greubel, May 09 2019

G.f.: 1/(1-x*y*(1+x)^2).

Number triangle T(n,k) = C(2*k, n-k) = C(n,k)*C(3*k,n)/C(3*k,k).

cp's OEIS Frontend

A005809 a(n) = binomial(3n,n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

Maxima

PARI

Sage

Formula

A047098 a(n) = 2*binomial(3*n, n) - Sum_{k=0..n} binomial(3*n, k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A116088 Riordan array (1, x*(1+x)^2).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

GAP

Magma

Mathematica

PARI

Sage

Formula

A109970 Riordan array (1,x(1-x)^2).

Views

Author

Keywords

Comments

Examples

Formula

A047098 a(n) = 2binomial(3n, n) - Sum_{k=0..n} binomial(3*n, k).