A049600 -id:A049600 - cp's OEIS frontend

A008288 Square array of Delannoy numbers D(i,j) (i >= 0, j >= 0) read by antidiagonals.

1, 1, 1, 1, 3, 1, 1, 5, 5, 1, 1, 7, 13, 7, 1, 1, 9, 25, 25, 9, 1, 1, 11, 41, 63, 41, 11, 1, 1, 13, 61, 129, 129, 61, 13, 1, 1, 15, 85, 231, 321, 231, 85, 15, 1, 1, 17, 113, 377, 681, 681, 377, 113, 17, 1, 1, 19, 145, 575, 1289, 1683, 1289, 575, 145, 19, 1, 1, 21, 181, 833, 2241, 3653, 3653

Offset: 0

N. J. A. Sloane

In the Formula section, some contributors use T(n,k) = D(n-k, k) (for 0 <= k <= n), which is the triangular version of the square array (D(n,k): n,k >= 0). Conversely, D(n,k) = T(n+k,k) for n,k >= 0. - Petros Hadjicostas, Aug 05 2020
Also called the tribonacci triangle [Alladi and Hoggatt (1977)]. - N. J. A. Sloane, Mar 23 2014
D(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (1,0), (0,1), (1,1). - Joerg Arndt, Jul 01 2011 [Corrected by N. J. A. Sloane, May 30 2020]
Or, triangle read by rows of coefficients of polynomials P[n](x) defined by P[0] = 1, P[1] = x+1; for n >= 2, P[n] = (x+1)*P[n-1] + x*P[n-2].
D(n, k) is the number of k-matchings of a comb-like graph with n+k teeth. Example: D(1, 3) = 7 because the graph consisting of a horizontal path ABCD and the teeth Aa, Bb, Cc, Dd has seven 3-matchings: four triples of three teeth and the three triples {Aa, Bb, CD}, {Aa, Dd, BC}, {Cc, Dd, AB}. Also D(3, 1)=7, the 1-matchings of the same graph being the seven edges: {AB}, {BC}, {CD}, {Aa}, {Bb}, {Cc}, {Dd}. - Emeric Deutsch, Jul 01 2002
Sum of n-th antidiagonal of the array D is A000129(n+1). - Reinhard Zumkeller, Dec 03 2004 [Edited by Petros Hadjicostas, Aug 05 2020 so that the counting of antidiagonals of D starts at n = 0. That is, the sum of the terms in the n-th row of the triangles T is  A000129(n+1).]
The A-sequence for this Riordan type triangle (see one of Paul Barry's comments under Formula) is A112478 and the Z-sequence the trivial: {1, 0, 0, 0, ...}. See the W. Lang link under A006232 for Sheffer a- and z-sequences where also Riordan A- and Z-sequences are explained. This leads to the recurrence for the triangle given below. - Wolfdieter Lang, Jan 21 2008
The triangle or chess sums, see A180662 for their definitions, link the Delannoy numbers with twelve different sequences, see the crossrefs. All sums come in pairs due to the symmetrical nature of this triangle. The knight sums Kn14 and Kn15 have been added. It is remarkable that all knight sums are related to the tribonacci numbers, that is, A000073 and A001590, but none of the others. - Johannes W. Meijer, Sep 22 2010
This sequence, A008288, is jointly generated with A035607 as an array of coefficients of polynomials u(n,x): initially, u(1,x) = v(1,x) = 1; for n > 1, u(n,x) = x*u(n-1,x) + v(n-1) and v(n,x) = 2*x*u(n-1,x) + v(n-1,x). See the Mathematica section. - Clark Kimberling, Mar 09 2012
Row n, for n > 0, of Roger L. Bagula's triangle in the Example section shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(n)*x^n which is the numerator of the n-th convergent of the continued fraction [k, k, k, ...], where k = sqrt(x) + 1/sqrt(x); see A230000. - Clark Kimberling, Nov 13 2013
In an n-dimensional hypercube lattice, D(n,k) gives the number of nodes situated at a Minkowski (Manhattan) distance of k from a given node. In cellular automata theory, the cells at Manhattan distance k are called the von Neumann neighborhood of radius k. For k=1, see A005843. - Dmitry Zaitsev, Dec 10 2015
These numbers appear as the coefficients of series relating spherical and bispherical harmonics, in the solutions of Laplace's equation in 3D. [Majic 2019, Eq. 22] - Matt Majic, Nov 24 2019
From Peter Bala, Feb 19 2020: (Start)
The following remarks assume an offset of 1 in the row and column indices of the triangle.
The sequence of row polynomials T(n,x), beginning with T(1,x) = x, T(2,x) = x + x^2, T(3,x) = x + 3*x^2 + x^3, ..., is a strong divisibility sequence of polynomials in the ring Z[x]; that is, for all positive integers n and m, poly_gcd(T(n,x), T(m,x)) = T(gcd(n, m), x) - apply Norfleet (2005), Theorem 3. Consequently, the sequence (T(n,x): n >= 1) is a divisibility sequence in the polynomial ring Z[x]; that is, if n divides m then T(n,x) divides T(m,x) in Z[x].
Let S(x) = 1 + 2*x + 6*x^2 + 22*x^3 + ... denote the o.g.f. for the large Schröder numbers A006318. The power series (x*S(x))^n, n = 2, 3, 4, ..., can be expressed as a linear combination with polynomial coefficients of S(x) and 1: (x*S(x))^n = T(n-1,-x) - T(n,-x)*S(x). The result can be extended to negative integer n if we define T(0,x) = 0 and T(-n,x) = (-1)^(n+1) * T(n,x)/x^n. Cf. A115139.
[In the previous two paragraphs, D(n,x) was replaced with T(n,x) because the contributor is referring to the rows of the triangle T(n,k), not the rows of the array D(n,k). - Petros Hadjicostas, Aug 05 2020] (End)
Named after the French amateur mathematician Henri-Auguste Delannoy (1833-1915). - Amiram Eldar, Apr 15 2021
D(i,j) = D(j,i). With this and Dmitry Zaitsev's Dec 10 2015 comment, D(i,j) can be considered the number of points at L1 distance <= i in Z^j or the number of points at L1 distance <= j in Z^i from any given point. The rows and columns of D(i,j) are the crystal ball sequences on cubic lattices. See the first example below. The n-th term in the k-th crystal ball sequence can be considered the number of points at distance <= n from any point in a k-dimensional cubic lattice, or the number of points at distance <= k from any point in an n-dimensional cubic lattice. - Shel Kaphan, Jan 01 2023 and Jan 07 2023
Dimensions of hom spaces Hom(R^{(i)}, R^{(j)}) in the Delannoy category attached to the oligomorphic group of order preserving self-bijections of the real line. - Noah Snyder, Mar 22 2023

			The square array D(i,j) (i >= 0, j >= 0) begins:
  1, 1,  1,   1,   1,   1,    1,    1,    1,    1, ... = A000012
  1, 3,  5,   7,   9,  11,   13,   15,   17,   19, ... = A005408
  1, 5, 13,  25,  41,  61,   85,  113,  145,  181, ... = A001844
  1, 7, 25,  63, 129, 231,  377,  575,  833, 1159, ... = A001845
  1, 9, 41, 129, 321, 681, 1289, 2241, 3649, 5641, ... = A001846
  ...
For D(2,5) = 61, which is seen above in the row labeled A001844, we calculate the sum (9 + 11 + 41) of the 3 nearest terms above and/or to the left. - _Peter Munn_, Jan 01 2023
D(2,5) = 61 can also be obtained from the row labeled A005408 using a recurrence mentioned in the formula section:  D(2,5) = D(1,5) + 2*Sum_{k=0..4} D(1,k), so D(2,5) = 11 + 2*(1+3+5+7+9) = 11 + 2*25. - _Shel Kaphan_, Jan 01 2023
As a triangular array (on its side) this begins:
   0,   0,   0,   0,   1,   0,  11,   0, ...
   0,   0,   0,   1,   0,   9,   0,  61, ...
   0,   0,   1,   0,   7,   0,  41,   0, ...
   0,   1,   0,   5,   0,  25,   0, 129, ...
   1,   0,   3,   0,  13,   0,  63,   0, ...
   0,   1,   0,   5,   0,  25,   0, 129, ...
   0,   0,   1,   0,   7,   0,  41,   0, ...
   0,   0,   0,   1,   0,   9,   0,  61, ...
   0,   0,   0,   0,   1,   0,  11,   0, ...
   [Edited by _Shel Kaphan_, Jan 01 2023]
From _Roger L. Bagula_, Dec 09 2008: (Start)
As a triangle T(n,k) (with rows n >= 0 and columns k = 0..n), this begins:
   1;
   1,  1;
   1,  3,   1;
   1,  5,   5,   1;
   1,  7,  13,   7,    1;
   1,  9,  25,  25,    9,    1;
   1, 11,  41,  63,   41,   11,    1;
   1, 13,  61, 129,  129,   61,   13,   1;
   1, 15,  85, 231,  321,  231,   85,  15,   1;
   1, 17, 113, 377,  681,  681,  377, 113,  17,  1;
   1, 19, 145, 575, 1289, 1683, 1289, 575, 145, 19, 1;
   ... (End)
Triangle T(n,k) recurrence: 63 = T(6,3) = 25 + 13 + 25 = T(5,2) + T(4,2) + T(5,3).
Triangle T(n,k) recurrence with A-sequence A112478: 63 = T(6,3) = 1*25 + 2*25 - 2*9 + 6*1 (T entries from row n = 5 only). [Here the formula T(n,k) = Sum_{j=0..n-k} A112478(j) * T(n-1, k-1+j) is used with n = 6 and k = 3; i.e., T(6,3) = Sum_{j=0..3} A111478(j) * T(5, 2+j). - _Petros Hadjicostas_, Aug 05 2020]
From _Philippe Deléham_, Mar 29 2012: (Start)
Subtriangle of the triangle given by (1, 0, 1, -1, 0, 0, 0, ...) DELTA (0, 1, 0, 0, 0, ...) where DELTA is the operator defined in A084938:
   1;
   1,  0;
   1,  1,  0;
   1,  3,  1,  0;
   1,  5,  5,  1,  0;
   1,  7, 13,  7,  1,  0;
   1,  9, 25, 25,  9,  1, 0;
   1, 11, 41, 63, 41, 11, 1, 0;
   ...
Subtriangle of the triangle given by (0, 1, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, ...) where DELTA is the operator defined in A084938:
   1;
   0, 1;
   0, 1,  1;
   0, 1,  3,  1;
   0, 1,  5,  5,  1;
   0, 1,  7, 13,  7,  1;
   0, 1,  9, 25, 25,  9,  1;
   0, 1, 11, 41, 63, 41, 11, 1;
   ... (End)

Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, page 593.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids (in Russian), FAN, Tashkent, 1990, ISBN 5-648-00738-8.
L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
L. Moser and W. Zayachkowski, Lattice paths with diagonal steps, Scripta Mathematica, 26 (1963), 223-229.
G. Picou, Note #2235, L'Intermédiaire des Mathématiciens, 8 (1901), page 281. - N. J. A. Sloane, Mar 02 2022
D. B. West, Combinatorial Mathematics, Cambridge, 2021, p. 28.

T. D. Noe, Table of n, a(n) for n = 0..5150
K. Alladi and V. E. Hoggatt Jr., On tribonacci numbers and related functions, Fibonacci Quart. 15 (1977), 42-45.
Said Amrouche and Hacène Belbachir, Asymmetric extension of Pascal-Dellanoy triangles, arXiv:2001.11665 [math.CO], 2020.
J.-M. Autebert et al., H.-A. Delannoy et les oeuvres posthumes d'Édouard Lucas, Gazette des Mathématiciens - no 95, Jan 2003 (in French).
Bela Bajnok, Additive Combinatorics: A Menu of Research Problems, arXiv:1705.07444 [math.NT], 2017. See Sect. 2.3.
C. Banderier and S. Schwer, Why Delannoy numbers?, arXiv:math/0411128 [math.CO], 2004.
Paul Barry, On Integer-Sequence-Based Constructions of Generalized Pascal Triangles, Journal of Integer Sequences, 9 (2006), #06.2.4.
Paul Barry, Riordan-Bernstein Polynomials, Hankel Transforms and Somos Sequences, Journal of Integer Sequences, 15 (2012), #12.8.2.
Paul Barry, On the Inverses of a Family of Pascal-Like Matrices Defined by Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.5.6.
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
Paul Barry, The Central Coefficients of a Family of Pascal-like Triangles and Colored Lattice Paths, J. Int. Seq., 22 (2019), #19.1.3.
Paul Barry, Generalized Catalan Numbers Associated with a Family of Pascal-like Triangles, J. Int. Seq., 22 (2019), #19.5.8.
Paul Barry, Notes on Riordan arrays and lattice paths, arXiv:2504.09719 [math.CO], 2025. See pp. 3, 29.
Frédéric Bihan, Francisco Santos, and Pierre-Jean Spaenlehauer, A Polyhedral Method for Sparse Systems with many Positive Solutions, arXiv:1804.05683 [math.CO], 2018.
Boris A. Bondarenko, Generalized Pascal Triangles and Pyramids, English translation published by Fibonacci Association, Santa Clara Univ., Santa Clara, CA, 1993; see p. 37.
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I, Math. Zeit. 233 (2000), 1-19.
C. Carré, N. Debroux, M. Deneufchatel, J.-P. Dubernard et al., Dirichlet convolution and enumeration of pyramid polycubes, hal-00905889, 2013.
C. Carre, N. Debroux, M. Deneufchatel, J.-Ph. Dubernard, C. Hillariet, J.-G. Luque, and O. Mallet, Enumeration of Polycubes and Dirichlet Convolutions, J. Int. Seq. 18 (2015), #15.11.4.
J. S. Caughman et al., A note on lattice chains and Delannoy numbers, Discrete Math., 308 (2008), 2623-2628.
Swee Hong Chan, Igor Pak, and Greta Panova, Log-concavity in planar random walks, arXiv:2106.10640 [math.CO], 2021.
H. Delannoy, Emploi de l'échiquier pour la résolution de certains problèmes de probabilités, Association Française pour l'Avancement des Sciences, 24th session, 1895, pp. 70-90 (see the table given on p. 76).
Jerry Ray Dias, Properties and relationships of conjugated polyenes having a reciprocal eigenvalue spectrum - dendralene and radialene hydrocarbons, Croatica Chem. Acta, 77 (2004), 325-330.
M. Dziemianczuk, Generalizing Delannoy numbers via counting weighted lattice paths, INTEGERS, 13 (2013), #A54.
James East and Nicholas Ham, Lattice paths and submonoids of Z^2, arXiv:1811.05735 [math.CO], 2018.
Steven Edwards and William Griffiths, Generalizations of Delannoy and cross polytope numbers, Fib. Q., Vol. 55, No. 4 (2017), pp. 356-366.
Steven Edwards and William Griffiths, On Generalized Delannoy Numbers, J. Int. Seq., 23 (2020), #20.3.6.
R. Feria-Puron, H. Perez-Roses, and J. Ryan, Searching for Large Circulant Graphs, arXiv:1503.07357 [math.CO], 2015.
R. Feria-Purón, J. Ryan, and H. Pérez-Rosés, Searching for Large Multi-Loop Networks, Electronic Notes in Discrete Mathematics, 46 (2014), 233-240.
Nate Harman, Andrew Snowden, and Noah Snyder, The Delannoy Category, arxiv:2211.15392 [math.RT], 2023.
Rebecca Hartman-Baker, The Diffusion Equation Method for Global Optimization and Its Application to Magnetotelluric Geoprospecting, University of Illinois, Urbana-Champaign, 2005.
G. Hetyei, Shifted Jacobi polynomials and Delannoy numbers, arXiv:0909.5512 [math.CO], 2009.
G. Hetyei, Links we almost missed between Delannoy numbers and Legendre polynomials.
V. E. Hoggatt, Jr., Letters to N. J. A. Sloane, 1974-1975.
Kentaro Ihara, Yayoi Nakamura, and Shuji Yamamoto, Interpolant of truncated multiple zeta functions, arXiv:2407.20509 [math.NT], 2024. See p. 21.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
Milan Janjić and B. Petkovic, A Counting Function, arXiv:1301.4550 [math.CO], 2013. - _N. J. A. Sloane_, Feb 13 2013
Milan Janjić and B. Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014), #14.3.5.
Svante Janson, Patterns in random permutations avoiding some sets of multiple patterns, arXiv:1804.06071 [math.PR], 2018.
Shel Kaphan, Tables of Sequences Related to Delannoy Numbers and Cubic Lattice Coordination Numbers
Shel Kaphan, Illustration of a recurrence relation on the Delannoy numbers and their connection with geometry.
G. Kreweras, Sur les hiérarchies de segments, Cahiers Bureau Universitaire Recherche Opérationnelle, # 20, Inst. Statistiques, Univ. Paris, 1973, pp. 4-10.
G. Kreweras, Sur les hiérarchies de segments, Cahiers du Bureau Universitaire de Recherche Opérationnelle, Institut de Statistique, Université de Paris, #20 (1973). (Annotated scanned copy)
G. Kreweras, Aires des chemins surdiagonaux et application à un problème économique, Cahiers du Bureau universitaire de recherche opérationnelle Série Recherche 24 (1976), 1-8. [Annotated scanned copy]
Eon Lee, Andrés R. Vindas-Meléndez, and Zhi Wang, Generalized snake posets, order polytopes, and lattice-point enumeration, arXiv:2411.18695 [math.CO], 2024. See p. 15.
Yi-Lin Lee, Off-diagonally symmetric domino tilings of the Aztec diamond of odd order, arXiv:2404.09057 [math.CO], 2024. See p. 20.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 1.
M. LLadser, Uniform formulas for coefficients of meromorphic functions, arXiv:math/0604152 [math.CO], 2006.
E. Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 174.
Matt Majic, Relationships between spherical and bispherical harmonics, and an electrostatic T-matrix for dimers, preprint, 2019, DOI:10.13140/RG.2.2.21203.12320.
J. W. Meijer, Famous numbers on a chessboard, Acta Nova, 4(4) (2010), 589-598.
Mirka Miller, Hebert Perez-Roses, and Joe Ryan, The Maximum Degree-and-Diameter-Bounded Subgraph in the Mesh, arXiv:1203.4069 [math.CO], 2012.
Alejandro H. Morales, Igor Pak, and Greta Panova, Hook formulas for skew shapes IV. Increasing tableaux and factorial Grothendieck polynomials, arXiv:2108.10140 [math.CO], 2021.
Lili Mu and Sai-nan Zheng, On the Total Positivity of Delannoy-Like Triangles, Journal of Integer Sequences, 20 (2017), #17.1.6.
M. Norfleet, Characterization of second-order strong divisibility sequences of polynomials, The Fibonacci Quarterly, 43(2) (2005), 166-169.
Richard L. Ollerton and Anthony G. Shannon, Some properties of generalized Pascal squares and triangles, Fib. Q., 36 (1998), 98-109. See Table 9.
L. Pachter and B. Sturmfels, The mathematics of phylogenomics, arXiv:math/0409132 [math.ST], 2004-2005.
R. Pemantle and M. C. Wilson, Asymptotics of multivariate sequences, I: smooth points of the singular variety, arXiv:math/0003192 [math.CO], 2000.
J. L. Ramirez and V. F. Sirvent, Incomplete Tribonacci Numbers and Polynomials, Journal of Integer Sequences, 17 (2014), #14.4.2. See Table 1. - _N. J. A. Sloane_, Mar 23 2014
Marko Razpet, A self-similarity structure generated by king's walk, Algebraic and topological methods in graph theory (Lake Bled, 1999). Discrete Math. 244(1-3) (2002), 423--433. MR1844050 (2002k:05022).
Shiva Samieinia, Digital straight line segments and curves, Licentiate Thesis. Stockholm University, Department of Mathematics, Report 2007:6.
Shiva Samieinia, The number of continuous curves in digital geometry, Port. Math. 67(1) (2010), 75-89.
Seunghyun Seo, The Catalan Threshold Arrangement, Journal of Integer Sequences, 20 (2017), #17.1.1.
Yuriy Shablya, Combinatorial Generation Algorithms for Some Lattice Paths Using the Method Based on AND/OR Trees, Algorithms (2023) Vol. 16, No. 6, 266.
M. Shattuck, Combinatorial identities for incomplete tribonacci polynomials, arXiv:1406.2755 [math.CO], 2014.
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Yi Wang, Zheng Sai-Nan, and Chen Xi, Analytic aspects of Delannoy numbers, Discrete Mathematics 342(8) (2019), 2270-2277.
Eric Weisstein's World of Mathematics, Delannoy Number.
Meng-Han Wu, Henryk A. Witek, Rafał Podeszwa, Clar Covers and Zhang-Zhang Polynomials of Zigzag and Armchair Carbon Nanotubes, MATCH Commun. Math. Comput. Chem. (2025) Vol. 93, 415-462. See p. 445.
Dmitry Zaitsev, k-neighborhood for Cellular Automata, arXiv:1605.08870 [cs.DM], 2016.
Liang Zhao and Fengyao Yan, Note on Total Positivity for a Class of Recursive Matrices, Journal of Integer Sequences, 19 (2016), #16.6.5.

Sums of antidiagonals: A000129 (Pell numbers).
Main diagonal: A001850 (central Delannoy numbers), which has further information and references.
A002002, A026002, and A190666 are +-k-diagonals for k=1, 2, 3 resp. - Shel Kaphan, Jan 01 2023
Rows 0..10: A000012, A005408, A001844, A001845, A001846, A001847, A001848, A001849, A008417, A008419, A008421.
See also A027618.
Cf. A059446.
Has same main diagonal as A064861. Different from A100936.
Cf. A101164, A101167, A128966, A131887, A131935.
Cf. A035607, A108625, A142979, A142992, A143007.
Read mod small primes: A211312, A211313, A211314, A211315.
Triangle sums (see the comments): A000129 (Row1); A056594 (Row2); A000073 (Kn11 & Kn21); A089068 (Kn12 & Kn22); A180668 (Kn13 & Kn23); A180669 (Kn14 & Kn24); A180670 (Kn15 & Kn25); A099463 (Kn3 & Kn4); A116404 (Fi1 & Fi2); A006498 (Ca1 & Ca2); A006498(3*n) (Ca3 & Ca4); A079972 (Gi1 & Gi2); A079972(4*n) (Gi3 & Gi4); A079973(3*n) (Ze1 & Ze2); A079973(2*n) (Ze3 & Ze4).
Cf. A102413, A128966. (D(n,1)) = A005843. Cf. A115139.

a008288 n k = a008288_tabl !! n !! k
a008288_row n = a008288_tabl !! n
a008288_tabl = map fst $ iterate
    (\(us, vs) -> (vs, zipWith (+) ([0] ++ us ++ [0]) $
                       zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [1, 1])
-- Reinhard Zumkeller, Jul 21 2013

A008288 := proc(n, k) option remember; if k = 0 then 1 elif n=k then 1 else procname(n-1, k-1) + procname(n-2, k-1) + procname(n-1, k) end if; end proc: seq(seq(A008288(n,k),k=0..n), n=0..10); # triangular indices n and k
P[0]:=1; P[1]:=x+1; for n from 2 to 12 do P[n]:=expand((x+1)*P[n-1]+x*P[n-2]); lprint(P[n]); lprint(seriestolist(series(P[n],x,200))); end do:

(* Next, A008288 jointly generated with A035607 *)
u[1, x_] := 1; v[1, x_] := 1; z = 16;
u[n_, x_] := x*u[n - 1, x] + v[n - 1, x];
v[n_, x_] := 2 x*u[n - 1, x] + v[n - 1, x];
Table[Expand[u[n, x]], {n, 1, z/2}]
Table[Expand[v[n, x]], {n, 1, z/2}]
cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
TableForm[cu]
Flatten[%]    (* A008288 *)
Table[Expand[v[n, x]], {n, 1, z}]
cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
TableForm[cv]
Flatten[%]    (* A035607 *)
(* Clark Kimberling, Mar 09 2012 *)
d[n_, k_] := Binomial[n+k, k]*Hypergeometric2F1[-k, -n, -n-k, -1]; A008288 = Flatten[Table[d[n-k, k], {n, 0, 12}, {k, 0, n}]] (* Jean-François Alcover, Apr 05 2012, after 3rd formula *)

from functools import cache
@cache
def delannoy_row(n: int) -> list[int]:
    if n == 0: return [1]
    if n == 1: return [1, 1]
    rov = delannoy_row(n - 2)
    row = delannoy_row(n - 1) + [1]
    for k in range(n - 1, 0, -1):
        row[k] += row[k - 1] + rov[k - 1]
    return row
for n in range(10): print(delannoy_row(n))  # Peter Luschny, Jul 30 2023

for k in range(8):  # seen as an array, read row by row
    a = lambda n: hypergeometric([-n, -k], [1], 2)
    print([simplify(a(n)) for n in range(11)]) # Peter Luschny, Nov 19 2014

D(n, 0) = 1 = D(0, n) for n >= 0; D(n, k) = D(n, k-1) + D(n-1, k-1) + D(n-1, k).
Bivariate o.g.f.: Sum_{n >= 0, k >= 0} D(n, k)*x^n*y^k = 1/(1 - x - y - x*y).
D(n, k) = Sum_{d = 0..min(n,k)} binomial(k, d)*binomial(n+k-d, k) = Sum_{d=0..min(n,k)} 2^d*binomial(n, d)*binomial(k, d). [Edited by Petros Hadjicostas, Aug 05 2020]
Seen as a triangle read by rows: T(n, 0) = T(n, n) = 1 for n >= 0 and T(n, k) = T(n-1, k-1) + T(n-2, k-1) + T(n-1, k), 0 < k < n and n > 1. - Reinhard Zumkeller, Dec 03 2004
Read as a number triangle, this is the Riordan array (1/(1-x), x(1+x)/(1-x)) with T(n, k) = Sum_{j=0..n-k} C(n-k, j) * C(k, j) * 2^j. - Paul Barry, Jul 18 2005
T(n,k) = Sum_{j=0..n-k} C(k,j)*C(n-j,k). - Paul Barry, May 21 2006
Let y^k(n) be the number of Khalimsky-continuous functions f from [0,n-1] to Z such that f(0) = 0 and f(n-1) = k. Then y^k(n) = D(i,j) for i = (1/2)*(n-1-k) and j = (1/2)*(n-1+k) where n-1+k belongs to 2Z. - Shiva Samieinia (shiva(AT)math.su.se), Oct 08 2007
Recurrence for triangle from A-sequence (see the Wolfdieter Lang comment above): T(n,k) = Sum_{j=0..n-k} A112478(j) * T(n-1, k-1+j), n >= 1, k >= 1. [For k > n, the sum is empty, in which case T(n,k) = 0.]
From Peter Bala, Jul 17 2008: (Start)
The n-th row of the square array is the crystal ball sequence for the product lattice A_1 x ... x A_1 (n copies). A035607 is the table of the associated coordination sequences for these lattices.
The polynomial p_n(x) := Sum {k = 0..n} 2^k * C(n,k) * C(x,k) = Sum_{k = 0..n} C(n,k) * C(x+k,n), whose values [p_n(0), p_n(1), p_n(2), ... ] give the n-th row of the square array, is the Ehrhart polynomial of the n-dimensional cross polytope (the hyperoctahedron) [Bump et al. (2000), Theorem 6].
The first few values are p_0(x) = 1, p_1(x) = 2*x + 1, p_2(x) = 2*x^2 + 2*x + 1 and p_3(x) = (4*x^3 + 6*x^2 + 8*x + 3)/3.
The reciprocity law p_n(m) = p_m(n) reflects the symmetry of the table.
The polynomial p_n(x) is the unique polynomial solution of the difference equation (x+1)*f(x+1) - x*f(x-1) = (2*n+1)*f(x), normalized so that f(0) = 1.
These polynomials have their zeros on the vertical line Re x = -1/2 in the complex plane; that is, the polynomials p_n(x-1), n = 1,2,3,..., satisfy a Riemann hypothesis [Bump et al. (2000), Theorem 4]. The o.g.f. for the p_n(x) is (1 + t)^x/(1 - t)^(x + 1) = 1 + (2*x + 1)*t + (2*x^2 + 2*x + 1)*t^2 + ... .
The square array of Delannoy numbers has a close connection with the constant log(2). The entries in the n-th row of the array occur in the series acceleration formula log(2) = (1 - 1/2 + 1/3 - ... + (-1)^(n+1)/n) + (-1)^n * Sum_{k>=1} (-1)^(k+1)/(k*D(n,k-1)*D(n,k)). [T(n,k) was replaced with D(n,k) in the formula to agree with the beginning of the paragraph. - Petros Hadjicostas, Aug 05 2020]
For example, the fourth row of the table (n = 3) gives the series log(2) = 1 - 1/2 + 1/3 - 1/(1*1*7) + 1/(2*7*25) - 1/(3*25*63) + 1/(4*63*129) - ... . See A142979 for further details.
Also the main diagonal entries (the central Delannoy numbers) give the series acceleration formula Sum_{n>=1} 1/(n*D(n-1,n-1)*D(n,n)) = (1/2)*log(2), a result due to Burnside. [T(n,n) was replaced here with D(n,n) to agree with the previous paragraphs. - Petros Hadjicostas, Aug 05 2020]
Similar relations hold between log(2) and the crystal ball sequences of the C_n lattices A142992. For corresponding results for the constants zeta(2) and zeta(3), involving the crystal ball sequences for root lattices of type A_n and A_n x A_n, see A108625 and A143007 respectively. (End)
From Peter Bala, Oct 28 2008: (Start)
Hilbert transform of Pascal's triangle A007318 (see A145905 for the definition of this term).
D(n+a,n) = P_n(a,0;3) for all integer a such that a >= -n, where P_n(a,0;x) is the Jacobi polynomial with parameters (a,0) [Hetyei]. The related formula A(n,k) = P_k(0,n-k;3) defines the table of asymmetric Delannoy numbers, essentially A049600. (End)
Seen as a triangle read by rows: T(n, k) = Hyper2F1([k-n, -k], [1], 2). - Peter Luschny, Aug 02 2014, Oct 13 2024.
From Peter Bala, Jun 25 2015: (Start)
O.g.f. for triangle T(n,k): A(z,t) = 1/(1 - (1 + t)*z - t*z^2) = 1 + (1 + t)*z + (1 + 3*t + t^2)*z^2 + (1 + 5*t + 5*t^2 + t^3)*z^3 + ....
1 + z*d/dz(A(z,t))/A(z,t) is the o.g.f. for A102413. (End)
E.g.f. for the n-th subdiagonal of T(n,k), n >= 0, equals exp(x)*P(n,x), where P(n,x) is the polynomial Sum_{k = 0..n} binomial(n,k)*(2*x)^k/k!. For example, the e.g.f. for the second subdiagonal is exp(x)*(1 + 4*x + 4*x^2/2) = 1 + 5*x + 13*x^2/2! + 25*x^3/3! + 41*x^4/4! + 61*x^5/5! + .... - Peter Bala, Mar 05 2017 [The n-th subdiagonal of triangle T(n,k) is the n-th row of array D(n,k).]
Let a_i(n) be multiplicative with a_i(p^e) = D(i, e), p prime and e >= 0, then Sum_{n > 0} a_i(n)/n^s = (zeta(s))^(2*i+1)/(zeta(2*s))^i for i >= 0. - Werner Schulte, Feb 14 2018
Seen as a triangle read by rows: T(n,k) = Sum_{i=0..k} binomial(n-i, i) * binomial(n-2*i, k-i) for 0 <= k <= n. - Werner Schulte, Jan 09 2019
Univariate generating function: Sum_{k >= 0} D(n,k)*z^k = (1 + z)^n/(1 - z)^(n+1). [Dziemianczuk (2013), Eq. 5.3] - Matt Majic, Nov 24 2019
(n+1)*D(n+1,k) = (2*k+1)*D(n,k) + n*D(n-1,k). [Majic (2019), Eq. 22] - Matt Majic, Nov 24 2019
For i, j >= 1, D(i,j) = D(i,j-1) + 2*Sum_{k=0..i-1} D(k,j-1), or, because D(i,j) = D(j,i), D(i,j) = D(i-1,j) + 2*Sum_{k=0..j-1} D(i-1,k). - Shel Kaphan, Jan 01 2023
Sum_{k=0..n} T(n,k)^2 = A026933(n). - R. J. Mathar, Nov 07 2023
Let S(x) = (1 - x - (1 - 6*x + x^2)^(1/2))/(2*x) denote the g.f. of the sequence of large Schröder numbers A006318. Read as a lower triangular array, the signed n-th row polynomial R(n, -x) = 1/sqrt(1 - 6*x + x^2) *( 1/S(x)^(n+1) + (x*S(x))^(n+1) ). For example, R(4, -x) = 1 - 7*x + 13*x^2 - 7*x^3 + x^4 = 1/sqrt(1 - 6*x + x^2) * ( 1/S(x)^5 + (x*S(x))^5 ). Cf. A102413. - Peter Bala, Aug 01 2024

Expanded description from Clark Kimberling, Jun 15 1997
Additional references from Sylviane R. Schwer (schwer(AT)lipn.univ-paris13.fr), Nov 28 2001
Changed the notation to make the formulas more precise. - N. J. A. Sloane, Jul 01 2002

A034856 a(n) = binomial(n+1, 2) + n - 1 = n*(n+3)/2 - 1.

Original entry on oeis.org

1, 4, 8, 13, 19, 26, 34, 43, 53, 64, 76, 89, 103, 118, 134, 151, 169, 188, 208, 229, 251, 274, 298, 323, 349, 376, 404, 433, 463, 494, 526, 559, 593, 628, 664, 701, 739, 778, 818, 859, 901, 944, 988, 1033, 1079, 1126, 1174, 1223, 1273, 1324, 1376, 1429, 1483

Offset: 1

N. J. A. Sloane

Number of 1's in the n X n lower Hessenberg (0,1)-matrix (i.e., the matrix having 1's on or below the superdiagonal and 0's above the superdiagonal).
If a 2-set Y and 2-set Z, having one element in common, are subsets of an n-set X then a(n-2) is the number of 3-subsets of X intersecting both Y and Z. - Milan Janjic, Oct 03 2007
Number of binary operations which have to be added to Moisil's algebras to obtain algebraic counterparts of n-valued Łukasiewicz logics. See the Wójcicki and Malinowski book, page 31. - Artur Jasinski, Feb 25 2010
Also (n + 1)!(-1)^(n + 1) times the determinant of the n X n matrix given by m(i,j) = i/(i+1) if i=j and otherwise 1. For example, (5+1)! * ((-1)^(5+1)) * Det[{{1/2, 1, 1, 1, 1}, {1, 2/3, 1, 1, 1}, {1, 1, 3/4, 1, 1}, {1, 1, 1, 4/5, 1}, {1, 1, 1, 1, 5/6}}] = 19 = a(5), and (6+1)! * ((-1)^(6+1)) * Det[{{1/2, 1, 1, 1, 1, 1}, {1, 2/3, 1, 1, 1, 1}, {1, 1, 3/4, 1, 1, 1}, {1, 1, 1, 4/5, 1, 1}, {1, 1, 1, 1, 5/6, 1}, {1, 1, 1, 1, 1, 6/7}}] = 26 = a(6). - John M. Campbell, May 20 2011
2*a(n-1) = n*(n+1) - 4, n>=0, with a(-1) = -2 and a(0) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 17 for b = 2*n + 1. In general D = b^2 - 4*a*c > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) is not divisible by 3, 5, 7, or 11. - Vladimir Shevelev, Feb 03 2014
With a(0) = 1 and a(1) = 2, a(n-1) is the number of distinct values of 1 +- 2 +- 3 +- ... +- n, for n > 0. - Derek Orr, Mar 11 2015
Also, numbers m such that 8*m+17 is a square. - Bruno Berselli, Sep 16 2015
Omar E. Pol's formula from Apr 23 2008 can be interpreted as the position of an element located on the third diagonal of an triangular array (read by rows) provided n > 1. - Enrique Pérez Herrero, Aug 29 2016
a(n) is the sum of the numerator and denominator of the fraction that is the sum of 2/(n-1) + 2/n; all fractions are reduced and n > 2. - J. M. Bergot, Jun 14 2017
a(n) is also the number of maximal irredundant sets in the (n+2)-path complement graph for n > 1. - Eric W. Weisstein, Apr 12 2018
From Klaus Purath, Dec 07 2020: (Start)
a(n) is not divisible by primes listed in A038890. The prime factors are given in A038889 and the prime terms of the sequence are listed in A124199.
Each odd prime factor p divides exactly 2 out of any p consecutive terms with the exception of 17, which appears only once in such an interval of terms. If a(i) and a(k) form such a pair that are divisible by p, then i + k == -3 (mod p), see examples.
If A is a sequence satisfying the recurrence t(n) = 5*t(n-1) - 2*t(n-2) with the initial values either A(0) = 1, A(1) = n + 4 or A(0) = -1, A(1) = n-1, then a(n) = (A(i)^2 - A(i-1)*A(i+1))/2^i for i>0. (End)
Mark each point on a 4^n grid with the number of points that are visible from the point; for n > 1, a(n) is the number of distinct values in the grid. - Torlach Rush, Mar 23 2021
The sequence gives the number of "ON" cells in the cellular automaton on a quadrant of a square grid after the n-th stage, where the "ON" cells lie only on the external perimeter and the perimeter of inscribed squares having the cell (1,1) as a unique common vertex. See Spezia link. - Stefano Spezia, May 28 2025

			From _Bruno Berselli_, Mar 09 2015: (Start)
By the definition (first formula):
----------------------------------------------------------------------
  1       4         8           13            19              26
----------------------------------------------------------------------
                                                              X
                                              X              X X
                                X            X X            X X X
                    X          X X          X X X          X X X X
          X        X X        X X X        X X X X        X X X X X
  X      X X      X X X      X X X X      X X X X X      X X X X X X
          X        X X        X X X        X X X X        X X X X X
----------------------------------------------------------------------
(End)
From _Klaus Purath_, Dec 07 2020: (Start)
Assuming a(i) is divisible by p with 0 < i < p and a(k) is the next term divisible by p, then from i + k == -3 (mod p) follows that k = min(p*m - i - 3) != i for any integer m.
(1) 17|a(7) => k = min(17*m - 10) != 7 => m = 2, k = 24 == 7 (mod 17). Thus every a(17*m + 7) is divisible by 17.
(2) a(9) = 53 => k = min(53*m - 12) != 9 => m = 1, k = 41. Thus every a(53*m + 9) and a(53*m + 41) is divisible by 53.
(3) 101|a(273) => 229 == 71 (mod 101) => k = min(101*m - 74) != 71 => m = 1, k = 27. Thus every a(101*m + 27) and a(101*m + 71) is divisible by 101. (End)
From _Omar E. Pol_, Aug 08 2021: (Start)
Illustration of initial terms:                             _ _
.                                           _ _           |_|_|_
.                              _ _         |_|_|_         |_|_|_|_
.                   _ _       |_|_|_       |_|_|_|_       |_|_|_|_|_
.          _ _     |_|_|_     |_|_|_|_     |_|_|_|_|_     |_|_|_|_|_|_
.   _     |_|_|    |_|_|_|    |_|_|_|_|    |_|_|_|_|_|    |_|_|_|_|_|_|
.  |_|    |_|_|    |_|_|_|    |_|_|_|_|    |_|_|_|_|_|    |_|_|_|_|_|_|
.
.   1       4         8          13            19              26
------------------------------------------------------------------------ (End)

A. S. Karpenko, Łukasiewicz's Logics and Prime Numbers, 2006 (English translation).
G. C. Moisil, Essais sur les logiques non-chrysippiennes, Ed. Academiei, Bucharest, 1972.
Wójcicki and Malinowski, eds., Łukasiewicz Sentential Calculi, Wrocław: Ossolineum, 1977.

Vincenzo Librandi, Table of n, a(n) for n = 1..10000
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 471.
Milan Janjic, Two Enumerative Functions.
W. F. Klostermeyer, M. E. Mays, L. Soltes, and G. Trapp, A Pascal rhombus, Fibonacci Quarterly, Vol. 35, No. 4 (1997), pp. 318-328.
Stavros Konstantinidis, António Machiavelo, Nelma Moreira, and Rogério Reis, On the size of partial derivatives and the word membership problem, Acta Informatica Vol. 58, 357-375.
G. C. Moisil, Recherches sur les logiques non-chrysippiennes, Ann. Sci. Univ. Jassy, 26 (1940), 431-466.
László Németh, Tetrahedron trinomial coefficient transform, arXiv:1905.13475 [math.CO], 2019.
D. D. Olesky, B. L. Shader, and P. van den Driessche, Permanents of Hessenberg (0,1)-matrices, Electronic Journal of Combinatorics, 12 (2005), #R70.
J. Riordan, Enumeration of trees by height and diameter, IBM J. Res. Dev., Vol. 4, No. 5 (1960), pp. 473-478.
Stefano Spezia, Illustrations for n = 1..9.
Renzo Sprugnoli, Alternating Weighted Sums of Inverses of Binomial Coefficients, J. Integer Sequences, 15 (2012), #12.6.3. - From _N. J. A. Sloane_, Nov 29 2012
Eric Weisstein's World of Mathematics, Maximal Irredundant Set.
Eric Weisstein's World of Mathematics, Path Complement Graph.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
Index entries for sequences related to Łukasiewicz.

Subsequence of A165157.
Triangular numbers (A000217) minus two.
Third diagonal of triangle in A059317.
Cf. A000096, A000124, A000217, A000290, A005408, A005843, A008277, A024206, A027379.
Cf. A028347, A038889, A038890, A049600, A091823, A101881, A113452, A113453, A113454.
Cf. A113455, A124199, A130296, A131818, A134225, A161680, A168244, A253909, A267633.

a034856 = subtract 1 . a000096 -- Reinhard Zumkeller, Feb 20 2015

[Binomial(n + 1, 2) + n - 1: n in [1..60]]; // Vincenzo Librandi, May 21 2011

a := n -> hypergeom([-2, n-1], [1], -1);
seq(simplify(a(n)), n=1..53); # Peter Luschny, Aug 02 2014

f[n_] := n (n + 3)/2 - 1; Array[f, 55] (* or *) k = 2; NestList[(k++; # + k) &, 1, 55] (* Robert G. Wilson v, Jun 11 2010 *)
Table[Binomial[n + 1, 2] + n - 1, {n, 53}] (* or *)
Rest@ CoefficientList[Series[x (1 + x - x^2)/(1 - x)^3, {x, 0, 53}], x] (* Michael De Vlieger, Aug 29 2016 *)

A034856(n) := block(
        n-1+(n+1)*n/2
)$ /* R. J. Mathar, Mar 19 2012 */

A034856(n)=(n+3)*n\2-1 \\ M. F. Hasler, Jan 21 2015

def A034856(n): return n*(n+3)//2 -1 # G. C. Greubel, Jun 15 2025

G.f.: A(x) = x*(1 + x - x^2)/(1 - x)^3.
a(n) = A049600(3, n-2).
a(n) = binomial(n+2, 2) - 2. - Paul Barry, Feb 27 2003
With offset 5, this is binomial(n, 0) - 2*binomial(n, 1) + binomial(n, 2), the binomial transform of (1, -2, 1, 0, 0, 0, ...). - Paul Barry, Jul 01 2003
Row sums of triangle A131818. - Gary W. Adamson, Jul 27 2007
Binomial transform of (1, 3, 1, 0, 0, 0, ...). Also equals A130296 * [1,2,3,...]. - Gary W. Adamson, Jul 27 2007
Row sums of triangle A134225. - Gary W. Adamson, Oct 14 2007
a(n) = A000217(n+1) - 2. - Omar E. Pol, Apr 23 2008
From Jaroslav Krizek, Sep 05 2009: (Start)
a(n) = a(n-1) + n + 1 for n >= 1.
a(n) = n*(n-1)/2 + 2*n - 1.
a(n) = A000217(n-1) + A005408(n-1) = A005843(n-1) + A000124(n-1). (End)
a(n) = Hyper2F1([-2, n-1], [1], -1). - Peter Luschny, Aug 02 2014
a(n) = floor[1/(-1 + Sum_{m >= n+1} 1/S2(m,n+1))], where S2 is A008277. - Richard R. Forberg, Jan 17 2015
a(n) = A101881(2*(n-1)). - Reinhard Zumkeller, Feb 20 2015
a(n) = A253909(n+3) - A000217(n+3). - David Neil McGrath, May 23 2015
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>3. - David Neil McGrath, May 23 2015
For n > 1, a(n) = 4*binomial(n-1,1) + binomial(n-2,2), comprising the third column of A267633. - Tom Copeland, Jan 25 2016
From Klaus Purath, Dec 07 2020: (Start)
a(n) = A024206(n) + A024206(n+1).
a(2*n-1) = -A168244(n+1).
a(2*n) = A091823(n). (End)
Sum_{n>=1} 1/a(n) = 3/2 + 2*Pi*tan(sqrt(17)*Pi/2)/sqrt(17). - Amiram Eldar, Jan 06 2021
a(n) + a(n+1) = A028347(n+2). - R. J. Mathar, Mar 13 2021
a(n) = A000290(n) - A161680(n-1). - Omar E. Pol, Mar 26 2021
E.g.f.: 1 + exp(x)*(x^2 + 4*x - 2)/2. - Stefano Spezia, Jun 05 2021
a(n) = A024916(n) - A244049(n). - Omar E. Pol, Aug 01 2021
a(n) = A000290(n) - A000217(n-2). - Omar E. Pol, Aug 05 2021

More terms from Zerinvary Lajos, May 12 2006

A100258 Triangle of coefficients of normalized Legendre polynomials, with increasing exponents.

Original entry on oeis.org

1, 0, 1, -1, 0, 3, 0, -3, 0, 5, 3, 0, -30, 0, 35, 0, 15, 0, -70, 0, 63, -5, 0, 105, 0, -315, 0, 231, 0, -35, 0, 315, 0, -693, 0, 429, 35, 0, -1260, 0, 6930, 0, -12012, 0, 6435, 0, 315, 0, -4620, 0, 18018, 0, -25740, 0, 12155, -63, 0, 3465, 0, -30030, 0, 90090, 0, -109395, 0, 46189

Offset: 0

Ralf Stephan, Nov 13 2004

For a relation to Jacobi quartic elliptic curves, see the MathOverflow link. For a self-convolution of the polynomials relating them to the Chebyshev and Fibonacci polynomials, see A049310 and A053117. For congruences and connections to other polynomials (Jacobi, Gegenbauer, and Chebyshev) see the Allouche et al. link. For relations to elliptic cohomology and modular forms, see references in Copeland link.- Tom Copeland, Feb 04 2016

			Triangle begins:
   1;
   0,   1;
  -1,   0,     3;
   0,  -3,     0,   5;
   3,   0,   -30,   0,   35;
   0,  15,     0, -70,    0,   63;
  -5,   0,   105,   0, -315,    0,    231;
   0, -35,     0, 315,    0, -693,      0, 429;
  35,   0, -1260,   0, 6930,    0, -12012,   0, 6435;
  ...

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 798.

T. D. Noe, Rows n=0..100 of triangle, flattened
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].
J. Allouche and G. Skordev, Schur congruences, Carlitz sequences of polynomials and automaticity, Discrete Mathematics, Vol. 214, Issue 1-3, 21 March 2000, p. 21-49.
Tom Copeland, The Elliptic Lie Triad: Riccati and KdV Equations, Infinigens, and Elliptic Genera
H. N. Laden, An historical, and critical development of the theory of Legendre polynomials before 1900, Master of Arts Thesis, University of Maryland 1938.
Shi-Mei Ma, On gamma-vectors and the derivatives of the tangent and secant functions, arXiv:1304.6654 [math.CO], 2013.
MathOverflow, Geometric picture of invariant differential of an elliptic curve, Dec 4 2011.

Without zeros: A008316. Row sums are A060818.
Columns (with interleaved zeros and signs) include A001790, A001803, A100259. Diagonals include A001790, A001800, A001801, A001802.
Cf. A049310, A005187, A049600, A053117.

row[n_] := CoefficientList[ LegendreP[n, x], x]*2^IntegerExponent[n!, 2]; Table[row[n], {n, 0, 10}] // Flatten (* Jean-François Alcover, Jan 15 2015 *)

a(k,n)=polcoeff(pollegendre(k,x),n)*2^valuation(k!,2)

from mpmath import *
mp.dps=20
def a007814(n):
    return 1 + bin(n - 1)[2:].count('1') - bin(n)[2:].count('1')
for n in range(11):
    y=2**sum(a007814(i) for i in range(2, n+1))
    l=chop(taylor(lambda x: legendre(n, x), 0, n))
    print([int(i*y) for i in l]) # Indranil Ghosh, Jul 02 2017

The n-th normalized Legendre polynomial is generated by 2^(-n-a(n)) (d/dx)^n (x^2-1)^n / n! with a(n) = A005187(n/2) for n even and a(n) = A005187((n-1)/2) for n odd. The non-normalized polynomials have the o.g.f. 1 / sqrt(1 - 2xz + z^2). - Tom Copeland, Feb 07 2016

The consecutive nonzero entries in the m-th row are, in order, (c+b)!/(c!(m-b)!(2b-m)!*A048896(m-1)) with sign (-1)^b where c = m/2-1, m/2, m/2+1, ..., (m-1) and b = c+1 if m is even and sign (-1)^c with c = (m-1)/2, (m-1)/2+1, (m-1)/2+2, ..., (m-1) with b = c+1 if m is odd. For the 9th row the 5 consecutive nonzero entries are 315, -4620, 18018, -25740, 12155 given by c = 4,5,6,7,8 and b = 5,6,7,8,9. - Richard Turk, Aug 22 2017

A002003 a(n) = 2 * Sum_{k=0..n-1} binomial(n-1, k)*binomial(n+k, k).

Original entry on oeis.org

0, 2, 8, 38, 192, 1002, 5336, 28814, 157184, 864146, 4780008, 26572086, 148321344, 830764794, 4666890936, 26283115038, 148348809216, 838944980514, 4752575891144, 26964373486406, 153196621856192, 871460014012682, 4962895187697048, 28292329581548718

Offset: 0

N. J. A. Sloane

nonn

a(n) is the number of order-preserving partial self maps of {1,...,n}. For example, a(2) = 8 because there are 8 order-preserving partial self maps of {1,2}: (1 2), (1 1), (2 2), (1 -), (2 -), (- 1), (- 2), (- -). Here for example (2 -) represents the partial map which maps 1 to 2 but does not include 2 in its domain. - James East, Oct 25 2005
From Peter Bala, Mar 02 2020: (Start)
For fixed m = 1,2,3,..., we conjecture that the sequence b(n) := a(m*n) satisfies a recurrence of the form P(2*m,n)*b(n+1) + P(2*m,-n)*b(n-1) = Q(2*m,n)*b(n), where the polynomials P(2*m,n) and Q(2*m,n) have degree 2*m. Conjecturally, the polynomial Q(2*m,n) is an even function of n; its 2*m zeros seem to belong to the interval [-1, 1] and 2*m - 2 of these zeros appear to lie close to the rational numbers of the form +-(2*k + 1)/(2*m), where 0 <= k <= m - 2. Cf. A103885. (End)
a(n), n>0, is the number of points at L1 distance = n from any given point in Z^n.  The sequence is also the difference between the central diagonal (A001850) and +-1 diagonal (A002002) of the Delannoy number triangle (A008288).  - Shel Kaphan, Feb 15 2023

			G.f. = 2*x + 8*x^2 + 38*x^3 + 192*x^4 + 1002*x^5 + 5336*x^6 + 28814*x^7 + ...

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).

Seiichi Manyama, Table of n, a(n) for n = 0..1000 (terms 0..200 from T. D. Noe)
Peter Bala, A supercongruence for A002003
J. Brzozowski, M. Szykula, Large Aperiodic Semigroups, arXiv preprint arXiv:1401.0157 [cs.FL], 2013-2014.
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
G. Rutledge and R. D. Douglass, Integral functions associated with certain binomial coefficient sums, Amer. Math. Monthly, 43 (1936), 27-32.

Cf. A002002, A006318, A027307, A108424, A103885.

A064861 := proc(n,k) option remember; if n = 1 then 1; elif k = 0 then 0; else A064861(n,k-1)+(3/2-1/2*(-1)^(n+k))*A064861(n-1,k); fi; end; seq(A064861(i,i-1),i=1..40);

Flatten[{0,Table[SeriesCoefficient[((1+x)/Sqrt[1-6*x+x^2]-1)/2,{x,0,n}],{n,1,20}]}] (* Vaclav Kotesovec, Oct 04 2012 *)
a[ n_] := If[ n < 1, 0, Hypergeometric2F1[ n, -n, 1, -1]]; (* Michael Somos, Aug 24 2014 *)
Table[2*Sum[Binomial[n-1,k]Binomial[n+k,k],{k,0,n-1}],{n,0,30}] (* Harvey P. Dale, Sep 18 2024 *)

{a(n) = if( n<1, 0, polcoeff( ((1 - x^2) / (1 - x)^2 + x * O(x^n))^n, n))} /* Michael Somos, Sep 24 2003 */

from math import comb
def A002003(n): return sum(comb(n,k)**2*k<Chai Wah Wu, Mar 22 2023

a(n) = 2*A047781(n).
From Vladeta Jovovic, Mar 28 2004: (Start)
G.f.: ((1+x)/sqrt(1-6*x+x^2)-1)/2.
E.g.f.: exp(3*x)*(2*BesselI(0, 2*sqrt(2)*x)+sqrt(2)*BesselI(1, 2*sqrt(2)*x)). (End)
a(n) = T(n, n-1), array T as in A064861.
a(n) = T(n, n-2), array T as in A049600.
a(n+1) = A110110(2n+1). - Tilman Neumann, Feb 05 2009
a(n) = 2 * JacobiP(n-1,0,1,3) = ((7*n+3)*LegendreP(n,3) - (n+1)*LegendreP(n+1,3)) /(2*n) for n > 0. - Mark van Hoeij, Jul 12 2010
Logarithmic derivative of A006318, the large Schroeder numbers. - Paul D. Hanna, Oct 25 2010
D-finite with recurrence: 4*(3*n^2-6*n+2)*a(n-1) - (n-2)*(2*n-1)*a(n-2) - n*(2*n-3)*a(n)=0. - Vaclav Kotesovec, Oct 04 2012
a(n) ~ (3+2*sqrt(2))^n/(2^(3/4)*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 04 2012
Recurrence (an alternative): n*a(n) = (6-n)*a(n-6) + 2*(5*n-27)*a(n-5) + (84-15*n)*a(n-4) + 52*(3-n)*a(n-3) + 3*(2-5*n)*a(n-2) + 2*(5*n-3)*a(n-1), n>=7. - Fung Lam, Feb 05 2014
a(n) = Hyper2F1([-n, n], [1], -1) for n > 0. - Peter Luschny, Aug 02 2014
a(n) = [x^n] ((1+x)/(1-x))^n for n > 0. - Seiichi Manyama, Jun 07 2018
From Peter Bala, Mar 13 2020: (Start)
a(n) = 2 * Sum_{k = 0..n-1}  2^k*C(n,k+1)*C(n-1,k).
a(n) = 2 * (-1)^(n+1) * Sum_{k = 0..n-1} (-2)^k*C(n+k,n-1)*C(n-1,k).
a(n) = Sum_{k = 0..n} C(n,k)*C(2*n-k-1,n-1).
Conjecture: a(n) = - [x^n] (1 - F(x))^n, where F(x) = 2*x + 6*x^2 + 34*x^3 + 238*x^4 + ... is the o.g.f. of A108424. Equivalently, a(n) = -[x^n](G(x))^(-n), where G(x) = 1 + 2*x + 10*x^2 + 66*x^3 + 498*x^4 + ... is the o.g.f. of A027307.
a(p) == 2 ( mod p^3 ) for prime p >= 5. (End)
a(n) = Sum_{k = 1..n} C(n, k) * C(n-1, k-1) * 2^k. - Michael Somos, May 23 2021
a(n) = A001850(n) - A002002(n), for n > 0. - Shel Kaphan, Feb 15 2023

More terms from Barbara Haas Margolius (b.margolius(AT)csuohio.edu), Oct 10 2001

A142978 Table of figurate numbers for the n-dimensional cross polytopes.

Original entry on oeis.org

1, 1, 2, 1, 4, 3, 1, 6, 9, 4, 1, 8, 19, 16, 5, 1, 10, 33, 44, 25, 6, 1, 12, 51, 96, 85, 36, 7, 1, 14, 73, 180, 225, 146, 49, 8, 1, 16, 99, 304, 501, 456, 231, 64, 9, 1, 18, 129, 476, 985, 1182, 833, 344, 81, 10

Offset: 1

Peter Bala, Jul 15 2008

The n-th row entries for this array are the regular polytope numbers for the n-dimensional cross polytope as defined by [Kim]. The rows are the partial sums of the rows of the square array of Delannoy numbers A008288.
The odd numbered rows of this array form A142977. For a triangular version of this table see A104698. Cf. also A101603.
The n-th row of the array is the binomial transform of n-th row of triangle A081277, followed by zeros. Example: row 4 (1, 6, 19, 44, 85, ...) = binomial transform of row 3 of A081277: (1, 5, 8, 4, 0, 0, 0, ...). - Gary W. Adamson, Jul 17 2008
The main diagonal of the array T(n,k) is A047781 Sum_{k=0..n-1} binomial(n-1,k)*binomial(n+k,k). Also a(n) = T(n,n), array T as in A049600. The link from A099193 to J. V. Post, Table of polytope numbers, Sorted, Through 1,000,000, includes all n-D Hyperoctahedron (n-Cross Polytope) Numbers through 10-Cross(20) = 1669752016. - Jonathan Vos Post, Jul 16 2008

			The square array A(n, k) begins:
  n\k| 1   2    3     4     5       6
  ---+-------------------------------
   1 | 1   2    3     4      5      6    A000027
   2 | 1   4    9    16     25     36    A000290
   3 | 1   6   19    44     85    146    A005900
   4 | 1   8   33    96    225    456    A014820
   5 | 1  10   51   180    501   1182    A069038
   6 | 1  12   73   304    985   2668    A069039
   7 | 1  14   99   476   1765   5418    A099193

Reinhard Zumkeller, Rows n = 1..100 of triangle, flattened
Peter Bala, A142978 and a family of continued fraction expansions for log(2)
Valentin Bonzom and Etera R. Livine, Self-duality of the 6j-symbol and Fisher zeros for the Tetrahedron, arXiv:1905.00348 [math-ph], 2019.
Steven Edwards and William Griffiths, Generalizations of Delannoy and cross polytope numbers, Fib. Q., Vol. 55, No. 4 (2017), pp. 356-366.
Steven Edwards and William Griffiths, On Generalized Delannoy Numbers, J. Int. Seq., Vol. 23 (2020), Article 20.3.6.
Milan Janjic and B. Petkovic, A Counting Function, arXiv 1301.4550 [math.CO], 2013.
Milan Janjic and B. Petkovic, A Counting Function Generalizing Binomial Coefficients and Some Other Classes of Integers, J. Int. Seq. 17 (2014) # 14.3.5.
Hyun Kwang Kim, On Regular Polytope Numbers, Proc. Amer. Math. Soc., 131 (2002), 65-75.
Eric W. Weisstein, Meixner polynomial of the first kind
Eric W. Weisstein, Mittag-Leffler polynomial
Yutong Yang, From Simplest Recursion to the Recursion of Generalizations of Cross Polytope Numbers, (2017), Honors College Capstone and Theses, 13.

Cf. A008288 (Delannoy numbers), A005900 (row 3), A014820 (row 4), A069038 (row 5), A069039 (row 6), A099193 (row 7), A099195 (row 8), A099196 (row 9), A099197 (row 10), A101603, A104698 (triangle version), A142977, A142983.

Cf. A049600, A047781, A081277.

a142978 n k = a142978_tabl !! (n-1) !! (k-1)
a142978_row n = a142978_tabl !! (n-1)
a142978_tabl = map reverse a104698_tabl
-- Reinhard Zumkeller, Jul 17 2015

with(combinat): T:=(n,k) -> add(binomial(n-1,i)*binomial(k+i,n),i = 0..n-1); for n from 1 to 10 do seq(T(n,k),k = 1..10) end do; # Program restored by Peter Bala, Oct 02 2008
A := (n, k) -> k*hypergeom([1 - n, 1 - k], [2], 2):
seq(print(seq(simplify(A(n, k)), k = 1..9)), n=1..7); # Peter Luschny, Mar 23 2023

t[n_, k_] := Sum[ Binomial[n-1, i]*Binomial[k+i, n], {i, 0, n-1}]; Table[t[n-k, k], {n, 1, 11}, {k, 1, n-1}] // Flatten (* Jean-François Alcover, Mar 06 2013 *)

T(n,k) = Sum_{i = 0..n-1} C(n-1,i)*C(k+i,n).
Reciprocity law: n*T(n,k) = k*T(k,n).
Recurrence relation: T(n,1) = 1, T(1,k) = k, T(n,k) = T(n,k-1) + T(n-1,k-1) + T(n-1,k), n,k > 1.
O.g.f. row n: x*(1 + x)^(n-1)/(1 - x)^(n+1).
O.g.f. for array: Sum_{n >= 1, k >= 1} T(n, k)*x^k*y^n = x*y/((1 - x)*(1 - x - y - x*y)).
The n-th row entries are the values [p_n(k)], k >= 1, of the polynomial function p_n(x) = Sum_{k = 1..n} 2^(k-1)*C(n-1,k-1)*C(x,k). The first few values are p_1(x) = x, p_2(x) = x^2, p_3(x) = (2*x^3 + x)/3 and p_4(x) = (x^4 + 2*x^2)/3.
The polynomial p_n(x) is the unique polynomial solution of the difference equation x*( f(x+1) - f(x-1) ) = 2*n*f(x), normalized so that f(1) = 1.
The o.g.f. for the p_n(x) is 1/2*((1 + t)/(1 - t))^x = 1/2 + x*t + x^2*t^2 + (2*x^3 + x)/3*t^3 + .... Thus p_n(x) is, apart from a constant factor, the Meixner polynomial of the first kind M_n(x;b,c) at b = 0, c = -1, also known as a Mittag-Leffler polynomial.
The entries in the n-th row appear in the series acceleration formula for the constant log(2): Sum_{k >= 1} (-1)^(k+1)/(T(n,k)*T(n,k+1)) = 1 + (-1)^(n+1) * (2*n)*(log(2) - (1 - 1/2 + 1/3 - ... + (-1)^(n+1)/n)). For example, n = 3 gives log(2) = 4/6 + (1/6)*(1/(1*6) - 1/(6*19) + 1/(19*44) - 1/(44*85) + ...). See A142983 for further details.
From Peter Bala, Oct 02 2008: (Start)
The odd-indexed columns of this array form the array A142992 of crystal ball sequences for lattices of type C_n.
Conjectural congruences for main diagonal entries: Put A(n) = T(n,n). Calculation suggests the following congruences: for prime p > 3 and m, r >= 1, A(m*p^r) == A(m*p^(r-1)) (mod p^(3*r));
Sum_{k = 0..p-1} A(k)^2 == 0 (mod p) if p is a prime of the form 8*n+1 or 8*n+7;
Sum_{k = 0..p-1} A(k)^2 == -1 (mod p) if p is a prime of the form 8*n+3 or 8*n+5.
(End)
From Peter Bala, Sep 27 2021: (Start)
T(n,k) = (1/2)*Sum_{i = 0..k} binomial(k,i)*binomial(n+k-1-i,k-1).
T(n,k) = (1/2)*[x^n] ((1+x)/(1-x))^k = (1/2)*(k/n)*[x^k] ((1+x)/(1-x))^n.
n*T(n,k) = 2*k*T(n-1,k) + (n - 2)*T(n-2,k). (End)
A(n,k) = k*hypergeom([1 - n, 1 - k], [2], 2). - Peter Luschny, Mar 23 2023
T(n,k) = 2*(Sum_{j=1..k-1} T(n-1,j)) + T(n-1,k) for n > 1. - Robert FERREOL, Jun 25 2024

A059576 Summatory Pascal triangle T(n,k) (0 <= k <= n) read by rows. Top entry is 1. Each entry is the sum of the parallelogram above it.

Original entry on oeis.org

1, 1, 1, 2, 3, 2, 4, 8, 8, 4, 8, 20, 26, 20, 8, 16, 48, 76, 76, 48, 16, 32, 112, 208, 252, 208, 112, 32, 64, 256, 544, 768, 768, 544, 256, 64, 128, 576, 1376, 2208, 2568, 2208, 1376, 576, 128, 256, 1280, 3392, 6080, 8016, 8016, 6080, 3392, 1280, 256

Offset: 0

Floor van Lamoen, Jan 23 2001

We may also relabel the entries as U(0,0), U(1,0), U(0,1), U(2,0), U(1,1), U(0,2), U(3,0), ... [That is, T(n,k) = U(n-k, k) for 0 <= k <= n and U(m,s) = T(m+s, s) for m,s >= 0.]
From Petros Hadjicostas, Jul 16 2020: (Start)
We explain the parallelogram definition of T(n,k).
    T(0,0) *
           |\
           | \
           |  * T(k,k)
  T(n-k,0) *  |
            \ |
             \|
              * T(n,k)
The definition implies that T(n,k) is the sum of all T(i,j) such that (i,j) has integer coordinates over the set
{(i,j): a(1,0) + b(1,1), 0 <= a <= n-k, 0 <= b <= k} - {(n,k)}.
The parallelogram can sometimes be degenerate; e.g., when k = 0 or n = k. (End)
T(n,k) is the number of 2-compositions of n having sum of the entries of the first row equal to k (0 <= k <= n). A 2-composition of n is a nonnegative matrix with two rows, such that each column has at least one nonzero entry and whose entries sum up to n. - Emeric Deutsch, Oct 12 2010
From Michel Marcus and Petros Hadjicostas, Jul 16 2020: (Start)
Robeva and Sun (2020) let A(m,n) = U(m-1, n-1) be the number of subdivisions of a 2-row grid with m points on the top and n points at the bottom (and such that the lower left point is the origin).
The authors proved that A(m,n) = 2*(A(m,n-1) + A(m-1,n) - A(m-1,n-1)) for m, n >= 2 (with (m,n) <> (2,2)), which is equivalent to a similar recurrence for U(n,k) given in the Formula section below. (They did not explicitly specify the value of A(1,1) = U(0,0) because they did not care about the number of subdivisions of a degenerate polygon with only one side.)
They also proved that, for (m,n) <> (1,1), A(m,n) = (2^(m-2)/(n-1)!) * Q_n(m) =
= (2^(m-2)/(n-1)!) * Sum_{k=1..n} A336244(n,k) * m^(n-k), where Q_n(m) is a polynomial in m of degree n-1. (End)
With the square array notation of Petros Hadjicostas, Jul 16 2020 below, U(i,j) is the number of lattice paths from (0,0) to (i,j) whose steps move north or east or have positive slope. For example, representing a path by its successive lattice points rather than its steps, U(1,2) = 8 counts {(0,0),(1,2)}, {(0,0),(0,1),(1,2)}, {(0,0),(0,2),(1,2)}, {(0,0),(1,0),(1,2)}, {(0,0),(1,1),(1,2)}, {(0,0),(0,1),(0,2),(1,2)}, {(0,0),(0,1),(1,1),(1,2)}, {(0,0),(1,0),(1,1),(1,2)}. If north (vertical) steps are excluded, the resulting paths are counted by A049600.  - David Callan, Nov 25 2021

			Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins
[0]   1;
[1]   1,   1;
[2]   2,   3,   2;
[3]   4,   8,   8,   4;
[4]   8,  20,  26,  20,   8;
[5]  16,  48,  76,  76,  48,  16;
[6]  32, 112, 208, 252, 208, 112, 32;
  ...
T(5,2) = 76 is the sum of the elements above it in the parallelogram bordered by T(0,0), T(5-2,0) = T(3,0), T(2,2) and T(5,2). We of course exclude T(5,2) from the summation. Thus
T(5,2) = Sum_{a=0..5-2, b=0..2, (a,b) <> (5-2,2)} T(a(1,0) + b(1,1)) =
= (1 + 1 + 2) + (1 + 3 + 8) + (2 + 8 + 26) + (4 + 20) = 76. [Edited by _Petros Hadjicostas_, Jul 16 2020]
From _Petros Hadjicostas_, Jul 16 2020: (Start)
Square array U(n,k) (with rows n >= 0 and columns k >= 0) begins
   1,   1,   2,    4,    8, ...
   1,   3,   8,   20,   48, ...
   2,   8,  26,   76,  208, ...
   4,  20,  76,  252,  768, ...
   8,  48, 208,  768, 2568, ...
  16, 112, 544, 2208, 8016, ...
  ...
Consider the following 2-row grid with n = 3 points at the top and k = 2 points at the bottom:
   A  B  C
   *--*--*
   |    /
   |   /
   *--*
   D  E
The sets of the dividing internal lines of the A(3,2) = U(3-1, 2-1) = 8 subdivisions of the above 2-row grid are as follows: { }, {DC}, {DB}, {EB}, {EA}, {DB, DC}, {DB, EB}, and {EA, EB}. See Robeva and Sun (2020).
These are the 2-compositions of n = 3 with sum of first row entries equal to k = 1:
[1; 2], [0,1; 2,0], [0,1; 1,1], [1,0; 0,2], [1,0; 1,1], [0,0,1; 1,1,0], [0,1,0; 1,0,1], and [1,0,0; 0,1,1]. We have T(3,2) = 8 such matrices. See _Emeric Deutsch_'s contribution above. See also Section 2 in Castiglione et al. (2007). (End)

Reinhard Zumkeller, Rows n = 0..120 of triangle, flattened
G. Castiglione, A. Frosini, E. Munarini, A. Restivo and S. Rinaldi, Combinatorial aspects of L-convex polyominoes, European Journal of Combinatorics, 28(6) (2007), 1724-1741; see Fig. 5, p. 1729.
Fang, E., Jenkins, J., Lee, Z., Li, D., Lu, E., Miller, S. J., ... & Siktar, J. (2019). Central Limit Theorems for Compound Paths on the 2-Dimensional Lattice, arXiv preprint arXiv:1906.10645. Also Fib. Q., 58:1 (2020), 208-225.
Elina Robeva and Melinda Sun, Bimonotone Subdivisions of Point Configurations in the Plane, arXiv:2007.00877 [math.CO], 2020.
Index entries for triangles and arrays related to Pascal's triangle

Cf. A000079, A001792, A003480 (row sums), A052141 (main diagonal).
Cf. A059474, A059473, A008288, A035002, A059226, A059283, A181292, A336244.
Cf. A011782, A049611, A049612, A055589, A055852, A055853, A181306.

a059576 n k = a059576_tabl !! n !! k
a059576_row n = a059576_tabl !! n
a059576_tabl = [1] : map fst (iterate f ([1,1], [2,3,2])) where
   f (us, vs) = (vs, map (* 2) ws) where
     ws = zipWith (-) (zipWith (+) ([0] ++ vs) (vs ++ [0]))
                      ([0] ++ us ++ [0])
-- Reinhard Zumkeller, Dec 03 2012

A011782:= func< n | n eq 0 select 1 else 2^(n-1) >;
function T(n,k) // T = A059576
  if k eq 0 or k eq n then return A011782(n);
  else return 2*T(n-1, k-1) + 2*T(n-1, k) - (2 - 0^(n-2))*T(n-2, k-1);
  end if; return T;
end function;
[T(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Sep 02 2022

A059576 := proc(n,k) local b,t1; t1 := min(n+k-2,n,k); add( (-1)^b * 2^(n+k-b-2) * (n+k-b-2)! * (1/(b! * (n-b)! * (k-b)!)) * (-2 * n-2 * k+2 * k^2+b^2-3 * k * b+2 * n^2+5 * n * k-3 * n * b), b=0..t1); end;
T := proc (n, k) if k <= n then add((-1)^j*2^(n-j-1)*binomial(k, j)*binomial(n-j, k), j = 0 .. min(k, n-k)) fi end proc: 1; for n to 10 do seq(T(n, k), k = 0 .. n) end do; # yields sequence in triangular form # Emeric Deutsch, Oct 12 2010
T := (n, k) -> `if`(n=0, 1, 2^(n-1)*binomial(n, k)*hypergeom([-k, k - n], [-n], 1/2)): seq(seq(simplify(T(n, k)), k=0..n), n=0..10); # Peter Luschny, Nov 26 2021

T[0, 0] = 1; T[n_, k_] := 2^(n-k-1)*n!*Hypergeometric2F1[ -k, -k, -n, -1 ] / (k!*(n-k)!); Flatten[ Table[ T[n, k], {n, 0, 9}, {k, 0, n}]] (* Jean-François Alcover, Feb 01 2012, after Robert Israel *)

def T(n,k): # T = A059576
    if (k==0 or k==n): return 1 if (n==0) else 2^(n-1) # A011782
    else: return 2*T(n-1, k-1) + 2*T(n-1, k) - (2 - 0^(n-2))*T(n-2, k-1)
flatten([[T(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Sep 02 2022

T(n, n-1) = A001792(n-1).
T(2*n, n) = A052141(n).
Sum_{k=0..n} T(n, k) = A003480(n).
G.f.: U(z, w) = Sum_{n >= 0, k >= 0} U(n, k)*z^n*w^k = Sum{n >= 0, k >= 0} T(n, k)*z^(n-k)*w^k = (1-z)*(1-w)/(1 - 2*w - 2*z + 2*z*w).
Maple code gives another explicit formula for U(n, k).
From Jon Stadler (jstadler(AT)capital.edu), Apr 30 2003: (Start)
U(n,k) is the number of ways of writing the vector (n,k) as an ordered sum of vectors, equivalently, the number of paths from (0,0) to (n,k) in which steps may be taken from (i,j) to (p,q) provided (p,q) is to the right or above (i,j).
2*U(n,k) = Sum_{i <= n, j <= k} U(i,j).
U(n,k) = 2*U(n-1,k) + Sum_{i < k} U(n,i).
U(n,k) = Sum_{j=0..n+k} C(n,j-k+1)*C(k,j-n+1)*2^j. (End)
T(n, k) = 2*(T(n-1, k-1) + T(n-1, k)) - (2 - 0^(n-2))*T(n-2, k-1) for n > 1 and 1 < k < n; T(n, 0) = T(n, n) = 2*T(n-1, 0) for n > 0; and T(0, 0) = 1. - Reinhard Zumkeller, Dec 03 2004
From Emeric Deutsch, Oct 12 2010: (Start)
Sum_{k=0..n} k*T(n,k) = A181292(n).
T(n,k) = Sum_{j=0..min(k, n-k)} (-1)^j*2^(n-j-1)*binomial(k, j)*binomial(n-j, k) for (n,k) != (0,0).
G.f.: G(t,z) = (1-z)*(1-t*z)/(1 - 2*z - 2*t*z + 2*t*z^2). (End)
U(n,k) = 0 if k < 0; else U(k,n) if k > n; else 1 if n <= 1; else 3 if n = 2 and k = 1; else 2*U(n,k-1) + 2*U(n-1,k) - 2*U(n-1,k-1). - David W. Wilson; corrected in the case k > n by Robert Israel, Jun 15 2011 [Corrected by Petros Hadjicostas, Jul 16 2020]
U(n,k) = binomial(n,k) * 2^(n-1) * hypergeom([-k,-k], [n+1-k], 2) if n >= k >= 0 with (n,k) <> (0,0). - Robert Israel, Jun 15 2011 [Corrected by Petros Hadjicostas, Jul 16 2020]
U(n,k) = Sum_{0 <= i+j <= n+k-1} (-1)^j*C(i+j+1, j)*C(n+i, n)*C(k+i, k). - Masato Maruoka, Dec 10 2019
T(n, k) = 2^(n - 1)*binomial(n, k)*hypergeom([-k, k - n], [-n], 1/2) = A059474(n, k)/2 for n >= 1. - Peter Luschny, Nov 26 2021
From G. C. Greubel, Sep 02 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = T(n, n) = A011782(n).
T(n, n-2) =  2*A049611(n-1), n >= 2.
T(n, n-3) =  4*A049612(n-2), n >= 3.
T(n, n-4) =  8*A055589(n-3), n >= 4.
T(n, n-5) = 16*A055852(n-4), n >= 5.
T(n, n-6) = 32*A055853(n-5), n >= 6.
Sum_{k=0..floor(n/2)} T(n, k) = A181306(n). (End)

A208341 Triangle read by rows, T(n,k) = hypergeometric_2F1([n-k+1, -k], [1], -1) for n>=0 and k>=0.

Original entry on oeis.org

1, 1, 2, 1, 3, 4, 1, 4, 8, 8, 1, 5, 13, 20, 16, 1, 6, 19, 38, 48, 32, 1, 7, 26, 63, 104, 112, 64, 1, 8, 34, 96, 192, 272, 256, 128, 1, 9, 43, 138, 321, 552, 688, 576, 256, 1, 10, 53, 190, 501, 1002, 1520, 1696, 1280, 512, 1, 11, 64, 253, 743, 1683, 2972, 4048

Offset: 0

Clark Kimberling, Feb 25 2012

Previous name was: Triangle of coefficients of polynomials v(n,x) jointly generated with A160232; see the Formula section.
Row sums: (1,3,8,...), even-indexed Fibonacci numbers.
Alt. row sums: (1,-1,2,-3,...), signed Fibonacci numbers.
v(n,2) = A107839(n), v(n,n) = 2^(n-1), v(n+1,n) = A001792(n),
v(n+2,n) = A049611, v(n+3,n) = A049612.
Subtriangle of the triangle T(n,k) given by (1, 0, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 12 2012
Essentially triangle in A049600. - Philippe Deléham, Mar 23 2012

			First five rows:
  1;
  1, 2;
  1, 3,  4;
  1, 4,  8,  8;
  1, 5, 13, 20, 16;
First five polynomials v(n,x):
  1
  1 + 2x
  1 + 3x +  4x^2
  1 + 4x +  8x^2 +  8x^3
  1 + 5x + 13x^2 + 20x^3 + 16x^4
(1, 0, -1/2, 1/2, 0, 0, ...) DELTA (0, 2, 0, 0, 0, ...) begins:
  1;
  1, 0;
  1, 2,  0;
  1, 3,  4,  0;
  1, 4,  8,  8,  0;
  1, 5, 13, 20, 16,  0;
  1, 6, 19, 38, 48, 32, 0;
Triangle in A049600 begins:
  0;
  0, 1;
  0, 1, 2;
  0, 1, 3,  4;
  0, 1, 4,  8,  8;
  0, 1, 5, 13, 20, 16;
  0, 1, 6, 19, 38, 48, 32;
  ... - _Philippe Deléham_, Mar 23 2012

Reinhard Zumkeller, Rows n = 0..124 of triangle, flattened

Cf. A160232, A000045, A049600, A106195.

a208341 n k = a208341_tabl !! (n-1) !! (k-1)
a208341_row n = a208341_tabl !! (n-1)
a208341_tabl = map reverse a106195_tabl
-- Reinhard Zumkeller, Dec 16 2013

T := (n,k) -> hypergeom([n-k+1, -k],[1],-1):
seq(lprint(seq(simplify(T(n,k)),k=0..n)),n=0..7); # Peter Luschny, May 20 2015

u[1, x_] := 1; v[1, x_] := 1; z = 13;
u[n_, x_] := u[n - 1, x] + x*v[n - 1, x];
v[n_, x_] := u[n - 1, x] + 2*x*v[n - 1, x];
Table[Expand[u[n, x]], {n, 1, z/2}]
Table[Expand[v[n, x]], {n, 1, z/2}]
cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
TableForm[cu]
Flatten[%]   (* A160232 *)
Table[Expand[v[n, x]], {n, 1, z}]
cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
TableForm[cv]
Flatten[%]   (* A208341 *)

T(n,k) = sum(i = 0, k, 2^(k-i)*binomial(n-k,i)*binomial(k,i));
tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n, k), ", ")); print();); \\ Michel Marcus, Aug 14 2015

u(n,x) = u(n-1,x) + x*v(n-1,x), v(n,x) = u(n-1,x) + 2x*v(n-1,x), where u(1,x) = 1, v(1,x) = 1.
As DELTA-triangle with 0 <= k <= n: T(n,k) = T(n-1,k) + 2*T(n-1,k-1) - T(n-2,k-1), T(0,0) = T(1,0) = T(2,0) = 1, T(1,1) = T(2,2) = 0, T(2,1) = 2 and T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Mar 12 2012
G.f.: (1-2*y*x+y*x^2)/(1-x-2*y*x+y*x^2). - Philippe Deléham, Mar 12 2012
T(n,k) = A106195(n-1,n-k), k = 1..n. - Reinhard Zumkeller, Dec 16 2013
From Peter Bala, Aug 11 2015: (Start)
The following remarks assume the row and column indexing start at 0.
T(n,k) = Sum_{i = 0..k} 2^(k-i)*binomial(n-k,i)*binomial(k,i) = Sum_{i = 0..k} binomial(n-k+i,i)*binomial(k,i).
Riordan array (1/(1 - x), x*(2 - x)/(1 - x)).
O.g.f. 1/(1 - (2*t + 1)*x + t*x^2) = 1 + (1 + 2*t)*x + (1 + 3*t + 4*t^2)*x^2 + ....
Read as a square array, this equals P * transpose(P^2), where P denotes Pascal's triangle A007318. (End)
For kGlen Whitney, Aug 17 2021

New name from Peter Luschny, May 20 2015

Offset corrected by Joerg Arndt, Aug 12 2015

A055584 Triangle of partial row sums (prs) of triangle A055252.

Original entry on oeis.org

1, 5, 1, 19, 6, 1, 63, 25, 7, 1, 192, 88, 32, 8, 1, 552, 280, 120, 40, 9, 1, 1520, 832, 400, 160, 49, 10, 1, 4048, 2352, 1232, 560, 209, 59, 11, 1, 10496, 6400, 3584, 1792, 769, 268, 70, 12, 1, 26624, 16896, 9984, 5376, 2561, 1037, 338, 82, 13, 1, 66304, 43520

Offset: 0

Wolfdieter Lang, May 26 2000

In the language of the Shapiro et al. reference (given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The G.f. for the row polynomials p(n,x) (increasing powers of x) is (((1-z)^3)/(1-2*z)^4)/(1-x*z/(1-z)).
This is the fourth member of the family of Riordan-type matrices obtained from A007318(n,m) (Pascal's triangle read as lower triangular matrix) by repeated application of the prs-procedure.
The column sequences appear as A049612(n+1), A055585, A001794, A001789(n+3), A027608, A055586 for m=0..5.

			[0] 1
[1] 5, 1
[2] 19, 6, 1
[3] 63, 25, 7, 1
[4] 192, 88, 32, 8, 1
[5] 552, 280, 120, 40, 9, 1
[6] 1520, 832, 400, 160, 49, 10, 1
[7] 4048, 2352, 1232, 560, 209, 59, 11, 1
Fourth row polynomial (n=3): p(3, x)= 63 + 25*x + 7*x^2 + x^3.

Cf. A007318, A055248, A055249, A055252. Row sums: A049600(n+1, 4).

T := (n, k) -> binomial(n, k)*hypergeom([4, k - n], [k + 1], -1):
for n from 0 to 7 do seq(simplify(T(n, k)), k = 0..n) od; # Peter Luschny, Sep 23 2024

a(n, m)=sum(A055252(n, k), k=m..n), n >= m >= 0, a(n, m) := 0 if n

Column m recursion: a(n, m)= sum(a(j, m), j=m..n-1)+ A055252(n, m), n >= m >= 0, a(n, m) := 0 if n

G.f. for column m: (((1-x)^3)/(1-2*x)^4)*(x/(1-x))^m, m >= 0.

T(n, k) = binomial(n, k)*hypergeom([4, k - n], [k + 1], -1). - Peter Luschny, Sep 23 2024

A055852 Convolution of A055589 with A011782.

Original entry on oeis.org

0, 1, 7, 34, 138, 501, 1683, 5336, 16172, 47264, 134048, 370688, 1003136, 2664192, 6960384, 17922048, 45552640, 114442240, 284508160, 700579840, 1710161920, 4141416448, 9955639296, 23770693632, 56400543744, 133041225728

Offset: 0

Wolfdieter Lang May 30 2000

Seventh column of triangle A055587.

T(n,5) of array T as in A049600.

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (12,-60,160,-240,192,-64).

Cf. A011782, A049600, A055587, A055589.

a:=[1,7,34,138,501,1683];; for n in [7..30] do a[n]:=12*a[n-1] -60*a[n-2] +160*a[n-3] -240*a[n-4] +192*a[n-5] -64*a[n-6]; od; Concatenation([0], a); # G. C. Greubel, Jan 16 2020

R:=PowerSeriesRing(Integers(), 30); [0] cat Coefficients(R!( x*(1-x)^5/(1-2*x)^6 )); // G. C. Greubel, Jan 16 2020

seq(coeff(series(x*(1-x)^5/(1-2*x)^6, x, n+1), x, n), n = 0..30); # G. C. Greubel, Jan 16 2020

CoefficientList[Series[x*(1-x)^5/(1-2*x)^6, {x,0,30}], x] (* G. C. Greubel, Jan 16 2020 *)

my(x='x+O('x^30)); concat([0], Vec(x*(1-x)^5/(1-2*x)^6)) \\ G. C. Greubel, Jan 16 2020

def A055852_list(prec):
    P. = PowerSeriesRing(ZZ, prec)
    return P( x*(1-x)^5/(1-2*x)^6 ).list()
A055852_list(30) # G. C. Greubel, Jan 16 2020

a(n) = T(n, 5) = A055587(n+5, 6).

G.f.: x*(1-x)^5/(1-2*x)^6.

A035029 a(n) = Sum_{k=0..n} (k+1) * Sum_{j=0..n} 2^jbinomial(n,j)binomial(n-k,j).

Original entry on oeis.org

1, 5, 26, 138, 743, 4043, 22180, 122468, 679757, 3789297, 21199998, 118973550, 669447123, 3775577367, 21336790152, 120795829128, 684962855705, 3889578815453, 22115533878178, 125892252068498, 717400693313471, 4092099111728355, 23362391663233196, 133488737662062188, 763310051648602213

Offset: 0

N. J. A. Sloane

nonn

Number of regions in all the dissections of a convex (n+3)-gon by non-intersecting diagonals. a(1)=5 because in the three dissections of a square we have altogether five regions: one in the "no-diagonals" dissection and two in each of the dissections by one of the two diagonals of the square. - Emeric Deutsch, Dec 28 2003

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Milan Janjic, Two Enumerative Functions
D. Merlini, R. Sprugnoli and M. C. Verri, Waiting patterns for a printer, FUN with algorithm'01, Isola d'Elba, 2001.

Cf. A001003, A001850, A002002, A049600.

I:=[1,5,26]; [n le 3 select I[n] else ( (7*n+1)*Self(n-1) - (7*n-15)*Self(n-2) + (n-3)*Self(n-3) )/(n+1): n in [1..30]]; // G. C. Greubel, Oct 20 2022

CoefficientList[Series[(1-x)^2/(8*x^2*Sqrt[1-6*x+x^2])-(1+x)/(8*x^2), {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 12 2014 *)
With[{P = LegendreP}, Table[((n+1)*(n+3)*P[n+3,3] -(6*n^2+22*n+17)*P[n+2,3] +(n+ 2)*(5*n+8)*P[n+1,3])/(8*(n+1)*(n+2)), {n,0,40}]] (* G. C. Greubel, Oct 20 2022 *)

def A001850(n): return gen_legendre_P(n,0,3)
def A035029(n): return ((n+1)*(n+3)*A001850(n+3) - (6*n^2 +22*n +17)*A001850(n+2) + (n+2)*(5*n+8)*A001850(n+1))/(8*(n+1)*(n+2))
[A035029(n) for n in range(40)] # G. C. Greubel, Oct 20 2022

a(n) = (1/4)*(A002002(n+2) - A002002(n+1)).
G.f.: (1-z)^2/(8*z^2*sqrt(1-6*z+z^2))-(1+z)/(8*z^2). - Emeric Deutsch, Dec 28 2003
a(n) = T(n+1, n+2), array T as in A049600.
Form an array with the m(n,1)=1 and m(1,n) = n*(n+1)/2 for n=1,2,3... The interior terms m(i,j) = m(i,j-1) + m(i-1,j-1) + m(i-1,j). The terms in the main diagonal = {a(n)}. - J. M. Bergot, Dec 01 2012
D-finite with recurrence: (n+2)*a(n) + (7*n+8)*a(n-1) - (7*n-8)*a(n-2) + (n-2)*a(n-3). - R. J. Mathar, Dec 03 2012
a(n) ~ (3+2*sqrt(2))^(n+3/2) / (2^(9/4)*sqrt(Pi*n)). - Vaclav Kotesovec, Feb 12 2014
a(n) = Jypergeometric2F1([-n, n+3]; [1]; -1), which satisfies the recurrence. - Benedict W. J. Irwin, Oct 14 2016

Previous Showing 11-20 of 28 results. Next

cp's OEIS Frontend

A008288 Square array of Delannoy numbers D(i,j) (i >= 0, j >= 0) read by antidiagonals.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Python

Sage

Formula

Extensions

A034856 a(n) = binomial(n+1, 2) + n - 1 = n*(n+3)/2 - 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Maxima

PARI

Python

Formula

Extensions

A100258 Triangle of coefficients of normalized Legendre polynomials, with increasing exponents.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A002003 a(n) = 2 * Sum_{k=0..n-1} binomial(n-1, k)*binomial(n+k, k).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A142978 Table of figurate numbers for the n-dimensional cross polytopes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Formula

A035029 a(n) = Sum_{k=0..n} (k+1) * Sum_{j=0..n} 2^jbinomial(n,j)binomial(n-k,j).