A008288 -id:A008288 - cp's OEIS frontend

A104698 Triangle read by rows: T(n,k) = Sum_{j=0..n-k} binomial(k, j)*binomial(n-j+1, k+1).

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

Offset: 0

Gary W. Adamson, Mar 19 2005

The n-th column of the triangle is the binomial transform of the n-th row of A081277, followed by zeros. Example: column 3, (1, 6, 19, 44, ...) = binomial transform of row 3 of A081277: (1, 5, 8, 4, 0, 0, 0, ...). A104698 = reversal by rows of A142978. - Gary W. Adamson, Jul 17 2008
This sequence is jointly generated with A210222 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) + 1 and v(n,x) = 2x*u(n-1,x) + v(n-1,x) + 1. See the Mathematica section at A210222. - Clark Kimberling, Mar 19 2012
This Riordan triangle T appears in a formula for A001100(n, 0) = A002464(n), for n >= 1. - Wolfdieter Lang, May 13 2025

			The Riordan triangle T begins:
  n\k  0   1   2    3    4    5    6   7   8  9 10 ...
  ----------------------------------------------------
  0:   1
  1:   2   1
  2:   3   4   1
  3:   4   9   6    1
  4:   5  16  19    8    1
  5:   6  25  44   33   10    1
  6:   7  36  85   96   51   12    1
  7:   8  49 146  225  180   73   14   1
  8:   9  64 231  456  501  304   99  16   1
  9:  10  81 344  833 1182  985  476 129  18  1
  10: 11 100 489 1408 2471 2668 1765 704 163 20  1
  ... reformatted and extended by _Wolfdieter Lang_, May 13 2025
From _Wolfdieter Lang_, May 13 2025: (Start)
Zumkeller recurrence (adapted for offset [0,0]): 19 = T(4, 2) = T(2, 1) + T(3, 1) + T(3,3) = 4 + 9 + 6 = 19.
A-sequence recurrence: 19 = T(4, 2) = 1*T(3. 1) + 2*T(3. 2) - 2*T(3, 3) = 9 + 12 - 2 = 19.
Z-sequence recurrence: 5 = T(4, 0) = 2*T(3, 0) - 1*T(3, 1) + 2*T(3, 2) - 6*T(3, 3) = 8 - 9 + 12 + 6 = 5.
Boas-Buck recurrence: 19 = T(4, 2) = (1/2)*((2 + 0)*T(2, 2) + (2 + 2*2)*T(3, 2)) = (1/2)*(2 + 36) = 19. (End)

Reinhard Zumkeller, First 100 rows of triangle, flattened

Diagonal sums are A008937(n+1).
Cf. A048739 (row sums), A008288, A005900 (column 3), A014820 (column 4)
Cf. A081277, A142978 by antidiagonals, A119328, A110271 (matrix inverse).
Cf. A001100, A002464, A010673, A040000, A112478, A133872, A383857.

a104698 n k = a104698_tabl !! (n-1) !! (k-1)
a104698_row n = a104698_tabl !! (n-1)
a104698_tabl = [1] : [2,1] : f [1] [2,1] where
   f us vs = ws : f vs ws where
     ws = zipWith (+) ([0] ++ us ++ [0]) $
          zipWith (+) ([1] ++ vs) (vs ++ [0])
-- Reinhard Zumkeller, Jul 17 2015

A104698 := proc(n, k) add(binomial(k, j)*binomial(n-j+1, n-k-j), j=0..n-k) ; end proc:
seq(seq(A104698(n, k), k=0..n), n=0..15); # R. J. Mathar, Sep 04 2011
T := (n, k) -> binomial(n + 1, k + 1)*hypergeom([-k, k - n], [-n - 1], -1):
for n from 0 to 9 do seq(simplify(T(n, k)), k = 0..n) od;
T := proc(n, k) option remember; if k = 0 then n + 1 elif k = n then 1 else T(n-2, k-1) + T(n-1, k-1) + T(n-1, k) fi end: # Peter Luschny, May 13 2025

u[1, ] = 1; v[1, ] = 1;
u[n_, x_] := u[n, x] = x u[n-1, x] + v[n-1, x] + 1;
v[n_, x_] := v[n, x] = 2 x u[n-1, x] + v[n-1, x] + 1;
Table[CoefficientList[u[n, x], x], {n, 1, 11}] // Flatten (* Jean-François Alcover, Mar 10 2019, after Clark Kimberling *)

T(n,k)=sum(j=0,n-k,binomial(k,j)*binomial(n-j+1,k+1)) \\ Charles R Greathouse IV, Jan 16 2012

The triangle is extracted from the product A * B; A = [1; 1, 1; 1, 1, 1; ...], B = [1; 1, 1; 1, 3, 1; 1, 5, 5, 1; ...] both infinite lower triangular matrices (rest of the terms are zeros). The triangle of matrix B by rows = A008288, Delannoy numbers.
From Paul Barry, Jul 18 2005: (Start)
Riordan array (1/(1-x)^2, x(1+x)/(1-x)) = (1/(1-x), x)*(1/(1-x), x(1+x)/(1-x)).
T(n, k) = Sum_{j=0..n} Sum_{i=0..j-k} C(j-k, i)*C(k, i)*2^i.
T(n, k) = Sum_{j=0..k} Sum_{i=0..n-k-j} (n-k-j-i+1)*C(k, j)*C(k+i-1, i). (End)
T(n, k) = binomial(n+1, k+1)*2F1([-k, k-n], [-n-1], -1) where 2F1 is a Gaussian hypergeometric function. - R. J. Mathar, Sep 04 2011
T(n, k) = T(n-2, k-1) + T(n-1, k-1) + T(n-1, k) for 1 < k < n; T(n, 0) = n + 1; T(n, n) = 1. - Reinhard Zumkeller, Jul 17 2015
From Wolfdieter Lang, May 13 2025: (Start)
The Riordan triangle T = (1/(1 - x)^2, x*(1 + x)/(1 - x)) has the o.g.f. G(x, y) = 1/((1 - x)*(1 - x - y*x*(1+x))) for the row polynomials R(n, y) = Sum_{k=0..n} T(n, k)*y^k.
The o.g.f. for column k is G(k, x) = (1/(1 - x)^2)*(x*(1 + x)/(1 - x))^k, for k >= 0.
The o.g.f. for the diagonal m is D(m, x) = N(m, x)/(1 - x)^(m+1), with the numerator polynomial N(m, x) = Sum_{k=0..floor(m/2)} A034867(m, k)*x^(2*k) for m >= 0.
The row sums with o.g.f. R(x) = 1/((1 -x)*(1 - 2*x -x^2) give A048739.
The alternating row sums with o.g.f. 1/((1 - x)(1 + x^2)) give A133872.
The A-sequence for this Riordan triangle has o.g.f. A(x) = 1 + x + sqrt(1 + 6*x + x^2))/2 giving A112478(n). Hence T(n, k) = Sum_{j=0..n-k} A112478(j)*T(n-1, k-1+j), for n >= 1, k >= 1, T(n, k) = 0 for n < k, and T(0, 0) = 1.
The Z-sequence has o.g.f. (5 + x - sqrt(1 + 6*x + x^2))/2 = 3 + x - A(x) giving Z(n) = {2, -1, -A112478(n >= 2)}. Hence T(n, 0) = Sum_{j=0..n-1} Z(j)*T(n-1, j), for n >= 1. For A- and Z-sequences of Riordan triangles see a W. Lang link at A006232 with references.
The Boas-Buck sequences alpha and beta for the Riordan triangle T (see A046521 for the Aug 10 2017 comment and reference) are alpha(n) = A040000(n+1) = repeat{2} and beta(n) = A010673(n+1) = repeat{2,0}. Hence the recurrence for column T(n, k){n>=k}, with input T(k, k) = 1, for k >= 0, is T(n, k) = (1/(n-k)) * Sum{j=k..n-1} (2 + k*(1 + (-1)^(n-1-j))) *T(j,k), for n >= k+1. (End)

A145905 Square array read by antidiagonals: Hilbert transform of triangle A060187.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 9, 5, 1, 1, 27, 25, 7, 1, 1, 81, 125, 49, 9, 1, 1, 243, 625, 343, 81, 11, 1, 1, 729, 3125, 2401, 729, 121, 13, 1, 1, 2187, 15625, 16807, 6561, 1331, 169, 15, 1, 1, 6561, 78125, 117649, 59049, 14641, 2197, 225, 17, 1, 1, 19683, 390625, 823543

Offset: 0

Peter Bala, Oct 27 2008

Definition of the Hilbert transform of a triangular array:
For many square arrays in the database the entries in a row are polynomial in the column index, of degree d say and hence the row generating function has the form P(x)/(1-x)^(d+1), where P is some polynomial function. Often the array whose rows are formed from the coefficients of these P polynomials is of independent interest. This suggests the following definition.
Let [L(n,k)]n,k>=0 be a lower triangular array and let R(n,x) := sum {k = 0 .. n} L(n,k)*x^k, denote the n-th row generating polynomial of L. Then we define the Hilbert transform of L, denoted Hilb(L), to be the square array whose n-th row, n >= 0, has the generating function R(n,x)/(1-x)^(n+1).
In this particular case, L is the array A060187, the array of Eulerian numbers of type B, whose row polynomials are the h-polynomials for permutohedra of type B. The Hilbert transform is an infinite Vandermonde matrix V(1,3,5,...).
We illustrate the Hilbert transform with a few examples:
(1) The Delannoy number array A008288 is the Hilbert transform of Pascal's triangle A007318 (view as the array of coefficients of h-polynomials of n-dimensional cross polytopes).
(2) The transpose of the array of nexus numbers A047969 is the Hilbert transform of the triangle of Eulerian numbers A008292 (best viewed in this context as the coefficients of h-polynomials of n-dimensional permutohedra of type A).
(3) The sequence of Eulerian polynomials begins [1, x, x + x^2, x + 4*x^2 + x^3, ...]. The coefficients of these polynomials are recorded in triangle A123125, whose Hilbert transform is A004248 read as square array.
(4) A108625, the array of crystal ball sequences for the A_n lattices, is the Hilbert transform of A008459 (viewed as the triangle of coefficients of h-polynomials of n-dimensional associahedra of type B).
(5) A142992, the array of crystal ball sequences for the C_n lattices, is the Hilbert transform of A086645, the array of h-vectors for type C root polytopes.
(6) A108553, the array of crystal ball sequences for the D_n lattices, is the Hilbert transform of A108558, the array of h-vectors for type D root polytopes.
(7) A086764, read as a square array, is the Hilbert transform of the rencontres numbers A008290.
(8) A143409 is the Hilbert transform of triangle A073107.

			Triangle A060187 (with an offset of 0) begins
1;
1, 1;
1, 6, 1;
so the entries in the first three rows of the Hilbert transform of
A060187 come from the expansions:
Row 0: 1/(1-x) = 1 + x + x^2 + x^3 + ...;
Row 1: (1+x)/(1-x)^2 = 1 + 3*x + 5*x^2 + 7*x^3 + ...;
Row 2: (1+6*x+x^2)/(1-x)^3 = 1 + 9*x + 25*x^2 + 49*x^3 + ...;
The array begins
n\k|..0....1.....2.....3......4
================================
0..|..1....1.....1.....1......1
1..|..1....3.....5.....7......9
2..|..1....9....25....49.....81
3..|..1...27...125...343....729
4..|..1...81...625..2401...6561
5..|..1..243..3125.16807..59049
...

Ghislain R. Franssens, On a Number Pyramid Related to the Binomial, Deleham, Eulerian, MacMahon and Stirling number triangles, Journal of Integer Sequences, Vol. 9 (2006), Article 06.4.1.
S. Parker, The Combinatorics of Functional Composition and Inversion, Ph.D. Dissertation, Brandeis Univ. (1993) [From _Tom Copeland_, Nov 09 2008]

Cf. A008292, A039755, A052750 (first superdiagonal), A060187, A114172, A145901.

T:=(n,k) -> (2*k + 1)^n: seq(seq(T(n-k,k),k = 0..n),n = 0..10);

T(n,k) = (2*k + 1)^n, (see equation 4.10 in [Franssens]). This array is the infinite Vandermonde matrix V(1,3,5,7, ....) having a LDU factorization equal to A039755 * diag(2^n*n!) * transpose(A007318).

A047665 Expansion of (1/sqrt(1-6*x+x^2)-1/(1-x))/2.

Original entry on oeis.org

0, 1, 6, 31, 160, 841, 4494, 24319, 132864, 731281, 4048726, 22523359, 125797984, 704966809, 3961924126, 22321190911, 126027618304, 712917362209, 4039658528934, 22924714957471, 130271906898720, 741188107113961, 4221707080583086, 24070622500965631, 137369104574280960, 784622537295845041

Offset: 0

Don Knuth, N. J. A. Sloane

nonn

Previous name was: Main diagonal of square array defined in A047662.

a(n) is the total number of weak plateaus in all Schroeder n-paths. A weak plateau is a subpath of the form UFF..FD where there are 0 or more Fs. For example, a(2)=6 counts the following weak plateaus (in parentheses) in the 6 Schroeder 2-paths: (UFD), U(UD)D, FF, (UD)F, F(UD), (UD)(UD). - David Callan, Aug 16 2006

Vincenzo Librandi, Table of n, a(n) for n = 0..200 (corrected by Sean A. Irvine, Jan 18 2019)
Y. Ding and R. R. X. Du, Counting Humps in Motzkin Paths, arXiv:1109.2661 [math.CO], 2011, Eq. (4.2).
D. E. Knuth and N. J. A. Sloane, Correspondence, December 1999
Matthew Roughan, Surreal Birthdays and Their Arithmetic, arXiv:1810.10373 [math.HO], 2018.

Cf. A001850, A002002 (Schroeder paths interpretation).

Cf. A008288 (Delannoy numbers triangle).

seq(add(multinomial(n+k,n-k,k,k)/2,k=1..n),n=1..22); # Zerinvary Lajos, Oct 18 2006
a:=n->add(add(binomial(n,j)*binomial(n,k)*binomial(k,j), j=0..n),k=1..n): seq(a(n)/2, n=1..22); # Zerinvary Lajos, Jun 02 2007

Table[SeriesCoefficient[(1/Sqrt[1-6*x+x^2]-1/(1-x))/2,{x,0,n}],{n,1,20}] (* Vaclav Kotesovec, Oct 08 2012 *)

x='x+O('x^66); Vec((1/sqrt(1-6*x+x^2)-1/(1-x))/2) \\ Joerg Arndt, May 04 2013

a = lambda n: (hypergeometric([-n, n+1], [1], -1)-1)/2
[simplify(a(n)) for n in (1..25)] # Peter Luschny, May 19 2015

2*a(n)+1 = A001850(n).
a(n)-a(n-1) = A002002(n).
a(n) = Sum_{k=0..n} Sum_{j=0..n} A008288(k, j).
a(n) = Sum_{j=1..n} C(2*j-1, j-1)*C(n+j, 2*j). - Stefan Hollos (stefan(AT)exstrom.com), Jul 21 2004
D-finite with recurrence: n*(2*n-3)*a(n) = (2*n-1)*(7*n-10)*a(n-1) - (2*n-3)*(7*n-4)*a(n-2) + (n-2)*(2*n-1)*a(n-3). - Vaclav Kotesovec, Oct 08 2012
a(n) ~ sqrt(8+6*sqrt(2))*(3+2*sqrt(2))^n/(8*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 08 2012
a(n) = (hyper2F1(-n,n+1,1,-1)-1)/2 = (hyper2F1(-n, -n, 1, 2)-1)/2. - Peter Luschny, May 19 2015
a(n) = Sum_{k=1..n} binomial(n,k)^2 * 2^(k-1). - Ilya Gutkovskiy, Nov 15 2021

Prepended 0, set offset to 0 and new name using a comment of Emeric Deutsch from Dec 25 2003 by Peter Luschny, May 20 2015

A050146 a(n) = T(n,n), array T as in A050143.

Original entry on oeis.org

1, 1, 4, 18, 88, 450, 2364, 12642, 68464, 374274, 2060980, 11414898, 63521352, 354870594, 1989102444, 11180805570, 63001648608, 355761664002, 2012724468324, 11406058224594, 64734486343480, 367891005738690, 2093292414443164, 11923933134635298, 67990160422313808

Offset: 0

Clark Kimberling

nonn

Also main diagonal of array : m(i,1)=1, i>=1; m(1,j)=2, j>1; m(i,j)=m(i,j-1)+m(i-1,j-1)+m(i-1,j): 1 2 2 2 ... / 1 4 8 12 ... / 1 6 18 38 ... / 1 8 32 88 ... / - Benoit Cloitre, Aug 05 2002
a(n) is also the number of order-preserving partial transformations (of an n-element chain) of waist n (waist(alpha) = max(Im(alpha))). - Abdullahi Umar, Aug 25 2008
Define a finite triangle T(r,c) with T(r,0) = binomial(n,r) for 0<=r<=n, and the other terms recursively with T(r,c) = T(r,c-1) + 2*T(r-1,c-1). The sum of the last terms in each row is Sum_{r=0..n} T(r,r)=a(n+1). For n=4 the triangle is 1; 4 6; 6 14 26; 4 16 44 96; 1 9 41 129 321 with the sum of the last terms being 1 + 6 + 26 + 96 + 321 = 450 = a(5). - J. M. Bergot, Jan 29 2013
It may be better to define a(0) = 0 for formulas without exceptions. - Michael Somos, Nov 25 2016
a(n) is the number of points at L1 distance n-1 from any point in Z^n, for n>=1. - Shel Kaphan, Mar 24 2023

			G.f. = 1 + x + 4*x^2 + 18*x^3 + 88*x^4 + 450*x^5 + 2364*x^6 + 12642*x^7 + ...

Vincenzo Librandi, Table of n, a(n) for n = 0..200
A. Laradji and A. Umar, A. Combinatorial results for semigroups of order-preserving partial transformations, Journal of Algebra, 278 (2004), 342-359.
A. Laradji and A. Umar, Combinatorial results for semigroups of order-decreasing partial transformations, J. Integer Seq., 7 (2004), 04.3.8.
Huyile Liang, Yanni Pei, and Yi Wang, Analytic combinatorics of coordination numbers of cubic lattices, arXiv:2302.11856 [math.CO], 2023. See p. 4.
Emanuele Munarini, Combinatorial properties of the antichains of a garland, Integers, 9 (2009), 353-374.

Cf. A002003, A006318, A050151, A008288, A002002, A026002.

-1-diagonal of A266213 for n>=1.

a050146 n = if n == 0 then 1 else a035607 (2 * n - 2) (n - 1)
-- Reinhard Zumkeller, Nov 05 2013, Jul 20 2013

Flatten[{1,RecurrenceTable[{(n-3)*(n-1)*a[n-2]-3*(n-2)*(2*n-3)*a[n-1]+(n-2)*(n-1)*a[n]==0,a[1]==1,a[2]==4},a,{n,20}]}] (* Vaclav Kotesovec, Oct 08 2012 *)
a[ n_] := If[ n == 0, 1, Sum[ Binomial[n, k] Binomial[n + k - 2, k - 1], {k, n}]]; (* Michael Somos, Nov 25 2016 *)
a[ n_] := If[ n == 0, 1, n Hypergeometric2F1[1 - n, n, 2, -1]]; (* Michael Somos, Nov 25 2016 *)

taylor(-(x^4+sqrt(x^2-6*x+1)*(x^3-5*x^2+5*x+1)-8*x^3+16*x^2-6*x+1)/(x^3+sqrt(x^2-6*x+1)*(x^2-4*x-1)-7*x^2+7*x-1),x,0,10); /* Vladimir Kruchinin, Nov 25 2016 */

a(n)=if(n==0, 1, sum(k=1,n, binomial(n, k)*binomial(n+k-2, k-1)) ); \\ Joerg Arndt, May 04 2013

A050146 = lambda n : n*hypergeometric([1-n, n], [2], -1) if n>0 else 1
[round(A050146(n).n(100)) for n in (0..24)] # Peter Luschny, Sep 17 2014

From Vladeta Jovovic, Mar 31 2004: (Start)
Coefficient of x^(n-1) in expansion of ((1+x)/(1-x))^n, n > 0.
a(n) = Sum_{k=1..n} binomial(n, k)*binomial(n+k-2, k-1), n > 0. (End)
D-finite with recurrence (n-1)*(n-2)*a(n) = 3*(2*n-3)*(n-2)*a(n-1) - (n-1)*(n-3)*a(n-2) for n > 2. - Vladeta Jovovic, Jul 16 2004
a(n+1) = Jacobi_P(n, 1, -1, 3); a(n+1) = Sum{k=0..n} C(n+1, k)*C(n-1, n-k)*2^k. - Paul Barry, Jan 23 2006
a(n) = n*A006318(n-1) - Abdullahi Umar, Aug 25 2008
a(n) ~ sqrt(3*sqrt(2)-4)*(3+2*sqrt(2))^n/(2*sqrt(Pi*n)). - Vaclav Kotesovec, Oct 08 2012
a(n+1) = A035607(2*n,n). - Reinhard Zumkeller, Jul 20 2013
a(n) = n*hypergeometric([1-n, n], [2], -1) for n >= 1. - Peter Luschny, Sep 17 2014
O.g.f.: -(x^4 + sqrt(x^2 - 6*x + 1)*(x^3 - 5*x^2 + 5*x + 1) - 8*x^3 + 16*x^2 - 6*x + 1)/(x^3 + sqrt(x^2 - 6*x + 1)*(x^2 - 4*x - 1)- 7*x^2 + 7*x - 1). - Vladimir Kruchinin, Nov 25 2016
0 = a(n)*(a(n+1) - 18*a(n+2) + 65*a(n+3) - 12*a(n+4)) + a(n+1)*(54*a(n+2) - 408*a(n+3) + 81*a(n+4)) + a(n+2)*(72*a(n+2) + 334*a(n+3) - 90*a(n+4)) + a(n+3)*(-24*a(n+3) + 9*a(n+4)) for all integer n if a(0) = 0 and a(n) = -2*A050151(-n) for n < 0. - Michael Somos, Nov 25 2016
O.g.f: (2 - x + x*(3 - x)/sqrt(x^2 - 6*x + 1))/2. - Petros Hadjicostas, Feb 14 2021
a(n) = A002002(n) - A026002(n-1) for n>=2. - Shel Kaphan, Mar 24 2023

A123562 Pascal-(1,-3,1) array, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, -1, 1, 1, -3, -3, 1, 1, -5, -3, -5, 1, 1, -7, 1, 1, -7, 1, 1, -9, 9, 11, 9, -9, 1, 1, -11, 21, 17, 17, 21, -11, 1, 1, -13, 37, 11, 1, 11, 37, -13, 1, 1, -15, 57, -15, -39, -39, -15, 57, -15, 1, 1, -17, 81, -69, -87, -81, -87, -69, 81, -17, 1

Offset: 0

Philippe Deléham, Nov 12 2006

Riordan array (1/(1-x), x*(1-3x)/(1-x)).

			Triangle begins:
  1;
  1,   1;
  1,  -1,   1;
  1,  -3,  -3,   1;
  1,  -5,  -3,  -5,   1;
  1,  -7,   1,   1,  -7,   1;
  1,  -9,   9,  11,   9,  -9,  1;
  1, -11,  21,  17,  17,  21, -11,   1;
  1, -13,  37,  11,   1,  11,  37, -13,   1;

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

Cf. Pascal (1,m,1) array: A098593 (m = -2), A000012 (m = -1), A007318 (m = 0), A008288 (m = 1), A081577 (m = 2), A081578 (m = 3), A081579 (m = 4), A081580 (m = 5), A081581 (m = 6), A081582 (m = 7), A143683 (m = 8).

T[n_, k_] := Sum[Binomial[n - j, k]*Binomial[k, j]*(-3)^j, {j, 0, n}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, Oct 15 2017 *)

for(n=0,10, for(k=0,n, print1(sum(j=0,n, binomial(n-j,k)*  binomial(k,j)*(-3)^j), ", "))) \\ G. C. Greubel, Oct 15 2017

Sum_{k=0..n} T(n,k) = A088137(n+1).
T(n,k) = T(n-1,k-1) + T(n-1,k) - 3*T(n-2,k-1), n>0.
From Paul Barry, Jan 24 2011: (Start)
T(n,k) = Sum_{j=0..n} binomial(n-j,k)*binomial(k,j)*(-3)^j.
T(n,k) = [k<=n]*Hypergeometric2F1(-k,k-n,1,-2). (End)
E.g.f. for the n-th subdiagonal: 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 - 3*x - 3*x^2/2! + x^3/3! + 9*x^4/4! + 21*x^5/5! + .... - Peter Bala, Mar 05 2017

A240876 Expansion of (1 + x)^11 / (1 - x)^12.

Original entry on oeis.org

1, 23, 265, 2047, 11969, 56695, 227305, 795455, 2485825, 7059735, 18474633, 45046719, 103274625, 224298231, 464387817, 921406335, 1759885185, 3248227095, 5812626185, 10113604735, 17152640321, 28418229623, 46082942185, 73265596607, 114375683009

Offset: 0

Bruno Berselli, Apr 16 2014

Also 11-dimensional centered hyperoctahedron numbers (see Deza in References) or Crystal ball sequence for 11-dimensional cubic lattice.

This is row/column 11 of the Delannoy numbers array, A008288, which is the main entry for these numbers, listing many more properties. - Peter Munn, Jan 05 2023

E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 230 (paragraph 3.6.6).

Bruno Berselli, Table of n, a(n) for n = 0..1000
D. Bump, K. Choi, P. Kurlberg, and J. Vaaler, A local Riemann hypothesis, I pages 16 and 17.
OEIS Wiki, Centered orthoplex numbers, see Table of formulas and values (row 11).
Index entries for linear recurrences with constant coefficients, signature (12,-66,220,-495,792,-924,792,-495,220,-66,12,-1).

Row/column 11 of A008288.
Cf. similar sequences with g.f. (1+x)^m/(1-x)^(m+1): A005408 (m=1), A001844 .. A001849 (m=2..7), A008417 (m=8), A008419 (m=9), A008421 (m=10), this sequence (m=11), A053805 (m=12).
Subsequence of the odd numbers, A005408.

m:=30; R:=PowerSeriesRing(Integers(),m); Coefficients(R!((1+x)^11/(1-x)^12));

CoefficientList[Series[(1 + x)^11/(1 - x)^12, {x, 0, 30}], x]
LinearRecurrence[{12,-66,220,-495,792,-924,792,-495,220,-66,12,-1},{1,23,265,2047,11969,56695,227305,795455,2485825,7059735,18474633,45046719},30] (* Harvey P. Dale, Apr 15 2018 *)

makelist(coeff(taylor((1+x)^11/(1-x)^12, x, 0, n), x, n), n, 0, 30);

PARI
```
Vec((1+x)^11/(1-x)^12+O(x^30))
```

m = 30; L. = PowerSeriesRing(ZZ, m)
f = (1+x)^11/(1-x)^12
print(f.coefficients())

G.f.: (1 + x)^11 / (1 - x)^12.
a(n) = 12*a(n-1) - 66*a(n-2) + 220*a(n-3) - 495*a(n-4) + 792*a(n-5) - 924*a(n-6) + 792*a(n-7) - 495*a(n-8) + 220*a(n-9) - 66*a(n-10) + 12*a(n-11) - a(n-12), with initial values as shown.
a(n) = (2*n + 1)*(2*n*(n + 1)*(n^2 + n + 5)*(2*n^2 + 2*n + 51)*(n^4 + 2*n^3 + 68*n^2 + 67*n + 537)/155925 + 1).
a(n) = A008421(n) + 2*Sum_{i=0..n-1} A008421(i) for n > 0, a(0) = 1.
Sum_{n >= 0} 1/a(n) = 1.047847848425287358769594801715758965260...
a(n) = Sum_{k = 0..min(11,n)} 2^k*binomial(11,k)*binomial(n,k). See Bump et al. - Tom Copeland, Sep 05 2014

Edited by M. F. Hasler, May 07 2018

A089068 a(n) = a(n-1)+a(n-2)+a(n-3)+2 with a(0)=0, a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 0, 1, 3, 6, 12, 23, 43, 80, 148, 273, 503, 926, 1704, 3135, 5767, 10608, 19512, 35889, 66011, 121414, 223316, 410743, 755475, 1389536, 2555756, 4700769, 8646063, 15902590, 29249424, 53798079, 98950095, 181997600, 334745776, 615693473

Offset: 0

Roger L. Bagula, Dec 03 2003

The a(n+2) represent the Kn12 and Kn22 sums of the square array of Delannoy numbers A008288. See A180662 for the definition of these knight and other chess sums. - Johannes W. Meijer, Sep 21 2010

Index entries for linear recurrences with constant coefficients, signature (2,0,0,-1).

Cf. A000931, A000073, A077939, A113300, A001057, A006054, A033505.

Cf. A000073 (Kn11 & Kn21), A089068 (Kn12 & Kn22), A180668 (Kn13 & Kn23), A180669 (Kn14 & Kn24), A180670 (Kn15 & Kn25). - Johannes W. Meijer, Sep 21 2010

Join[{a=0,b=0,c=1},Table[d=a+b+c+2;a=b;b=c;c=d,{n,50}]] (* Vladimir Joseph Stephan Orlovsky, Apr 19 2011 *)
RecurrenceTable[{a[0]==a[1]==0,a[2]==1,a[n]==a[n-1]+a[n-2]+a[n-3]+2}, a[n],{n,40}] (* or *) LinearRecurrence[{2,0,0,-1},{0,0,1,3},40] (* Harvey P. Dale, Sep 19 2011 *)

a(n) = A008937(n-2)+A008937(n-1). - Johannes W. Meijer, Sep 21 2010
a(n) = A018921(n-5)+A018921(n-4), n>4. - Johannes W. Meijer, Sep 21 2010
a(n) = A000073(n+2)-1. - R. J. Mathar, Sep 22 2010
From Johannes W. Meijer, Sep 22 2010: (Start)
a(n) = a(n-1)+A001590(n+1).
a(n) = Sum_{m=0..n} A040000(m)*A000073(n-m).
a(n+2) = Sum_{k=0..floor(n/2)} A008288(n-k+1,k+1).
G.f. = x^2*(1+x)/((1-x)*(1-x-x^2-x^3)). (End)
a(n) = 2*a(n-1)-a(n-4), a(0)=0, a(1)=0, a(2)=1, a(3)=3. - Bruno Berselli, Sep 23 2010

Corrected and information added by Johannes W. Meijer, Sep 22 2010, Oct 22 2010

Definition based on arbitrarily set floating-point precision removed by R. J. Mathar, Sep 30 2010

A143409 Square array read by antidiagonals: form the Euler-Seidel matrix for the sequence {k!} and then divide column k by k!.

Original entry on oeis.org

1, 2, 1, 5, 3, 1, 16, 11, 4, 1, 65, 49, 19, 5, 1, 326, 261, 106, 29, 6, 1, 1957, 1631, 685, 193, 41, 7, 1, 13700, 11743, 5056, 1457, 316, 55, 8, 1, 109601, 95901, 42079, 12341, 2721, 481, 71, 9, 1, 986410, 876809, 390454, 116125, 25946, 4645, 694, 89, 10, 1

Offset: 0

Peter Bala, Aug 14 2008

The Euler-Seidel matrix for the sequence {k!} is array A076571 read as a square, whose k-th column entries have a common factor of k!. Removing these common factors gives the current table.
This table is closely connected to the constant 1/e. The row, column and diagonal entries of this table occur in series acceleration formulas for 1/e.
For a similar table based on the differences of the sequence {k!} and related to the constant e, see A086764. For other arrays similarly related to constants see A143410 (for sqrt(e)), A143411 (for 1/sqrt(e)), A008288 (for log(2)), A108625 (for zeta(2)) and A143007 (for zeta(3)).

			The Euler-Seidel matrix for the sequence {k!} begins
==============================================
n\k|.....0.....1.....2.....3.....4.....5.....6
==============================================
0..|.....1.....1.....2.....6....24...120...720
1..|.....2.....3.....8....30...144...840
2..|.....5....11....38...174...984
3..|....16....49...212..1158
4..|....65...261..1370
5..|...326..1631
6..|..1957
...
Dividing the k-th column by k! gives
==============================================
n\k|.....0.....1.....2.....3.....4.....5.....6
==============================================
0..|.....1.....1.....1.....1.....1.....1.....1
1..|.....2.....3.....4.....5.....6.....7
2..|.....5....11....19....29....41
3..|....16....49...106...193
4..|....65...261...685
5..|...326..1631
6..|..1957
...
Examples of series formula for 1/e:
Row 2: 1/e = 2*(1/5 - 1/(1!*5*11) + 1/(2!*11*19) - 1/(3!*19*29) + ...).
Column 4: 24/e = 9 - (0!/(1*6) + 1!/(6*41) + 2!/(41*316) + ...).
...
Displayed as a triangle:
0 |     1
1 |     2,     1
2 |     5,     3,    1
3 |    16,    11,    4,    1
4 |    65,    49,   19,    5,   1
5 |   326,   261,  106,   29,   6,  1
6 |  1957,  1631,  685,  193,  41,  7, 1
7 | 13700, 11743, 5056, 1457, 316, 55, 8, 1

Robert Israel, Table of n, a(n) for n = 0..10010 (antidiagonals 0 to 140, flattened)
D. Dumont, Matrices d'Euler-Seidel, Sem. Loth. Comb. B05c (1981) 59-78.
Eric Weisstein's World of Mathematics Poisson-Charlier polynomial

Cf. A008288, A076571, A086764, A108625, A143007, A143410, A143411, A143413, A001517 (main diagonal), A028387 (row 2), A000522 (column 0), A001339 (column 1), A082030 (column 2), A095000 (column 3), A095177 (column 4).

Cf. A273596, A003470.

T := (n, k) -> 1/k!*add(binomial(n,j)*(k+j)!, j = 0..n):
for n from 0 to 9 do seq(T(n, k), k = 0..9) end do;
# Alternate:
T:= proc(n,k) option remember;
  if n = 0 then return 1 fi;
  (n+k)*procname(n-1,k) + procname(n-1,k-1);
end proc:
seq(seq(T(s-n,n),n=0..s),s=0..10); # Robert Israel, Jul 07 2017
# Or:
A143409 := (n,k) -> hypergeom([k+1, k-n], [], -1):
seq(seq(simplify(A143409(n,k)),k=0..n),n=0..9); # Peter Luschny, Oct 05 2017

T[n_, k_] := HypergeometricPFQ[{k+1,k-n}, {}, -1];
Table[T[n,k], {n,0,9}, {k,0,n}] // Flatten (* Peter Luschny, Oct 05 2017 *)

T(n,k) = (1/k!)*Sum_{j = 0..n} binomial(n,j)*(k+j)!.
T(n,k) = ((n+k)!/k!)*Num_Pade(n,k), where Num_Pade(n,k) denotes the numerator of the Padé approximation for the function exp(x) of degree (n,k) evaluated at x = 1.
Recurrence relations:
T(n,k) = T(n-1,k) + (k+1)*T(n-1,k+1);
T(n,k) = (n+k)*T(n-1,k) + T(n-1,k-1).
E.g.f. for column k: exp(y)/(1-y)^(k+1).
E.g.f. for array: exp(y)/(1-x-y) = (1 + x + x^2 + ...) + (2 + 3*x + 4*x^2 + ...)*y + (5 + 11*x + 19*x^2 + ...)*y^2/2! + ... .
Row n lists the values of the Poisson-Charlier polynomial x^(n) + C(n,1)*x^(n-1) + C(n,2)*x^(n-2) + ... + C(n,n) for x = 1,2,3,..., where x^(m) denotes the rising factorial x*(x+1)*...*(x+m-1).
Main diagonal is A001517.
Series formulas for 1/e:
Row n: 1/e = n!*[1/T(n,0) - 1/(1!*T(n,0)*T(n,1)) + 1/(2!*T(n,1)*T(n,2)) - 1/(3!*T(n,2)*T(n,3)) + ...].
Column k: k!/e = A000166(k) + (-1)^(k+1)*[0!/(T(0,k)*T(1,k)) + 1!/(T(1,k)*T(2,k)) + 2!/(T(2,k)*T(3,k)) + ...].
Main diagonal: 1/e = 1 - 2*Sum_{n>=0} (-1)^n/(T(n,n)*T(n+1,n+1)) = 1 - 2*[1/(1*3) - 1/(3*19) + 1/(19*193) - ...].
Second subdiagonal: 1/e = 2*(1^2/(1*5) - 2^2/(5*49) + 3^2/(49*685) - ...).
Compare with A143413.
From Peter Luschny, Oct 05 2017: (Start)
T(n, k) = hypergeom([k+1, k-n], [], -1).
When seen as a triangular array then the row sums are A273596 and the alternating row sums are A003470. (End)

A143683 Pascal-(1,8,1) array.

Original entry on oeis.org

1, 1, 1, 1, 10, 1, 1, 19, 19, 1, 1, 28, 118, 28, 1, 1, 37, 298, 298, 37, 1, 1, 46, 559, 1540, 559, 46, 1, 1, 55, 901, 4483, 4483, 901, 55, 1, 1, 64, 1324, 9856, 21286, 9856, 1324, 64, 1, 1, 73, 1828, 18388, 67006, 67006, 18388, 1828, 73, 1, 1, 82, 2413, 30808, 164242, 304300, 164242, 30808, 2413, 82, 1

Offset: 0

Paul Barry, Aug 28 2008

			Square array begins as:
  1,  1,    1,     1,      1,       1,        1, ... A000012;
  1, 10,   19,    28,     37,      46,       55, ... A017173;
  1, 19,  118,   298,    559,     901,     1324, ...
  1, 28,  298,  1540,   4483,    9856,    18388, ...
  1, 37,  559,  4483,  21286,   67006,   164242, ...
  1, 46,  901,  9856,  67006,  304300,  1004590, ...
  1, 55, 1324, 18388, 164242, 1004590,  4443580, ...
Antidiagonal triangle begins as:
  1;
  1,  1;
  1, 10,   1;
  1, 19,  19,    1;
  1, 28, 118,   28,    1;
  1, 37, 298,  298,   37,   1;
  1, 46, 559, 1540,  559,  46,  1;
  1, 55, 901, 4483, 4483, 901, 55, 1;

Reinhard Zumkeller, Rows n = 0..125 of table, flattened

Cf.Pascal (1,m,1) array: A123562 (m = -3), A098593 (m = -2), A000012 (m = -1), A007318 (m = 0), A008288 (m = 1), A081577 (m = 2), A081578 (m = 3), A081579 (m = 4), A081580 (m = 5), A081581 (m = 6), A081582 (m = 7).

Cf. A003683, A017173, A143680.

a143683 n k = a143683_tabl !! n !! k
a143683_row n = a143683_tabl !! n
a143683_tabl = map fst $ iterate
   (\(us, vs) -> (vs, zipWith (+) (map (* 8) ([0] ++ us ++ [0])) $
                      zipWith (+) ([0] ++ vs) (vs ++ [0]))) ([1], [1, 1])
-- Reinhard Zumkeller, Mar 16 2014

A143683:= func< n,k,q | (&+[Binomial(k, j)*Binomial(n-j, k)*q^j: j in [0..n-k]]) >;
[A143683(n,k,8): k in [0..n], n in [0..12]]; // G. C. Greubel, May 27 2021

Table[Hypergeometric2F1[-k, k-n, 1, 9], {n,0,12}, {k,0,n}]//Flatten (* Jean-François Alcover, May 24 2013 *)

flatten([[hypergeometric([-k, k-n], [1], 9).simplify() for k in (0..n)] for n in (0..12)]) # G. C. Greubel, May 27 2021

Square array: T(n, 0) = T(0, k) = 1, T(n, k) = T(n, k-1) + 8*T(n-1, k-1) + T(n-1, k).
Number triangle: T(n,k) = Sum_{j=0..n-k} binomial(n-k,j)*binomial(k,j)*9^j.
Rows are the expansions of (1+8*x)^k/(1-x)^(k+1).
Riordan array (1/(1-x), x*(1+8*x)/(1-x)).
T(n, k) = Hypergeometric2F1([-k, k-n], [1], 9). - Jean-François Alcover, May 24 2013
E.g.f. for the n-th subdiagonal, n = 0,1,2,..., equals exp(x)*P(n,x), where P(n,x) is the polynomial Sum_{k = 0..n} binomial(n,k)*(9*x)^k/k!. For example, the e.g.f. for the second subdiagonal is exp(x)*(1 + 18*x + 81*x^2/2) = 1 + 19*x + 118*x^2/2! + 298*x^3/3! + 559*x^4/4! + 901*x^5/5! + .... - Peter Bala, Mar 05 2017
Sum_{k=0..n} T(n,k) = A003683(n+1). - G. C. Greubel, May 27 2021

A001849 Crystal ball sequence for 7-dimensional cubic lattice.

Original entry on oeis.org

1, 15, 113, 575, 2241, 7183, 19825, 48639, 108545, 224143, 433905, 795455, 1392065, 2340495, 3800305, 5984767, 9173505, 13726991, 20103025, 28875327, 40754369, 56610575, 77500017, 104692735, 139703809, 184327311, 240673265, 311207743, 398796225, 506750351

Offset: 0

N. J. A. Sloane

This is row/column 7 of the Delannoy numbers array, A008288, which is the main entry for these numbers, listing many more properties. - Shel Kaphan, Jan 06 2023

L. Comtet, Advanced Combinatorics, Reidel, 1974, p. 81.
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).

T. D. Noe, Table of n, a(n) for n = 0..1000
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
R. G. Stanton and D. D. Cowan, Note on a "square" functional equation, SIAM Rev., 12 (1970), 277-279.
Index entries for crystal ball sequences
Index entries for linear recurrences with constant coefficients, signature (8,-28,56,-70,56,-28,8,-1).

Cf. A001848, A001849.
Cf. A240876.
Row/column 7 of A008288.

A001849:=(z+1)**7/(z-1)**8; # conjectured by Simon Plouffe in his 1992 dissertation

CoefficientList[Series[(z + 1)^7/(z - 1)^8, {z, 0, 200}], z] (* Vladimir Joseph Stephan Orlovsky, Jun 19 2011 *)

G.f.: (1+x)^7 /(1-x)^8.
a(n) = (8*n^7 + 28*n^6 + 224*n^5 + 490*n^4 + 1232*n^3 + 1372*n^2 + 1056*n + 315)/315. - Johannes W. Meijer, Jul 14 2013
Sum_{n >= 1} (-1)^(n+1)/(n*a(n-1)*a(n)) =  319/420 - log(2) = (1 - 1/2 + 1/3 - 1/4 + 1/5 - 1/6 + 1/7) - log(2). - Peter Bala, Mar 23 2024

cp's OEIS Frontend

A104698 Triangle read by rows: T(n,k) = Sum_{j=0..n-k} binomial(k, j)*binomial(n-j+1, k+1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A145905 Square array read by antidiagonals: Hilbert transform of triangle A060187.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Formula

A047665 Expansion of (1/sqrt(1-6*x+x^2)-1/(1-x))/2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

Extensions

A050146 a(n) = T(n,n), array T as in A050143.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

Maxima

PARI

Sage

Formula

A123562 Pascal-(1,-3,1) array, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A240876 Expansion of (1 + x)^11 / (1 - x)^12.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Mathematica

Maxima

PARI