cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-9 of 9 results.

A131053 A000012 * A131047.

Original entry on oeis.org

1, 1, 2, 2, 2, 3, 2, 6, 3, 4, 3, 6, 13, 4, 5, 3, 12, 13, 24, 5, 6, 4, 12, 34, 24, 40, 6, 7, 4, 20, 34, 80, 40, 62, 7, 8, 5, 20, 70, 80, 166, 62, 91, 8, 9
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Row sums = (2^n - 1): (1, 3, 7, 15, ...).
Left border = A004526: (1, 1, 2, 2, 3, 3, ...).

Examples

			First few rows of the triangle:
  1;
  1,  2;
  2,  2,  3;
  2,  6,  3,  4;
  3,  6, 13,  4,  5;
  3, 12, 13, 24,  5,  6;
  4, 12, 34, 24, 40,  6,  7;
  ...
		

Crossrefs

Formula

A000012 * A131047 as infinite lower triangular matrices.

A131052 A131047 * A000012.

Original entry on oeis.org

1, 2, 2, 4, 3, 3, 8, 8, 4, 4, 16, 15, 15, 5, 5, 32, 32, 26, 26, 6, 6, 64, 63, 63, 42, 42, 7, 7, 128, 128, 120, 120, 64, 64, 8, 8
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Row sums = A087447: (1, 4, 10, 24, 56, ...). A131053 = A000012 * A131047.

Examples

			First few rows of the triangle:
   1;
   2,  2;
   4,  3,  3;
   8,  8,  4,  4;
  16, 15, 15,  5,  5;
  32, 32, 26, 26,  6,  6;
  64, 63, 63, 42, 42,  7,  7;
  ...
		

Crossrefs

Formula

A131047 * A000012 as infinite lower triangular matrices.

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

Original entry on oeis.org

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

Views

Author

Paul D. Hanna, Jan 14 2008

Keywords

Comments

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

Examples

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

References

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

Crossrefs

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

Programs

  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    {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)}

Formula

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)

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

Views

Author

Paul Barry, May 21 2006

Keywords

Comments

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

Examples

			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
		

References

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

Crossrefs

From Peter Luschny, Jul 14 2009: (Start)
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)

Programs

  • Haskell
    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
    
  • Magma
    /* As triangle */ [[Binomial(n, k)*(1 + (-1)^(n - k))/2: k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 26 2015
  • Maple
    # 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
  • Mathematica
    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 *)
  • Sage
    @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
    

Formula

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

Extensions

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

Views

Author

Peter Bala, Oct 06 2011

Keywords

Comments

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

Examples

			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.
		

Crossrefs

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

Formula

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

A162590 Polynomials with e.g.f. exp(x*t)/csch(t), triangle of coefficients read by rows.

Original entry on oeis.org

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

Views

Author

Peter Luschny, Jul 07 2009

Keywords

Comments

Comment from Peter Bala (Dec 06 2011): "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" A196776(n,k), which gives the number of ordered partitions of an n set into k odd-sized blocks. - Peter Luschny, Dec 06 2011
The n-th row of the triangle is formed 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 probability 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 in the n-th row of this triangle is 2^(n-1). Furthermore, by the properties of Markov chains, we can interpret P^(2k) as the (2k)-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). - Luca Onnis, Oct 29 2023

Examples

			Triangle begins:
  0
  1,  0
  0,  2,  0
  1,  0,  3,  0
  0,  4,  0,  4,  0
  1,  0, 10,  0,  5,  0
  0,  6,  0, 20,  0,  6,  0
  1,  0, 21,  0, 35,  0,  7,  0
  ...
  p[0](x) = 0;
  p[1](x) = 1
  p[2](x) = 2*x
  p[3](x) = 3*x^2 +  1
  p[4](x) = 4*x^3 +  4*x
  p[5](x) = 5*x^4 + 10*x^2 +  1
  p[6](x) = 6*x^5 + 20*x^3 +  6*x
  p[7](x) = 7*x^6 + 35*x^4 + 21*x^2 + 1
  p[8](x) = 8*x^7 + 56*x^5 + 56*x^3 + 8*x
.
Cf. the triangle of odd-numbered terms in rows of Pascal's triangle (A034867).
p[n] (k), n=0,1,...
k=0:  0, 1,  0,   1,    0,     1, ... A000035, (A059841)
k=1:  0, 1,  2,   4,    8,    16, ... A131577, (A000079)
k=2:  0, 1,  4,  13,   40,   121, ... A003462
k=3:  0, 1,  6,  28,  120,   496, ... A006516
k=4:  0, 1,  8,  49,  272,  1441, ... A005059
k=5:  0, 1, 10,  76,  520,  3376, ... A081199, (A016149)
k=6:  0, 1, 12, 109,  888,  6841, ... A081200, (A016161)
k=7:  0, 1, 14, 148, 1400, 12496, ... A081201, (A016170)
k=8:  0, 1, 16, 193, 2080, 21121, ... A081202, (A016178)
k=9:  0, 1, 18, 244, 2952, 33616, ... A081203, (A016186)
k=10: 0, 1, 20, 301, 4040, 51001, ... ......., (A016190)
.
p[n] (k), k=0,1,...
p[0]: 0,  0,   0,    0,    0,     0, ... A000004
p[1]: 1,  1,   1,    1,    1,     1, ... A000012
p[2]: 0,  2,   4,    6,    8,    10, ... A005843
p[3]: 1,  4,  13,   28,   49,    76, ... A056107
p[4]: 0,  8,  40,  120,  272,   520, ... A105374
p[5]: 1, 16, 121,  496, 1441,  3376, ...
p[6]: 0, 32, 364, 2016, 7448, 21280, ...
		

Crossrefs

Cf. A119467.

Programs

  • Maple
    # Polynomials: p_n(x)
    p := proc(n,x) local k;
    pow := (n,k) -> `if`(n=0 and k=0,1,n^k);
    add((k 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)/csch(t), t,16),t,i),x,n), n=0..i)), i=0..8);
  • Mathematica
    p[n_, x_] := Sum[Binomial[n, 2*k-1]*x^(n-2*k+1), {k, 0, n+2}]; row[n_] := CoefficientList[p[n, x], x] // Append[#, 0]&; Table[row[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Jun 28 2013 *)
    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}]
    even = R1*2^(n - 1) (* Luca Onnis, Oct 29 2023 *)

Formula

p_n(x) = Sum_{k=0..n} (k mod 2)*binomial(n,k)*x^(n-k).
E.g.f.: exp(x*t)/csch(t) = 0*(t^0/0!) + 1*(t^1/1!) + (2*x)*(t^2/2!) + (3*x^2+1)*(t^3/3!) + ...
The 'co'-polynomials with generating function exp(x*t)*sech(t) are the Swiss-Knife polynomials (A153641).

A131048 (1/3) * (A007318^2 - A007318^(-1)).

Original entry on oeis.org

1, 1, 2, 3, 3, 3, 5, 12, 6, 4, 11, 25, 30, 10, 5, 21, 66, 75, 60, 15, 6, 43, 147, 231, 175, 105, 21, 7, 85, 344, 588, 616, 350, 168, 28, 8, 171, 765, 1548, 1764, 1386, 630, 252, 36, 9
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Left border = A001045: (1, 1, 3, 5, 11, 21, 43, 85, ...).
Row sums = (1, 3, 9, 27, ...).
Analogous triangles for other powers of P are: A131047, A131049, A131050 and A131051.

Examples

			First few rows of the triangle:
   1;
   1,   2;
   3,   3,   3;
   5,  12,   6,   4;
  11,  25,  30,  10,   5;
  21,  66,  75,  60,  15,  6;
  43, 147, 231, 175, 105, 21, 7;
  ...
		

Crossrefs

Formula

Let A007318 (Pascal's triangle) = P. then A131048 = (1/3) * (P^2 - 1/P). Delete right border of zeros.
From Peter Bala, Oct 24 2007: (Start)
O.g.f.: 1/(1 - (2*x + 1)*t + (x^2 + x - 2)*t^2) = 1 + (1 + 2*x)*t + (3 + 3*x + 3*x^2)*t^2 + ....
T(n,n-k) = (1/3)*C(n,k)*(2^k - (-1)^k) = C(n,k)*A001045(k).
The row polynomials R(n,x) := Sum_{k = 0..n} T(n,n-k)*x^(n-k) = (1/3)*((x + 2)^n - (x - 1)^n) and have the divisibility property R(n,x) divides R(m,x) in the polynomial ring Z[x] if n divides m.
The polynomials R(n,-x), n >= 2, satisfy a Riemann hypothesis: their zeros lie on the vertical line Re x = 1/2 in the complex plane. Compare with A094440. (End)

A131050 (1/5) * (A007318^4 - A007318^(-1)).

Original entry on oeis.org

1, 3, 2, 13, 9, 3, 51, 52, 18, 4, 205, 255, 130, 30, 5, 819, 1230, 765, 260, 45, 6, 3277, 5733, 4305, 1785, 455, 63, 7, 13107, 26216, 22932, 11480, 3570, 728, 84, 8
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Row sums = powers of 5: (1, 5, 25, 125, ...).
Left border = A015521: (1, 3, 13, 51, 205, 819, ...).

Examples

			First few rows of the triangle:
    1;
    3,   2;
   13,   9,   3;
   51,  52,  18,  4;
  205, 255, 130, 30, 5;
  ...
		

Crossrefs

Formula

Let P = Pascal's triangle, A007318. Then A131050 = (1/5) * (P^4 - 1/P); deleting the right border of zeros.

A131049 (1/4) * (A007318^3 - A007318^(-1)).

Original entry on oeis.org

1, 2, 2, 7, 6, 3, 20, 28, 12, 4, 61, 100, 70, 20, 5, 182, 366, 300, 140, 30, 6, 547, 1274, 1281, 700, 245, 42, 7, 1640, 4376, 5096, 3416, 1400, 392, 56, 8, 4921, 14760, 19692, 15288, 7686, 2520, 588, 72, 9
Offset: 1

Views

Author

Gary W. Adamson, Jun 12 2007

Keywords

Comments

Row sums = powers of 4: (1, 4, 16, 64, ...).
Left border = A015518: (1, 2, 7, 20, 61, 182, ...).

Examples

			First few rows of the triangle:
    1;
    2,    2;
    7,    6,    3;
   20,   28,   12,   4;
   61,  100,   70,  20,   5;
  182,  366,  300, 140,  30,  6;
  547, 1274, 1281, 700, 245, 42, 7;
  ...
		

Crossrefs

Formula

(1/4) * (P^3 - 1/P), where P = Pascal's triangle, A007318. Delete right border of zeros.
Showing 1-9 of 9 results.