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

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

Original entry on oeis.org

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

Views

Author

Keywords

Comments

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

Examples

			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)
		

References

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

Crossrefs

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
See also A027618.
Cf. A059446.
Has same main diagonal as A064861. Different from A100936.
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.

Programs

  • Haskell
    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
    
  • Maple
    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:
  • Mathematica
    (* 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 *)
  • Python
    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
  • Sage
    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
    

Formula

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

Extensions

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

A037027 Skew Fibonacci-Pascal triangle read by rows.

Original entry on oeis.org

1, 1, 1, 2, 2, 1, 3, 5, 3, 1, 5, 10, 9, 4, 1, 8, 20, 22, 14, 5, 1, 13, 38, 51, 40, 20, 6, 1, 21, 71, 111, 105, 65, 27, 7, 1, 34, 130, 233, 256, 190, 98, 35, 8, 1, 55, 235, 474, 594, 511, 315, 140, 44, 9, 1, 89, 420, 942, 1324, 1295, 924, 490, 192, 54, 10, 1, 144, 744, 1836
Offset: 0

Views

Author

Floor van Lamoen, Jan 01 1999

Keywords

Comments

T(n,k) is the number of lattice paths from (0,0) to (n,k) using steps (0,1), (1,0), (2,0). - Joerg Arndt, Jun 30 2011
T(n,k) is the number of lattice paths of length n, starting from the origin and ending at (n,k), using horizontal steps H=(1,0), up steps U=(1,1) and down steps D=(1,-1), never containing UUU, DD, HD. For instance, for n=4 and k=2, we have the paths; HHUU, HUHU, HUUH, UHHU, UHUH, UUHH, UUDU, UDUU, UUUD. - Emanuele Munarini, Mar 15 2011
Row sums form Pell numbers A000129, T(n,0) forms Fibonacci numbers A000045, T(n,1) forms A001629. T(n+k,n-k) is polynomial sequence of degree k.
T(n,k) gives a convolved Fibonacci sequence (A001629, A001872, etc.).
As a Riordan array, this is (1/(1-x-x^2),x/(1-x-x^2)). An interesting factorization is (1/(1-x^2),x/(1-x^2))*(1/(1-x),x/(1-x)) [abs(A049310) times A007318]. Diagonal sums are the Jacobsthal numbers A001045(n+1). - Paul Barry, Jul 28 2005
T(n,k) = T'(n+1,k+1), T' given by [0, 1, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 19 2005
Equals A049310 * A007318 as infinite lower triangular matrices. - Gary W. Adamson, Oct 28 2007
This triangle may also be obtained from the coefficients of the Morgan-Voyce polynomials defined by: Mv(x, n) = (x + 1)*Mv(x, n - 1) + Mv(x, n - 2). - Roger L. Bagula, Apr 09 2008
Row sums are A000129. - Roger L. Bagula, Apr 09 2008
Absolute value of coefficients of the characteristic polynomial of tridiagonal matrices with 1's along the main diagonal, and i's along the superdiagonal and the subdiagonal (where i=sqrt(-1), see Mathematica program). - John M. Campbell, Aug 23 2011
A037027 is jointly generated with A122075 as an array of coefficients of polynomials v(n,x): initially, u(1,x)=v(1,x)=1; for n>1, u(n,x)=u(n-1,x)+(x+1)*v(n-1)x and v(n,x)=u(n-1,x)+x*v(n-1,x). See the Mathematica section at A122075. - Clark Kimberling, Mar 05 2012
For a closed-form formula for arbitrary left and right borders of Pascal like triangle see A228196. - Boris Putievskiy, Aug 18 2013
For a closed-form formula for generalized Pascal's triangle see A228576. - Boris Putievskiy, Sep 09 2013
Row n, for n>=0, shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(n)*x^n which is the denominator of the n-th convergent of the continued fraction [x+1, x+1, x+1, ...]; see A230000. - Clark Kimberling, Nov 13 2013
T(n,k) is the number of ternary words of length n having k letters 2 and avoiding a runs of odd length for the letter 0. - Milan Janjic, Jan 14 2017
Let T(m, n, k) be an m-bonacci Pascal's triangle, where T(m, n, 0) gives the values of F(m, n), the n-th m-bonacci number, and T(m, n, k) gives the values for the k-th convolution of F(m, n). Then the classic Pascal triangle is T(1, n, k) and this sequence is T(2, n, k). T(m, n, k) is the number of compositions of n using only the positive integers 1, 1' and 2 through m, with the part 1' used exactly k times. G.f. for k-th column of T(m, n, k): x/(1 - x - x^2 - ... - x^m)^k. The row sum for T(m, n, k) is the number of compositions of n using only the positive integers 1, 1' and 2 through m. G.f. for row sum of T(m, n, k): 1/(1 - 2x - x^2 - ... - x^m). - Gregory L. Simay, Jul 24 2021

Examples

			Ratio of row polynomials R(3)/R(2) = (3 + 5*x + 3*x^2 + x^3)/(2 + 2*x + x^2) = [1+x; 1+x, 1+x].
Triangle begins:
                                 1;
                              1,    1;
                           2,    2,    1;
                        3,    5,    3,    1;
                     5,   10,    9,    4,    1;
                  8,   20,   22,   14,    5,    1;
              13,   38,   51,   40,   20,    6,    1;
           21,   71,  111,  105,   65,   27,    7,    1;
        34,  130,  233,  256,  190,   98,   35,    8,    1;
     55,  235,  474,  594,  511,  315,  140,   44,    9,    1;
  89,  420,  942, 1324, 1295,  924,  490,  192,   54,   10,    1;
		

Crossrefs

A038112(n) = T(2n, n). A038137 is reflected version. Maximal row entries: A038149.
Diagonal differences are in A055830. Vertical sums are in A091186.
Some other Fibonacci-Pascal triangles: A027926, A036355, A074829, A105809, A109906, A111006, A114197, A162741, A228074.

Programs

  • Haskell
    a037027 n k = a037027_tabl !! n !! k
    a037027_row n = a037027_tabl !! n
    a037027_tabl = [1] : [1,1] : f [1] [1,1] where
       f xs ys = ys' : f ys ys' where
         ys' = zipWith3 (\u v w -> u + v + w) (ys ++ [0]) (xs ++ [0,0]) ([0] ++ ys)
    -- Reinhard Zumkeller, Jul 07 2012
  • Maple
    T := (n,k) -> `if`(n=0,1,binomial(n,k)*hypergeom([(k-n)/2, (k-n+1)/2], [-n], -4)): seq(seq(simplify(T(n,k)), k=0..n), n=0..10); # Peter Luschny, Apr 25 2016
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> combinat:-fibonacci(n)); # Peter Luschny, Oct 07 2022
  • Mathematica
    Mv[x, -1] = 0; Mv[x, 0] = 1; Mv[x, 1] = 1 + x; Mv[x_, n_] := Mv[x, n] = ExpandAll[(x + 1)*Mv[x, n - 1] + Mv[x, n - 2]]; Table[ CoefficientList[ Mv[x, n], x], {n, 0, 10}] // Flatten (* Roger L. Bagula, Apr 09 2008 *)
    Abs[Flatten[Table[CoefficientList[CharacteristicPolynomial[Array[KroneckerDelta[#1,#2]+KroneckerDelta[#1,#2+1]*I+KroneckerDelta[#1,#2-1]*I&,{n,n}],x],x],{n,1,20}]]] (* John M. Campbell, Aug 23 2011 *)
    T[n_, k_] := Binomial[n, k] Hypergeometric2F1[(k-n)/2, (k-n+1)/2, -n, -4];
    Table[T[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Feb 16 2019, after Peter Luschny *)
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, if( n==0 && k==0, 1, T(n-1, k) + T(n-1, k-1) + T(n-2, k)))}; /* Michael Somos, Sep 29 2003 */
    
  • PARI
    T(n,k)=if(nPaul D. Hanna, Feb 27 2004
    

Formula

T(n, m) = T'(n-1, m) + T'(n-2, m) + T'(n-1, m-1), where T'(n, m) = T(n, m) for n >= 0 and 0< = m <= n and T'(n, m) = 0 otherwise.
G.f.: 1/(1 - y - y*z - y^2).
G.f. for k-th column: x/(1-x-x^2)^k.
T(n, m) = Sum_{k=0..n-m} binomial(m+k, m)*binomial(k, n-k-m), n >= m >= 0, otherwise 0. - Wolfdieter Lang, Jun 17 2002
T(n, m) = ((n-m+1)*T(n, m-1) + 2*(n+m)*T(n-1, m-1))/(5*m), n >= m >= 1; T(n, 0)= A000045(n+1); T(n, m)= 0 if n < m. - Wolfdieter Lang, Apr 12 2000
Chebyshev coefficient triangle (abs(A049310)) times Pascal's triangle (A007318) as product of lower triangular matrices. T(n, k) = Sum_{j=0..n} binomial((n+j)/2, j)*(1+(-1)^(n+j))*binomial(j, k)/2. - Paul Barry, Dec 22 2004
Let R(n) = n-th row polynomial in x, with R(0)=1, then R(n+1)/R(n) equals the continued fraction [1+x;1+x, ...(1+x) occurring (n+1) times ..., 1+x] for n >= 0. - Paul D. Hanna, Feb 27 2004
T(n,k) = Sum_{j=0..n} binomial(n-j,j)*binomial(n-2*j,k); in Egorychev notation, T(n,k) = res_w(1-w-w^2)^(-k-1)*w^(-n+k+1). - Paul Barry, Sep 13 2006
Sum_{k=0..n} T(n,k)*x^k = A000045(n+1), A000129(n+1), A006190(n+1), A001076(n+1), A052918(n), A005668(n+1), A054413(n), A041025(n), A099371(n+1), A041041(n), A049666(n+1), A041061(n), A140455(n+1), A041085(n), A154597(n+1), A041113(n) for x = 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 respectively. - Philippe Deléham, Nov 29 2009
T((m+1)*n+r-1, m*n+r-1)*r/(m*n+r) = Sum_{k=1..n} k/n*T((m+1)*n-k-1, m*n-1)*(r+k,r), n >= m > 1.
T(n-1,m-1) = (m/n)*Sum_{k=1..n-m+1} k*A000045(k)*T(n-k-1,m-2), n >= m > 1. - Vladimir Kruchinin, Mar 17 2011
T(n,k) = binomial(n,k)*hypergeom([(k-n)/2, (k-n+1)/2], [-n], -4) for n >= 1. - Peter Luschny, Apr 25 2016

Extensions

Examples from Paul D. Hanna, Feb 27 2004

A230001 Triangular array read by rows: row n shows the coefficients of the polynomial v(n) = d(0) + d(1)*x + ... + d(k)*x^k which is the denominator of the n-th convergent of the continued fraction [1, 1/x, 1/x^2, ... ,1/x^n].

Original entry on oeis.org

1, 1, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 0, 2, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 2, 0, 2, 0, 2, 0, 1, 1, 1, 1, 0, 1
Offset: 0

Views

Author

Clark Kimberling, Oct 11 2013

Keywords

Comments

In the Name section, k = n(n+1)/2 if n is even, and k = (n-1)(n+2)/2 if n is odd. For the numerator polynomials, see A230000. Conjecture: every nonnegative integer occurs infinitely many times.

Examples

			The first 7 rows:
1 . . . . . . . . . . . . polynomial v(0) = 1
1 . . . . . . . . . . . . polynomial v(1) = 1
1 0 0 1 . . . . . . . . . v(2) = 1 + x^3
1 0 0 1 0 1 . . . . . . . v(3) = 1 + x^3 + x^5
1 0 0 1 0 1 0 1 0 0 1
1 0 0 1 0 1 0 1 0 1 1 0 1 0 1
		

Crossrefs

Cf. A230000.

Programs

  • Mathematica
    t[n_] := t[n] = Table[1/x^k, {k, 0, n}];
    b = Table[Factor[Denominator[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]
    Flatten[v]

Formula

Write the numerator polynomials as u(0), u(1), u(2), ... and the denominator polynomials as v(0), v(1), v(2),... Let p(0) = 1, q(0) = 1; p(1) = (1 + x)/x; q(1) = 1/x; p(n ) = p(n-1)/x^n + p(n-2), q(n) = q(n-1)/x^n + q(n-2). Then u(n)/v(n) = p(n)/q(n) for n>=0.

A231154 Triangular array read by rows: row n 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 = (x + 1)/(x - 1).

Original entry on oeis.org

1, 1, 2, 0, 2, 3, 1, 1, 3, 5, 0, 6, 0, 5, 8, 0, 8, 8, 0, 8, 13, -2, 19, 4, 19, -2, 13, 21, -5, 33, 15, 15, 33, -5, 21, 34, -12, 64, 12, 60, 12, 64, -12, 34, 55, -25, 116, 20, 90, 90, 20, 116, -25, 55, 89, -50, 213, 8, 210, 84, 210, 8, 213, -50, 89
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: 2^n. Left and right edges: A000045 (Fibonacci numbers).

Examples

			First 5 rows:
1 . . . 1
2 . . . 0 . . . 2
3 . . . 1 . . . 1 . . . 3
5 . . . 0 . . . 6 . . . 0 . . . 5
8 . . . 0 . . . 8 . . . 8 . . . 0 . . . 8
First 3 polynomials:  1 + x, 2 + 2*x^2, 3 + x + x^2 + 3*x^3.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(x + 1)/(x - 1), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231727 Triangular array read by rows: row n shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(n)*x^n which is the denominator of the n-th convergent of the continued fraction [k, k, k, ... ], where k = (x + 1)/(x - 1).

Original entry on oeis.org

-1, 1, -1, 0, 1, -2, 2, -2, 2, -3, 2, 0, -2, 3, -5, 5, -6, 6, -5, 5, -8, 8, -8, 0, 8, -8, 8, -13, 15, -21, 15, -15, 21, -15, 13, -21, 26, -38, 18, 0, -18, 38, -26, 21, -34, 46, -76, 52, -48, 48, -52, 76, -46, 34, -55, 80, -141, 96, -70, 0, 70, -96, 141, -80
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: 0. Left and right edges: A000045 (Fibonacci numbers).

Examples

			First 5 rows:
-1 . . . 1
-1 . . . 0 . . . 1
-2 . . . 2 . . . -2 . . . 2
-3 . . . 2 . . . 0 . . . -2 . . . 3
-5 . . . 5 . . . -6 . . . 6 . . . -5 . . . 5
First 3 polynomials:  -1 + x, -1 + x^2, -2 + 2*x - 2*x^2 + 2*x^3.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(x + 1)/(x - 1), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Denominator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231728 Triangular array read by rows: row n shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(2n)*x^(2n) which is the numerator of the n-th convergent of the continued fraction [k, k, k, ... ], where k = (x^2 + 1)/(x + 1).

Original entry on oeis.org

1, 0, 1, 2, 2, 3, 0, 1, 3, 4, 7, 4, 5, 0, 1, 5, 10, 19, 16, 16, 6, 7, 0, 1, 8, 20, 42, 48, 55, 36, 29, 8, 9, 0, 1, 13, 40, 94, 132, 164, 138, 119, 64, 46, 10, 11, 0, 1, 21, 76, 197, 324, 451, 464, 439, 304, 219, 100, 67, 12, 13, 0, 1, 34, 142, 405, 760, 1170
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: 2*A087206(n). Left edge: A000045 (Fibonacci numbers).

Examples

			First 3 rows:
1 . . . 0 . . . 1
2 . . . 2 . . . 3 . . . 0 . . . 1
3 . . . 4 . . . 7 . . . 4 . . . 5 . . . 0 . . . 1
First 3 polynomials:  1 + x^2, 2 + 2*x + 3*x^2 + x^4.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(1 + x^2)/(1 + x), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231729 Triangular array read by rows: row n shows the coefficients of the polynomial u(n) = c(0) + c(1)*x + ... + c(2n-1)*x^(2n-1) which is the denominator of the n-th convergent of the continued fraction [k, k, k, ... ], where k = (x^2 + 1)/(x + 1).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 4, 5, 3, 1, 1, 3, 7, 11, 11, 9, 5, 1, 1, 5, 15, 29, 35, 32, 22, 13, 7, 1, 1, 8, 28, 62, 90, 103, 91, 65, 37, 17, 9, 1, 1, 13, 53, 134, 226, 296, 302, 257, 183, 110, 56, 21, 11, 1, 1, 21, 97, 273, 521, 775, 915, 903, 743, 523, 319, 167
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: 2*A063727(n). Left edge: A000045 (Fibonacci numbers).

Examples

			First 3 rows:
1 . . . 1
1 . . . 1 . . . 1 . . . 1
2 . . . 4 . . . 5 . . . 3 . . . 1 . . . 1 . . . 1
First 3 polynomials:  1 + x, 1 + x + x^2 + x^3, 2 + 4*x + 5*x^2 + 3*x^3 + x^4 + x^5.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(1 + x^2)/(1 + x), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Denominator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231730 Triangular array read by rows: row n 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 = x + 1/2.

Original entry on oeis.org

1, 2, 5, 4, 4, 9, 22, 12, 8, 29, 56, 72, 32, 16, 65, 202, 232, 208, 80, 32, 181, 556, 924, 800, 560, 192, 64, 441, 1726, 2964, 3480, 2480, 1440, 448, 128, 1165, 4832, 10112, 12608, 11680, 7168, 3584, 1024, 256, 2929, 14066, 31632, 46752, 46816, 36288, 19712
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: A015521(n). Left edge: A006131. Right edge: powers of 2

Examples

			First 3 rows:
1 .... 2
5 .... 4 .... 4
9 .... 22 ... 12 ... 8
First 3 polynomials:  1 + 2*x, 5 + 4*x + 4*x^2, 9 + 22*x + 12*x^2 + 8*x^3.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[x + 1/2, {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231731 Triangular array read by rows: row n 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 = 2*x + 1.

Original entry on oeis.org

1, 2, 2, 4, 4, 3, 10, 12, 8, 5, 20, 36, 32, 16, 8, 40, 88, 112, 80, 32, 13, 76, 204, 320, 320, 192, 64, 21, 142, 444, 840, 1040, 864, 448, 128, 34, 260, 932, 2048, 3040, 3136, 2240, 1024, 256, 55, 470, 1896, 4752, 8176, 10080, 8960, 5632, 2304, 512, 89, 840
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: A006190(n). Left edge: A000045 (Fibonacci numbers). Right edge: powers of 2

Examples

			First 3 rows:
1 ... 2
2 ... 4 .... 4
3 ... 10 ... 12 ... 8
First 3 polynomials:  1 + 2*x, 2 + 4*x + 4*x^2, 3 + 10*x + 12*x^2 + 8*x^3.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[2 x + 1, {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]

A231732 Triangular array read by rows: row n 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 = (x + 2)/(x + 1).

Original entry on oeis.org

2, 1, 5, 6, 2, 12, 22, 14, 3, 29, 72, 69, 30, 5, 70, 219, 280, 182, 60, 8, 169, 638, 1021, 884, 436, 116, 13, 408, 1804, 3468, 3750, 2460, 978, 218, 21, 985, 4992, 11206, 14532, 11895, 6288, 2095, 402, 34, 2378, 13589, 34888, 52760, 51750, 34119, 15112, 4334
Offset: 1

Views

Author

Clark Kimberling, Nov 13 2013

Keywords

Comments

Sum of numbers in row n: A015521(n). Left edge: A000129. Right edge: A000045 (Fibonacci numbers).

Examples

			First 3 rows:
2 .... 1
5 .... 6 .... 2
12 ... 22 ... 14 ... 3
First 3 polynomials:  2 + x, 5 + 6*x + 2*x^2, 12 + 22*x + 14*x^2 + 3*x^3.
		

Crossrefs

Programs

  • Mathematica
    t[n_] := t[n] = Table[(x + 2)/(x + 1), {k, 0, n}];
    b = Table[Factor[Convergents[t[n]]], {n, 0, 10}];
    p[x_, n_] := p[x, n] = Last[Expand[Numerator[b]]][[n]];
    u = Table[p[x, n], {n, 1, 10}]
    v = CoefficientList[u, x]; Flatten[v]
Showing 1-10 of 13 results. Next