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

A059090 Triangle T(n,m) giving number of m-element intersecting antichains on a labeled n-set or n-variable Boolean functions with m nonzero values in the Post class F(7,2), m=0,.., A037952(n).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 7, 3, 1, 1, 15, 30, 30, 5, 1, 31, 195, 605, 780, 543, 300, 135, 45, 10, 1, 1, 63, 1050, 9030, 41545, 118629, 233821, 329205, 327915, 224280, 100716, 29337, 5950, 910, 105, 7
Offset: 0

Views

Author

Vladeta Jovovic, Goran Kilibarda, Dec 28 2000

Keywords

Comments

An antichain is called intersecting (or proper) antichain if every two members have a nonempty intersection. Row sums give the number of intersecting antichains on a labeled n-set or n-variable Boolean functions in the Post class F(7,2) or self-dual monotone Boolean functions of n+1 variables. Cf. A001206.

Examples

			1;
1, 1;
1, 3;
1, 7, 3, 1;
1, 15, 30, 30, 5;
1, 31, 195, 605, 780, 543, 300, 135, 45, 10, 1;
1, 63, 1050, 9030, 41545, 118629, 233821, 329205, 327915, 224280, 100716, 29337, 5950, 910, 105, 7;
		

References

  • Jovovic V., Kilibarda G., The number of n-variable Boolean functions in the Post class F(7,2), Belgrade, 2001, in preparation.
  • Pogosyan G., Miyakawa M., Nozaki A., Rosenberg I., The Number of Clique Boolean Functions, IEICE Trans. Fundamentals, Vol. E80-A, No. 8, pp. 1502-1507, 1997/8.

Crossrefs

Formula

T(n, 0)=1, T(n, 1)=2^n-1, T(n, 2)=A032263(n), T(n, 3)=A051303(n), T(n, 4)=A051304(n), T(n, 5)=A051305(n), T(n, 6)=A051306(n), T(n, 7)=A051307(n).

A348012 Invert transform of A037952.

Original entry on oeis.org

1, 1, 2, 6, 15, 41, 107, 287, 758, 2022, 5362, 14276, 37916, 100874, 268084, 713006, 1895415, 5040433, 13400863, 35634419, 94745705, 251932417, 669862581, 1781165785, 4736000199, 12592950131, 33484024019, 89033206579, 236735594539
Offset: 0

Views

Author

R. J. Mathar, Sep 24 2021

Keywords

Comments

Essentially the number of paths starting at the origin, consisting of n U or D steps or starting at the horizontal axis afresh after any number of steps, and never descending below the horizontal axis.

Crossrefs

Row sums of A348013.

Formula

G.f.: 1/(1-g037952(x)) where g037952(x) = x +x^2 + 3*x^3+.. is the g.f. of A037952.
D-finite with recurrence n*a(n) +(-3*n+2)*a(n-1) +(-3*n+8)*a(n-2) +(11*n-46)*a(n-3) +2*(-n+12)*a(n-4) +4*(n-1)*a(n-5) +8*(-n+5)*a(n-6)=0.

A001405 a(n) = binomial(n, floor(n/2)).

Original entry on oeis.org

1, 1, 2, 3, 6, 10, 20, 35, 70, 126, 252, 462, 924, 1716, 3432, 6435, 12870, 24310, 48620, 92378, 184756, 352716, 705432, 1352078, 2704156, 5200300, 10400600, 20058300, 40116600, 77558760, 155117520, 300540195, 601080390, 1166803110
Offset: 0

Views

Author

Keywords

Comments

Sperner's theorem says that this is the maximal number of subsets of an n-set such that no one contains another.
When computed from index -1, [seq(binomial(n,floor(n/2)), n = -1..30)]; -> [1,1,1,2,3,6,10,20,35,70,126,...] and convolved with aerated Catalan numbers [seq(((n+1) mod 2)*binomial(n,n/2)/((n/2)+1), n = 0..30)]; -> [1,0,1,0,2,0,5,0,14,0,42,0,132,0,...] shifts left by one: [1,1,2,3,6,10,20,35,70,126,252,...] and if again convolved with aerated Catalan numbers, gives A037952 apart from the initial term. - Antti Karttunen, Jun 05 2001 [This is correct because the g.f.'s satisfy (1+x*g001405(x))*g126120(x) = g001405(x) and g001405(x)*g126120(x) = g037952(x)/x. - R. J. Mathar, Sep 23 2021]
Number of ordered trees with n+1 edges, having nonroot nodes of outdegree 0 or 2. - Emeric Deutsch, Aug 02 2002
Gives for n >= 1 the maximum absolute column sum norm of the inverse of the Vandermonde matrix (a_ij) i=0..n-1, j=0..n-1 with a_00=1 and a_ij=i^j for (i,j) != (0,0). - Torsten Muetze, Feb 06 2004
Image of Catalan numbers A000108 under the Riordan array (1/(1-2x),-x/(1-2x)) or A065109. - Paul Barry, Jan 27 2005
Number of left factors of Dyck paths, consisting of n steps. Example: a(4)=6 because we have UDUD, UDUU, UUDD, UUDU, UUUD and UUUU, where U=(1,1) and D=(1,-1). - Emeric Deutsch, Apr 23 2005
Number of dispersed Dyck paths of length n; they are defined as concatenations of Dyck paths and (1,0)-steps on the x-axis; equivalently, Motzkin paths with no (1,0)-steps at positive height. Example: a(4)=6 because we have HHHH, HHUD, HUDH, UDHH, UDUD, and UUDD, where U=(1,1), H=(1,0), and D=(1,-1). - Emeric Deutsch, Jun 04 2011
a(n) is odd iff n=2^k-1. - Jon Perry, May 05 2005
An inverse Chebyshev transform of binomial(1,n)=(1,1,0,0,0,...) where g(x)->(1/sqrt(1-4*x^2))*g(x*c(x^2)), with c(x) the g.f. of A000108. - Paul Barry, May 13 2005
In a random walk on the number line, starting at 0 and with 0 absorbing after the first step, number of ways of ending up at a positive integer after n steps. - Joshua Zucker, Jul 31 2005
Maximum number of sums of the form Sum_{i=1..n} e(i)*a(i) that are congruent to 0 mod q, where e_i=0 or 1 and gcd(a_i,q)=1, provided that q > ceiling(n/2). - Ralf Stephan, Apr 27 2003
Also the number of standard tableaux of height <= 2. - Mike Zabrocki, Mar 24 2007
Hankel transform of this sequence forms A000012 = [1,1,1,1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
A001263 * [1, -2, 3, -4, 5, ...] = [1, -1, -2, 3, 6, -10, -20, 35, 70, -126, ...]. - Gary W. Adamson, Jan 02 2008
Equals right border of triangle A153585. - Gary W. Adamson, Dec 28 2008
Second binomial transform of A168491. - Philippe Deléham, Nov 27 2009
a(n) is also the number of distinct strings of length n, each of which is a prefix of a string of balanced parentheses; see example. - Lee A. Newberg, Apr 26 2010
Number of symmetric balanced strings of n pairs of parentheses; see example. - Joerg Arndt, Jul 25 2011
a(n) is the number of permutation patterns modulo 2. - Olivier Gérard, Feb 25 2011
For n >= 2, a(n-1) is the number of incongruent two-color bracelets of 2*n-1 beads, n of which are black (A007123), having a diameter of symmetry. - Vladimir Shevelev, May 03 2011
The number of permutations of n elements where p(k-2) < p(k) for all k. - Joerg Arndt, Jul 23 2011
Also size of the equivalence class of S_{n+1} containing the identity permutation under transformations of positionally adjacent elements of the form abc <--> cba where a < b < c, cf. A210668. - Tom Roby, May 15 2012
a(n) is the number of symmetric Dyck paths of length 2n. - Matt Watson, Sep 26 2012
a(n) is divisible by A000108(floor(n/2)) = abs(A129996(n-2)). - Paul Curtz, Oct 23 2012
a(n) is the number of permutations of length n avoiding both 213 and 231 in the classical sense which are breadth-first search reading words of increasing unary-binary trees. For more details, see the entry for permutations avoiding 231 at A245898. - Manda Riehl, Aug 05 2014
Number of symmetric standard Young tableaux of shape (n,n). - Ran Pan, Apr 10 2015
From Luciano Ancora, May 09 2015: (Start)
Also "stepped path" in the array formed by partial sums of the all 1's sequence (or a Pascal's triangle displayed as a square). Example:
[1], [1], 1, 1, 1, 1, 1, ... A000012
1, [2], [3], 4, 5, 6, 7, ...
1, 3, [6], [10], 15, 21, 28, ...
1, 4, 10, [20], [35], 56, 84, ...
1, 5, 15, 35, [70], [126], 210, ...
Sequences in second formula are the mixed diagonals shown in this array. (End)
a(n) = A265848(n,n). - Reinhard Zumkeller, Dec 24 2015
The constant Sum_{n >= 0} a(n)/n! is 1 + A130820. - Peter Bala, Jul 02 2016
Number of meanders (walks starting at the origin and ending at any altitude >= 0 that may touch but never go below the x-axis) with n steps from {-1,1}. - David Nguyen, Dec 20 2016
a(n) is also the number of paths of n steps (either up or down by 1) that end at the maximal value achieved along the path. - Winston Luo, Jun 01 2017
Number of binary n-tuples such that the number of 1's in the even positions is the same as the number of 1's in the odd positions. - Juan A. Olmos, Dec 21 2017
Equivalently, a(n) is the number of subsets of {1,...,n} containing as many even numbers as odd numbers. - Gus Wiseman, Mar 17 2018
a(n) is the number of Dyck paths with semilength = n+1, returns to the x-axis = floor((n+3)/2) and up movements in odd positions = floor((n+3)/2). Example: a(4)=6, U=up movement in odd position, u=up movement in even position, d=down movement, -=return to x-axis: Uududd-Ud-Ud-, Ud-Uudd-Uudd-, Uudd-Uudd-Ud-, Ud-Ud-Uududd-, Uudd-Ud-Uudd-, Ud-Uududd-Ud-. - Roger Ford, Dec 29 2017
Let C_n(R, H) denote the transition matrix from the ribbon basis to the homogeneous basis of the graded component of the algebra of noncommutative symmetric functions of order n. Letting I(2^(n-1)) denote the identity matrix of order 2^(n-1), it has been conjectured that the dimension of the kernel of C_n(R, H) - I(2^(n-1)) is always equal to a(n-1). - John M. Campbell, Mar 30 2018
The number of U-equivalence classes of Łukasiewicz paths. Łukasiewicz paths are U-equivalent iff the positions of pattern U are identical in these paths. - Sergey Kirgizov, Apr 08 2018
All binary self-dual codes of length 2n, for n > 0, must contain at least a(n) codewords of weight n. More to the point, there will always be at least one, perhaps unique, binary self-dual code of length 2n that will contain exactly a(n) codewords that have a hamming weight equal to half the length of the code (n). This code can be constructed by direct summing the unique binary self-dual code of length 2 (up to permutation equivalence) to itself n times. A permutation equivalent code can be constructed by augmenting two identity matrices of length n together. - Nathan J. Russell, Nov 25 2018
Closed under addition. - Torlach Rush, Apr 18 2019
The sequence starting (1, 2, 3, 6, ...) is the invert transform of A097331: (1, 1, 0, 1, 0, 2, 0, 5, 0, 14, 0, 42, ...). - Gary W. Adamson, Feb 22 2020
From Gary W. Adamson, Feb 24 2020: (Start)
The sequence is the culminating limit of an infinite set of sequences with convergents of 2*cos(Pi/N), N = (3, 5, 7, 9, ...).
The first few such sequences are:
N = 3: (1, 1, 1, 1, 1, 1, 1, 1, ...)
N = 5: (1, 1, 2, 3, 5, 8, 13, 21, ...) = A000045
N = 7: (1, 1, 2, 3, 6, 10, 19, 33, ...) = A028495, a(n)/a(n-1) tends to 1.801937...
N = 9 (1, 1, 2, 3, 6, 10, 20, 35, ...) = A061551, a(n)/a(n_1) tends to 1.879385...
...
In the limit one gets the current sequence with ratio 2. (End)
a(n) is also the number of monotone lattice paths from (0,0) to (floor(n/2),ceiling(n/2)). These are the number of Grand Dyck paths when n is even. - Nachum Dershowitz, Aug 12 2020
The maximum number of preimages that a permutation of length n+1 can have under the consecutive-132-avoiding stack-sorting map. - Colin Defant, Aug 28 2020
Counts faro permutations of length n. Faro permutations are permutations avoiding the three consecutive patterns 231, 321 and 312. They are obtained by a perfect faro shuffle of two nondecreasing words of lengths differing by at most one. - Sergey Kirgizov, Jan 12 2021
Per "Sperner's Theorem", the largest possible familes of finite sets none of which contain any other sets in the family. - Renzo Benedetti, May 26 2021
a(n-1) are the incomplete, primitive Dyck paths of n steps without a first return: paths of U and D steps starting at the origin, never touching the horizontal axis later on, and ending above the horizontal axis. n=1: {U}, n=2: {UU}, n=3: {UUU, UUD}, n=4: {UUUU, UUUD, UUDU}, n=5: {UUUUU, UUUUD, UUUDD, UUDUU, UUUDU, UUDUD}. For comparison: A037952 counts incomplete Dyck paths with n steps with any number of intermediate returns to the horizontal axis, ending above the horizontal axis. - R. J. Mathar, Sep 24 2021
a(n) is the number of noncrossing partitions of [n] whose nontrivial blocks are of type {a,b}, with a <= n/2, b > n/2. - Francesca Aicardi, May 29 2022
Maximal coefficient of (1+x)^n. - Vaclav Kotesovec, Dec 30 2022
Sums of lower-left-to-upper-right diagonals of the Catalan Triangle A001263. - Howard A. Landman, Sep 16 2024

Examples

			For n = 4, the a(4) = 6 distinct strings of length 4, each of which is a prefix of a string of balanced parentheses, are ((((, (((), (()(, ()((, ()(), and (()). - _Lee A. Newberg_, Apr 26 2010
There are a(5)=10 symmetric balanced strings of 5 pairs of parentheses:
[ 1] ((((()))))
[ 2] (((()())))
[ 3] ((()()()))
[ 4] ((())(()))
[ 5] (()()()())
[ 6] (()(())())
[ 7] (())()(())
[ 8] ()()()()()
[ 9] ()((()))()
[10] ()(()())() - _Joerg Arndt_, Jul 25 2011
G.f. = 1 + x + 2*x^2 + 3*x^3 + 6*x^4 + 10*x^5 + 20*x^6 + 35*x^7 + 70*x^8 + ...
The a(4)=6 binary 4-tuples such that the number of 1's in the even positions is the same as the number of 1's in the odd positions are 0000, 1100, 1001, 0110, 0011, 1111. - _Juan A. Olmos_, Dec 21 2017
		

References

  • 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.
  • M. Aigner and G. M. Ziegler, Proofs from The Book, Springer-Verlag, Berlin, 1999; see p. 135.
  • K. Engel, Sperner Theory, Camb. Univ. Press, 1997; Theorem 1.1.1.
  • P. Frankl, Extremal sets systems, Chap. 24 of R. L. Graham et al., eds, Handbook of Combinatorics, North-Holland.
  • J. C. P. Miller, editor, Table of Binomial Coefficients. Royal Society Mathematical Tables, Vol. 3, Cambridge Univ. Press, 1954.
  • 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).
  • R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 7.16(b), p. 452.

Crossrefs

Row sums of Catalan triangle A053121 and of symmetric Dyck paths A088855.
Enumerates the structures encoded by A061854 and A061855.
First differences are in A037952.
Apparently a(n) = lim_{k->infinity} A094718(k, n).
Partial sums are in A036256. Column k=2 of A182172. Column k=1 of A335570.
Bisections: A000984 (even part), A001700 (odd part).
Cf. A097331.
Cf. A107373, A340567, A340568, A340569 (popularity of certain patterns in faro permutations).

Programs

  • GAP
    List([0..40],n->Binomial(n,Int(n/2))); # Muniru A Asiru, Apr 08 2018
    
  • Haskell
    a001405 n = a007318_row n !! (n `div` 2) -- Reinhard Zumkeller, Nov 09 2011
    
  • Magma
    [Binomial(n, Floor(n/2)): n in [0..40]]; // Vincenzo Librandi, Nov 16 2014
    
  • Maple
    A001405 := n->binomial(n, floor(n/2)): seq(A001405(n), n=0..33);
  • Mathematica
    Table[Binomial[n, Floor[n/2]], {n, 0, 40}] (* Stefan Steinerberger, Apr 08 2006 *)
    Table[DifferenceRoot[Function[{a,n},{-4 n a[n]-2 a[1+n]+(2+n) a[2+n] == 0,a[1] == 1,a[2] == 1}]][n], {n, 30}] (* Luciano Ancora, Jul 08 2015 *)
    Array[Binomial[#,Floor[#/2]]&,40,0] (* Harvey P. Dale, Mar 05 2018 *)
  • Maxima
    A001405(n):=binomial(n,floor(n/2))$
    makelist(A001405(n),n,0,30); /* Martin Ettl, Nov 01 2012 */
    
  • PARI
    a(n) = binomial(n, n\2);
    
  • PARI
    first(n) = x='x+O('x^n); Vec((-1+2*x+sqrt(1-4*x^2))/(2*x-4*x^2)) \\ Iain Fox, Dec 20 2017 (edited by Iain Fox, May 07 2018)
    
  • Python
    from math import comb
    def A001405(n): return comb(n,n//2) # Chai Wah Wu, Jun 07 2022

Formula

a(n) = max_{k=0..n} binomial(n, k).
a(2*n) = A000984(n), a(2*n+1) = A001700(n).
By symmetry, a(n) = binomial(n, ceiling(n/2)). - Labos Elemer, Mar 20 2003
P-recursive with recurrence: a(0) = 1, a(1) = 1, and for n >= 2, (n+1)*a(n) = 2*a(n-1) + 4*(n-1)*a(n-2). - Peter Bala, Feb 28 2011
G.f.: (1+x*c(x^2))/sqrt(1-4*x^2) = 1/(1 - x - x^2*c(x^2)); where c(x) = g.f. for Catalan numbers A000108.
G.f.: (-1 + 2*x + sqrt(1-4*x^2))/(2*x - 4*x^2). - Lee A. Newberg, Apr 26 2010
G.f.: 1/(1 - x - x^2/(1 - x^2/(1 - x^2/(1 - x^2/(1 - ... (continued fraction). - Paul Barry, Aug 12 2009
a(0) = 1; a(2*m+2) = 2*a(2*m+1); a(2*m+1) = Sum_{k = 0..2*m} (-1)^k*a(k)*a(2*m-k). - Len Smiley, Dec 09 2001
G.f.: (sqrt((1+2*x)/(1-2*x)) - 1)/(2*x). - Vladeta Jovovic, Apr 28 2003
The o.g.f. A(x) satisfies A(x) + x*A^2(x) = 1/(1-2*x). - Peter Bala, Feb 28 2011
E.g.f.: BesselI(0, 2*x) + BesselI(1, 2*x). - Vladeta Jovovic, Apr 28 2003
a(0) = 1; a(2*m+2) = 2*a(2*m+1); a(2*m+1) = 2*a(2*m) - c(m), where c(m)=A000108(m) are the Catalan numbers. - Christopher Hanusa (chanusa(AT)washington.edu), Nov 25 2003
a(n) = Sum_{k=0..n} (-1)^k*2^(n-k)*binomial(n, k)*A000108(k). - Paul Barry, Jan 27 2005
a(n) = Sum_{k=0..floor(n/2)} binomial(n, k)*binomial(1, n-2*k). - Paul Barry, May 13 2005
From Paul Barry, Nov 02 2004: (Start)
a(n) = Sum_{k=0..floor((n+1)/2)} (binomial(n+1, k)*(cos((n-2*k+1)*Pi/2) + sin((n-2*k+1)*Pi/2))).
a(n) = Sum_{k=0..n+1}, (binomial(n+1, (n-k+1)/2)*(1-(-1)^(n-k))*(cos(k*Pi/2) + sin(k*Pi))/2). (End)
a(n) = Sum_{k=floor(n/2)..n} (binomial(n,n-k) - binomial(n,n-k-1)). - Paul Barry, Sep 06 2007
Inverse binomial transform of A005773 starting (1, 2, 5, 13, 35, 96, ...) and double inverse binomial transform of A001700. Row sums of triangle A132815. - Gary W. Adamson, Aug 31 2007
a(n) = Sum_{k=0..n} A120730(n,k). - Philippe Deléham, Oct 16 2008
a(n) = Sum_{k = 0..floor(n/2)} (binomial(n,k) - binomial(n,k-1)). - Nishant Doshi (doshinikki2004(AT)gmail.com), Apr 06 2009
Sum_{n>=0} a(n)/10^(n+1) = 0.1123724... = (sqrt(3)-sqrt(2))/(2*sqrt(2)); Sum_{n>=0} a(n)/100^(n+1) = 0.0101020306102035... = (sqrt(51)-sqrt(49))/(2*sqrt(49)). - Mark Dols, Jul 15 2010
Conjectured: a(n) = 2^n*2F1(1/2,-n;2;2), useful for number of paths in 1-d for which the coordinate is never negative. - Benjamin Phillabaum, Feb 20 2011
a(2*m+1) = (2*m+1)*a(2*m)/(m+1), e.g., a(7) = (7/4)*a(6) = (7/4)*20 = 35. - Jon Perry, Jan 20 2011
From Peter Bala, Feb 28 2011: (Start)
Let F(x) be the logarithmic derivative of the o.g.f. A(x). Then 1+x*F(x) is the o.g.f. for A027306.
Let G(x) be the logarithmic derivative of 1+x*A(x). Then x*G(x) is the o.g.f. for A058622. (End)
Let M = an infinite tridiagonal matrix with 1's in the super and subdiagonals and [1,0,0,0,...] in the main diagonal; and V = the vector [1,0,0,0,...]. a(n) = M^n*V, leftmost term. - Gary W. Adamson, Jun 13 2011
Let M = an infinite tridiagonal matrix with 1's in the super and subdiagonals and [1,0,0,0,...] in the main diagonal. a(n) = M^n_{1,1}. - Corrected by Gary W. Adamson, Jan 30 2012
a(n) = A007318(n, floor(n/2)). - Reinhard Zumkeller, Nov 09 2011
a(n+1) = Sum_{k=0..n} a(n-k)*A097331(k) = a(n) + Sum_{k=0..(n-1)/2} A000108(k)*a(n-2*k-1). - Philippe Deléham, Nov 27 2011
a(n) = A214282(n) - A214283(n), for n > 0. - Reinhard Zumkeller, Jul 14 2012
a(n) = Sum_{k=0..n} A168511(n,k)*(-1)^(n-k). - Philippe Deléham, Mar 19 2013
a(n+2*p-2) = Sum_{k=0..floor(n/2)} A009766(n-k+p-1, k+p-1) + binomial(n+2*p-2, p-2), for p >= 1. - Johannes W. Meijer, Aug 02 2013
O.g.f.: (1-x*c(x^2))/(1-2*x), with the o.g.f. c(x) of Catalan numbers A000108. See the rewritten formula given by Lee A. Newberg above. This is the o.g.f. for the row sums the Riordan triangle A053121. - Wolfdieter Lang, Sep 22 2013
a(n) ~ 2^n / sqrt(Pi * n/2). - Charles R Greathouse IV, Oct 23 2015
a(n) = 2^n*hypergeom([1/2,-n], [2], 2). - Vladimir Reshetnikov, Nov 02 2015
a(2*k) = Sum_{i=0..k} binomial(k, i)*binomial(k, i), a(2*k+1) = Sum_{i=0..k} binomial(k+1, i)*binomial(k, i). - Juan A. Olmos, Dec 21 2017
a(0) = 1, a(n) = 2 * a(n-1) for even n, a(n) = (2*n/(n+1)) * a(n-1) for odd n. - James East, Sep 25 2019
a(n) = A037952(n) + A000108(n/2) where A(.)=0 for non-integer argument. - R. J. Mathar, Sep 23 2021
From Amiram Eldar, Mar 10 2022: (Start)
Sum_{n>=0} 1/a(n) = 2*Pi/(3*sqrt(3)) + 2.
Sum_{n>=0} (-1)^n/a(n) = 2/3 - 2*Pi/(9*sqrt(3)). (End)
For k>2, Sum_{n>=0} a(n)/k^n = (sqrt((k+2)/(k-2)) - 1)*k/2. - Vaclav Kotesovec, May 13 2022
From Peter Bala, Mar 24 2023: (Start)
a(n) = Sum_{k = 0..n+1} (-1)^(k+binomial(n+2,2)) * k/(n+1) * binomial(n+1,k)^2.
(n + 1)*(2*n - 1)*a(n) = (-1)^(n+1)*2*a(n-1) + 4*(n - 1)*(2*n + 1)*a(n-2) with a(0) = a(1) = 1. (End)
a(n) = Integral_{x=-2..2} x^n*W(x)*dx, n>=0, where W(x) = sqrt((2+x)/(2-x))/(2*Pi) is a positive function on x=(-2,2) and is singular at x = 2. Therefore a(n) is a positive definite sequence. - Karol A. Penson, May 12 2025

A053121 Catalan triangle (with 0's) read by rows.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 2, 0, 1, 2, 0, 3, 0, 1, 0, 5, 0, 4, 0, 1, 5, 0, 9, 0, 5, 0, 1, 0, 14, 0, 14, 0, 6, 0, 1, 14, 0, 28, 0, 20, 0, 7, 0, 1, 0, 42, 0, 48, 0, 27, 0, 8, 0, 1, 42, 0, 90, 0, 75, 0, 35, 0, 9, 0, 1, 0, 132, 0, 165, 0, 110, 0, 44, 0, 10, 0, 1, 132, 0, 297, 0, 275, 0, 154, 0, 54, 0, 11, 0
Offset: 0

Views

Author

Keywords

Comments

Inverse lower triangular matrix of A049310(n,m) (coefficients of Chebyshev's S polynomials).
Walks with a wall: triangle of number of n-step walks from (0,0) to (n,m) where each step goes from (a,b) to (a+1,b+1) or (a+1,b-1) and the path stays in the nonnegative quadrant.
T(n,m) is the number of left factors of Dyck paths of length n ending at height m. Example: T(4,2)=3 because we have UDUU, UUDU, and UUUD, where U=(1,1) and D=(1,-1). (This is basically a different formulation of the previous - walks with a wall - property.) - Emeric Deutsch, Jun 16 2011
"The Catalan triangle is formed in the same manner as Pascal's triangle, except that no number may appear on the left of the vertical bar." [Conway and Smith]
G.f. for row polynomials p(n,x) := Sum_{m=0..n} (a(n,m)*x^m): c(z^2)/(1-x*z*c(z^2)). Row sums (x=1): A001405 (central binomial).
In the language of the Shapiro et al. reference such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Bell-subgroup of the Riordan-group. The g.f. Ginv(x) of the m=0 column of the inverse of a given Bell-matrix (here A049310) is obtained from its g.f. of the m=0 column (here G(x)=1/(1+x^2)) by Ginv(x)=(f^{(-1)}(x))/x, with f(x) := x*G(x) and f^{(-1)}is the compositional inverse function of f (here one finds, with Ginv(0)=1, c(x^2)). See the Shapiro et al. reference.
Number of involutions of {1,2,...,n} that avoid the patterns 132 and have exactly k fixed points. Example: T(4,2)=3 because we have 2134, 4231 and 3214. Number of involutions of {1,2,...,n} that avoid the patterns 321 and have exactly k fixed points. Example: T(4,2)=3 because we have 1243, 1324 and 2134. Number of involutions of {1,2,...,n} that avoid the patterns 213 and have exactly k fixed points. Example: T(4,2)=3 because we have 1243, 1432 and 4231. - Emeric Deutsch, Oct 12 2006
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0) -> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Riordan array (c(x^2),xc(x^2)), where c(x) is the g.f. of Catalan numbers A000108. - Philippe Deléham, Nov 25 2007
A053121^2 = triangle A145973. Convolved with A001405 = triangle A153585. - Gary W. Adamson, Dec 28 2008
By columns without the zeros, n-th row = A000108 convolved with itself n times; equivalent to A = (1 + x + 2x^2 + 5x^3 + 14x^4 + ...), then n-th row = coefficients of A^(n+1). - Gary W. Adamson, May 13 2009
Triangle read by rows,product of A130595 and A064189 considered as infinite lower triangular arrays; A053121 = A130595*A064189 = B^(-1)*A097609*B where B = A007318. - Philippe Deléham, Dec 07 2009
From Mark Dols, Aug 17 2010: (Start)
As an upper right triangle, rows represent powers of 5-sqrt(24):
5 - sqrt(24)^1 = 0.101020514...
5 - sqrt(24)^2 = 0.010205144...
5 - sqrt(24)^3 = 0.001030928...
(Divided by sqrt(96) these powers give a decimal representation of the columns of A007318, with 1/sqrt(96) being the middle column.) (End)
T(n,k) is the number of dispersed Dyck paths of length n (i.e., Motzkin paths of length n with no (1,0) steps at positive heights) having k (1,0)-steps. Example: T(5,3)=4 because, denoting U=(1,1), D=(1,-1), H=(1,0), we have HHHUD, HHUDH, HUDHH, and UDHHH. - Emeric Deutsch, Jun 01 2011
Let S(N,x) denote the N-th Chebyshev S-polynomial in x (see A049310, cf. [W. Lang]). Then x^n = sum_{k=0..n} T(n,k)*S(k,x). - L. Edson Jeffery, Sep 06 2012
This triangle a(n,m) appears also in the (unreduced) formula for the powers rho(N)^n for the algebraic number over the rationals rho(N) = 2*cos(Pi/N) = R(N, 2), the smallest diagonal/side ratio R in the regular N-gon:
rho(N)^n = sum(a(n,m)*R(N,m+1),m=0..n), n>=0, identical in N >= 1. R(N,j) = S(j-1, x=rho(N)) (Chebyshev S (A049310)). See a comment on this under A039599 (even powers) and A039598 (odd powers). Proof: see the Sep 06 2012 comment by L. Edson Jeffery, which follows from T(n,k) (called here a(n,k)) being the inverse of the Riordan triangle A049310. - Wolfdieter Lang, Sep 21 2013
The so-called A-sequence for this Riordan triangle of the Bell type (c(x^2), x*c(x^2)) (see comments above) is A(x) = 1 + x^2. This proves the recurrence given in the formula section by Henry Bottomley for a(n, m) = a(n-1, m-1) + a(n-1, m+1) for n>=1 and m>=1, with inputs. The Z-sequence for this Riordan triangle is Z(x) = x which proves the recurrence a(n,0) = a(n-1,1), n>=1, a(0,0) = 1. For A- and Z-sequences for Riordan triangles see the W. Lang link under A006232. - Wolfdieter Lang, Sep 22 2013
Rows of the triangle describe decompositions of tensor powers of the standard (2-dimensional) representation of the Lie algebra sl(2) into irreducibles. Thus a(n,m) is the multiplicity of the m-th ((m+1)-dimensional) irreducible representation in the n-th tensor power of the standard one. - Mamuka Jibladze, May 26 2015
The Riordan row polynomials p(n, x) belong to the Boas-Buck class (see a comment and references in A046521), hence they satisfy the Boas-Buck identity: (E_x - n*1)*p(n, x) = (E_x + 1)*Sum_{j=0..n-1} (1/2)*(1 - (-1)^j)*binomial(j+1, (j+1)/2)*p(n-1-j, x), for n >= 0, where E_x = x*d/dx (Euler operator). For the triangle a(n, m) this entails a recurrence for the sequence of column m, given in the formula section. - Wolfdieter Lang, Aug 11 2017
From Roger Ford, Jan 22 2018: (Start)
For row n, the nonzero values represent the odd components (loops) formed by n+1 nonintersecting arches above and below the x-axis with the following constraints: The top has floor((n+3)/2) starting arches at position 1 and the next consecutive odd positions. All other starting top arches are in even positions. The bottom arches are a rainbow of arches. If the component=1 then the arch configuration is a semimeander solution.
Examples: For row 3 {0, 2, 0, 1} there are 3 arch configurations: 2 arch configurations have a component=1; 1 has a component=3. c=components, U=top arch starting in odd position, u=top arch starting in an even position, d=ending top arch:
.
top UuUdUddd c=3 top UdUuUddd c=1 top UdUdUudd c=1
/\ /\
//\\ / \
// \\ / /\ \ /\
// \\ / / \ \ / \
///\ /\\\ /\ / / /\ \ \ /\ /\ / /\ \
\\\ \/ /// \ \ \ \/ / / / \ \ \ \/ / / /
\\\ /// \ \ \ / / / \ \ \ / / /
\\\/// \ \ \/ / / \ \ \/ / /
\\// \ \ / / \ \ / /
\/ \ \/ / \ \/ /
\ / \ /
\/ \/
For row 4 {2, 0, 3, 0, 1} there are 6 arch configurations: 2 have a component=1; 3 have a component=3: 1 has a component=1. (End)

Examples

			Triangle a(n,m) begins:
  n\m  0   1   2   3   4   5   6  7  8  9 10 ...
  0:   1
  1:   0   1
  2:   1   0   1
  3:   0   2   0   1
  4:   2   0   3   0   1
  5:   0   5   0   4   0   1
  6:   5   0   9   0   5   0   1
  7:   0  14   0  14   0   6   0  1
  8:  14   0  28   0  20   0   7  0  1
  9:   0  42   0  48   0  27   0  8  0  1
  10: 42   0  90   0  75   0  35  0  9  0  1
  ... (Reformatted by _Wolfdieter Lang_, Sep 20 2013)
E.g., the fourth row corresponds to the polynomial p(3,x)= 2*x + x^3.
From _Paul Barry_, May 29 2009: (Start)
Production matrix is
  0, 1,
  1, 0, 1,
  0, 1, 0, 1,
  0, 0, 1, 0, 1,
  0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 0, 0, 0, 1, 0, 1 (End)
Boas-Buck recurrence for column k = 2, n = 6: a(6, 2) = (3/4)*(0 + 2*a(4 ,2) + 0 + 6*a(2, 2)) = (3/4)*(2*3 + 6) = 9. - _Wolfdieter Lang_, Aug 11 2017
		

References

  • J. H. Conway and D. A. Smith, On Quaternions and Octonions, A K Peters, Ltd., Natick, MA, 2003. See p. 60. MR1957212 (2004a:17002)
  • A. Nkwanta, Lattice paths and RNA secondary structures, in African Americans in Mathematics, ed. N. Dean, Amer. Math. Soc., 1997, pp. 137-147.

Crossrefs

Cf. A008315, A049310, A000108, A001405 (row sums), A145973, A153585, A108786, A037952. Another version: A008313. A039598 and A039599 without zeros, and odd and even numbered rows.
Variant without zero-diagonals: A033184 and with rows reversed: A009766.

Programs

  • Haskell
    a053121 n k = a053121_tabl !! n !! k
    a053121_row n = a053121_tabl !! n
    a053121_tabl = iterate
       (\row -> zipWith (+) ([0] ++ row) (tail row ++ [0,0])) [1]
    -- Reinhard Zumkeller, Feb 24 2012
    
  • Maple
    T:=proc(n,k): if n+k mod 2 = 0 then (k+1)*binomial(n+1,(n-k)/2)/(n+1) else 0 fi end: for n from 0 to 13 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form; Emeric Deutsch, Oct 12 2006
    F:=proc(l,p) if ((l-p) mod 2) = 1 then 0 else (p+1)*l!/( ( (l-p)/2 )! * ( (l+p)/2 +1)! ); fi; end;
    r:=n->[seq( F(n,p),p=0..n)]; [seq(r(n),n=0..15)]; # N. J. A. Sloane, Jan 29 2011
    A053121 := proc(n,k) option remember; `if`(k>n or k<0,0,`if`(n=k,1,
    procname(n-1,k-1)+procname(n-1,k+1))) end proc:
    seq(print(seq(A053121(n,k), k=0..n)),n=0..12); # Peter Luschny, May 01 2011
  • Mathematica
    a[n_, m_] /; n < m || OddQ[n-m] = 0; a[n_, m_] = (m+1) Binomial[n+1, (n-m)/2]/(n+1); Flatten[Table[a[n, m], {n, 0, 12}, {m, 0, n}]] [[1 ;; 90]] (* Jean-François Alcover, May 18 2011 *)
    T[0, 0] := 1; T[n_, k_]/;0<=k<=n := T[n, k] = T[n-1, k-1]+T[n-1, k+1]; T[n_, k_] := 0; Flatten@Table[T[n, k], {n, 0, 12}, {k, 0, n}] (* Oliver Seipel, Dec 31 2024 *)
  • PARI
    T(n, m)=if(nCharles R Greathouse IV, Mar 09 2016
  • Sage
    def A053121_triangle(dim):
        M = matrix(ZZ,dim,dim)
        for n in (0..dim-1): M[n,n] = 1
        for n in (1..dim-1):
            for k in (0..n-1):
                M[n,k] = M[n-1,k-1] + M[n-1,k+1]
        return M
    A053121_triangle(13) # Peter Luschny, Sep 19 2012
    

Formula

a(n, m) := 0 if n
a(n, m) = (4*(n-1)*a(n-2, m) + 2*(m+1)*a(n-1, m-1))/(n+m+2), a(n, m)=0 if n
G.f. for m-th column: c(x^2)*(x*c(x^2))^m, where c(x) = g.f. for Catalan numbers A000108.
G.f.: G(t,z) = c(z^2)/(1 - t*z*c(z^2)), where c(z) = (1 - sqrt(1-4*z))/(2*z) is the g.f. for the Catalan numbers (A000108). - Emeric Deutsch, Jun 16 2011
a(n, m) = a(n-1, m-1) + a(n-1, m+1) if n > 0 and m >= 0, a(0, 0)=1, a(0, m)=0 if m > 0, a(n, m)=0 if m < 0. - Henry Bottomley, Jan 25 2001
Sum_{k>=0} T(m,k)^2 = A000108(m). - Paul D. Hanna, Apr 23 2005
Sum_{k>=0} T(m, k)*T(n, k) = 0 if m+n is odd; Sum_{k>=0} T(m, k)*T(n, k) = A000108((m+n)/2) if m+n is even. - Philippe Deléham, May 26 2005
T(n,k)=sum{i=0..n, (-1)^(n-i)*C(n,i)*sum{j=0..i, C(i,j)*(C(i-j,j+k)-C(i-j,j+k+2))}}; Column k has e.g.f. BesselI(k,2x)-BesselI(k+2,2x). - Paul Barry, Feb 16 2006
Sum_{k=0..n} T(n,k)*(k+1) = 2^n. - Philippe Deléham, Mar 22 2007
Sum_{j>=0} T(n,j)*binomial(j,k) = A054336(n,k). - Philippe Deléham, Mar 30 2007
T(2*n+1,2*k+1) = A039598(n,k), T(2*n,2*k) = A039599(n,k). - Philippe Deléham, Apr 16 2007
Sum_{k=0..n} T(n,k)^x = A000027(n+1), A001405(n), A000108(n), A003161(n), A129123(n) for x = 0,1,2,3,4 respectively. - Philippe Deléham, Nov 22 2009
Sum_{k=0..n} T(n,k)*x^k = A126930(n), A126120(n), A001405(n), A054341(n), A126931(n) for x = -1, 0, 1, 2, 3 respectively. - Philippe Deléham, Nov 28 2009
Sum_{k=0..n} T(n,k)*A000045(k+1) = A098615(n). - Philippe Deléham, Feb 03 2012
Recurrence for row polynomials C(n, x) := Sum_{m=0..n} a(n, m)*x^m = x*Sum_{k=0..n} Chat(k)*C(n-1-k, x), n >= 0, with C(-1, 1/x) = 1/x and Chat(k) = A000108(k/2) if n is even and 0 otherwise. From the o.g.f. of the row polynomials: G(z; x) := Sum_{n >= 0} C(n, x)*z^n = c(z^2)*(1 + x*z*G(z, x)), with the o.g.f. c of A000108. - Ahmet Zahid KÜÇÜK and Wolfdieter Lang, Aug 23 2015
The Boas-Buck recurrence (see a comment above) for the sequence of column m is: a(n, m) = ((m+1)/(n-m))*Sum_{j=0..n-1-m} (1/2)*(1 - (-1)^j)*binomial(j+1, (j+1)/2)* a(n-1-j, k), for n > m >= 0 and input a(m, m) = 1. - Wolfdieter Lang, Aug 11 2017
Sum_{m=1..n} a(n,m) = A037952(n). - R. J. Mathar, Sep 23 2021

Extensions

Edited by N. J. A. Sloane, Jan 29 2011

A061554 Square table read by antidiagonals: a(n,k) = binomial(n+k, floor(k/2)).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 6, 4, 4, 1, 1, 10, 10, 5, 5, 1, 1, 20, 15, 15, 6, 6, 1, 1, 35, 35, 21, 21, 7, 7, 1, 1, 70, 56, 56, 28, 28, 8, 8, 1, 1, 126, 126, 84, 84, 36, 36, 9, 9, 1, 1, 252, 210, 210, 120, 120, 45, 45, 10, 10, 1, 1, 462, 462, 330, 330, 165, 165, 55, 55, 11, 11, 1, 1
Offset: 0

Author

Henry Bottomley, May 17 2001

Keywords

Comments

Equivalently, a triangle read by rows, where the rows are obtained by sorting the elements of the rows of Pascal's triangle (A007318) into descending order. - Philippe Deléham, May 21 2005
Equivalently, as a triangle read by rows, this is T(n,k)=binomial(n,floor((n-k)/2)); column k then has e.g.f. Bessel_I(k,2x)+Bessel_I(k+1,2x). - Paul Barry, Feb 28 2006
Antidiagonal sums are A037952(n+1) = C(n+1,[n/2]). Matrix inverse is the row reversal of triangle A066170. Eigensequence is A125094(n) = Sum_{k=0..n-1} A125093(n-1,k)*A125094(k). - Paul D. Hanna, Nov 20 2006
Riordan array (1/(1-x-x^2*c(x^2)),x*c(x^2)); where c(x)=g.f.for Catalan numbers A000108. - Philippe Deléham, Mar 17 2007
Triangle T(n,k), 0<=k<=n, read by rows given by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 27 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
T(n,k) is the number of paths from (0,k) to some (n,m) which never dip below y=0, touch y=0 at least once and are made up only of the steps (1,1) and (1,-1). This can be proved using the recurrence supplied by Deléham. - Gerald McGarvey, Oct 15 2008
Triangle read by rows = partial sums of A053121 terms starting from the right. - Gary W. Adamson, Oct 24 2008
As a subset of the "family of triangles" (Deleham comment of Sep 25 2007), beginning with a signed variant of A061554, M = (-1,0) = (1; -1, 1; 2, -1, 1; -3, 3, -1, 1; ...) successive binomial transforms of M yield (0,1) - A089942; (1,2) - A039599; (2,3) - A124733; (3,4) - A124574; (4,5) - A126331; ... such that the binomial transform of the triangle generated from (n,n+1) = the triangle generated from (n+1,n+2). Similarly, another subset beginning with A053121 - (0,0), and taking successive binomial transforms yields (1,1) - A064189; (2,2) - A039598; (3,3) - A091965, ... By rows, the triangle generated from (n,n) can be obtained by taking pairwise sums from the (n-1,n) triangle starting from the right. For example, row 2 of (1,2) - A039599 = (2, 3, 1); and taking pairwise sums from the right we obtain (5, 4, 1) = row 2 of (2,2) - A039598. - Gary W. Adamson, Aug 04 2011
The triangle by rows (n) with alternating signs (+-+...) from the top as a set of simultaneous equations solves for diagonal lengths of odd N (N = 2n+1) regular polygons. The constants in each case are powers of c = 2*cos(2*Pi/N). By way of example, the first 3 rows relate to the heptagon and the simultaneous equations are (1,0,0) = 1; (-1,1,0) = c = 1.24697...; and (2,-1,1) = c^2. The answers are 1, 2.24697..., and 1.801...; the 3 distinct diagonal lengths of the heptagon with edge = 1. - Gary W. Adamson, Sep 07 2011

Examples

			The array starts:
   1, 1,  2,  3,  6, 10,  20,  35,   70,  126, ...
   1, 1,  3,  4, 10, 15,  35,  56,  126,  210, ...
   1, 1,  4,  5, 15, 21,  56,  84,  210,  330, ...
   1, 1,  5,  6, 21, 28,  84, 120,  330,  495, ...
   1, 1,  6,  7, 28, 36, 120, 165,  495,  715, ...
   1, 1,  7,  8, 36, 45, 165, 220,  715, 1001, ...
   1, 1,  8,  9, 45, 55, 220, 286, 1001, 1365, ...
   1, 1,  9, 10, 55, 66, 286, 364, 1365, 1820, ...
   1, 1, 10, 11, 66, 78, 364, 455, 1820, 2380, ...
   1, 1, 11, 12, 78, 91, 455, 560, 2380, 3060, ...
Triangle (antidiagonal) version begins:
    1;
    1,   1;
    2,   1,   1;
    3,   3,   1,   1;
    6,   4,   4,   1,   1;
   10,  10,   5,   5,   1,   1;
   20,  15,  15,   6,   6,   1,  1;
   35,  35,  21,  21,   7,   7,  1,  1;
   70,  56,  56,  28,  28,   8,  8,  1,  1;
  126, 126,  84,  84,  36,  36,  9,  9,  1,  1;
  252, 210, 210, 120, 120,  45, 45, 10, 10,  1, 1;
  462, 462, 330, 330, 165, 165, 55, 55, 11, 11, 1, 1; ...
Matrix inverse begins:
   1;
  -1,  1;
  -1, -1,   1;
   1, -2,  -1,   1;
   1,  2,  -3,  -1,  1;
  -1,  3,   3,  -4, -1,  1;
  -1, -3,   6,   4, -5, -1,  1;
   1, -4,  -6,  10,  5, -6, -1,  1;
   1,  4, -10, -10, 15,  6, -7, -1, 1; ...
From _Paul Barry_, May 21 2009: (Start)
Production matrix is
  1, 1,
  1, 0, 1,
  0, 1, 0, 1,
  0, 0, 1, 0, 1,
  0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 1, 0, 1,
  0, 0, 0, 0, 0, 1, 0, 1 (End)
		

Crossrefs

Rows are A001405, A037952, A037955, A037951, A037956, A037953, A037957 etc. Columns are truncated pairs of A000012, A000027, A000217, A000292, A000332, A000389, A000579, etc. Main diagonal is alternate values of A051036.

Programs

  • Maple
    T := proc(n, k) option remember;
    if n = k then 1 elif k < 0 or n < 0 or k > n then 0
    elif k = 0 then T(n-1, 0) + T(n-1, 1) else T(n-1, k-1) + T(n-1, k+1) fi end:
    for n from 0 to 9 do seq(T(n, k), k = 0..n) od; # Peter Luschny, May 25 2021
  • Mathematica
    t[n_, k_] = Binomial[n, Floor[(n+1)/2 - (-1)^(n-k)*(k+1)/2]]; Flatten[Table[t[n, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, May 31 2011 *)
  • PARI
    T(n,k)=binomial(n,(n+1)\2-(-1)^(n-k)*((k+1)\2))

Formula

As a triangle: T(n,k) = binomial(n,m) where m = floor((n+1)/2 - (-1)^(n-k)*(k+1)/2).
a(0, k) = binomial(k, floor(k/2)) = A001405(k); for n>0 T(n, k) = T(n+1, k-2) + T(n-1, k).
n-th row = M^n * V, where M = the infinite tridiagonal matrix with all 1's in the super and subdiagonals and (1,0,0,0,...) in the main diagonal. V = the infinite vector [1,0,0,0,...]. Example: (3,3,1,1,0,0,0,...) = M^3 * V. - Gary W. Adamson, Nov 04 2006
Sum_{k=0..n} T(m,k)*T(n,k) = T(m+n,0) = A001405(m+n). - Philippe Deléham, Feb 26 2007
Sum_{k=0..n} T(n,k)=2^n. - Philippe Deléham, Mar 27 2007
Sum_{k=0..n} T(n,k)*x^k = A127361(n), A126869(n), A001405(n), A000079(n), A127358(n), A127359(n), A127360(n) for x = -2, -1, 0, 1, 2, 3, 4 respectively. - Philippe Deléham, Dec 04 2009

Extensions

Entry revised by N. J. A. Sloane, Nov 22 2006

A001206 Number of self-dual monotone Boolean functions of n variables.

Original entry on oeis.org

0, 1, 2, 4, 12, 81, 2646, 1422564, 229809982112, 423295099074735261880
Offset: 0

Keywords

Comments

Sometimes called Hosten-Morris numbers (or HM numbers).
Also the number of simplicial complexes on the set {1, ..., n-1} such that no pair of faces covers all of {1, ..., n-1}. [Miller-Sturmfels] - N. J. A. Sloane, Feb 18 2008
Also the maximal number of generators of a neighborly monomial ideal in n variables. [Miller-Sturmfels]. - N. J. A. Sloane, Feb 18 2008
Also the number of intersecting antichains on a labeled (n-1)-set or (n-1)-variable Boolean functions in the Post class F(7,2). Cf. A059090. - Vladeta Jovovic, Goran Kilibarda, Dec 28 2000
Also the number of nondominated coteries on n members. - Don Knuth, Sep 01 2005
The number of maximal families of intersecting subsets of an n-element set. - Bridget Tenner, Nov 16 2006
Rivière gives a(n) for n <= 5. - N. J. A. Sloane, May 12 2012

Examples

			a(2) = 1 + 1 = 2;
a(3) = 1 + 3 = 4;
a(4) = 1 + 7 + 3 + 1 = 12;
a(5) = 1 + 15 + 30 + 30 + 5 = 81;
a(6) = 1 + 31 + 195 + 605 + 780 + 543 + 300 + 135 + 45 + 10 + 1 = 2646;
a(7) = 1 + 63 + 1050 + 9030 + 41545 + 118629 + 233821 + 329205 + 327915 + 224280 + 100716 + 29337 + 5950 + 910 + 105 + 1 = 1422564.
Cf. A059090.
From _Gus Wiseman_, Jul 03 2019: (Start)
The a(1) = 1 through a(4) = 12 intersecting antichains of nonempty sets (see Jovovic and Kilibarda's comment):
  {}  {}     {}       {}
      {{1}}  {{1}}    {{1}}
             {{2}}    {{2}}
             {{1,2}}  {{3}}
                      {{1,2}}
                      {{1,3}}
                      {{2,3}}
                      {{1,2,3}}
                      {{1,2},{1,3}}
                      {{1,2},{2,3}}
                      {{1,3},{2,3}}
                      {{1,2},{1,3},{2,3}}
(End)
		

References

  • Martin Aigner and Günter M. Ziegler, Proofs from THE BOOK, Third Edition, Springer-Verlag, 2004. See chapter 22.
  • V. Jovovic and G. Kilibarda, The number of n-variable Boolean functions in the Post class F(7,2), Belgrade, 2001, in preparation.
  • D. E. Knuth, The Art of Computer Programming, Vol. 4A, Section 7.1.1, p. 79.
  • W. F. Lunnon, The IU function: the size of a free distributive lattice, pp. 173-181 of D. J. A. Welsh, editor, Combinatorial Mathematics and Its Applications. Academic Press, NY, 1971.
  • Charles F. Mills and W. M. Mills, The calculation of λ(8), preprint, 1979. Gives a(8).
  • E. Miller and B. Sturmfels, Combinatorial Commutative Algebra, Springer, 2005.
  • 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

The case with empty edges allowed is A326372.
The maximal case is A007363, or A326363 with empty edges allowed.
The case with empty intersection is A326366.
The inverse binomial transform is the covering case A305844.

Programs

  • Mathematica
    stableSets[u_,Q_]:=If[Length[u]==0,{{}},With[{w=First[u]},Join[stableSets[DeleteCases[u,w],Q],Prepend[#,w]&/@stableSets[DeleteCases[u,r_/;r==w||Q[r,w]||Q[w,r]],Q]]]];
    Table[Length[stableSets[Subsets[Range[n],{1,n}],Or[Intersection[#1,#2]=={},SubsetQ[#1,#2]]&]],{n,0,5}] (* Gus Wiseman, Jul 03 2019 *)

Formula

a(n+1) = Sum_{m=0..A037952(n)} A059090(n, m).
For n > 0, a(n) = A326372(n - 1) - 1. - Gus Wiseman, Jul 03 2019

Extensions

a(8) due to C. F. Mills & W. H. Mills, 1979
a(8) from Daniel E. Loeb, Jan 04 1996
a(8) confirmed by Don Knuth, Feb 08 2008
a(9) from Andries E. Brouwer, Aug 25 2012

A057979 a(n) = 1 for even n and (n-1)/2 for odd n.

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 1, 3, 1, 4, 1, 5, 1, 6, 1, 7, 1, 8, 1, 9, 1, 10, 1, 11, 1, 12, 1, 13, 1, 14, 1, 15, 1, 16, 1, 17, 1, 18, 1, 19, 1, 20, 1, 21, 1, 22, 1, 23, 1, 24, 1, 25, 1, 26, 1, 27, 1, 28, 1, 29, 1, 30, 1, 31, 1, 32, 1, 33, 1, 34, 1, 35, 1, 36, 1, 37, 1, 38, 1, 39, 1, 40, 1, 41, 1, 42, 1, 43, 1
Offset: 0

Author

Labos Elemer, Nov 13 2000

Keywords

Comments

a(n) = b(n)/c(n) where b(n) = A001405(n+1) - A001405(n), c(n) = gcd(A001405(n+1), A001405(n)).
Also the minimal number of disjoint edge-paths into which the complete graph on n edges can be partitioned - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jan 19 2001
For n >= 2, number of partitions of n-2 into parts that are distinct mod 2. - Giovanni Resta, Feb 06 2006
Sequence starting with a(3) obeys the rule "smallest positive value such that the ordered pair (a(n-1),a(n)) has not occurred previously", or the rule "smallest positive value such that the ratio a(n-1)/a(n) has not occurred previously". The same subsequence has its ordinal transform equal to itself, shifted left. (The ordinal transform has as its n-th term the number of values in a(1),...,a(n) that are equal to a(n).) - Franklin T. Adams-Watters, Dec 13 2006
Numerators of floor(n/2)/n, n > 0. - Wesley Ivan Hurt, Jun 14 2013
Number of nonisomorphic outer planar graphs of order n >= 3, maximum degree 3, and largest possible size. The size is (3n-2)/2 when n is even and (3n-3)/2 when n is odd. - Christian Barrientos and Sarah Minion, Feb 27 2018

Examples

			For n=12, C(12,6) - C(11,5) = 924 - 462 = 462, gcd(C(12,6), C(11,5)) = 462, and the quotient is 1.
For n=13, C(13,6) - C(12,6) = 792, gcd(C(13,6), C(12,6)) = 132, and the quotient is 6.
		

Programs

  • Haskell
    import Data.List (transpose)
    a057979 n = 1 - rest * (1 - n') where (n', rest) = divMod n 2
    a057979_list = concat $ transpose [repeat 1, [0..]]
    -- Reinhard Zumkeller, Aug 11 2014
    
  • Magma
    [Floor(n/2)^(n mod 2): n in [0..100]]; // Vincenzo Librandi, Mar 17 2015
    
  • Maple
    A057979:=n->(n+1)/4+(3-n)*(-1)^n/4; seq(A057979(k), k=0..100); # Wesley Ivan Hurt, Oct 14 2013
  • Mathematica
    With[{no=45},Riffle[Table[1,{no}],Range[0,no-1]]] (* Harvey P. Dale, Feb 18 2011 *)
  • PARI
    a(n)=if(n%2,n-1,2)/2 \\ Charles R Greathouse IV, Sep 02 2015
    
  • Python
    def A057979(n): return n>>1 if n&1 else 1 # Chai Wah Wu, Jan 04 2024

Formula

a(n) = (n+1)/4+(3-n)*(-1)^n/4. - Paul Barry, Mar 21 2003, corrected by Hieronymus Fischer, Sep 25 2007
a(n) = (a(n-2) + a(n-3)) / a(n-1).
From Paul Barry, Oct 21 2004: (Start)
G.f.: (1-x^2+x^3)/((1+x)^2(1-x)^2);
a(n) = 2*a(n-2) - a(n-4);
a(n) = 0^n + Sum_{k=0..floor((n-2)/2)} C(n-k-2,k) * C(1,n-2k-2). (End)
a(n) = gcd(n-1, floor((n-1)/2)). - Paul Barry, May 02 2005
a(n) = binomial((2n-3)/4-(-1)^n/4, (1-(-1)^n)/2). - Paul Barry, Jun 29 2006
G.f.: (x^3-x^2+1)/(1-x^2)^2 = 1 + x^2*G(0) where G(k) = 1 + x*(k+1)/(1 - x/(x + (k+1)/G(k+1) )); (continued fraction, 3-step). - Sergei N. Gladkovskii, Nov 29 2012
a(n) = binomial(floor(n/2), n mod 2). - Wesley Ivan Hurt, Oct 14 2013
a(n) = 1 - n mod 2 * (1 - floor(n/2)). - Reinhard Zumkeller, Aug 11 2014
a(n) = floor(n/2)^(n mod 2). - Wesley Ivan Hurt, Mar 16 2015
E.g.f.: ((2 + x)*cosh(x) - sinh(x))/2. - Stefano Spezia, Mar 26 2022

A026010 a(n) = number of (s(0), s(1), ..., s(n)) such that s(i) is a nonnegative integer and |s(i) - s(i-1)| = 1 for i = 1,2,...,n and s(0) = 2. Also a(n) = sum of numbers in row n+1 of array T defined in A026009.

Original entry on oeis.org

1, 2, 4, 7, 14, 25, 50, 91, 182, 336, 672, 1254, 2508, 4719, 9438, 17875, 35750, 68068, 136136, 260338, 520676, 999362, 1998724, 3848222, 7696444, 14858000, 29716000, 57500460, 115000920, 222981435, 445962870, 866262915, 1732525830, 3370764540
Offset: 0

Keywords

Comments

Conjecture: a(n) is the number of integer compositions of n + 2 in which the even parts appear as often at even positions as at odd positions (confirmed up to n = 19). - Gus Wiseman, Mar 17 2018

Examples

			The a(3) = 7 compositions of 5 in which the even parts appear as often at even positions as at odd positions are (5), (311), (131), (113), (221), (122), (11111). Missing are (41), (14), (32), (23), (212), (2111), (1211), (1121), (1112). - _Gus Wiseman_, Mar 17 2018
		

Crossrefs

Programs

  • Magma
    [(&+[Binomial(Floor((n+k)/2), Floor(k/2)): k in [0..n]]): n in [0..40]]; // G. C. Greubel, Nov 08 2018
  • Mathematica
    Array[Sum[Binomial[Floor[(# + k)/2], Floor[k/2]], {k, 0, #}] &, 34, 0] (* Michael De Vlieger, May 16 2018 *)
    Table[2^(-1 + n)*(((2 + 3*#)*Gamma[(1 + #)/2])/(Sqrt[Pi]*Gamma[2 + #/2]) &[n + Mod[n, 2]]), {n,0,40}] (* Peter Pein, Nov 08 2018 *)
    Table[(1/2)^((5 - (-1)^n)/2)*(6*n + 7 - 3*(-1)^n)*CatalanNumber[(2*n + 1 - (-1)^n)/4], {n, 0, 40}] (* G. C. Greubel, Nov 08 2018 *)
  • PARI
    vector(40, n, n--; sum(k=0,n, binomial(floor((n+k)/2), floor(k/2)))) \\ G. C. Greubel, Nov 08 2018
    

Formula

a(2*n) = ((3*n + 1)/(2*n + 1))*C(2*n + 1, n)= A051924(1+n), n>=0, a(2*n-1) = a(2*n)/2 = A097613(1+n), n >= 1. - Herbert Kociemba, May 08 2004
a(n) = Sum_{k=0..n} binomial(floor((n+k)/2), floor(k/2)). - Paul Barry, Jul 15 2004
Inverse binomial transform of A005774: (1, 3, 9, 26, 75, 216, ...). - Gary W. Adamson, Oct 22 2007
Conjecture: (n+3)*a(n) - 2*a(n-1) + (-5*n-3)*a(n-2) + 2*a(n-3) + 4*(n-3)*a(n-4) = 0. - R. J. Mathar, Jun 20 2013
a(n) = (1/2)^((5 - (-1)^n)/2)*(6*n + 7 - 3*(-1)^n)*Catalan((2*n + 1 - (-1)^n)/4), where Catalan is the Catalan number = A000108. - G. C. Greubel, Nov 08 2018

A014495 Central binomial coefficient - 1.

Original entry on oeis.org

0, 0, 1, 2, 5, 9, 19, 34, 69, 125, 251, 461, 923, 1715, 3431, 6434, 12869, 24309, 48619, 92377, 184755, 352715, 705431, 1352077, 2704155, 5200299, 10400599, 20058299, 40116599, 77558759, 155117519, 300540194, 601080389, 1166803109, 2333606219, 4537567649
Offset: 0

Author

N. J. A. Sloane, Denis Pochuev (denis(AT)cdc.Informatik.TH-Darmstadt.de)

Keywords

Comments

For n > 0: sum of positive elements in row (n-1) of triangle A214292. - Reinhard Zumkeller, Jul 12 2012
Number of Young tableaux with n cells and 2 rows. Also number of self-inverse permutations in S_n with longest increasing subsequence of length 2. The a(4) = 5 permutations are 1432, 2143, 3214, 3412, 4231 and the a(5) = 9 permutations are 15432, 21543, 32154, 35142, 42513, 43215, 45312, 52431, 53241. - Alois P. Heinz, Oct 03 2012
Number of nonempty subsets of {1,2,...,n} that contain the same number of even and odd numbers. For example, a(5)=9 and the 9 subsets are {1,2}, {1,4}, {2,3}, {2,5}, {3,4}, {4,5}, {1,2,3,4}, {1,2,4,5}, {2,3,4,5}. - Enrique Navarrete, Feb 10 2018

Crossrefs

Cf. A001405, A037952 (first differences).
a(n) = A094718(n, n) = A094718(n-1, n)+1.
a(n) = A047884(n, 2) for n>=2. - Alois P. Heinz, Oct 03 2012
Cf. A214292.

Programs

  • Magma
    [Binomial(n, Floor(n/2))-1: n in [0..50]]; // Vincenzo Librandi, Feb 11 2018
  • Maple
    a:= n-> binomial(n, iquo(n, 2))-1:
    seq(a(n), n=0..40);  # Alois P. Heinz, Oct 03 2012
  • Mathematica
    Table[Binomial[n, Floor[n/2]] - 1, {n, 0, 50}] (* Bruno Berselli, Oct 03 2012 *)
  • Maxima
    A014495(n):=binomial(n,floor(n/2))-1$
    makelist(A014495(n),n,0,30); /* Martin Ettl, Nov 01 2012 */
    

Formula

a(n) = A001405(n)-1.
a(n) = C(n,floor(n/2))-1. - Alois P. Heinz, Oct 03 2012
(n+1)*a(n)-2*a(n-1)-4*(n-1)*a(n-2) = 3*n-3 with n>1, a(0)=a(1)=0. - Bruno Berselli, Oct 03 2012
D-finite with recurrence: -(n+1)*(n-2)*a(n) +(n^2+n-4)*a(n-1) +2*(n-1)*(2*n-5)*a(n-2) -4*(n-1)*(n-2)*a(n-3)=0. - Conjectured by R. J. Mathar, Jan 04 2017, confirmed by Robert Israel, Feb 11 2018
G.f.: (x+1)/(2*x*(x-1)) - sqrt(1-4*x^2)/(2*x*(2*x-1)). - Robert Israel, Feb 11 2018

Extensions

Edited by Andrey Zabolotskiy, Apr 14 2025

A116382 Riordan array (1/sqrt(1-4*x^2), (1-2*x^2*c(x^2))*(x^2*c(x^2))/(x*(1-x-x^2*c(x^2)))) where c(x) is the g.f. of A000108.

Original entry on oeis.org

1, 0, 1, 2, 1, 1, 0, 3, 2, 1, 6, 4, 5, 3, 1, 0, 10, 10, 8, 4, 1, 20, 15, 21, 19, 12, 5, 1, 0, 35, 42, 42, 32, 17, 6, 1, 70, 56, 84, 92, 77, 50, 23, 7, 1, 0, 126, 168, 192, 180, 131, 74, 30, 8, 1, 252, 210, 330, 405, 400, 326, 210, 105, 38, 9, 1
Offset: 0

Author

Paul Barry, Feb 12 2006

Keywords

Comments

Row sums are A116383. Diagonal sums are A116384.
First column has e.g.f. Bessel_I(0,2*x) (A000984 with interpolated zeros).
Second column has e.g.f. Bessel_I(1,2*x) + Bessel_I(2,2*x) (A037952).
Third column has e.g.f. Bessel_I(2,2*x) + 2*Bessel_I(3,2*x) + Bessel_I(4,2*x) (A116385).
A binomial-Bessel triangle: column k has e.g.f. Sum_{j=0..k} C(k,j) * Bessel_I(k+j,2*x).

Examples

			Triangle begins
    1;
    0,   1;
    2,   1,   1;
    0,   3,   2,   1;
    6,   4,   5,   3,   1;
    0,  10,  10,   8,   4,   1;
   20,  15,  21,  19,  12,   5,   1;
    0,  35,  42,  42,  32,  17,   6,   1;
   70,  56,  84,  92,  77,  50,  23,   7,  1;
    0, 126, 168, 192, 180, 131,  74,  30,  8, 1;
  252, 210, 330, 405, 400, 326, 210, 105, 38, 9, 1;
		

Programs

  • GAP
    Flat(List([0..10], n-> List([0..n], k-> Sum([0..n], j-> (-1)^(n-j)*Binomial(n,j)*Sum([0..j], m-> Binomial(j,m-k)*Binomial(m,j-m)  ))))); # G. C. Greubel, May 22 2019
  • Magma
    T:= func< n,k | (&+[(-1)^(n-j)*Binomial(n,j)*(&+[Binomial(j,m-k)* Binomial(m,j-m): m in [0..j]]): j in [0..n]]) >;
    [[T(n,k): k in [0..n]]: n in [0..10]]; // G. C. Greubel, May 22 2019
    
  • Mathematica
    T[n_, k_] := Sum[(-1)^(n-j)*Binomial[n, j]*Sum[Binomial[j, i-k]* Binomial[i, j-i], {i, 0, j}], {j, 0, n}];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 24 2018 *)
  • PARI
    {T(n,k) = sum(j=0,n, (-1)^(n-j)*binomial(n,j)*sum(m=0,j, binomial(j,m-k)*binomial(m,j-m) ))}; \\ G. C. Greubel, May 22 2019
    
  • Sage
    def T(n, k): return sum((-1)^(n-j)*binomial(n,j)*sum(binomial(j,m-k)*binomial(m,j-m) for m in (0..j)) for j in (0..n))
    [[T(n, k) for k in (0..n)] for n in (0..10)] # G. C. Greubel, May 22 2019
    

Formula

Riordan array (1/sqrt(1-4*x^2), sqrt(1-4*x^2)*(1-sqrt(1-4*x^2))/(x-2*x^2 + x*sqrt(1-4*x^2))).
Number triangle T(n,k) = Sum{j=0..n} (-1)^(n-j)* C(n,j)*Sum_{i=0..j} C(j,i-k)*C(i,j-i).
Showing 1-10 of 29 results. Next