A162590 -id:A162590 - cp's OEIS frontend

A136630 Triangular array: T(n,k) counts the partitions of the set [n] into k odd sized blocks.

1, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 4, 0, 1, 0, 1, 0, 10, 0, 1, 0, 0, 16, 0, 20, 0, 1, 0, 1, 0, 91, 0, 35, 0, 1, 0, 0, 64, 0, 336, 0, 56, 0, 1, 0, 1, 0, 820, 0, 966, 0, 84, 0, 1, 0, 0, 256, 0, 5440, 0, 2352, 0, 120, 0, 1, 0, 1, 0, 7381, 0, 24970, 0, 5082, 0, 165, 0, 1, 0, 0, 1024, 0, 87296, 0

Offset: 0

Paul D. Hanna, Jan 14 2008

For partitions into blocks of even size see A156289.
Essentially the unsigned matrix inverse of triangle A121408.
From Peter Bala, Jul 28 2014: (Start)
Define a polynomial sequence x_(n) by setting x_(0) = 1 and for n = 1,2,... setting x_(n) = x*(x + n - 2)*(x + n - 4)*...*(x + n - 2*(n - 1)). Then this table is the triangle of connection constants for expressing the monomial polynomials x^n in terms of the basis x_(k), that is, x^n = sum {k = 0..n} T(n,k)*x_(k) for n = 0,1,2,.... An example is given below.
Let M denote the lower unit triangular array A119467 and for k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0  M/ having the k x k identity matrix I_k as the upper left block; in particular, M(0) = M. Then the present triangle, omitting the first row and column, equals the infinite matrix product M(0)*M(1)*M(2)*.... (End)
Also the Bell transform of A000035(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016

			Triangle begins:
  1;
  0, 1;
  0, 0,   1;
  0, 1,   0,    1;
  0, 0,   4,    0,    1;
  0, 1,   0,   10,    0,     1;
  0, 0,  16,    0,   20,     0,    1;
  0, 1,   0,   91,    0,    35,    0,    1;
  0, 0,  64,    0,  336,     0,   56,    0,   1;
  0, 1,   0,  820,    0,   966,    0,   84,   0,   1;
  0, 0, 256,    0, 5440,     0, 2352,    0, 120,   0, 1;
  0, 1,   0, 7381,    0, 24970,    0, 5082,   0, 165, 0, 1;
T(5,3) = 10. The ten partitions of the set [5] into 3 odd-sized blocks are
(1)(2)(345), (1)(3)(245), (1)(4)(235), (1)(5)(234), (2)(3)(145),
(2)(4)(135), (2)(5)(134), (3)(4)(125), (3)(5)(124), (4)(5)(123).
Connection constants: Row 5 = [0,1,0,10,0,1]. Hence, with the polynomial sequence x_(n) as defined in the Comments section we have x^5 = x_(1) + 10*x_(3) + x_(5) = x + 10*x*(x+1)*(x-1) + x*(x+3)*(x+1)*(x-1)*(x-3).

L. Comtet, Analyse Combinatoire, Presses Univ. de France, 1970, Vol. II, pages 61-62.
L. Comtet, Advanced Combinatorics, Reidel, 1974, pp. 225-226.

Ch. A. Charalambides, Central factorial numbers and related expansions, Fib. Quarterly, Vol. 19, No 5, Dec 1981, pp. 451-456.
Feng Qi and Peter Taylor, Series expansions for powers of sinc function and closed-form expressions for specific partial Bell polynomials, Appl. Anal. Disc. Math. (2024) Vol. 18, No. 1, 1-24. See p. 13.

Cf. A121408; A136631 (antidiagonal sums), A003724 (row sums), A136632; A002452 (column 3), A002453 (column 5); A008958 (central factorial triangle), A156289. A185690, A196776.

A136630 := proc (n, k) option remember; if k < 0 or n < k then 0 elif k = n then 1 else procname(n-2, k-2) + k^2*procname(n-2, k) end if end proc: seq(seq(A136630(n, k), k = 1 .. n), n = 1 .. 12); # Peter Bala, Jul 27 2014
# The function BellMatrix is defined in A264428.
BellMatrix(n -> (n+1) mod 2, 9); # Peter Luschny, Jan 27 2016

t[n_, k_] := Coefficient[ x^k/Product[ 1 - (2*j + k - 2*Quotient[k, 2])^2*x^2, {j, 0, k/2}] + x*O[x]^n, x, n]; Table[t[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Nov 22 2013, after Pari *)
BellMatrix[f_Function, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len-1}, {k, 0, len-1}]];
rows = 13;
M = BellMatrix[Mod[#+1, 2]&, rows];
Table[M[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)

{T(n,k)=polcoeff(x^k/prod(j=0,k\2,1-(2*j+k-2*(k\2))^2*x^2 +x*O(x^n)),n)}

G.f. for column k: x^k/Product_{j=0..floor(k/2)} (1 - (2*j + k-2*floor(k/2))^2 * x^2).
G.f. for column 2*k: x^(2*k)/Product_{j=0..k} (1 - (2*j)^2*x^2).
G.f. for column 2*k+1: x^(2*k+1)/Product_{j=0..k} (1 - (2*j+1)^2*x^2).
From Peter Bala, Feb 21 2011 (Start)
T(n,k) = 1/(2^k*k!)*Sum_{j = 0..k} (-1)^(k-j)*binomial(k,j)*(2*j-k)^n,
Recurrence relation T(n+2,k) = T(n,k-2) + k^2*T(n,k).
E.g.f.: F(x,z) = exp(x*sinh(z)) = Sum_{n>=0} R(n,x)*z^n/n! = 1 + x*z + x^2*z^2/2! + (x+x^3)*z^3/3! + ....
The row polynomials R(n,x) begin
R(1,x) = x
R(2,x) = x^2
R(3,x) = x+x^3.
The e.g.f. F(x,z) satisfies the partial differential equation d^2/dz^2(F) = x^2*F + x*F' + x^2*F'' where ' denotes differentiation w.r.t. x.
Hence the row polynomials satisfy the recurrence relation R(n+2,x) = x^2*R(n,x) + x*R'(n,x) + x^2*R''(n,x) with R(0,x) = 1.
The recurrence relation for T(n,k) given above follows from this.
(End)
For the corresponding triangle of ordered partitions into odd-sized blocks see A196776. Let P denote Pascal's triangle A070318 and put M = 1/2*(P-P^-1). M is A162590 (see also A131047). Then the first column of exp(t*M) lists the row polynomials for the present triangle. - Peter Bala, Oct 06 2011
Row generating polynomials equal D^n(exp(x*t)) evaluated at x = 0, where D is the operator sqrt(1+x^2)*d/dx. Cf. A196776. - Peter Bala, Dec 06 2011
From Peter Bala, Jul 28 2014: (Start)
E.g.f.: exp(t*sinh(x)) = 1 + t*x + t^2*x^2/2! + (t + t^3)*x^3/3! + ....
Hockey-stick recurrence: T(n+1,k+1) = Sum_{i = 0..floor((n-k)/2)} binomial(n,2*i)*T(n-2*i,k).
Recurrence equation for the row polynomials R(n,t):
R(n+1,t) = t*Sum_{k = 0..floor(n/2)} binomial(n,2*k)*R(n-2*k,t) with R(0,t) = 1. (End)

A153641 Nonzero coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent, and Bernoulli numbers (triangle read by rows).

Original entry on oeis.org

1, 1, 1, -1, 1, -3, 1, -6, 5, 1, -10, 25, 1, -15, 75, -61, 1, -21, 175, -427, 1, -28, 350, -1708, 1385, 1, -36, 630, -5124, 12465, 1, -45, 1050, -12810, 62325, -50521, 1, -55, 1650, -28182, 228525, -555731, 1, -66, 2475, -56364, 685575, -3334386, 2702765, 1

Offset: 0

Peter Luschny, Dec 29 2008

In the following the expression [n odd] is 1 if n is odd, 0 otherwise.
(+) W_n(0) = E_n are the Euler (or secant) numbers A122045.
(+) W_n(1) = T_n are the signed tangent numbers, see A009006.
(+) W_{n-1}(1) n / (4^n - 2^n) = B_n gives for n > 1 the Bernoulli number A027641/A027642.
(+) W_n(-1) 2^{-n}(n+1) = G_n the Genocchi number A036968.
(+) W_n(1/2) 2^{n} are the signed generalized Euler (Springer) number, see A001586.
(+) | W_n([n odd]) | the number of alternating permutations A000111.
(+) | W_n([n odd]) / n! | for 0<=n the Euler zeta number A099612/A099617 (see Wikipedia on Bernoulli number). - Peter Luschny, Dec 29 2008
The diagonals in the full triangle (with zero coefficients) of the polynomials have the general form E(k)*binomial(n+k,k) (k>=0 fixed, n=0,1,...) where E(n) are the Euler numbers in the enumeration A122045. For k=2 we find the triangular numbers A000217 and for k=4 A154286. - Peter Luschny, Jan 06 2009
From Peter Bala, Jun 10 2009: (Start)
The Swiss-Knife polynomials W_n(x) may be expressed in terms of the Bernoulli polynomials B(n,x) as
... W_n(x) = 4^(n+1)/(2*n+2)*[B(n+1,(x+3)/4) - B(n+1,(x+1)/4)].
The Swiss-Knife polynomials are, apart from a multiplying factor, examples of generalized Bernoulli polynomials.
Let X be the Dirichlet character modulus 4 defined by X(4*n+1) = 1, X(4*n+3) = -1 and X(2*n) = 0. The generalized Bernoulli polynomials B(X;n,x), n = 1,2,..., associated with the character X are defined by means of the generating function
... t*exp(x*t)*(exp(t)-exp(3*t))/(exp(4*t)-1) = sum {n = 1..inf} B(X;n,x)*t^n/n!.
The first few values are B(X;1,x) = -1/2, B(X;2,x) = -x, B(X,3,x) = -3/2*(x^2-1) and B(X;4,x) = -2*(x^3-3*x).
In general, W_n(x) = -2/(n+1)*B(X;n+1,x).
For the theory of generalized Bernoulli polynomials associated to a periodic arithmetical function see [Cohen, Section 9.4].
The generalized Bernoulli polynomials may be used to evaluate twisted sums of k-th powers. For the present case the result is
sum{n = 0..4*N-1} X(n)*n^k = 1^k - 3^k + 5^k - 7^k + ... - (4*N-1)^k
= [B(X;k+1,4*N) - B(X;k+1,0)]/(k+1) = [W_k(0) - W_k(4*N)]/2.
For the proof apply [Cohen, Corollary 9.4.17 with m = 4 and x = 0].
The generalized Bernoulli polynomials and the Swiss-Knife polynomials are also related to infinite sums of powers through their Fourier series - see the formula section below. For a table of the coefficients of generalized Bernoulli polynomials attached to a Dirichlet character modulus 8 see A151751.
(End)
The Swiss-Knife polynomials provide a general formula for alternating sums of powers similar to the formula which are provided by the Bernoulli polynomials for non-alternating sums of powers (see the Luschny link). Sequences covered by this formula include A001057, A062393, A062392, A011934, A144129, A077221, A137501, A046092. - Peter Luschny, Jul 12 2009
The greatest common divisor of the nonzero coefficients of the decapitated Swiss-Knife polynomials is exp(Lambda(n)), where Lambda(n) is the von Mangoldt function for odd primes, symbolically:
gcd(coeffs(SKP_{n}(x) - x^n)) = A155457(n) (n>1). - Peter Luschny, Dec 16 2009
Another version is at A119879. - Philippe Deléham, Oct 26 2013

			1
x
x^2  -1
x^3  -3x
x^4  -6x^2   +5
x^5 -10x^3  +25x
x^6 -15x^4  +75x^2  -61
x^7 -21x^5 +175x^3 -427x

H. Cohen, Number Theory - Volume II: Analytic and Modern Tools, Graduate Texts in Mathematics. Springer-Verlag. [From Peter Bala, Jun 10 2009]

G. C. Greubel, Table of n, a(n) for the first 76 rows, flattened
Ayse Yilmaz Ceylan and Yilmaz Simsek, Formulae for Generalization of Touchard Polynomials with Their Generating Functions, Symmetry (2025) Vol. 17, Issue 7, Art. No. 1126. See Eq. 28 and after.
Kwang-Wu Chen, Algorithms for Bernoulli numbers and Euler numbers, J. Integer Sequences, 4 (2001), #01.1.6.
Suyoung Choi and Hanchul Park, A new graph invariant arises in toric topology, arXiv preprint arXiv:1210.3776 [math.AT], 2012.
Leonhard Euler (1735), De summis serierum reciprocarum, Opera Omnia I.14, E 41, 73-86; On the sums of series of reciprocals, arXiv:math/0506415 [math.HO], 2005-2008.
A. Hodges and C. V. Sukumar, Bernoulli, Euler, permutations and quantum algebras, Proc. R. Soc. A Oct. 2007 vol 463 no. 463 2086 2401-2414. [Added by Tom Copeland, Aug 31 2015]
Peter Luschny, The Swiss-Knife polynomials.
Peter Luschny, Swiss-Knife polynomials and Euler numbers
Wikipedia, Bernoulli number
J. Worpitzky, Studien über die Bernoullischen und Eulerschen Zahlen, Journal für die reine und angewandte Mathematik, 94 (1883), 203-232.

Cf. A151751, A162590, A000111, A001586, A009006, A027641/A027642, A036968, A099612/A099617, A119879, A104033.
W_n(k), k=0,1,...
W_0:  1,  1,  1,  1,   1,   1, ........ A000012
W_1:  0,  1,  2,  3,   4,   5, ........ A001477
W_2: -1,  0,  3,  8,  15,  24, ........ A067998
W_3:  0, -2,  2, 18,  52, 110, ........ A121670
W_4:  5,  0, -3, 32, 165, 480, ........
W_n(k), n=0,1,...
k=0:  1,  0, -1,  0,   5,   0, -61, ... A122045
k=1:  1,  1,  0, -2,   0,  16,   0, ... A155585
k=2:  1,  2,  3,  2,  -3,   2,  63, ... A119880
k=3:  1,  3,  8, 18,  32,  48, 128, ... A119881
k=4:  1,  4, 15, 52, 165, 484, ........         [Peter Luschny, Jul 07 2009]

w := proc(n,x) local v,k,pow,chen; pow := (a,b) -> if a = 0 and b = 0 then 1 else a^b fi; chen := proc(m) if irem(m+1,4) = 0 then RETURN(0) fi; 1/((-1)^iquo(m+1,4) *2^iquo(m,2)) end; add(add((-1)^v*binomial(k,v)*pow(v+x+1,n)*chen(k),v=0..k), k=0..n) end:
# Coefficients with zeros:
seq(print(seq(coeff(i!*coeff(series(exp(x*t)*sech(t),t,16),t,i),x,i-n),n=0..i)), i=0..8);
# Recursion
W := proc(n,z) option remember; local k,p;
if n = 0 then 1 else p := irem(n+1,2);
z^n - p + add(`if`(irem(k,2)=1,0,
W(k,0)*binomial(n,k)*(power(z,n-k)-p)),k=2..n-1) fi end:
# Peter Luschny, edited and additions Jul 07 2009, May 13 2010, Oct 24 2011

max = 9; rows = (Reverse[ CoefficientList[ #, x]] & ) /@ CoefficientList[ Series[ Exp[x*t]*Sech[t], {t, 0, max}], t]*Range[0, max]!; par[coefs_] := (p = Partition[ coefs, 2][[All, 1]]; If[ EvenQ[ Length[ coefs]], p, Append[ p, Last[ coefs]]]); Flatten[ par /@ rows] (* Jean-François Alcover, Oct 03 2011, after g.f. *)
sk[n_, x_] := Sum[Binomial[n, k]*EulerE[k]*x^(n-k), {k, 0, n}]; Table[CoefficientList[sk[n, x], x] // Reverse // Select[#, # =!= 0 &] &, {n, 0, 13}] // Flatten (* Jean-François Alcover, May 21 2013 *)
Flatten@Table[Binomial[n, 2k] EulerE[2k], {n, 0, 12}, {k, 0, n/2}](* Oliver Seipel, Jan 14 2025 *)

def A046978(k):
    if k % 4 == 0:
        return 0
    return (-1)**(k // 4)
def A153641_poly(n, x):
    return expand(add(2**(-(k // 2))*A046978(k+1)*add((-1)**v*binomial(k,v)*(v+x+1)**n for v in (0..k)) for k in (0..n)))
for n in (0..7): print(A153641_poly(n, x))  # Peter Luschny, Oct 24 2011

W_n(x) = Sum_{k=0..n}{v=0..k} (-1)^v binomial(k,v)*c_k*(x+v+1)^n where c_k = frac((-1)^(floor(k/4))/2^(floor(k/2))) [4 not div k] (Iverson notation).
From Peter Bala, Jun 10 2009: (Start)
E.g.f.: 2*exp(x*t)*(exp(t)-exp(3*t))/(1-exp(4*t))= 1 + x*t + (x^2-1)*t^2/2! + (x^3-3*x)*t^3/3! + ....
W_n(x) = 1/(2*n+2)*Sum_{k=0..n+1} 1/(k+1)*Sum_{i=0..k} (-1)^i*binomial(k,i)*((x+4*i+3)^(n+1) - (x+4*i+1)^(n+1)).
Fourier series expansion for the generalized Bernoulli polynomials:
B(X;2*n,x) = (-1)^n*(2/Pi)^(2*n)*(2*n)! * {sin(Pi*x/2)/1^(2*n) - sin(3*Pi*x/2)/3^(2*n) + sin(5*Pi*x/2)/5^(2*n) - ...}, valid for 0 <= x <= 1 when n >= 1.
B(X;2*n+1,x) = (-1)^(n+1)*(2/Pi)^(2*n+1)*(2*n+1)! * {cos(Pi*x/2)/1^(2*n+1) - cos(3*Pi*x/2)/3^(2*n+1) + cos(5*Pi*x/2)/5^(2*n+1) - ...}, valid for 0 <= x <= 1 when n >= 1 and for 0 <= x < 1 when n = 0.
(End)
E.g.f.: exp(x*t) * sech(t). - Peter Luschny, Jul 07 2009
O.g.f. as a J-fraction: z/(1-x*z+z^2/(1-x*z+4*z^2/(1-x*z+9*z^2/(1-x*z+...)))) = z + x*z^2 + (x^2-1)*z^3 + (x^3-3*x)*z^4 + .... - Peter Bala, Mar 11 2012
Conjectural o.g.f.: Sum_{n >= 0} (1/2^((n-1)/2))*cos((n+1)*Pi/4)*( Sum_{k = 0..n} (-1)^k*binomial(n,k)/(1 - (k + x)*t) ) = 1 + x*t + (x^2 - 1)*t^2 + (x^3 - 3*x)*t^3 + ... (checked up to O(t^13)), which leads to W_n(x) = Sum_{k = 0..n} 1/2^((k - 1)/2)*cos((k + 1)*Pi/4)*( Sum_{j = 0..k} (-1)^j*binomial(k, j)*(j + x)^n ). - Peter Bala, Oct 03 2016

A119467 A masked Pascal triangle.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 0, 3, 0, 1, 1, 0, 6, 0, 1, 0, 5, 0, 10, 0, 1, 1, 0, 15, 0, 15, 0, 1, 0, 7, 0, 35, 0, 21, 0, 1, 1, 0, 28, 0, 70, 0, 28, 0, 1, 0, 9, 0, 84, 0, 126, 0, 36, 0, 1, 1, 0, 45, 0, 210, 0, 210, 0, 45, 0, 1, 0, 11, 0, 165, 0, 462, 0, 330, 0, 55, 0, 1, 1, 0, 66, 0, 495, 0, 924

Offset: 0

Paul Barry, May 21 2006

Row sums are A011782. Diagonal sums are F(n+1)*(1+(-1)^n)/2 (aerated version of A001519). Product by Pascal's triangle A007318 is A119468. Schur product of (1/(1-x),x/(1-x)) and (1/(1-x^2),x).
Exponential Riordan array (cosh(x),x). Inverse is (sech(x),x) or A119879. - Paul Barry, May 26 2006
Rows give coefficients of polynomials p_n(x) = Sum_{k=0..n} (k+1 mod 2)*binomial(n,k)*x^(n-k) having e.g.f. exp(x*t)*cosh(t)= 1*(t^0/0!) + x*(t^1/1!) + (1+x^2)*(t^2/2!) + ... - Peter Luschny, Jul 14 2009
Inverse of the coefficient matrix of the Swiss-Knife polynomials in ascending order of x^i (reversed and aerated rows of A153641). - Peter Luschny, Jul 16 2012
Call this array M and for k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0  M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. The infinite matrix product M(0)*M(1)*M(2)*... is equal to A136630 but with the first row and column omitted. - Peter Bala, Jul 28 2014
The row polynomials SKv(n,x) = [(x+1)^n + (x-1)^n]/2 , with e.g.f. cosh(t)*exp(xt), are the umbral compositional inverses of the row polynomials of A119879 (basically the Swiss Knife polynomials SK(n,x) of A153641); i.e., umbrally SKv(n,SK(.,x)) = x^n = SK(n,SKv(.,x)). Therefore, this entry's matrix and A119879 are an inverse pair. Both sequences of polynomials are Appell sequences, i.e., d/dx P(n,x) = n * P(n-1,x) and (P(.,x)+y)^n = P(n,x+y). In particular, (SKv(.,0)+x)^n = SKv(n,x), reflecting that the first column has the e.g.f. cosh(t). The raising operator is R = x + tanh(d/dx); i.e., R SKv(n,x) = SKv(n+1,x). The coefficients of this operator are basically the signed and aerated zag numbers A000182, which can be expressed as normalized Bernoulli numbers. The triangle is formed by multiplying the n-th diagonal of the lower triangular Pascal matrix by the Taylor series coefficient a(n) of cosh(x). More relations for this type of triangle and its inverse are given by the formalism of A133314. - Tom Copeland, Sep 05 2015
The signed version of this matrix has the e.g.f. cos(t) e^{xt}, generating Appell polynomials that have only real, simple zeros and whose extrema are maxima above the x-axis and minima below and situated above and below the zeros of the next lower degree polynomial. The bivariate versions appear on p. 27 of Dimitrov and Rusev in conditions for entire functions that are cosine transforms of a class of functions to have only real zeros. - Tom Copeland, May 21 2020
The n-th row of the triangle is obtained by multiplying by 2^(n-1) the elements of the first row of the limit as k approaches infinity of the stochastic matrix P^(2k-1) where P is the stochastic matrix associated with the Ehrenfest model with n balls. The elements of a stochastic matrix P give the probabilities of arriving in a state j given the previous state i. In particular the sum of every row of the matrix must be 1, and so the sum of the terms of the n-th row of this triangle is 2^(n-1). Furthermore, by the properties of Markov chains, we can interpret P^(2k-1) as the (2k-1)-step transition matrix of the Ehrenfest model and its limit exists and it is again a stochastic matrix. The rows of the triangle divided by 2^(n-1) are the even rows (second, fourth, ...) and the odd rows (first, third, ...) of the limit matrix P^(2k-1). - Luca Onnis, Oct 29 2023

			Triangle begins
  1,
  0, 1,
  1, 0,  1,
  0, 3,  0,  1,
  1, 0,  6,  0,   1,
  0, 5,  0, 10,   0,   1,
  1, 0, 15,  0,  15,   0,   1,
  0, 7,  0, 35,   0,  21,   0,  1,
  1, 0, 28,  0,  70,   0,  28,  0,  1,
  0, 9,  0, 84,   0, 126,   0, 36,  0, 1,
  1, 0, 45,  0, 210,   0, 210,  0, 45, 0, 1
p[0](x) = 1
p[1](x) = x
p[2](x) = 1 + x^2
p[3](x) = 3*x + x^3
p[4](x) = 1 + 6*x^2 + x^4
p[5](x) = 5*x + 10*x^3 + x^5
Connection with A136630: With the arrays M(k) as defined in the Comments section, the infinite product M(0)*M(1)*M(2)*... begins
/1        \/1        \/1        \      /1         \
|0 1      ||0 1      ||0 1      |      |0 1       |
|1 0 1    ||0 0 1    ||0 0 1    |... = |1 0  1    |
|0 3 0 1  ||0 1 0 1  ||0 0 0 1  |      |0 4  0 1  |
|1 0 6 0 1||0 0 3 0 1||0 0 1 0 1|      |1 0 10 0 1|
|...      ||...      ||...      |      |...       |
- _Peter Bala_, Jul 28 2014

Paul and Tatjana Ehrenfest, Über zwei bekannte Einwände gegen das Boltzmannsche H-Theorem, Physikalische Zeitschrift, vol. 8 (1907), pp. 311-314.

Reinhard Zumkeller, Rows n = 0..125 of table, flattened
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 28.
Paul Barry, On the inversion of Riordan arrays, arXiv:2101.06713 [math.CO], 2021.
Tom Copeland, Skipping over Dimensions, Juggling Zeros in the Matrix, 2020.
D. Dimitrov and P. Rusev, Zeros of entire Fourier transforms, East Journal on Approximations, Vol. 17, No. 1, p. 1-108, 2011.
Miguel Méndez and Rafael Sánchez, On the combinatorics of Riordan arrays and Sheffer polynomials: monoids, operads and monops, arXiv:1707.00336 [math.CO], 2017, Section 4.3, Example 4.
Miguel A. Méndez and Rafael Sánchez Lamoneda, Monops, Monoids and Operads: The Combinatorics of Sheffer Polynomials, The Electronic Journal of Combinatorics 25(3) (2018), #P3.25.
Luca Onnis, Animation of the Ehrenfest model.
Wikipedia, Ehrenfest model.
Index entries for triangles and arrays related to Pascal's triangle

Cf. A131047, A153641, A162590.
From Peter Luschny, Jul 14 2009: (Start)
Cf. A034839, A162590.
p[n](k), n=0,1,...
k= 0: 1,  0,   1,    0,    1,   0, ... A128174
k= 1: 1,  1,   2,    4,    8,  16, ... A011782
k= 2: 1,  2,   5,   14,   41, 122, ... A007051
k= 3: 1,  3,  10,   36,  136,      ... A007582
k= 4: 1,  4,  17,   76,  353,      ... A081186
k= 5: 1,  5,  26,  140,  776,      ... A081187
k= 6: 1,  6,  37,  234, 1513,      ... A081188
k= 7: 1,  7,  50,  364, 2696,      ... A081189
k= 8: 1,  8,  65,  536, 4481,      ... A081190
k= 9: 1,  9,  82,  756, 7048,      ... A060531
k=10: 1, 10, 101, 1030,            ... A081192
p[n](k), k=0,1,...
p[0]: 1,1,1,1,1,1, ....... A000012
p[1]: 0,1,2,3,4,5, ....... A001477
p[2]: 1,2,5,10,17,26, .... A002522
p[3]: 0,4,14,36,76,140, .. A079908 (End)
Cf. A000182, A133314, A153641.

a119467 n k = a119467_tabl !! n !! k
a119467_row n = a119467_tabl !! n
a119467_tabl = map (map (flip div 2)) $
               zipWith (zipWith (+)) a007318_tabl a130595_tabl
-- Reinhard Zumkeller, Mar 23 2014

/* As triangle */ [[Binomial(n, k)*(1 + (-1)^(n - k))/2: k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 26 2015

# Polynomials: p_n(x)
p := proc(n,x) local k, pow; pow := (n,k) -> `if`(n=0 and k=0,1,n^k);
add((k+1 mod 2)*binomial(n,k)*pow(x,n-k),k=0..n) end;
# Coefficients: a(n)
seq(print(seq(coeff(i!*coeff(series(exp(x*t)*cosh(t),t,16),t,i),x,n),n=0..i)),i=0..8); # Peter Luschny, Jul 14 2009

Table[Binomial[n, k] (1 + (-1)^(n - k))/2, {n, 0, 12}, {k, 0, n}] // Flatten (* Michael De Vlieger, Sep 06 2015 *)
n = 15; "n-th row"
mat = Table[Table[0, {j, 1, n + 1}], {i, 1, n + 1}];
mat[[1, 2]] = 1;
mat[[n + 1, n]] = 1;
For[i = 2, i <= n, i++, mat[[i, i - 1]] = (i - 1)/n ];
For[i = 2, i <= n, i++, mat[[i, i + 1]] = (n - i + 1)/n];
mat // MatrixForm;
P2 = Dot[mat, mat];
R1 = Simplify[
  Eigenvectors[Transpose[P2]][[1]]/
   Total[Eigenvectors[Transpose[P2]][[1]]]]
R2 = Table[Dot[R1, Transpose[mat][[k]]], {k, 1, n + 1}]
odd = R2*2^(n - 1) (* _Luca Onnis *)

@CachedFunction
def A119467_poly(n):
    R = PolynomialRing(ZZ, 'x')
    x = R.gen()
    return R.one() if n==0 else R.sum(binomial(n,k)*x^(n-k) for k in range(0,n+1,2))
def A119467_row(n):
    return list(A119467_poly(n))
for n in (0..10) : print(A119467_row(n)) # Peter Luschny, Jul 16 2012

G.f.: (1-x*y)/(1-2*x*y-x^2+x^2*y^2);
T(n,k) = C(n,k)*(1+(-1)^(n-k))/2;
Column k has g.f. (1/(1-x^2))*(x/(1-x^2))^k*Sum_{j=0..k+1} binomial(k+1,j)*sin((j+1)*Pi/2)^2*x^j.
Column k has e.g.f. cosh(x)*x^k/k!. - Paul Barry, May 26 2006
Let Pascal's triangle, A007318 = P; then this triangle = (1/2) * (P + 1/P). Also A131047 = (1/2) * (P - 1/P). - Gary W. Adamson, Jun 12 2007
Equals A007318 - A131047 since the zeros of the triangle are masks for the terms of A131047. Thus A119467 + A131047 = Pascal's triangle. - Gary W. Adamson, Jun 12 2007
T(n,k) = (A007318(n,k) + A130595(n,k))/2, 0<=k<=n. - Reinhard Zumkeller, Mar 23 2014

Edited by N. J. A. Sloane, Jul 14 2009

A196776 Triangle T(n,k) gives the number of ordered partitions of an n set into k odd-sized blocks.

Original entry on oeis.org

1, 0, 2, 1, 0, 6, 0, 8, 0, 24, 1, 0, 60, 0, 120, 0, 32, 0, 480, 0, 720, 1, 0, 546, 0, 4200, 0, 5040, 0, 128, 0, 8064, 0, 40320, 0, 40320, 1, 0, 4920, 0, 115920, 0, 423360, 0, 362880, 0, 512, 0, 130560, 0, 1693440, 0, 4838400, 0, 3628800

Offset: 1

Peter Bala, Oct 06 2011

See A136630 for the case of unordered partitions into odd-sized blocks. See A193474 for this triangle in row reverse form (but with an offset of 0).

			Triangle begins
.n\k.|..1....2....3....4.....5....6.....7
= = = = = = = = = = = = = = = = = = = = =
..1..|..1
..2..|..0....2
..3..|..1....0....6
..4..|..0....8....0...24
..5..|..1....0...60....0...120
..6..|..0...32....0..480.....0..720
..7..|..1....0..546....0..4200....0..5040
...
T(4,2) = 8: The 8 ordered partitions of the set {1,2,3,4} into 2 odd-sized blocks are {1}{2,3,4}, {2,3,4}{1}, {2}{1,3,4}, {1,3,4}{2}, {3}{1,2,4}, {1,2,4}{3}, {4}{1,2,3} and {1,2,3}{4}.
Example of recurrence relation: T(7,3) = 3^2*T(5,3) + 3*(3-1)*T(5,1) = 9*60 + 6*1 = 546.

Cf. A006154 (row sums), A136630, A162590, A193474 (row reverse).

T(n,k) = 1/(2^k)*sum {j = 0..k}(-1)^(k-j)*binomial(k,j)*(2*j-k)^n.
Recurrence: T(n+2,k) = k^2*T(n,k) + k*(k-1)*T(n,k-2).
E.g.f.: x*sinh(t)/(1-x*sinh(t)) = x*t + 2*x^2*t^2/2! + (x+6*x^3)*t^3/3! + (8*x^2+24*x^4)*t^4/4! + (x+60*x^3+120*x^5)*t^5/5! + ....
O.g.f. for column 2*k: (2*k)!*x^(2*k)/Product {j = 0..k} (1 - (2*j)^2*x^2).
O.g.f. for column 2*k+1: (2*k+1)!*x^(2*k+1)/Product {j = 0..k} (1 - (2*j+1)^2*x^2).
Let P denote Pascal's triangle A070318 and put M = 1/2*(P-P^-1). M is A162590 (see also A131047). Then the first column of (I-t*M)^-1 (apart from the initial 1) lists the row polynomials for the present triangle.
n-th row sum = A006154(n).
Row generating polynomials equal D^n(1/(1-x*t)) evaluated at x = 0, where D is the operator sqrt(1+x^2)*d/dx. Cf. A136630. - Peter Bala, Dec 06 2011

A135685 Triangular sequence of the coefficients of the numerator of the rational recursive sequence for tan(n*y) with x = tan(y).

Original entry on oeis.org

0, 0, 1, 0, -2, 0, -3, 0, 1, 0, 4, 0, -4, 0, 5, 0, -10, 0, 1, 0, -6, 0, 20, 0, -6, 0, -7, 0, 35, 0, -21, 0, 1, 0, 8, 0, -56, 0, 56, 0, -8, 0, 9, 0, -84, 0, 126, 0, -36, 0, 1, 0, -10, 0, 120, 0, -252, 0, 120, 0, -10, 0, -11, 0, 165, 0, -462, 0, 330, 0, -55, 0, 1

Offset: 0

Roger L. Bagula, Feb 17 2008

Signed version of A034867 with interlaced zeros. - Joerg Arndt, Sep 14 2014

The negatives of these terms gives the coefficients for the numerators for when n is negative (i.e. tan(-n*y) = -tan(n*y)). - James Burling, Sep 14 2014

			Triangle starts:
  0;
  0,   1;
  0,  -2;
  0,  -3,  0,   1;
  0,   4,  0,  -4;
  0,   5,  0, -10,  0,    1;
  0,  -6,  0,  20,  0,   -6;
  0,  -7,  0,  35,  0,  -21,  0,   1;
  0,   8,  0, -56,  0,   56,  0,  -8;
  0,   9,  0, -84,  0,  126,  0, -36,  0,   1;
  0, -10,  0, 120,  0, -252,  0, 120,  0, -10;
  0, -11,  0, 165,  0, -462,  0, 330,  0, -55,  0,  1;

Robert Israel, Table of n, a(n) for n = 0..10082
Clark Kimberling, Polynomials associated with reciprocation, JIS 12 (2009) 09.3.4, section 5.

Cf. A095704, A162590.

g[0]:= 0:
g[1]:= x;
for n from 2 to 20 do
g[n]:= expand(-2*(-1)^n*g[n-1]+(x^2+1)*g[n-2])
od:
0, seq(seq(coeff(g[n],x,j),j=0..degree(g[n])),n=1..20); # Robert Israel, Sep 14 2014

p[n_, x_]:= p[n, x]= If[n<2, n*x, (p[n-1, x] + x)/(1 - x*p[n-1, x])];
Table[CoefficientList[Numerator[FullSimplify[p[n, x]]], x], {n,0,12}]//Flatten

def p(n, x): return n*x if (n<2) else 2*(-1)^(n+1)*p(n-1,x) + (1+x^2)*p(n-2,x)
def A135685(n,k): return ( p(n,x) ).series(x,n+1).list()[k]
flatten([[A135685(n,k) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Nov 26 2021

p(n, x) = (p(n-1, x) + x)/(1 - x*p(n-1, x)), with p(0, x) = 0, p(1, x) = x.

Sum_{j} T(n,j)*x^j = g(n,x) where g(0,x) = 0, g(1,x) = x, g(n,x) = -2*(-1)^n*g(n-1,x) + (x^2+1)*g(n-2,x). - Robert Israel, Sep 14 2014

Prepended first term and offset corrected by James Burling, Sep 14 2014

cp's OEIS Frontend

A136630 Triangular array: T(n,k) counts the partitions of the set [n] into k odd sized blocks.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A153641 Nonzero coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent, and Bernoulli numbers (triangle read by rows).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

A119467 A masked Pascal triangle.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

Sage

Formula

Extensions

A196776 Triangle T(n,k) gives the number of ordered partitions of an n set into k odd-sized blocks.

Views

Author

Keywords

Comments

Examples

Crossrefs

Formula

A135685 Triangular sequence of the coefficients of the numerator of the rational recursive sequence for tan(n*y) with x = tan(y).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Sage

Formula

Extensions