A091965 -id:A091965 - cp's OEIS frontend

A111418 Right-hand side of odd-numbered rows of Pascal's triangle.

1, 3, 1, 10, 5, 1, 35, 21, 7, 1, 126, 84, 36, 9, 1, 462, 330, 165, 55, 11, 1, 1716, 1287, 715, 286, 78, 13, 1, 6435, 5005, 3003, 1365, 455, 105, 15, 1, 24310, 19448, 12376, 6188, 2380, 680, 136, 17, 1, 92378, 75582, 50388

Offset: 0

Philippe Deléham, Nov 13 2005

Riordan array (c(x)/sqrt(1-4*x),x*c(x)^2) where c(x) is g.f. of A000108. Unsigned version of A113187. Diagonal sums are A014301(n+1).
Triangle T(n,k),0<=k<=n, read by rows defined by :T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=3*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+2*T(n-1,k)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 22 2007
Reversal of A122366. - Philippe Deléham, Mar 22 2007
Column k has e.g.f. exp(2x)(Bessel_I(k,2x)+Bessel_I(k+1,2x)). - Paul Barry, Jun 06 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Diagonal sums are A014301(n+1). - Paul Barry, Mar 08 2011
This triangle T(n,k) appears in the expansion of odd powers of Fibonacci numbers F=A000045 in terms of F-numbers with multiples of odd numbers as indices. See the Ozeki reference, p. 108, Lemma 2. The formula is: F_l^(2*n+1) = sum(T(n,k)*(-1)^((n-k)*(l+1))* F_{(2*k+1)*l}, k=0..n)/5^n, n >= 0, l >= 0. - Wolfdieter Lang, Aug 24 2012
Central terms give A052203. - Reinhard Zumkeller, Mar 14 2014
This triangle appears in the expansion of (4*x)^n in terms of the polynomials Todd(n, x):= T(2*n+1, sqrt(x))/sqrt(x) = sum(A084930(n,m)*x^m), n >= 0. This follows from the inversion of the lower triangular Riordan matrix built from A084930 and comparing the g.f. of the row polynomials. - Wolfdieter Lang, Aug 05 2014
From Wolfdieter Lang, Aug 15 2014: (Start)
This triangle is the inverse of the signed Riordan triangle (-1)^(n-m)*A111125(n,m).
This triangle T(n,k) appears in the expansion of x^n in terms of the polynomials todd(k, x):= T(2*k+1, sqrt(x)/2)/(sqrt(x)/2) = S(k, x-2) - S(k-1, x-2) with the row polynomials T and S for the triangles A053120 and A049310, respectively: x^n = sum(T(n,k)*todd(k, x), k=0..n). Compare this with the preceding comment.
The A- and Z-sequences for this Riordan triangle are [1, 2, 1, repeated 0] and [3, 1, repeated 0]. For A- and Z-sequences for Riordan triangles see the W. Lang link under A006232. This corresponds to the recurrences given in the Philippe Deléham, Mar 22 2007 comment above. (End)

			From _Wolfdieter Lang_, Aug 05 2014: (Start)
The triangle T(n,k) begins:
n\k      0      1      2      3     4     5    6    7   8  9  10 ...
0:       1
1:       3      1
2:      10      5      1
3:      35     21      7      1
4:     126     84     36      9     1
5:     462    330    165     55    11     1
6:    1716   1287    715    286    78    13    1
7:    6435   5005   3003   1365   455   105   15    1
8:   24310  19448  12376   6188  2380   680  136   17   1
9:   92378  75582  50388  27132 11628  3876  969  171  19  1
10: 352716 293930 203490 116280 54264 20349 5985 1330 210 21   1
...
Expansion examples (for the Todd polynomials see A084930 and a comment above):
(4*x)^2 = 10*Todd(n,  0) + 5*Todd(n, 1) + 1*Todd(n, 2) = 10*1 + 5*(-3 + 4*x) + 1*(5 - 20*x + 16*x^2).
(4*x)^3 =  35*1 + 21*(-3 + 4*x) + 7*(5 - 20*x + 16*x^2) + (-7 + 56*x - 112*x^2 +64*x^3)*1. (End)
---------------------------------------------------------------------
Production matrix is
3, 1,
1, 2, 1,
0, 1, 2, 1,
0, 0, 1, 2, 1,
0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 0, 1, 2, 1,
0, 0, 0, 0, 0, 0, 0, 1, 2, 1
- _Paul Barry_, Mar 08 2011
Application to odd powers of Fibonacci numbers F, row n=2:
F_l^5 = (10*(-1)^(2*(l+1))*F_l + 5*(-1)^(1*(l+1))*F_{3*l} + 1*F_{5*l})/5^2, l >= 0. - _Wolfdieter Lang_, Aug 24 2012

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
Paul Barry, On the Hurwitz Transform of Sequences, Journal of Integer Sequences, Vol. 15 (2012), #12.8.7.
E. Deutsch, L. Ferrari and S. Rinaldi, Production Matrices, Advances in Applied Mathematics, 34 (2005) pp. 101-122.
Asamoah Nkwanta and Earl R. Barnes, Two Catalan-type Riordan Arrays and their Connections to the Chebyshev Polynomials of the First Kind, Journal of Integer Sequences, Article 12.3.3, 2012.
A. Nkwanta, A. Tefera, Curious Relations and Identities Involving the Catalan Generating Function and Numbers, Journal of Integer Sequences, 16 (2013), #13.9.5.
K. Ozeki, On Melham's sum, The Fibonacci Quart. 46/47 (2008/2009), no. 2, 107-110.
Sun, Yidong; Ma, Luping Minors of a class of Riordan arrays related to weighted partial Motzkin paths. Eur. J. Comb. 39, 157-169 (2014), Table 2.2.

Cf. A000108, A113187.

Columns are : A001700, A002054, A003516, A030053, A030054, A030055, A030056.

a111418 n k = a111418_tabl !! n !! k
a111418_row n = a111418_tabl !! n
a111418_tabl = map reverse a122366_tabl
-- Reinhard Zumkeller, Mar 14 2014

Table[Binomial[2*n+1, n-k], {n,0,10}, {k,0,n}] (* G. C. Greubel, May 22 2017 *)
T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0,
T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]];
Table[T[n, k, 3, 2], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)

T(n, k) = C(2*n+1, n-k).
Sum_{k=0..n} T(n, k) = 4^n.
Sum_{k, 0<=k<=n}(-1)^k *T(n,k) = binomial(2*n,n) = A000984(n). - Philippe Deléham, Mar 22 2007
T(n,k) = sum{j=k..n, C(n,j)*2^(n-j)*C(j,floor((j-k)/2))}. - Paul Barry, Jun 06 2007
Sum_{k, k>=0} T(m,k)*T(n,k) = T(m+n,0)= A001700(m+n). - Philippe Deléham, Nov 22 2009
G.f. row polynomials: ((1+x) - (1-x)/sqrt(1-4*z))/(2*(x - (1+x)^2*z))
(see the Riordan property mentioned in a comment above). - Wolfdieter Lang, Aug 05 2014

A033321 Binomial transform of Fine's sequence A000957: 1, 0, 1, 2, 6, 18, 57, 186, ...

Original entry on oeis.org

1, 1, 2, 6, 21, 79, 311, 1265, 5275, 22431, 96900, 424068, 1876143, 8377299, 37704042, 170870106, 779058843, 3571051579, 16447100702, 76073821946, 353224531663, 1645807790529, 7692793487307, 36061795278341, 169498231169821

Offset: 0

Emeric Deutsch

nonn

Number of permutations avoiding the patterns {2431,4231,4321}; number of weak sorting class based on 2431. - Len Smiley, Nov 01 2005
Number of permutations avoiding the patterns {2413, 3142, 2143}. - Vincent Vatter, Aug 16 2006
Number of permutations avoiding the patterns {2143, 3142, 4132}. - Alexander Burstein and Jonathan Bloom, Aug 03 2013
Number of unimodal Lehmer codes. Those are exactly the inversion sequences for permutations avoiding the patterns {2143, 3142, 4132}. - Alexander Burstein, Jun 16 2015
Number of skew Dyck paths of semilength n ending with a down step (1,-1). A skew Dyck path is a path in the first quadrant which begins at the origin, ends on the x-axis, consists of steps U=(1,1)(up), D=(1,-1)(down) and L=(-1,-1)(left) so that up and left steps do not overlap. The length of the path is defined to be the number of its steps. Number of skew Dyck paths of semilength n and ending with a left step is A128714(n). - Emeric Deutsch, May 11 2007
Number of permutations sortable by a pop stack followed directly by a stack. Equivalently, the number of permutations avoiding {2431, 3142, 3241}. - Vincent Vatter, Mar 06 2013
Hankel transform of this sequence gives A000012 = [1,1,1,1,1,1,...]. - Philippe Deléham, Oct 24 2007
Starting with offset 1, Hankel transform = odd-indexed Fibonacci numbers. - Gary W. Adamson, Dec 27 2008
Starting with offset 1 = INVERT transform of A002212: (1, 1, 3, 10, 36, 137, ...). - Gary W. Adamson, May 19 2009
Equals INVERTi transform of A007317: (1, 2, 5, 15, 51, 188, ...). - Gary W. Adamson, May 17 2009
Number of sequences (e(1), ..., e(n)), 0 <= e(i) < i, such that there is no triple i < j < k with e(i) > e(j) < e(k). [Martinez and Savage, 2.20] - Eric M. Schmidt, Jul 17 2017
From David Callan, Jul 21 2017: (Start)
a(n) is the number of permutations of [n] in which the excedances and subcedances are both increasing. (For example, the 3 permutations of [4] NOT counted by a(4)=21 are 3421, 4312, 4321 with excedances/subcedances 34/21, 43/12, 43/21 respectively.)
Proof. It suffices to show that (*) the number of such permutations of [n] containing k fixed points is binomial(n,k)*F(n-k), where F is the Fine number A000957. Since F(n) is the number of 321-avoiding derangements of [n] and because inserting or deleting a fixed point in a permutation does not change the excedance/fixed point/subcedance status of any other entry, (*) is an immediate consequence of the following claim: The excedances and subcedances of a permutation p are both increasing if and only if p avoids 321. The claim is a nice exercise utilizing the cycles of p for the "if" direction and the pigeonhole principle for the "only if" direction. (End)
Conjectured to be the number of permutations of length n that are sorted to the identity by a consecutive-231-avoiding stack followed by a classical-21-avoiding stack. - Colin Defant, Aug 30 2020
a(n) is the number of permutations of length n avoiding the partially ordered pattern (POP) {3>1, 3>4, 1>2} of length 4. That is, the number of length n permutations having no subsequences of length 4 in which the third element is the largest and the first element is larger than the second element. - Sergey Kitaev, Dec 10 2020

G. C. Greubel, Table of n, a(n) for n = 0..1000 (terms 0 to 200 by T. D. Noe)
M. Albert, R. Aldred, M. Atkinson, C Handley, D. Holton, D. McCaughan, and H. van Ditmarsch, Sorting Classes, Elec. J. of Comb. 12 (2005) R31.
Andrei Asinowski and Cyril Banderier, From geometry to generating functions: rectangulations and permutations, arXiv:2401.05558 [cs.DM], 2024. See page 2.
Jean-Luc Baril, José L. Ramírez, and Lina M. Simbaqueba, Counting prefixes of skew Dyck paths, J. Int. Seq., Vol. 24 (2021), Article 21.8.2.
Paul Barry, Series reversion with Jacobi and Thron continued fractions, arXiv:2107.14278 [math.NT], 2021.
Christian Bean, Finding structure in permutation sets, Ph.D. Dissertation, Reykjavík University, School of Computer Science, 2018.
Christian Bean, Émile Nadeau, and Henning Ulfarsson, Enumeration of Permutation Classes and Weighted Labelled Independent Sets, arXiv:1912.07503 [math.CO], 2019.
David Bevan, The permutation class Av(4213,2143), arXiv:1510.06328 [math.CO], 2014.
Jonathan Bloom and Alex Burstein, Egge triples and unbalanced Wilf-equivalence, arXiv:1410.0230 [math.CO], 2014.
Miklos Bona, Long increasing subsequences and non-algebraicity, arXiv:2310.13649 [math.CO], 2023.
Robert Brignall, Sophie Huczynska, and Vincent Vatter, Simple permutations and algebraic generating functions, arXiv:math/0608391 [math.CO], 2006.
Robert Brignall and Jakub Sliacan, Juxtaposing Catalan permutation classes with monotone ones, arXiv:1611.05370 [math.CO], 2016.
Colin Defant and Kai Zheng, Stack-Sorting with Consecutive-Pattern-Avoiding Stacks, arXiv:2008.12297 [math.CO], 2020.
Isaac DeJager, Madeleine Naquin, and Frank Seidl, Colored Motzkin Paths of Higher Order, VERUM 2019.
Emeric Deutsch, Emanuele Munarini, and Simone Rinaldi, Skew Dyck paths, J. Stat. Plann. Infer. 140 (8) (2010) 2191-2203.
Rigoberto Flórez, Leandro Junes, Luisa M. Montoya, and José L. Ramírez, Counting Subwords in Non-Decreasing Dyck Paths, J. Int. Seq. (2025) Vol. 28, Art. No. 25.1.6. See p. 19.
Alice L. L. Gao and Sergey Kitaev, On partially ordered patterns of length 4 and 5 in permutations, arXiv:1903.08946 [math.CO], 2019.
Alice L. L. Gao and Sergey Kitaev, On partially ordered patterns of length 4 and 5 in permutations, The Electronic Journal of Combinatorics 26(3) (2019), P3.26.
Elizabeth Hartung, Hung Phuc Hoang, Torsten Mütze, and Aaron Williams, Combinatorial generation via permutation languages. I. Fundamentals, arXiv:1906.06069 [cs.DM], 2019.
John W. Layman, The Hankel Transform and Some of its Properties, J. Integer Sequences, 4 (2001), #01.1.5.
Toufik Mansour and Mark Shattuck, Nine classes of permutations enumerated by binomial transform of Fine's sequence, Discrete Applied Mathematics, Vol. 226, 31 July 2017, p. 94-105.
Megan A. Martinez and Carla D. Savage, Patterns in Inversion Sequences II: Inversion Sequences Avoiding Triples of Relations, arXiv:1609.08106 [math.CO], 2016.
Arturo Merino and Torsten Mütze, Combinatorial generation via permutation languages. III. Rectangulations, arXiv:2103.09333 [math.CO], 2021.
Sam Miner, Enumeration of several two-by-four classes, arXiv preprint arXiv:1610.01908 [math.CO], 2016.
N. J. A. Sloane, Transforms
Rebecca Smith and Vincent Vatter, A stack and a pop stack in series, arXiv:1303.1395 [math.CO], 2013.
Index entries for reversions of series

Cf. A000957, A002212, A007317, A128714, A214611.

a[0] := 1: a[1] := 1: a[2] := 2: for n from 3 to 23 do a[n] := ((13*n-5)*a[n-1]-(16*n-23)*a[n-2]+5*(n-2)*a[n-3])/2/(n+1) od;

f[n_] := Sum[Binomial[n, k]*g[n - k], {k, 0, n}]; g[n_] := Sum[(-1)^(m + n)(n + m)!/n!/m!(n - m + 1)/(n + 1), {m, 0, n}]; Table[ f[n], {n, 24}] (* Robert G. Wilson v, Nov 04 2005 *)

a(n):=sum(sum(binomial(n-m-1,k-1)*m/(k+m)*binomial(2*k+m-1,k+m-1),k,1,n-m),m,1,n-1)+1; /* Vladimir Kruchinin, May 12 2011 */

a(n)=1+sum(m=1,n-1,sum(k=1,n-m,binomial(n-m-1,k-1)/(k+m)* binomial(2*k+m-1,k+m-1)*m)) \\ Charles R Greathouse IV, Mar 06 2013

x='x+O('x^50); Vec(2/(1+x+sqrt(1-6*x+5*x^2))) \\ Altug Alkan, Oct 22 2015

Also REVERT transform of x*(2*x-1)/(x^2+x-1). - Olivier Gérard
G.f.: 2/(1 + x + sqrt(1 - 6*x + 5*x^2)).
D-finite with recurrence a(n) = ((13*n-5)*a(n-1) - (16*n-23)*a(n-2) + 5*(n-2)*a(n-3))/(2*(n+1)) (n>=3); a(0)=a(1)=1, a(2)=2. - Emeric Deutsch, Mar 21 2004
Binomial transform of Fine's sequence: a(n) = Sum_{k=0..n} binomial(n, k)*A000957(n-k).
G.f.: 1/(1-x-x^2/(1-3x-x^2/(1-3x-x^2/(1-3x-x^2/(1-... (continued fraction). - Paul Barry, Jun 15 2009
a(n) = Sum_{k=0..n} A091965(n,k)*(-2)^k. - Philippe Deléham, Nov 28 2009
a(n) = 1 + Sum_{m=1..n-1} Sum_{k=1..n-m} binomial(n-m-1, k-1)*(m/(k+m))*binomial(2*k+m-1, k+m-1). - Vladimir Kruchinin, May 12 2011
a(n) = upper left term in M^n, M = the production matrix:
  1, 1, 0, 0, 0, 0, 0, ...
  1, 2, 1, 0, 0, 0, 0, ...
  1, 2, 1, 1, 0, 0, 0, ...
  1, 2, 1, 2, 1, 0, 0, ...
  1, 2, 1, 2, 1, 1, 0, ...
  1, 2, 1, 2, 1, 2, 1, ...
  ...
- Gary W. Adamson, Jul 08 2011
a(n) ~ 5^(n+3/2)/(18*sqrt(Pi)*n^(3/2)). - Vaclav Kotesovec, Aug 09 2013
G.f.: 1/(1-x*C(x/(1-x))), where C(x) = g.f. for A000108(n). - Alexander Burstein, Oct 05 2014

More terms from Robert G. Wilson v, Nov 04 2005

Entry revised by N. J. A. Sloane, Aug 07 2006

A124575 Triangle read by rows: row n is the first row of the matrix M[n]^(n-1), where M[n] is the n X n tridiagonal matrix with main diagonal (2,4,4,...) and super- and subdiagonals (1,1,1,...).

Original entry on oeis.org

1, 2, 1, 5, 6, 1, 16, 30, 10, 1, 62, 146, 71, 14, 1, 270, 717, 444, 128, 18, 1, 1257, 3582, 2621, 974, 201, 22, 1, 6096, 18206, 15040, 6718, 1800, 290, 26, 1, 30398, 93960, 85084, 43712, 14208, 2986, 395, 30, 1, 154756, 491322, 478008, 274140, 103530

Offset: 0

Gary W. Adamson & Roger L. Bagula, Nov 05 2006

Column k=0 yields A033543 (2nd binomial transform of the sequence A000957(n+1)). Row sums yield A133158. [Corrected by Philippe Deléham, Oct 24 2007, Dec 05 2009]
Triangle T(n,k), 0 <= k <= n, read by rows given by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + 4*T(n-1,k) + T(n-1,k+1) for k >= 1. - Philippe Deléham, Mar 27 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1)  +y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise from choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007

			Row 2 is (5,6,1) because M[3]= [2,1,0;1,4,1;0,1,4] and M[3]^2=[5,6,1;6,18,8;1,8,17].
Triangle starts:
    1;
    2,   1;
    5,   6,   1;
   16,  30,  10,   1;
   62, 146,  71,  14,  1;
  270, 717, 444, 128, 18, 1;

G. C. Greubel, Table of n, a(n) for the first 100 rows, flattened

Cf. A124576, A124574, A052179, A064613, A133158.

with(linalg): m:=proc(i,j) if i=1 and j=1 then 2 elif i=j then 4 elif abs(i-j)=1 then 1 else 0 fi end: for n from 3 to 11 do A[n]:=matrix(n,n,m): B[n]:=multiply(seq(A[n],i=1..n-1)) od: 1; 2,1; for n from 3 to 11 do seq(B[n][1,j],j=1..n) od; # yields sequence in triangular form

M[n_] := SparseArray[{{1, 1} -> 2, Band[{2, 2}] -> 4, Band[{1, 2}] -> 1, Band[{2, 1}] -> 1}, {n, n}]; row[1] = {1}; row[n_] := MatrixPower[M[n], n-1] // First // Normal; Table[row[n], {n, 1, 10}] // Flatten (* Jean-François Alcover, Jan 09 2014 *)

T(n,k) = T(n-1,k-1) + 4*T(n-1,k) + T(n-1,k-1) for k >= 2.
Sum_{k=0..n} T(n,k)*(3*k+1) = 6^n. - Philippe Deléham, Mar 27 2007
Sum_{k>=0} T(m,k)*T(n,k) = T(m+n,0) = A033543(m+n). - Philippe Deléham, Nov 22 2009

Edited by N. J. A. Sloane, Dec 04 2006

A126075 Triangle T(n,k), 0 <= k <= n, read by rows, defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k+1) for k >= 1.

Original entry on oeis.org

1, 2, 1, 5, 2, 1, 12, 6, 2, 1, 30, 14, 7, 2, 1, 74, 37, 16, 8, 2, 1, 185, 90, 45, 18, 9, 2, 1, 460, 230, 108, 54, 20, 10, 2, 1, 1150, 568, 284, 128, 64, 22, 11, 2, 1, 2868, 1434, 696, 348, 150, 75, 24, 12, 2, 1

Offset: 0

Philippe Deléham, Mar 02 2007

Riordan array (c(x^2)/(1-2xc(x^2)),xc(x^2)) where c(x)=g.f. of Catalan numbers A000108. - Philippe Deléham, Mar 18 2007

This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007

			Triangle begins:
     1;
     2,    1;
     5,    2,   1;
    12,    6,   2,   1;
    30,   14,   7,   2,   1;
    74,   37,  16,   8,   2,  1;
   185,   90,  45,  18,   9,  2,  1;
   460,  230, 108,  54,  20, 10,  2,  1;
  1150,  568, 284, 128,  64, 22, 11,  2, 1;
  2868, 1434, 696, 348, 150, 75, 24, 12, 2, 1;

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
P. Bala, A 4-parameter family of embedded Riordan arrays

Cf. A000108, A054341, A127358.

A126075 := proc (n, k)
add( 2^(n-k-2*j)*binomial(n, j), j = 0..floor((n-k)/2) ) - add( 2^(n-k-2-2*j)*binomial(n, j), j = 0..floor((n-k-2)/2) )
end proc:
# display sequence in triangular form
for n from 0 to 10 do seq(A126075(n, k), k = 0..n) end do;
# Peter Bala, Feb 20 2018

T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0, T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]]; Table[T[n, k, 2, 0], {n, 0, 49}, {k, 0, n}] // Flatten  (* G. C. Greubel, Apr 21 2017 *)

Sum_{k=0..n} T(n,k) = A127358(n). T(n,0)=A054341(n).
Sum_{k=0..n} T(n,k)*(-k+1) = 2^n. - Philippe Deléham, Mar 25 2007
From Peter Bala, Feb 20 2018: (Start)
T(n,k) = Sum_{j = 0..floor((n-k)/2)} 2^(n-k-2*j)*binomial(n, j) - Sum_{j = 0..floor((n-k-2)/2)} 2^(n-k-2-2*j)*binomial(n, j), 0 <= k <= n. - Peter Bala, Feb 20 2018
The n-th row polynomial in descending powers of x is the n-th Taylor polynomial of the rational function (1 - x^2)/(1 - 2*x) * (1 + x^2)^n about 0. For example, for n = 4, (1 - x^2)/(1 - 2*x) * (1 + x^2)^4 = (30*x^4 + 14*x*3 + 7*x^2 + 2*x + 1) + O(x^5). (End)

A110877 Triangle T(n,k), 0 <= k <= n, read by rows, defined by: T(0,0) = 1, T(n,k) = 0 if n= 1: T(n,k) = T(n-1,k-1) + x*T(n-1,k) + T(n-1,k+1) with x = 3.

Original entry on oeis.org

1, 1, 1, 2, 4, 1, 6, 15, 7, 1, 21, 58, 37, 10, 1, 79, 232, 179, 68, 13, 1, 311, 954, 837, 396, 108, 16, 1, 1265, 4010, 3861, 2133, 736, 157, 19, 1, 5275, 17156, 17726, 10996, 4498, 1226, 215, 22, 1, 22431, 74469, 81330, 55212, 25716, 8391, 1893

Offset: 0

Philippe Deléham, Sep 19 2005

Similar to A064189 (x = 1) and to A039599 (x = 2).
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
Row sums yield A126568. - Philippe Deléham, Oct 10 2007
5^n = (n-th row terms) dot (first n+1 terms in the series (1, 4, 7, 10, ...)). Example for row 4: 5^4 = 625 = (21, 58, 37, 10, 1) dot (1, 4, 7, 10, 13) = (21 + 232 + 259 + 100 + 13). - Gary W. Adamson, Jun 15 2011
Riordan array (2/(1+x+sqrt(1-6*x+5*x^2)), (1-3*x-sqrt(1-6*x+5*x^2))/(2*x)). - Philippe Deléham, Mar 04 2013

			Triangle begins:
      1;
      1,     1;
      2,     4,     1;
      6,    15,     7,     1;
     21,    58,    37,    10,     1;
     79,   232,   179,    68,    13,    1;
    311,   954,   837,   396,   108,   16,    1;
   1265,  4010,  3861,  2133,   736,  157,   19,   1;
   5275, 17156, 17726, 10996,  4498, 1226,  215,  22,  1;
  22431, 74469, 81330, 55212, 25716, 8391, 1893, 282, 25, 1;
  ...
From _Philippe Deléham_, Nov 07 2011: (Start)
Production matrix begins:
  1, 1;
  1, 3, 1;
  0, 1, 3, 1;
  0, 0, 1, 3, 1;
  0, 0, 0, 1, 3, 1;
  0, 0, 0, 0, 1, 3, 1;
  0, 0, 0, 0, 0, 1, 3, 1;
  0, 0, 0, 0, 0, 0, 1, 3, 1;
  0, 0, 0, 0, 0, 0, 0, 1, 3, 1;
  ... (End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

Cf. A002212, A033321, A039599, A064189.

The inverse of A126126.

A110877 := proc(n,k)
    if k > n then
        0;
    elif n= 0 then
        1;
    elif k = 0 then
        procname(n-1,0)+procname(n-1,1) ;
    else
        procname(n-1,k-1)+3*procname(n-1,k)+procname(n-1,k+1) ;
    end if;
end proc: # R. J. Mathar, Sep 06 2013

T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0, T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]]; Table[T[n, k, 1, 3], {n, 0, 49}, {k, 0, n}] // Flatten (* G. C. Greubel, Apr 21 2017 *)

T(n, 0) = A033321(n) and for k >= 1: T(n, k) = Sum_{j>=1} T(n-j, k-1)*A002212(j).
Sum_{k=0..n} T(m, k)*T(n, k) = T(m+n, 0) = A033321(m+n).
The triangle may also be generated from M^n * [1,0,0,0,...], where M = an infinite tridiagonal matrix with 1's in the super and subdiagonals and (1,3,3,3,...) in the main diagonal. - Gary W. Adamson, Dec 17 2006
Sum_{k=0..n} T(n,k)*(3*k+1) = 5^n. - Philippe Deléham, Feb 26 2007
Sum_{k=0..n} T(n,k) = A126568(n). - Philippe Deléham, Oct 10 2007

A125906 Riordan array (1/(1 + 5x + x^2), x/(1 + 5x + x^2))^(-1); inverse of Riordan array A123967.

Original entry on oeis.org

1, 5, 1, 26, 10, 1, 140, 77, 15, 1, 777, 540, 153, 20, 1, 4425, 3630, 1325, 254, 25, 1, 25755, 23900, 10509, 2620, 380, 30, 1, 152675, 155764, 79065, 23989, 4550, 531, 35, 1, 919139, 1010560, 575078, 203560, 47270, 7240, 707, 40, 1

Offset: 0

Philippe Deléham, Feb 04 2007

T(0)=A053121, T(1)=A064189, T(2)=A039598, T(3)=A091965, T(4)=A052179.
Triangle read by rows: T(n,k) = number of lattice paths from (0,0) to (n,k) that do not go below the line y=0 and consist of steps U=(1,1), D=(1,-1) and five types of steps H=(1,0); example: T(3,1)=77 because we have UDU, UUD, 25 HHU paths, 25 HUH paths and 25 UHH paths. - Philippe Deléham, Sep 25 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise from choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
7^n = (n-th row terms) dot (first n+1 terms in 1,2,3,...). Example: 7^3 = 343 = (140, 77, 15, 1) dot (1, 2, 3, 4) = (140 + 154 + 45 + 4) = 343. - Gary W. Adamson, Jun 17 2011
A subset of the "family of triangles" (Deleham comment of Sep 25 2007) is the succession of binomial transforms beginning with triangle A053121, (0,0); giving -> A064189, (1,1); -> A039598, (2,2); -> A091965, (3,3); -> A052179, (4,4); -> A125906, (5,5) ->, etc; generally the binomial transform of the triangle generated from (n,n) = that generated from ((n+1),(n+1)). - Gary W. Adamson, Aug 03 2011
Riordan array (f(x), x*f(x)) where f(x) is the o.g.f. of A182401. - Philippe Deléham, Mar 04 2013

			Triangle begins
       1;
       5,       1;
      26,      10,      1;
     140,      77,     15,      1;
     777,     540,    153,     20,     1;
    4425,    3630,   1325,    254,    25,    1;
   25755,   23900,  10509,   2620,   380,   30,   1;
  152675,  155764,  79065,  23989,  4550,  531,  35,  1;
  919139, 1010560, 575078, 203560, 47270, 7240, 707, 40, 1;
From _Philippe Deléham_, Nov 07 2011: (Start)
Production matrix begins
  5, 1;
  1, 5, 1,;
  0, 1, 5, 1;
  0, 0, 1, 5, 1;
  0, 0, 0, 1, 5, 1;
  0, 0, 0, 0, 1, 5, 1;
  0, 0, 0, 0, 0, 1, 5, 1;
  0, 0, 0, 0, 0, 0, 1, 5, 1;
  0, 0, 0, 0, 0, 0, 0, 1, 5, 1; (End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

Cf. A182401.

T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0,  T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]];
Table[T[n, k, 5, 5], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)

Triangle T(5) where T(x) is defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,k) = T(n-1,k-1) + x*T(n-1,k) + T(n-1,k+1). Sum_{k=0..n} T(m,k)*T(n,k) = T(m+n,0). Sum_{k=0..n} T(n,k) = A122898(n).
Sum_{k=0..n} T(n,k)*(k+1) = 7^n. - Philippe Deléham, Mar 26 2007
T(n,0) = A182401(n). - Philippe Deléham, Mar 04 2013
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 - x^2)*(1 + 5*x + x^2)^n expanded about the point x = 0. - Peter Bala, Sep 06 2022

A126093 Inverse binomial matrix applied to A110877.

Original entry on oeis.org

1, 0, 1, 1, 2, 1, 2, 6, 4, 1, 6, 18, 15, 6, 1, 18, 57, 54, 28, 8, 1, 57, 186, 193, 118, 45, 10, 1, 186, 622, 690, 474, 218, 66, 12, 1, 622, 2120, 2476, 1856, 976, 362, 91, 14, 1, 2120, 7338, 8928, 7164, 4170, 1791, 558, 120, 16, 1

Offset: 0

Philippe Deléham, Mar 03 2007

Diagonal sums are A065601. - Philippe Deléham, Mar 05 2007

This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007

			Triangle begins:
     1;
     0,    1;
     1,    2,    1;
     2,    6,    4,    1;
     6,   18,   15,    6,    1;
    18,   57,   54,   28,    8,    1;
    57,  186,  193,  118,   45,   10,   1;
   186,  622,  690,  474,  218,   66,  12,   1;
   622, 2120, 2476, 1856,  976,  362,  91,  14,  1;
  2120, 7338, 8928, 7164, 4170, 1791, 558, 120, 16, 1;
Production matrix begins
  0, 1;
  1, 2, 1;
  0, 1, 2, 1;
  0, 0, 1, 2, 1;
  0, 0, 0, 1, 2, 1;
  0, 0, 0, 0, 1, 2, 1;
  0, 0, 0, 0, 0, 1, 2, 1;
  0, 0, 0, 0, 0, 0, 1, 2, 1;
  0, 0, 0, 0, 0, 0, 0, 1, 2, 1;
- _Philippe Deléham_, Nov 07 2011

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Yidong Sun and Luping Ma, Minors of a class of Riordan arrays related to weighted partial Motzkin paths, Eur. J. Comb. 39, 157-169 (2014), Table 2.2.

T[0, 0, x_, y_]:= 1; T[n_, 0, x_, y_]:= x*T[n-1,0,x,y] + T[n-1,1,x,y]; T[n_, k_, x_, y_]:= T[n, k, x, y]= If[k<0 || k>n, 0, T[n-1,k-1,x,y] + y*T[n-1,k,x,y] + T[n-1,k+1,x,y]]; Table[T[n,k,0,2], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Apr 21 2017 *)

@CachedFunction
def T(n, k, x, y):
    if (k<0 or k>n): return 0
    elif (n==0 and k==0): return 1
    elif (k==0): return x*T(n-1,0,x,y) + T(n-1,1,x,y)
    else: return T(n-1,k-1,x,y) + y*T(n-1,k,x,y) + T(n-1,k+1,x,y)
[[T(n,k,0,2) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Jan 27 2020

Triangle T(n,k), 0<=k<=n, read by rows defined by : T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0) = T(n-1,1), T(n,k) = T(n-1,k-1) + 2*T(n-1,k) + T(n-1,k+1) for k>=1.
Sum_{k=0..n} T(m,k)*T(n,k) = T(m+n,0) = A000957(m+n+1).
Sum_{k=0..n-1} T(n,k) = A026641(n), for n>=1. - Philippe Deléham, Mar 05 2007
Sum_{k=0..n} T(n,k)*(3k+1) = 4^n. - Philippe Deléham, Mar 22 2007

A126954 Triangle T(n,k), 0 <= k <= n, read by rows given by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 3*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k) + T(n-1,k+1) for k >= 1.

Original entry on oeis.org

1, 3, 1, 10, 4, 1, 34, 15, 5, 1, 117, 54, 21, 6, 1, 405, 192, 81, 28, 7, 1, 1407, 678, 301, 116, 36, 8, 1, 4899, 2386, 1095, 453, 160, 45, 9, 1, 17083, 8380, 3934, 1708, 658, 214, 55, 10, 1, 59629, 29397, 14022, 6300, 2580, 927, 279, 66, 11, 1

Offset: 0

Philippe Deléham, Mar 19 2007

This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = x*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + y*T(n-1,k) + T(n-1,k+1) for k >= 1. Other triangles arise from choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007

			Triangle begins:
     1;
     3,    1;
    10,    4,    1;
    34,   15,    5,   1;
   117,   54,   21,   6,   1;
   405,  192,   81,  28,   7,  1;
  1407,  678,  301, 116,  36,  8, 1;
  4899, 2386, 1095, 453, 160, 45, 9, 1;

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0, T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]]; Table[T[n, k, 3, 1], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, May 22 2017 *)

Sum_{k=0..n} T(n,k) = A126932(n).
Sum_{k>=0} T(m,k)*T(n,k) = T(m+n,0) = A059738(m+n).
Sum_{k=0..n} T(n,k)*(-k+1) = 3^n. - Philippe Deléham, Mar 26 2007

A124733 Triangle read by rows: row n is the first row of the matrix M[n]^(n-1), where M[n] is the n X n tridiagonal matrix with main diagonal (2,3,3,...) and super- and subdiagonals (1,1,1,...).

Original entry on oeis.org

1, 2, 1, 5, 5, 1, 15, 21, 8, 1, 51, 86, 46, 11, 1, 188, 355, 235, 80, 14, 1, 731, 1488, 1140, 489, 123, 17, 1, 2950, 6335, 5397, 2730, 875, 175, 20, 1, 12235, 27352, 25256, 14462, 5530, 1420, 236, 23, 1, 51822, 119547, 117582, 74172, 32472, 10026, 2151, 306, 26, 1

Offset: 1

Gary W. Adamson & Roger L. Bagula, Nov 05 2006

With a different offset: Triangle T(n,k), 0<=k<=n, read by rows given by : T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=2*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+3*T(n-1,k)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 27 2007
Equals A007318*A039599 (when written as lower triangular matrix). - Philippe Deléham, Jun 16 2007
This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007
5^n = (n-th row terms) dot (first n+1 odd integers). Example: 5^4 = 625 = (51, 86, 46, 11, 1) dot (1, 3, 5, 7, 9) = (51 + 258 + 230 + 77 + 9) = 625. [Gary W. Adamson, Jun 13 2011]

			Row 3 is (5,5,1) because M[3]=[2,1,0;1,3,1;0,1,3] and M[3]^2=[5,5,1;5,11,6;1,6,10].
Triangle starts:
1;
2, 1;
5, 5, 1;
15, 21, 8, 1;
51, 86, 46, 11, 1;
188, 355, 235, 80, 14, 1;

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Shu-Chiuan Chang and Robert Shrock, Structure of the Partition Function and Transfer Matrices for the Potts Model in a Magnetic Field on Lattice Strips, Journal of Statistical Physics 137 (2009) 667.

Cf. A110877, A091965, A002212, A007317, A026375 (row sums).

with(linalg): m:=proc(i,j) if i=1 and j=1 then 2 elif i=j then 3 elif abs(i-j)=1 then 1 else 0 fi end: for n from 3 to 11 do A[n]:=matrix(n,n,m): B[n]:=multiply(seq(A[n],i=1..n-1)) od: 1; 2,1; for n from 3 to 11 do seq(B[n][1,j],j=1..n) od; # yields sequence in triangular form
T := (n,k) -> (-1)^(n-k)*simplify(GegenbauerC(n-k,-n+1,3/2) + GegenbauerC(n-k-1,-n+1,3/2)): seq(seq(T(n,k), k=1..n), n=1..10); # Peter Luschny, May 13 2016

T[0, 0, x_, y_] := 1; T[n_, 0, x_, y_] := x*T[n - 1, 0, x, y] + T[n - 1, 1, x, y]; T[n_, k_, x_, y_] := T[n, k, x, y] = If[k < 0 || k > n, 0,  T[n - 1, k - 1, x, y] + y*T[n - 1, k, x, y] + T[n - 1, k + 1, x, y]];
Table[T[n, k, 2, 3], {n, 0, 49}, {k, 0, n}] // Flatten (* G. C. Greubel, Apr 21 2017 *)

Sum_{k=0..n} (-1)^(n-k)*T(n,k) = (-1)^n. - Philippe Deléham, Feb 27 2007
Sum_{k=0..n} T(n,k)*(2*k+1) = 5^n. - Philippe Deléham, Mar 27 2007
T(n,k) = (-1)^(n-k)*(GegenbauerC(n-k,-n+1,3/2) + GegenbauerC(n-k-1,-n+1,3/2)). - Peter Luschny, May 13 2016
From Peter Bala, Sep 06 2022: (Start)
The following assume the row and column indexing start at 0.
Riordan array (f(x), x*g(x)), where f(x) = ( 1 - sqrt((1 - 5*x)/(1 - x)) )/(2*x) = 1 + 2*x + 5*x^2 + 15*x^3 + 51*x^4 + ... is the o.g.f. of A007317 and g(x) = ( 1 - 3*x - sqrt(1 - 6*x + 5*x^2) )/(2*x^2) = 1 + 3*x + 10*x^2 + 36*x^3 + 137*x^4 + .... See A002212.
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 - x)*(1 + 3*x + x^2)^n expanded about the point x = 0.
T(n,k) = a(n,k) - a(n,k+1), where a(n,k) = Sum_{j = 0..n} binomial(n,j)* binomial(j,n-k-j)*3^(2*j+k-n). (End)

Edited by N. J. A. Sloane, Dec 04 2006

A124576 Triangle read by rows: row n is the first row of the matrix M[n]^(n-1), where M[n] is the n X n tridiagonal matrix with main diagonal (1,4,4,...) and super- and subdiagonals (1,1,1,...).

Original entry on oeis.org

1, 1, 1, 2, 5, 1, 7, 23, 9, 1, 30, 108, 60, 13, 1, 138, 522, 361, 113, 17, 1, 660, 2587, 2079, 830, 182, 21, 1, 3247, 13087, 11733, 5581, 1579, 267, 25, 1, 16334, 67328, 65600, 35636, 12164, 2672, 368, 29, 1, 83662, 351246, 365364, 220308, 86964, 23220, 4173

Offset: 1

Gary W. Adamson & Roger L. Bagula, Nov 05 2006

Triangle T(n,k), 0<=k<=n, read by rows given by : T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+4*T(n-1,k)+T(n-1,k+1) for k>=1. - Philippe Deléham, Mar 27 2007

This triangle belongs to the family of triangles defined by: T(0,0)=1, T(n,k)=0 if k<0 or if k>n, T(n,0)=x*T(n-1,0)+T(n-1,1), T(n,k)=T(n-1,k-1)+y*T(n-1,k)+T(n-1,k+1) for k>=1 . Other triangles arise by choosing different values for (x,y): (0,0) -> A053121; (0,1) -> A089942; (0,2) -> A126093; (0,3) -> A126970; (1,0)-> A061554; (1,1) -> A064189; (1,2) -> A039599; (1,3) -> A110877; (1,4) -> A124576; (2,0) -> A126075; (2,1) -> A038622; (2,2) -> A039598; (2,3) -> A124733; (2,4) -> A124575; (3,0) -> A126953; (3,1) -> A126954; (3,2) -> A111418; (3,3) -> A091965; (3,4) -> A124574; (4,3) -> A126791; (4,4) -> A052179; (4,5) -> A126331; (5,5) -> A125906. - Philippe Deléham, Sep 25 2007

			Row 3 is (2,5,1) because M[3]=[1,1,0;1,4,1;0,1,4] and M[3]^2=[2,5,1;5,18,8;1,8,17].
Triangle starts:
1;
1, 1;
2, 5, 1;
7, 23, 9, 1;
30, 108, 60, 13, 1;
138, 522, 361, 113, 17, 1;

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened

Cf. A124575, A124574, A052179, A227081 (row sums).

with(linalg): m:=proc(i,j) if i=1 and j=1 then 1 elif i=j then 4 elif abs(i-j)=1 then 1 else 0 fi end: for n from 3 to 11 do A[n]:=matrix(n,n,m): B[n]:=multiply(seq(A[n],i=1..n-1)) od: 1; 1,1; for n from 3 to 11 do seq(B[n][1,j],j=1..n) od; # yields sequence in triangular form
# alternative
A124576_row := proc(n)
    if n = 0 then
        return [1] ;
    else
        M := Matrix(n,n) ;
        M[1,1] := 1;
        for c from 2 to n do
            if c = 2 then
                M[1,c] := 1;
            else
                M[1,c] := 0;
            end if;
        end do:
        for r from 2 to n do
            for c from 1 to n do
                if r = c then
                    M[r,c] := 4;
                elif abs(r-c) = 1 then
                    M[r,c] := 1;
                else
                    M[r,c] := 0;
                end if;
            end do:
        end do:
        LinearAlgebra[MatrixPower](M,n-1) ;
        return [seq(%[1,r],r=1..n)] ;
    end if;
end proc:
for n from 0 to 10 do
    A124576_row(n) ;
    print(%) ;
end do: # R. J. Mathar, May 20 2025

M[n_] := SparseArray[{{1, 1} -> 1, Band[{2, 2}] -> 4, Band[{1, 2}] -> 1, Band[{2, 1}] -> 1}, {n, n}]; row[1] = {1}; row[n_] := MatrixPower[M[n], n-1] // First // Normal; Table[row[n], {n, 1, 10}] // Flatten (* Jean-François Alcover, Jan 09 2014 *)

Sum_{k=0..n} T(n,k)*(4*k+1) = 6^n. - Philippe Deléham, Mar 27 2007

Edited by N. J. A. Sloane, Dec 04 2006

cp's OEIS Frontend

A111418 Right-hand side of odd-numbered rows of Pascal's triangle.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Formula

A033321 Binomial transform of Fine's sequence A000957: 1, 0, 1, 2, 6, 18, 57, 186, ...

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Maxima

PARI

PARI

Formula

Extensions

A124575 Triangle read by rows: row n is the first row of the matrix M[n]^(n-1), where M[n] is the n X n tridiagonal matrix with main diagonal (2,4,4,...) and super- and subdiagonals (1,1,1,...).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A126075 Triangle T(n,k), 0 <= k <= n, read by rows, defined by: T(0,0)=1, T(n,k)=0 if k < 0 or if k > n, T(n,0) = 2*T(n-1,0) + T(n-1,1), T(n,k) = T(n-1,k-1) + T(n-1,k+1) for k >= 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A110877 Triangle T(n,k), 0 <= k <= n, read by rows, defined by: T(0,0) = 1, T(n,k) = 0 if n= 1: T(n,k) = T(n-1,k-1) + x*T(n-1,k) + T(n-1,k+1) with x = 3.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A125906 Riordan array (1/(1 + 5*x + x^2), x/(1 + 5*x + x^2))^(-1); inverse of Riordan array A123967.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A126093 Inverse binomial matrix applied to A110877.

Views

Author

Keywords

A125906 Riordan array (1/(1 + 5x + x^2), x/(1 + 5x + x^2))^(-1); inverse of Riordan array A123967.