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.

Previous Showing 11-20 of 58 results. Next

A002057 Fourth convolution of Catalan numbers: a(n) = 4*binomial(2*n+3,n)/(n+4).

Original entry on oeis.org

1, 4, 14, 48, 165, 572, 2002, 7072, 25194, 90440, 326876, 1188640, 4345965, 15967980, 58929450, 218349120, 811985790, 3029594040, 11338026180, 42550029600, 160094486370, 603784920024, 2282138106804, 8643460269248, 32798844771700, 124680849918352
Offset: 0

Views

Author

Keywords

Comments

a(n) is sum of the (flattened) list obtained by the iteration of: replace each integer k with the list 0,...,k+1 on the starting value 0. Length of this list is Catalan(n) or A000108. - Wouter Meeussen, Nov 11 2001
a(n-2) is the number of n-th generation vertices in the tree of sequences with unit increase labeled by 3 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of standard tableaux of shape (n+2,n-1). - Emeric Deutsch, May 30 2004
a(n) = CatalanNumber(n+3) - 2*CatalanNumber(n+2). Proof. From its definition as a convolution of Catalan numbers, a(n) counts lists of 4 Dyck paths of total size (semilength) = n. Connect the 4 paths by 3 upsteps (U) and append 3 downsteps (D). This is a reversible procedure. So a(n) is also the number of Dyck (n+3)-paths that end DDD (D for downstep). Let C(n) denote CatalanNumber(n) (A000108). Since C(n+3) is the total number of Dyck (n+3)-paths and C(n+2) is the number that end UD, we have (*) C(n+3) - C(n+2) is the number of Dyck (n+3)-paths that end DD. Also, (**) C(n+2) is the number of Dyck (n+3)-paths that end UDD (change the last D in a Dyck (n+2)-path to UDD). Subtracting (**) from (*) yields a(n) = C(n+3) - 2C(n+2) as claimed. - David Callan, Nov 21 2006
Convolution square of the Catalan sequence without one of the initial "1"'s: (1 + 4x + 14x^2 + 48x^3 + ...) = (1/x^2) * square(x + 2x^2 + 5x^3 + 14x^4 + ...)
a(n) is the number of binary trees with n+3 internal nodes in which both subtrees of the root are nonempty. Cf. A068875 [Sedgewick and Flajolet]. - Geoffrey Critzer, Jan 05 2013
With offset 4, a(n) is the number of permutations on {1,2,...,n} that are 123-avoiding, i.e., do not contain a three-term monotone subsequence, for which the first ascent is at positions (4,5); for example, there are 48 123-avoiding permutations on n=7 for which the first ascent is at spots (4,5). See Connolly link. There it is shown in general that the k-th Catalan Convolution is the number of 123-avoiding permutations for which the first ascent is at (k, k+1). (For n=k, the first ascent is defined to be at positions (k,k+1) if the permutation is the decreasing permutation with no ascents.) - Anant Godbole, Jan 17 2014
With offset 4, a(n) is the number of permutations on {1,2,...,n} that are 123-avoiding and for which the integer n is in the 4th spot; see Connolly link. - Anant Godbole, Jan 17 2014
a(n) is the number of North-East lattice paths from (0,0) to (n+2,n+2) that have exactly one east step below the subdiagonal y = x-1. Details can be found in Section 3.1 in Pan and Remmel's link. - Ran Pan, Feb 04 2016
a(n) is the number of North-East lattice paths from (0,0) to (n+2,n+2) that bounce off the diagonal y = x to the right exactly once but do not bounce off y = x to the left. Details can be found in Section 4.2 in Pan and Remmel's link. - Ran Pan, Feb 04 2016
a(n) is the number of North-East lattice paths from (0,0) to (n+2,n+2) that horizontally cross the diagonal y = x exactly once but do not cross the diagonal vertically. Details can be found in Section 4.3 in Pan and Remmel's link. - Ran Pan, Feb 04 2016
Apparently also Young tableaux of (non-partition) shape [n+1, 1, 1, n+1], see example file. - Joerg Arndt, Dec 30 2023

Examples

			From _Peter Bala_, Apr 14 2017: (Start)
This sequence appears on the main diagonal of a generalized Catalan triangle. Construct a lower triangular array (T(n,k)), n,k >= 0 by placing the sequence [0,0,0,1,1,1,1,...] in the first column and then filling in the remaining entries in the array using the rule T(n,k) = T(n,k-1) + T(n-1,k). The resulting array begins
  n\k| 0 1  2  3  4   5   6   7  ...
  ---+-------------------------------
   0 | 0
   1 | 0 0
   2 | 0 0  0
   3 | 1 1  1  1
   4 | 1 2  3  4  4
   5 | 1 3  6 10 14  14
   6 | 1 4 10 20 34  48  48
   7 | 1 5 15 35 69 117 165 165
   ...
(see Tedford 2011; this is essentially the array C_4(n,k) in the notation of Lee and Oh). Compare with A279004. (End)
		

References

  • Pierre de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 11, coefficients of P_4(z).
  • C. Krishnamachary and M. Bheemasena Rao, Determinants whose elements are Eulerian, prepared Bernoullian and other numbers, J. Indian Math. Soc., Vol. 14 (1922), pp. 55-62, 122-138 and 143-146.
  • Robert Sedgewick and Phillipe Flajolet, Analysis of Algorithms, Addison Wesley, 1996, page 225.
  • 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

T(n, n+4) for n=0, 1, 2, ..., array T as in A047072. Also a diagonal of A059365 and of A009766.
Cf. A001003.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.
Cf. A145596 (row sums), A279004.

Programs

  • GAP
    List([0..25],n->4*Binomial(2*n+3,n)/(n+4)); # Muniru A Asiru, Mar 05 2018
    
  • Magma
    [4*Binomial(2*n+3,n)/(n+4): n in [0..30]]; // Vincenzo Librandi, Feb 04 2016
    
  • Maple
    a := n -> 32*4^n*GAMMA(5/2+n)*(1+n)/(sqrt(Pi)*GAMMA(5+n)):
    seq(a(n),n=0..23); # Peter Luschny, Dec 14 2015
    A002057List := proc(m) local A, P, n; A := [1]; P := [1,1,1];
    for n from 1 to m - 2 do P := ListTools:-PartialSums([op(P), P[-1]]);
    A := [op(A), P[-1]] od; A end: A002057List(27); # Peter Luschny, Mar 26 2022
  • Mathematica
    Table[Plus@@Flatten[Nest[ #/.a_Integer:> Range[0, a+1]&, {0}, n]], {n, 0, 10}]
    Table[4 Binomial[2n+3,n]/(n+4),{n,0,30}] (* or *) CoefficientList[ Series[ (1-Sqrt[1-4 x]+2 x (-2+Sqrt[1-4 x]+x))/(2 x^4),{x,0,30}],x] (* Harvey P. Dale, May 05 2011 *)
  • PARI
    {a(n) = if( n<0, 0, n+=2; 2*binomial(2*n, n-2) / n)}; /* Michael Somos, Jul 31 2005 */
    
  • PARI
    x='x+O('x^100); Vec((1-(1-4*x)^(1/2)+2*x*(-2+(1-4*x)^(1/2)+x))/(2*x^4)) \\ Altug Alkan, Dec 14 2015
    
  • SageMath
    [2*(n+1)*catalan_number(n+2)/(n+4) for n in (0..30)] # G. C. Greubel, May 27 2022

Formula

a(n) = A033184(n+4, 4) = 4*binomial(2*n+3, n)/(n+4) = 2*(n+1)*A000108(n+2)/(n+4).
G.f.: c(x)^4 with c(x) g.f. of A000108 (Catalan).
Row sums of A145596. Column 4 of A033184. By specializing the identities for the row polynomials given in A145596 we obtain the results a(n) = Sum_{k = 0..n} (-1)^k*binomial(n+1,k+1)*a(k)*4^(n-k) and a(n) = Sum_{k = 0..floor(n/2)} binomial(n+1,2*k+1) * Catalan(k+1) * 2^(n-2*k). From the latter identity we can derive the congruences a(2n+1) == 0 (mod 4) and a(2n) == Catalan(n+1) (mod 4). It follows that a(n) is odd if and only if n = (2^m - 4) for some m >= 2. - Peter Bala, Oct 14 2008
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=3, a(n-3) = (-1)^(n-3) * coeff(charpoly(A,x), x^3). - Milan Janjic, Jul 08 2010
G.f.: (1-sqrt(1-4*x) + 2*x*(-2+sqrt(1-4*x) + x))/(2*x^4). - Harvey P. Dale, May 05 2011
a(n+1) = A214292(2*n+4,n). - Reinhard Zumkeller, Jul 12 2012
D-finite with recurrence: (n+4)a(n) = 8*(2*n-1)*a(n-3) - 20*(n+1)*a(n-2) + 4*(2*n+5)*a(n-1). - Fung Lam, Jan 29 2014
D-finite with recurrence: (n+4)*a(n) - 2*(3*n+7)*a(n-1) + 4*(2*n+1)*a(n-2) = 0. - R. J. Mathar, Jun 03 2014
Asymptotics: a(n) ~ 4^(n+3)/sqrt(4*Pi*n^3). - Fung Lam, Mar 31 2014
a(n) = 32*4^n*Gamma(5/2+n)*(1+n)/(sqrt(Pi)*Gamma(5+n)). - Peter Luschny, Dec 14 2015
a(n) = C(n+1) - 2*C(n) where C is Catalan number A000108. Yuchun Ji, Oct 18 2017 [Note: Offset is off by 2]
E.g.f.: d/dx ( 2*exp(2*x)*BesselI(2,2*x)/x ). - Ilya Gutkovskiy, Nov 01 2017
From Bradley Klee, Mar 05 2018: (Start)
With F(x) = 16/(1+sqrt(1-4*x))^4 g.f. of A002057, xi(x) = F(x/4)*(x/4)^2, K(16*x) = 2F1(1/2,1/2;1;16*x) g.f. of A002894, q(x) g.f. of A005797, and q'(x) g.f. of A274344:
K(x) = (1+sqrt(xi(x)))*K(xi(x)).
2*K(1-x) = (1+sqrt(xi(x)))*K(1-xi(x)).
q(x) = sqrt(q(xi(16*x)/16)) = q'(xi(16*x)/16)/sqrt(xi(16*x)/16). (End)
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=0} 1/a(n) = 5/4 + Pi/(18*sqrt(3)).
Sum_{n>=0} (-1)^n/a(n) = 183*log(phi)/(25*sqrt(5)) - 77/100, where phi is the golden ratio (A001622). (End)
a(n) = Integral_{x=0..4} x^n*W(x) dx where W(x) = -x^(3/2)*(1 - x/2)*sqrt(4 - x)/Pi, defined on the open interval (0,4). - Karol A. Penson, Nov 13 2022

A109466 Riordan array (1, x(1-x)).

Original entry on oeis.org

1, 0, 1, 0, -1, 1, 0, 0, -2, 1, 0, 0, 1, -3, 1, 0, 0, 0, 3, -4, 1, 0, 0, 0, -1, 6, -5, 1, 0, 0, 0, 0, -4, 10, -6, 1, 0, 0, 0, 0, 1, -10, 15, -7, 1, 0, 0, 0, 0, 0, 5, -20, 21, -8, 1, 0, 0, 0, 0, 0, -1, 15, -35, 28, -9, 1, 0, 0, 0, 0, 0, 0, -6, 35, -56, 36, -10, 1, 0, 0, 0, 0, 0, 0, 1, -21, 70, -84, 45, -11, 1, 0, 0, 0, 0
Offset: 0

Views

Author

Philippe Deléham, Aug 28 2005

Keywords

Comments

Inverse is Riordan array (1, xc(x)) (A106566).
Triangle T(n,k), 0 <= k <= n, read by rows, given by [0, -1, 1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.
Modulo 2, this sequence gives A106344. - Philippe Deléham, Dec 18 2008
Coefficient array of the polynomials Chebyshev_U(n, sqrt(x)/2)*(sqrt(x))^n. - Paul Barry, Sep 28 2009

Examples

			Rows begin:
  1;
  0,  1;
  0, -1,  1;
  0,  0, -2,  1;
  0,  0,  1, -3,  1;
  0,  0,  0,  3, -4,   1;
  0,  0,  0, -1,  6,  -5,   1;
  0,  0,  0,  0, -4,  10,  -6,   1;
  0,  0,  0,  0,  1, -10,  15,  -7,  1;
  0,  0,  0,  0,  0,   5, -20,  21, -8,  1;
  0,  0,  0,  0,  0,  -1,  15, -35, 28, -9, 1;
From _Paul Barry_, Sep 28 2009: (Start)
Production array is
  0,    1,
  0,   -1,    1,
  0,   -1,   -1,   1,
  0,   -2,   -1,  -1,   1,
  0,   -5,   -2,  -1,  -1,  1,
  0,  -14,   -5,  -2,  -1, -1,  1,
  0,  -42,  -14,  -5,  -2, -1, -1,  1,
  0, -132,  -42, -14,  -5, -2, -1, -1,  1,
  0, -429, -132, -42, -14, -5, -2, -1, -1, 1 (End)
		

Crossrefs

Cf. A026729 (unsigned version), A000108, A030528, A124644.

Programs

  • Magma
    /* As triangle */ [[(-1)^(n-k)*Binomial(k, n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Jan 14 2016
  • Mathematica
    (* The function RiordanArray is defined in A256893. *)
    RiordanArray[1&, #(1-#)&, 13] // Flatten (* Jean-François Alcover, Jul 16 2019 *)

Formula

Number triangle T(n, k) = (-1)^(n-k)*binomial(k, n-k).
T(n, k)*2^(n-k) = A110509(n, k); T(n, k)*3^(n-k) = A110517(n, k).
Sum_{k=0..n} T(n,k)*A000108(k)=1. - Philippe Deléham, Jun 11 2007
From Philippe Deléham, Oct 30 2008: (Start)
Sum_{k=0..n} T(n,k)*A144706(k) = A082505(n+1).
Sum_{k=0..n} T(n,k)*A002450(k) = A100335(n).
Sum_{k=0..n} T(n,k)*A001906(k) = A100334(n).
Sum_{k=0..n} T(n,k)*A015565(k) = A099322(n).
Sum_{k=0..n} T(n,k)*A003462(k) = A106233(n). (End)
Sum_{k=0..n} T(n,k)*x^(n-k) = A053404(n), A015447(n), A015446(n), A015445(n), A015443(n), A015442(n), A015441(n), A015440(n), A006131(n), A006130(n), A001045(n+1), A000045(n+1), A000012(n), A010892(n), A107920(n+1), A106852(n), A106853(n), A106854(n), A145934(n), A145976(n), A145978(n), A146078(n), A146080(n), A146083(n), A146084(n) for x = -12,-11,-10,-9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12 respectively. - Philippe Deléham, Oct 27 2008
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A010892(n), A099087(n), A057083(n), A001787(n+1), A030191(n), A030192(n), A030240(n), A057084(n), A057085(n+1), A057086(n) for x = 0,1,2,3,4,5,6,7,8,9,10 respectively. - Philippe Deléham, Oct 28 2008
G.f.: 1/(1-y*x+y*x^2). - Philippe Deléham, Dec 15 2011
T(n,k) = T(n-1,k-1) - T(n-2,k-1), T(n,0) = 0^n. - Philippe Deléham, Feb 15 2012
Sum_{k=0..n} T(n,k)*x^(n-k) = F(n+1,-x) where F(n,x)is the n-th Fibonacci polynomial in x defined in A011973. - Philippe Deléham, Feb 22 2013
Sum_{k=0..n} T(n,k)^2 = A051286(n). - Philippe Deléham, Feb 26 2013
Sum_{k=0..n} T(n,k)*T(n+1,k) = -A110320(n). - Philippe Deléham, Feb 26 2013
For T(0,0) = 0, the signed triangle below has the o.g.f. G(x,t) = [t*x(1-x)]/[1-t*x(1-x)] = L[t*Cinv(x)] where L(x) = x/(1-x) and Cinv(x)=x(1-x) with the inverses Linv(x) = x/(1+x) and C(x)= [1-sqrt(1-4*x)]/2, an o.g.f. for the shifted Catalan numbers A000108, so the inverse o.g.f. is Ginv(x,t) = C[Linv(x)/t] = [1-sqrt[1-4*x/(t(1+x))]]/2 (cf. A124644 and A030528). - Tom Copeland, Jan 19 2016

A003517 Number of permutations of [n+1] with exactly 1 increasing subsequence of length 3.

Original entry on oeis.org

1, 6, 27, 110, 429, 1638, 6188, 23256, 87210, 326876, 1225785, 4601610, 17298645, 65132550, 245642760, 927983760, 3511574910, 13309856820, 50528160150, 192113383644, 731508653106, 2789279908316, 10649977831752, 40715807302800, 155851062397940, 597261490737912
Offset: 2

Views

Author

Keywords

Comments

a(n-4) = number of n-th generation vertices in the tree of sequences with unit increase labeled by 5 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of standard tableaux of shape (n+3,n-2). - Emeric Deutsch, May 30 2004
a(n) = A214292(2*n,n-3) for n > 2. - Reinhard Zumkeller, Jul 12 2012
a(n) is the number of North-East paths from (0,0) to (n,n) that cross the diagonal y = x horizontally exactly once. By symmetry, it is also the number of North-East paths from (0,0) to (n,n) that cross the diagonal y = x vertically exactly once. Details can be found in Section 3.3 in Pan and Remmel's link. - Ran Pan, Feb 02 2016
a(n) is the number of permutations pi of [n+3] such that s(pi)=321456...(n+3), where s denotes West's stack-sorting map. - Colin Defant, Jan 14 2019
a(n) is also the number of permutations of [n+1] avoiding the pattern 321. For permutations avoiding any of the other permutations of [3] (that is, any of 132, 213, 231, or 312) see A002054. - N. J. A. Sloane, Nov 26 2022

Examples

			a(3)=6 because the only permutations of 1234 with exactly 1 increasing subsequence of length 3 are 1423, 4123, 1342, 2314, 2341, 3124.
		

References

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

Crossrefs

T(n, n+6) for n=0, 1, 2, ..., array T as in A047072.
See also A002054.
First differences are in A026017.
A diagonal of any of the essentially equivalent arrays: A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Maple
    A003517List := proc(m) local A, P, n; A := [1]; P := [1,1,1,1,1];
    for n from 1 to m - 2 do P := ListTools:-PartialSums([op(P), P[-1]]);
    A := [op(A), P[-1]] od; A end: A003517List(25); # Peter Luschny, Mar 26 2022
  • Mathematica
    f[x_] = (Sqrt[1 - 4 x] - 1)^6/(64 x^4); CoefficientList[Series[f[x], {x, 0, 25}], x][[3 ;; 26]] (* Jean-François Alcover, Jul 13 2011, after g.f. *)
    Table[6 Binomial[2n+1,n-2]/(n+4),{n,2,30}] (* Harvey P. Dale, Feb 27 2012 *)
  • PARI
    a(n)=6*binomial(2*n+1,n-2)/(n+4) \\ Charles R Greathouse IV, May 18 2015
    
  • PARI
    x='x+O('x^50); Vec(x^2*((1-(1-4*x)^(1/2))/(2*x))^6) \\ Altug Alkan, Nov 01 2015

Formula

a(n) = 6*binomial(2*n+1, n-2)/(n+4).
G.f.: x^2*C(x)^6, where C(x) is g.f. for the Catalan numbers (A000108). - Emeric Deutsch, May 30 2004
E.g.f.: exp(2*x)*(Bessel_I(2,2*x) - Bessel_I(4,2*x)). - Paul Barry, Jun 04 2007
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n >= 5, a(n-3) = (-1)^(n-5)*coeff(charpoly(A,x),x^5). - Milan Janjic, Jul 08 2010
a(n) = Sum_{i>=1, j>=1, k>=1, i+j+k=n+1} Catalan(i)*Catalan(j)*Catalan(k). T. D. Noe, Dec 22 2010
D-finite with recurrence -(n+4)*(n-2)*a(n) + 2*n*(2*n+1)*a(n-1) = 0. - R. J. Mathar, Dec 04 2012
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=2} 1/a(n) = 7/2 - 34*Pi/(27*sqrt(3)).
Sum_{n>=2} (-1)^n/a(n) = 828*log(phi)/(25*sqrt(5)) - 2819/450, where phi is the golden ratio (A001622). (End)
a(n) ~ 3*4^(n+1)/(n^(3/2)*sqrt(Pi)). - Stefano Spezia, Apr 17 2024
a(n) = A000108(n+3) - 4*A000108(n+2) + 3*A000108(n+1). - Taras Goy, Jul 15 2024
a(n) = 6*(2*n+1)!*(n-1)!/((2*n-4)!*(n+4)!)*A000108(n-2). - Taras Goy, Dec 21 2024

A026729 Square array of binomial coefficients T(n,k) = binomial(n,k), n >= 0, k >= 0, read by downward antidiagonals.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 2, 1, 0, 0, 1, 3, 1, 0, 0, 0, 3, 4, 1, 0, 0, 0, 1, 6, 5, 1, 0, 0, 0, 0, 4, 10, 6, 1, 0, 0, 0, 0, 1, 10, 15, 7, 1, 0, 0, 0, 0, 0, 5, 20, 21, 8, 1, 0, 0, 0, 0, 0, 1, 15, 35, 28, 9, 1, 0, 0, 0, 0, 0, 0, 6, 35, 56, 36, 10, 1, 0, 0, 0, 0, 0, 0, 1, 21, 70, 84, 45, 11, 1, 0, 0, 0, 0
Offset: 0

Views

Author

N. J. A. Sloane, Jan 19 2003

Keywords

Comments

The signed triangular matrix T(n,k)*(-1)^(n-k) is the inverse matrix of the triangular Catalan convolution matrix A106566(n,k), n=k>=0, with A106566(n,k) = 0 if nPhilippe Deléham, Aug 01 2005
As a number triangle: unsigned version of A109466. - Philippe Deléham, Oct 26 2008
A063967*A130595 as infinite lower triangular matrices. - Philippe Deléham, Dec 11 2008
Modulo 2, this sequence becomes A106344. - Philippe Deléham, Dec 18 2008
Let {a_(k,i)}, k>=1, i=0,...,k, be the k-th antidiagonal of the array. Then s_k(n) = Sum_{i=0..k}a_(k,i)* binomial(n,k) is the n-th element of the k-th column of A111808. For example, s_1(n) = binomial(n,1) = n is the first column of A111808 for n>1, s_2(n) = binomial(n,1) + binomial(n,2) is the second column of A111808 for n>1, etc. Therefore, in cases k=3,4,5,6,7,8, s_k(n) is A005581(n), A005712(n), A000574(n), A005714(n), A005715(n), A005716(n), respectively. Besides, s_k(n+5) = A064054(n). - Vladimir Shevelev and Peter J. C. Moses, Jun 22 2012
As a triangle, T(n,k) = binomial(k,n-k). - Peter Bala, Nov 27 2015
For all n >= 0, k >= 0, the k-th homology group of the n-torus H_k(T^n) is the free abelian group of rank T(n,k) = binomial(n,k). See the Math Stack Exchange link below. - Jianing Song, Mar 13 2023

Examples

			Array begins
  1 0 0 0 0 0 ...
  1 1 0 0 0 0 ...
  1 2 1 0 0 0 ...
  1 3 3 1 0 0 ...
  1 4 6 4 1 0 ...
As a triangle, this begins
  1
  0 1
  0 1 1
  0 0 2 1
  0 0 1 3 1
  0 0 0 3 4 1
  0 0 0 1 6 5 1
  ...
Production array is
  0    1
  0    1   1
  0   -1   1   1
  0    2  -1   1  1
  0   -5   2  -1  1  1
  0   14  -5   2 -1  1  1
  0  -42  14  -5  2 -1  1  1
  0  132 -42  14 -5  2 -1  1  1
  0 -429 132 -42 14 -5  2 -1  1  1
  ... (Cf. A000108)
		

Crossrefs

The official entry for Pascal's triangle is A007318. See also A052553 (the same array read by upward antidiagonals).
Cf. A030528 (subtriangle for 1<=k<=n).

Programs

  • GAP
    nmax:=15;; T:=List([0..nmax],n->List([0..nmax],k->Binomial(n,k)));;
    b:=List([2..nmax],n->OrderedPartitions(n,2));;
    a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][1]][b[i][j][2]]))); # Muniru A Asiru, Jul 17 2018
  • Magma
    /* As triangle: */ [[Binomial(k, n-k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Nov 29 2015
    
  • Maple
    seq(seq(binomial(k,n-k),k=0..n),n=0..12); # Peter Luschny, May 31 2014
  • Mathematica
    Table[Binomial[k, n - k], {n, 0, 12}, {k, 0, n}] // Flatten (* Michael De Vlieger, Nov 28 2015 *)

Formula

As a number triangle, this is defined by T(n,0) = 0^n, T(0,k) = 0^k, T(n,k) = T(n-1,k-1) + Sum_{j, j>=0} (-1)^j*T(n-1,k+j)*A000108(j) for n>0 and k>0. - Philippe Deléham, Nov 07 2005
As a triangle read by rows, it is [0, 1, -1, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 22 2006
As a number triangle, this is defined by T(n, k) = Sum_{i=0..n} (-1)^(n+i)*binomial(n, i)*binomial(i+k, i-k) and is the Riordan array ( 1, x*(1+x) ). The row sums of this triangle are F(n+1). - Paul Barry, Jun 21 2004
Sum_{k=0..n} x^k*T(n,k) = A000007(n), A000045(n+1), A002605(n), A030195(n+1), A057087(n), A057088(n), A057089(n), A057090(n), A057091(n), A057092(n), A057093(n) for n=0,1,2,3,4,5,6,7,8,9,10. - Philippe Deléham, Oct 16 2006
T(n,k) = A109466(n,k)*(-1)^(n-k). - Philippe Deléham, Dec 11 2008
G.f. for the triangular interpretation: -1/(-1+x*y+x^2*y). - R. J. Mathar, Aug 11 2015
For T(0,0) = 0, the triangle below has the o.g.f. G(x,t) = [t*x(1+x)]/[1-t*x(1+x)]. See A109466 for a signed version and inverse, A030528 for reverse and A102426 for a shifted version. - Tom Copeland, Jan 19 2016

A000344 a(n) = 5*binomial(2n, n-2)/(n+3).

Original entry on oeis.org

1, 5, 20, 75, 275, 1001, 3640, 13260, 48450, 177650, 653752, 2414425, 8947575, 33266625, 124062000, 463991880, 1739969550, 6541168950, 24647883000, 93078189750, 352207870014, 1335293573130, 5071418015120, 19293438101000, 73514652074500, 280531912316292
Offset: 2

Views

Author

Keywords

Comments

a(n-3) is the number of n-th generation vertices in the tree of sequences with unit increase labeled by 4 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of lattice paths from (0,0) to (n,n) with steps E=(1,0) and N=(0,1) which touch but do not cross the line x-y=2. Example: For n=3 there are the 5 paths EENENN, EENNEN, EENNNE, ENEENN, NEEENN. - Herbert Kociemba, May 24 2004
Number of standard tableaux of shape (n+2,n-2). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^2 + 5*x^3 + 20*x^4 + 75*x^5 + 275*x^6 + 1001*x^7 + 3640*x^8 + ...
		

References

  • C. Krishnamachary and M. Bheemasena Rao, Determinants whose elements are Eulerian, prepared Bernoullian and other numbers, J. Indian Math. Soc., Vol. 14 (1922), pp. 55-62, 122-138 and 143-146.
  • 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

T(n, n+5) for n=0, 1, 2, ..., array T as in A047072.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • GAP
    List([2..30],n->5*Binomial(2*n,n-2)/(n+3)); # Muniru A Asiru, Aug 09 2018
  • Magma
    [5*Binomial(2*n,n-2)/(n+3): n in [2..30]]; // Vincenzo Librandi, May 03 2011
    
  • Maple
    A000344List := proc(m) local A, P, n; A := [1]; P := [1,1,1,1];
    for n from 1 to m - 2 do P := ListTools:-PartialSums([op(P), P[-1]]);
    A := [op(A), P[-1]] od; A end: A000344List(27); # Peter Luschny, Mar 26 2022
  • Mathematica
    Table[5 Binomial[2n,n-2]/(n+3),{n,2,40}] (* or *) CoefficientList[Series[ (1-Sqrt[1-4 x]+x (-5+3 Sqrt[1-4 x]-(-5+Sqrt[1-4 x]) x))/(2 x^5), {x,0,38}],x]  (* Harvey P. Dale, May 01 2011 *)
    a[ n_] := If[ n < 0, 0, 5 Binomial[2 n, n - 2] / (n + 3)]; (* Michael Somos, May 28 2014 *)
  • PARI
    a(n)=5*binomial(2*n,n-2)/(n+3) \\ Charles R Greathouse IV, Jul 25 2011
    

Formula

Integral representation as n-th moment of a function on [0, 4]: a(n) = Integral_{x=0..4} x^n*((1/2)/Pi*x^(3/2)*(x^2-3*x+1)*(4-x)^(1/2)) dx, n >= 0, for which offset=0. - Karol A. Penson, Oct 11 2001
Expansion of x^2*C^5, where C = (1-(1-4*x)^(1/2))/(2*x) is g.f. for Catalan numbers (A000108). - Herbert Kociemba, May 02 2004
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=4, a(n-2)=(-1)^(n-4)*coeff(charpoly(A,x),x^4). - Milan Janjic, Jul 08 2010
a(n) = A000108(n+2) - 3*A000108(n+1)+ A000108(n). - David Scambler, May 20 2012
D-finite with recurrence: (n+3)*(n-2)*a(n) = 2*n*(2n-1)*a(n-1). - R. J. Mathar, Jun 27 2012
a(n) = A214292(2*n-1,n-3) for n > 2. - Reinhard Zumkeller, Jul 12 2012
0 = a(n)*(-528*a(n+1) + 9162*a(n+2) - 9295*a(n+3) + 1859*a(n+4)) + a(n+1)*(-1650*a(n+1) - 762*a(n+2) + 4188*a(n+3) - 946*a(n+4)) + a(n+2)*(-1050*a(n+2) - 126*a(n+3) + 84*a(n+4)) for all n in Z. - Michael Somos, May 28 2014
0 = a(n)*(a(n)*(+16*a(n+1) + 6*a(n+2)) + a(n+1)*(+66*a(n+1) - 105*a(n+2) + 40*a(n+3)) + a(n+2)*(-69*a(n+2) + 15*a(n+3))) +a(n+1)*(a(n+1)*(50*a(n+1) + 42*a(n+2) - 28*a(n+3)) +a(n+2)*(+12*a(n+2))) for all n in Z. - Michael Somos, May 28 2014
0 = a(n)^2*(-16*a(n+1)^2 - 38*a(n+1)*a(n+2) - 12*a(n+2)^2) + a(n)*a(n+1)*(-66*a(n+1)^2 + 149*a(n+1)*a(n+2) - 23*a(n+2)^2) + a(n+1)^2*(-50*a(n+1)^2 + 2*a(n+2)^2) for all n in Z. - Michael Somos, May 28 2014
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: (x*(2 + x) * BesselI(0, 2*x) - (2+x+x^2) * BesselI(1, 2*x)) * exp(2*x)/x^2.
a(n) ~ 5*4^n/(sqrt(Pi)*n^(3/2)). (End)
a(n) = (1/(n+1))*Sum_{i=0..n-2} (-1)^(n+i)*(n-i+1)*binomial(2n+2,i), n >= 2. - Taras Goy, Aug 09 2018
G.f.: x^2* 2F1(5/2,3;6;4*x) . - R. J. Mathar, Jan 27 2020
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=2} 1/a(n) = 14/5 - 38*Pi/(45*sqrt(3)).
Sum_{n>=2} (-1)^n/a(n) = 1956*log(phi)/(125*sqrt(5)) - 316/125, where phi is the golden ratio (A001622). (End)
a(n) = 5*(2*n)!*(n-1)!/((2*n-4)!*(n+3)!)*A000108(n-2). - Taras Goy, Jul 15 2024
a(n) = Sum_{i+j+k+l+m = n-2} C(i)C(j)C(k)C(l)C(m), where C(s) = A000108(s). (Fifth convolution of Catalan numbers). - Taras Goy, Dec 21 2024

A059365 Another version of the Catalan triangle: T(r,s) = binomial(2*r-s-1,r-1) - binomial(2*r-s-1,r), r>=0, 0 <= s <= r.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 5, 5, 3, 1, 0, 14, 14, 9, 4, 1, 0, 42, 42, 28, 14, 5, 1, 0, 132, 132, 90, 48, 20, 6, 1, 0, 429, 429, 297, 165, 75, 27, 7, 1, 0, 1430, 1430, 1001, 572, 275, 110, 35, 8, 1, 0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44
Offset: 0

Views

Author

N. J. A. Sloane, Jan 28 2001

Keywords

Examples

			Triangle starts
  0;
  0,    1;
  0,    1,    1;
  0,    2,    2,    1;
  0,    5,    5,    3,    1;
  0,   14,   14,    9,    4,    1;
  0,   42,   42,   28,   14,    5,   1;
  0,  132,  132,   90,   48,   20,   6,   1;
  0,  429,  429,  297,  165,   75,  27,   7,  1;
  0, 1430, 1430, 1001,  572,  275, 110,  35,  8, 1;
  0, 4862, 4862, 3432, 2002, 1001, 429, 154, 44, 9, 1;
  ...
		

Crossrefs

See also the triangle in A009766. First 2 diagonals both give A000108, next give A000245, A002057.
The three triangles A059365, A106566 and A099039 are the same except for signs and the leading term.
Essentially the same as A033184.
The following are all versions of (essentially) the same Catalan triangle: A009766, A030237, A033184, A053121, A059365, A099039, A106566, A130020, A047072, A171567, A181645.

Programs

  • Magma
    /* as triangle */ [[[0] cat [Binomial(2*r-s-1, r-1)- Binomial(2*r-s-1, r): s in [1..r]]: r in [0..10]]]; // Vincenzo Librandi, Jan 09 2017
  • Mathematica
    Table[Binomial[2*r - s - 1, r - 1] - Binomial[2*r - s - 1, r], {r, 0, 10}, {s, 0, r}] // Flatten (* G. C. Greubel, Jan 08 2017 *)
  • PARI
    tabl(nn) = { print(0); for (r=1, nn, for (s=0, r, print1(binomial(2*r-s-1,r-1)-binomial(2*r-s-1,r), ", ");); print(););}  \\ Michel Marcus, Nov 01 2013
    

Formula

Essentially the same triangle as [0, 1, 1, 1, 1, 1, 1, ...] DELTA A000007, where DELTA is Deléham's operator defined in A084938, but the first term is T(0,0) = 0.

A000588 a(n) = 7*binomial(2n,n-3)/(n+4).

Original entry on oeis.org

0, 0, 0, 1, 7, 35, 154, 637, 2548, 9996, 38760, 149226, 572033, 2187185, 8351070, 31865925, 121580760, 463991880, 1771605360, 6768687870, 25880277150, 99035193894, 379300783092, 1453986335186, 5578559816632, 21422369201800, 82336410323440, 316729578421620
Offset: 0

Views

Author

Keywords

Comments

a(n-5) is the number of n-th generation vertices in the tree of sequences with unit increase labeled by 6 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of lattice paths from (0,0) to (n,n) with steps E=(1,0) and N=(0,1) which touch but do not cross the line x-y=3. Example: For n=3 there is only one path EEENNN. - Herbert Kociemba, May 24 2004
Number of standard tableaux of shape (n+3,n-3). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^3 + 7*x^4 + 35*x^5 + 154*x^6 + 637*x^7 + 2548*x^8 + 9996*x^9 + ...
		

References

  • 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

First differences are in A026014.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Mathematica
    a[n_] := 7*Binomial[2n, n-3]/(n + 4); Table[a[n],{n,0,27}] (* James C. McMahon, Dec 05 2023 *)
  • PARI
    A000588(n)=7*binomial(2*n,n-3)/(n+4) \\ M. F. Hasler, Aug 25 2012
    
  • PARI
    my(x='x+O('x^50)); concat([0, 0, 0], Vec(x^3*((1-(1-4*x)^(1/2))/(2*x))^7)) \\ Altug Alkan, Nov 01 2015

Formula

Expansion of x^3*C^7, where C = (1-(1-4*x)^(1/2))/(2*x) is the g.f. for the Catalan numbers, A000108. - Philippe Deléham, Feb 03 2004
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=6, a(n-3)=(-1)^(n-6)*coeff(charpoly(A,x),x^6). - Milan Janjic, Jul 08 2010
a(n) = A214292(2*n-1,n-4) for n > 3. - Reinhard Zumkeller, Jul 12 2012
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: (1/6)*x^3*1F1(7/2; 8; 4*x).
a(n) ~ 7*4^n/(sqrt(Pi)*n^(3/2)). (End)
0 = a(n)*(+1456*a(n+1) - 87310*a(n+2) + 132834*a(n+3) - 68068*a(n+4) + 9724*a(n+5)) + a(n+1)*(+8918*a(n+1) - 39623*a(n+2) + 51726*a(n+3) - 299*a(n+4) - 1573*a(n+5)) + a(n+2)*(-24696*a(n+2) - 1512*a(n+3) + 1008*a(n+4)) for all n in Z. - Michael Somos, Jan 22 2017
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=3} 1/a(n) = 27/14 - 26*Pi/(63*sqrt(3)).
Sum_{n>=3} (-1)^(n+1)/a(n) = 11364*log(phi)/(175*sqrt(5)) - 4583/350, where phi is the golden ratio (A001622). (End)
a(n) = Integral_{x=0..4} x^(n)*W(x)dx, n>=0, where W(x) = sqrt(4/x - 1)*(x^3 - 5*x^2 + 6*x - 1)/(2*Pi). The function W(x) for x->0 tends to -infinity (which is its absolute minimum), and W(4) = 0. W(x) is a signed function on the interval x = (0, 4) where it has two maxima separated by one local minimum. - Karol A. Penson, Jun 17 2024
D-finite with recurrence -(n+4)*(n-3)*a(n) +2*n*(2*n-1)*a(n-1)=0. - R. J. Mathar, Jul 30 2024
a(n) = A000108(n+3) - 5*A000108(n+2) + 6*A000108(n+1) - A000108(n). - Taras Goy, Dec 21 2024

Extensions

More terms from N. J. A. Sloane, Jul 13 2010

A015565 a(n) = 7*a(n-1) + 8*a(n-2), a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 7, 57, 455, 3641, 29127, 233017, 1864135, 14913081, 119304647, 954437177, 7635497415, 61083979321, 488671834567, 3909374676537, 31274997412295, 250199979298361, 2001599834386887, 16012798675095097, 128102389400760775, 1024819115206086201, 8198552921648689607
Offset: 0

Views

Author

Keywords

Comments

A linear 2nd order recurrence. A Jacobsthal number sequence.
Binomial transform of A053573 (preceded by zero). - Paul Barry, Apr 09 2003
Second binomial transform of A080424. Binomial transform of A053573, with leading zero. Binomial transform is 0,1,9,81,729,....(9^n - 0^n)/9. Second binomial transform is 0,1,11,111,1111,... (A002275: repunits). - Paul Barry, Mar 14 2004
Number of walks of length n between any two distinct nodes of the complete graph K_9. Example: a(2)=7 because the walks of length 2 between the nodes A and B of the complete graph ABCDEFGHI are: ACB, ADB, AEB, AFB, AGB, AHB and AIB. - Emeric Deutsch, Apr 01 2004
Unsigned version of A014990. - Philippe Deléham, Feb 13 2007
The ratio a(n+1)/a(n) converges to 8 as n approaches infinity. - Felix P. Muga II, Mar 09 2014

Examples

			G.f. = x + 7*x^2 + 57*x^3 + 455*x^4 + 3641*x^5 + 29127*x^6 + 233017*x^7 + ...
		

Crossrefs

Programs

Formula

From Paul Barry, Apr 09 2003: (Start)
a(n) = (8^n - (-1)^n)/9.
a(n) = J(3*n)/3 = A001045(3*n)/3. (End)
From Emeric Deutsch, Apr 01 2004: (Start)
a(n) = 8^(n-1) - a(n-1).
G.f.: x/(1-7*x-8*x^2). (End)
a(n) = Sum_{k = 0..n} A106566(n,k)*A099322(k). - Philippe Deléham, Oct 30 2008
a(n) = round(8^n/9). - Mircea Merca, Dec 28 2010
From Peter Bala, May 31 2024: (Start)
G.f: A(x) = x/(1 - x^2) o x/(1 - x^2), where o denotes the black diamond product of power series as defined by Dukes and White. Cf. A054878.
The black diamond product A(x) o A(x) is the g.f. for the number of walks of length n between any two distinct nodes of the complete graph K_81.
Row 8 of A062160. (End)
E.g.f.: exp(-x)*(exp(9*x) - 1)/9. - Elmo R. Oliveira, Aug 17 2024

A001392 a(n) = 9*binomial(2n,n-4)/(n+5).

Original entry on oeis.org

1, 9, 54, 273, 1260, 5508, 23256, 95931, 389367, 1562275, 6216210, 24582285, 96768360, 379629720, 1485507600, 5801732460, 22626756594, 88152205554, 343176898988, 1335293573130, 5193831553416, 20198233818840, 78542105700240, 305417807763705
Offset: 4

Views

Author

Keywords

Comments

Number of n-th generation vertices in the tree of sequences with unit increase labeled by 8 (cf. Zoran Sunic reference) - Benoit Cloitre, Oct 07 2003
Number of lattice paths from (0,0) to (n,n) with steps E=(1,0) and N=(0,1) which touch but do not cross the line x-y=4. - Herbert Kociemba, May 24 2004
Number of standard tableaux of shape (n+4,n-4). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^4 + 9*x^5 + 54*x^6 + 273*x^7 + 1260*x^8 + 5508*x^9 + 23256*x^10 + ...
		

References

  • 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

First differences are in A026015.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

Formula

Expansion of x^4*C^9, where C = (1-(1-4*x)^(1/2))/(2*x) is g.f. for Catalan numbers, A000108. - Philippe Deléham, Feb 03 2004
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=8, a(n-4)=(-1)^(n-8)*coeff(charpoly(A,x),x^8). - Milan Janjic, Jul 08 2010
a(n) = A214292(2*n-1,n-5) for n > 4. - Reinhard Zumkeller, Jul 12 2012
D-finite with recurrence -(n+5)*(n-4)*a(n) +2*n*(2*n-1)*a(n-1)=0. - R. J. Mathar, Jun 20 2013
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: (1/24)*x^4*1F1(9/2; 10; 4*x).
a(n) ~ 9*4^n/(sqrt(Pi)*n^(3/2)). (End)
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=4} 1/a(n) = 158*Pi/(81*sqrt(3)) - 649/270.
Sum_{n>=4} (-1)^n/a(n) = 52076*log(phi)/(225*sqrt(5)) - 22007/450, where phi is the golden ratio (A001622). (End)

Extensions

More terms from Harvey P. Dale, Mar 03 2011

A003518 a(n) = 8*binomial(2*n+1,n-3)/(n+5).

Original entry on oeis.org

1, 8, 44, 208, 910, 3808, 15504, 62016, 245157, 961400, 3749460, 14567280, 56448210, 218349120, 843621600, 3257112960, 12570420330, 48507033744, 187187399448, 722477682080, 2789279908316, 10772391370048, 41620603020640, 160878516023680, 622147386185325
Offset: 3

Views

Author

Keywords

Comments

a(n-6) is the number of n-th generation nodes in the tree of sequences with unit increase labeled by 7 (cf. Zoran Sunic reference). - Benoit Cloitre, Oct 07 2003
Number of standard tableaux of shape (n+4,n-3). - Emeric Deutsch, May 30 2004

Examples

			G.f. = x^3 + 8*x^4 + 44*x^5 + 208*x^6 + 910*x^7 + 3808*x^8 + 15504*x^9 + ...
		

References

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

Crossrefs

Cf. A002057.
First differences are in A026018.
A diagonal of any of the essentially equivalent arrays A009766, A030237, A033184, A059365, A099039, A106566, A130020, A047072.

Programs

  • Magma
    [8*Binomial(2*n+1,n-3)/(n+5): n in [3..30]]; // Vincenzo Librandi, Jan 23 2017
  • Mathematica
    Table[8 Binomial[2 n + 1, n - 3]/(n + 5), {n, 3, 25}] (* Michael De Vlieger, Oct 26 2016 *)
    CoefficientList[Series[((1 - Sqrt[1 - 4 x])/(2 x))^8, {x, 0, 30}], x] (* Vincenzo Librandi, Jan 23 2017 *)
  • PARI
    {a(n) = if( n<3, 0, 8 * binomial(2*n + 1, n-3) / (n + 5))}; /* Michael Somos, Mar 14 2011 */
    
  • PARI
    my(x='x+O('x^50)); Vec(x^3*((1-(1-4*x)^(1/2))/(2*x))^8) \\ Altug Alkan, Nov 01 2015
    

Formula

G.f.: x^3*C(x)^8, where C(x)=(1-sqrt(1-4*x))/(2*x) is g.f. for the Catalan numbers (A000108). - Emeric Deutsch, May 30 2004
The convolution of A002057 with itself. - Gerald McGarvey, Nov 08 2007
Let A be the Toeplitz matrix of order n defined by: A[i,i-1]=-1, A[i,j]=Catalan(j-i), (i<=j), and A[i,j]=0, otherwise. Then, for n>=7, a(n-4)=(-1)^(n-7)*coeff(charpoly(A,x),x^7). - Milan Janjic, Jul 08 2010
a(n) = A214292(2*n,n-4) for n > 3. - Reinhard Zumkeller, Jul 12 2012
Integral representation as the n-th moment of the signed weight function W(x) on (0,4), i.e.: a(n+3) = Integral_{x=0..4} x^n*W(x) dx, n >= 0, with W(x) = (1/2)*x^(7/2)*(x-2)*(x^2-4*x+2)*sqrt(4-x)/Pi. - Karol A. Penson, Oct 26 2016
From Ilya Gutkovskiy, Jan 22 2017: (Start)
E.g.f.: 4*BesselI(4,2*x)*exp(2*x)/x.
a(n) ~ 4^(n+2)/(sqrt(Pi)*n^(3/2)). (End)
D-finite with recurrence: -(n+5)*(n-3)*a(n) +2*n*(2*n+1)*a(n-1)=0. - R. J. Mathar, Feb 20 2020
From Amiram Eldar, Jan 02 2022: (Start)
Sum_{n>=3} 1/a(n) = 43*Pi/(36*sqrt(3)) - 81/80.
Sum_{n>=3} (-1)^(n+1)/a(n) = 6213*log(phi)/(50*sqrt(5)) - 10339/400, where phi is the golden ratio (A001622). (End)

Extensions

More terms from Jon E. Schoenfield, May 06 2010
Previous Showing 11-20 of 58 results. Next