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-8 of 8 results.

A097805 Number of compositions of n with k parts, T(n, k) = binomial(n-1, k-1) for n, k >= 1 and T(n, 0) = 0^n, triangle read by rows for n >= 0 and 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 3, 1, 0, 1, 4, 6, 4, 1, 0, 1, 5, 10, 10, 5, 1, 0, 1, 6, 15, 20, 15, 6, 1, 0, 1, 7, 21, 35, 35, 21, 7, 1, 0, 1, 8, 28, 56, 70, 56, 28, 8, 1, 0, 1, 9, 36, 84, 126, 126, 84, 36, 9, 1, 0, 1, 10, 45, 120, 210, 252, 210, 120, 45, 10, 1, 0, 1, 11, 55, 165, 330, 462, 462, 330, 165, 55, 11, 1
Offset: 0

Views

Author

Paul Barry, Aug 25 2004

Keywords

Comments

Previous name was: Riordan array (1, 1/(1-x)) read by rows.
Note this Riordan array would be denoted (1, x/(1-x)) by some authors.
Columns have g.f. (x/(1-x))^k. Reverse of A071919. Row sums are A011782. Antidiagonal sums are Fibonacci(n-1). Inverse as Riordan array is (1, 1/(1+x)). A097805=B*A059260*B^(-1), where B is the binomial matrix.
(0,1)-Pascal triangle. - Philippe Deléham, Nov 21 2006
(n+1) * each term of row n generates triangle A127952: (1; 0, 2; 0, 3, 3; 0, 4, 8, 4; ...). - Gary W. Adamson, Feb 09 2007
Triangle T(n,k), 0<=k<=n, read by rows, given by [0,1,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 12 2008
From Paul Weisenhorn, Feb 09 2011: (Start)
Triangle read by rows: T(r,c) is the number of unordered partitions of n=r*(r+1)/2+c into (r+1) parts < (r+1) and at most pairs of equal parts and parts in neighboring pairs have difference 2.
Triangle read by rows: T(r,c) is the number of unordered partitions of the number n=r*(r+1)/2+(c-1) into r parts < (r+1) and at most pairs of equal parts and parts in neighboring pairs have difference 2. (End)
Triangle read by rows: T(r,c) is the number of ordered partitions (compositions) of r into c parts. - Juergen Will, Jan 04 2016
From Tom Copeland, Oct 25 2012: (Start)
Given a basis composed of a sequence of polynomials p_n(x) characterized by ladder (creation / annihilation, or raising / lowering) operators defined by R p_n(x) = p_(n+1)(x) and L p_n(x) = n p_(n-1)(x) with p_0(x)=1, giving the number operator # p_n(x) = RL p_n(x) = n p_n(x), the lower triangular padded Pascal matrix Pd (A097805) serves as a matrix representation of the operator exp(R^2*L) = exp(R#) =
1) exp(x^2D) for the set x^n and
2) D^(-1) exp(t*x)D for the set x^n/n! (see A218234).
(End)
From James East, Apr 11 2014: (Start)
Square array a(m,n) with m,n=0,1,2,... read by off-diagonals.
a(m,n) gives the number of order-preserving functions f:{1,...,m}->{1,...,n}. Order-preserving means that x
a(n,n)=A088218(n) is the size of the semigroup O_n of all order-preserving transformations of {1,...,n}.
Read as a triangle, this sequence may be obtained by augmenting Pascal's triangle by appending the column 1,0,0,0,... on the left.
(End)
A formula based on the partitions of n with largest part k is given as a Sage program below. The 'conjugate' formula leads to A048004. - Peter Luschny, Jul 13 2015
From Wolfdieter Lang, Feb 17 2017: (Start)
The transposed of this lower triangular Riordan matrix of the associated type T provides the transition matrix between the monomial basis {x^n}, n >= 0, and the basis {y^n}, n >= 0, with y = x/(1-x): x^0 = 1 = y^0, x^n = Sum_{m >= n} Ttrans(n,m) y^m, for n >= 1, with Ttrans(n,m) = binomial(m-1,n-1).
Therefore, if a transformation with this Riordan matrix from a sequence {a} to the sequence {b} is given by b(n) = Sum_{m=0..n} T(n, m)*a(m), with T(n, m) = binomial(n-1, m-1), for n >= 1, then Sum_{n >= 0} a(n)*x^n = Sum_{n >= 0} b(n)*y^n, with y = x/(1-x) and vice versa. This is a modified binomial transformation; the usual one belongs to the Pascal Riordan matrix A007318. (End)
From Gus Wiseman, Jan 23 2022: (Start)
Also the number of compositions of n with alternating sum k, with k ranging from -n to n in steps of 2. For example, row n = 6 counts the following compositions (empty column indicated by dot):
. (15) (24) (33) (42) (51) (6)
(141) (132) (123) (114)
(1113) (231) (222) (213)
(1212) (1122) (321) (312)
(1311) (1221) (1131) (411)
(2112) (2121)
(2211) (3111)
(11121) (11112)
(12111) (11211)
(111111) (21111)
The reverse-alternating version is the same. Counting compositions by all three parameters (sum, length, alternating sum) gives A345197. Compositions of 2n with alternating sum 2k with k ranging from -n + 1 to n are A034871. (End)
Also the convolution triangle of A000012. - Peter Luschny, Oct 07 2022
From Sergey Kitaev, Nov 18 2023: (Start)
Number of permutations of length n avoiding simultaneously the patterns 123 and 132 with k right-to-left maxima. A right-to-left maximum in a permutation a(1)a(2)...a(n) is position i such that a(j) < a(i) for all i < j.
Number of permutations of length n avoiding simultaneously the patterns 231 and 312 with k right-to-left minima (resp., left-to-right maxima). A right-to-left minimum (resp., left-to-right maximum) in a permutation a(1)a(2)...a(n) is position i such that a(j) > a(i) for all j > i (resp., a(j) < a(i) for all j < i).
Number of permutations of length n avoiding simultaneously the patterns 213 and 312 with k right-to-left maxima (resp., left-to-right maxima).
Number of permutations of length n avoiding simultaneously the patterns 213 and 231 with k right-to-left maxima (resp., right-to-left minima). (End)

Examples

			G.f. = 1 + x * (x + x^3 * (1 + x) + x^6 * (1 + x)^2 + x^10 * (1 + x)^3 + ...). - _Michael Somos_, Aug 20 2006
The triangle T(n, k) begins:
n\k  0 1 2  3  4   5   6  7  8 9 10 ...
0:   1
1:   0 1
2:   0 1 1
3:   0 1 2  1
4:   0 1 3  3  1
5:   0 1 4  6  4   1
6:   0 1 5 10 10   5   1
7:   0 1 6 15 20  15   6  1
8:   0 1 7 21 35  35  21  7  1
9:   0 1 8 28 56  70  56 28  8 1
10:  0 1 9 36 84 126 126 84 36 9  1
... reformatted _Wolfdieter Lang_, Jul 31 2017
From _Paul Weisenhorn_, Feb 09 2011: (Start)
T(r=5,c=3) = binomial(4,2) = 6 unordered partitions of the number n = r*(r+1)/2+c = 18 with (r+1)=6 summands: (5+5+4+2+1+1), (5+5+3+3+1+1), (5+4+4+3+1+1), (5+5+3+2+2+1), (5+4+4+2+2+1), (5+4+3+3+2+1).
T(r=5,c=3) = binomial(4,2) = 6 unordered partitions of the number n = r*(r+1)/2+(c-1) = 17 with r=5 summands: (5+5+4+2+1), (5+5+3+3+1), (5+5+3+2+2), (5+4+4+3+1), (5+4+4+2+2), (5+4+3+3+2).  (End)
From _James East_, Apr 11 2014: (Start)
a(0,0)=1 since there is a unique (order-preserving) function {}->{}.
a(m,0)=0 for m>0 since there is no function from a nonempty set to the empty set.
a(3,2)=4 because there are four order-preserving functions {1,2,3}->{1,2}: these are [1,1,1], [2,2,2], [1,1,2], [1,2,2]. Here f=[a,b,c] denotes the function defined by f(1)=a, f(2)=b, f(3)=c.
a(2,3)=6 because there are six order-preserving functions {1,2}->{1,2,3}: these are [1,1], [1,2], [1,3], [2,2], [2,3], [3,3].
(End)
		

References

  • D. E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, Part 1, Section 7.2.1.3, 2011.

Crossrefs

Case m=0 of the polynomials defined in A278073.
Cf. A000012 (diagonal), A011782 (row sums), A088218 (central terms).
The terms just left of center in odd-indexed rows are A001791, even A002054.
The odd-indexed rows are A034871.
Row sums without the center are A058622.
The unordered version is A072233, without zeros A008284.
Right half without center has row sums A027306(n-1).
Right half with center has row sums A116406(n).
Left half without center has row sums A294175(n-1).
Left half with center has row sums A058622(n-1).
A025047 counts alternating compositions.
A098124 counts balanced compositions, unordered A047993.
A106356 counts compositions by number of maximal anti-runs.
A344651 counts partitions by sum and alternating sum.
A345197 counts compositions by sum, length, and alternating sum.

Programs

  • Maple
    b:= proc(n, i, p) option remember; `if`(n=0, p!, `if`(i<1, 0,
          expand(add(b(n-i*j, i-1, p+j)/j!*x^j, j=0..n/i))))
        end:
    T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2, 0)):
    seq(T(n), n=0..20);  # Alois P. Heinz, May 25 2014
    # Alternatively:
    T := proc(k,n) option remember;
    if k=n then 1 elif k=0 then 0 else
    add(T(k-1,n-i), i=1..n-k+1) fi end:
    A097805 := (n,k) -> T(k,n):
    for n from 0 to 12 do seq(A097805(n,k), k=0..n) od; # Peter Luschny, Mar 12 2016
    # Uses function PMatrix from A357368.
    PMatrix(10, n -> 1);  # Peter Luschny, Oct 07 2022
  • Mathematica
    T[0, 0] = 1; T[n_, k_] := Binomial[n-1, k-1]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Sep 03 2014, after Paul Weisenhorn *)
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Length[#]==k&]],{n,0,10},{k,0,n}] (* Gus Wiseman, Jan 23 2022 *)
  • PARI
    {a(n) = my(m); if( n<2, n==0, n--; m = (sqrtint(8*n + 1) - 1)\2; binomial(m-1, n - m*(m + 1)/2))}; /* Michael Somos, Aug 20 2006 */
    
  • PARI
    T(n,k) = if (k==0, 0^n, binomial(n-1, k-1)); \\ Michel Marcus, May 06 2022
    
  • PARI
    row(n) = vector(n+1, k, k--; if (k==0, 0^n, binomial(n-1, k-1))); \\ Michel Marcus, May 06 2022
    
  • Python
    from math import comb
    def T(n, k): return comb(n-1, k-1) if k != 0 else k**n  # Peter Luschny, May 06 2022
  • Sage
    # Illustrates a basic partition formula, is not efficient as a program for large n.
    def A097805_row(n):
        r = []
        for k in (0..n):
            s = 0
            for q in Partitions(n, max_part=k, inner=[k]):
                s += mul(binomial(q[j],q[j+1]) for j in range(len(q)-1))
            r.append(s)
        return r
    [A097805_row(n) for n in (0..9)] # Peter Luschny, Jul 13 2015
    

Formula

Number triangle T(n, k) defined by T(n,k) = Sum_{j=0..n} binomial(n, j)*if(k<=j, (-1)^(j-k), 0).
T(r,c) = binomial(r-1,c-1), 0 <= c <= r. - Paul Weisenhorn, Feb 09 2011
G.f.: (-1+x)/(-1+x+x*y). - R. J. Mathar, Aug 11 2015
a(0,0) = 1, a(n,k) = binomial(n-1,n-k) = binomial(n-1,k-1) Juergen Will, Jan 04 2016
G.f.: (x^1 + x^2 + x^3 + ...)^k = (x/(1-x))^k. - Juergen Will, Jan 04 2016
From Tom Copeland, Nov 15 2016: (Start)
E.g.f.: 1 + x*[e^((x+1)t)-1]/(x+1).
This padded Pascal matrix with the odd columns negated is NpdP = M*S = S^(-1)*M^(-1) = S^(-1)*M, where M(n,k) = (-1)^n A130595(n,k), the inverse Pascal matrix with the odd rows negated, S is the summation matrix A000012, the lower triangular matrix with all elements unity, and S^(-1) = A167374, a finite difference matrix. NpdP is self-inverse, i.e., (M*S)^2 = the identity matrix, and has the e.g.f. 1 - x*[e^((1-x)t)-1]/(1-x).
M = NpdP*S^(-1) follows from the well-known recursion property of the Pascal matrix, implying NpdP = M*S.
The self-inverse property of -NpdP is implied by the self-inverse relation of its embedded signed Pascal submatrix M (cf. A130595). Also see A118800 for another proof.
Let P^(-1) be A130595, the inverse Pascal matrix. Then T = A200139*P^(-1) and T^(-1) = padded P^(-1) = P*A097808*P^(-1). (End)
The center (n>0) is T(2n+1,n+1) = A000984(n) = 2*A001700(n-1) = 2*A088218(n) = A126869(2n) = 2*A138364(2n-1). - Gus Wiseman, Jan 25 2022

Extensions

Corrected by Philippe Deléham, Oct 05 2005
New name using classical terminology by Peter Luschny, Feb 05 2019

A214292 Triangle read by rows: T(n,k) = T(n-1,k-1) + T(n-1,k), 0 < k < n with T(n,0) = n and T(n,n) = -n.

Original entry on oeis.org

0, 1, -1, 2, 0, -2, 3, 2, -2, -3, 4, 5, 0, -5, -4, 5, 9, 5, -5, -9, -5, 6, 14, 14, 0, -14, -14, -6, 7, 20, 28, 14, -14, -28, -20, -7, 8, 27, 48, 42, 0, -42, -48, -27, -8, 9, 35, 75, 90, 42, -42, -90, -75, -35, -9, 10, 44, 110, 165, 132, 0, -132, -165, -110, -44, -10
Offset: 0

Author

Reinhard Zumkeller, Jul 12 2012

Keywords

Examples

			The triangle begins:
    0:                              0
    1:                            1   -1
    2:                          2   0   -2
    3:                       3    2   -2   -3
    4:                     4    5   0   -5   -4
    5:                  5    9    5   -5   -9   -5
    6:                6   14   14   0  -14  -14   -6
    7:             7   20   28   14  -14  -28  -20   -7
    8:           8   27   48   42   0  -42  -48  -27   -8
    9:        9   35   75   90   42  -42  -90  -75  -35   -9
   10:     10   44  110  165  132   0 -132 -165 -110  -44  -10
   11:  11   54  154  275  297  132 -132 -297 -275 -154  -54  -11  .
		

Programs

  • Haskell
    a214292 n k = a214292_tabl !! n !! k
    a214292_row n = a214292_tabl !! n
    a214292_tabl = map diff $ tail a007318_tabl
       where diff row = zipWith (-) (tail row) row
  • Mathematica
    row[n_] := Table[Binomial[n, k], {k, 0, n}] // Differences;
    T[n_, k_] := row[n + 1][[k + 1]];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jul 31 2018 *)

Formula

T(n,k) = A007318(n+1,k+1) - A007318(n+1,k), 0<=k<=n, i.e. first differences of rows in Pascal's triangle;
T(n,k) = -T(n,k);
row sums and central terms equal 0, cf. A000004;
sum of positive elements of n-th row = A014495(n+1);
T(n,0) = n;
T(n,1) = A000096(n-2) for n > 1; T(n,1) = - A080956(n) for n > 0;
T(n,2) = A005586(n-4) for n > 3; T(n,2) = A129936(n-2);
T(n,3) = A005587(n-6) for n > 5;
T(n,4) = A005557(n-9) for n > 8;
T(n,5) = A064059(n-11) for n > 10;
T(n,6) = A064061(n-13) for n > 12;
T(n,7) = A124087(n) for n > 14;
T(n,8) = A124088(n) for n > 16;
T(2*n+1,n) = T(2*n+2,n) = A000108(n+1), Catalan numbers;
T(2*n+3,n) = A000245(n+2);
T(2*n+4,n) = A002057(n+1);
T(2*n+5,n) = A000344(n+3);
T(2*n+6,n) = A003517(n+3);
T(2*n+7,n) = A000588(n+4);
T(2*n+8,n) = A003518(n+4);
T(2*n+9,n) = A001392(n+5);
T(2*n+10,n) = A003519(n+5);
T(2*n+11,n) = A000589(n+6);
T(2*n+12,n) = A090749(n+6);
T(2*n+13,n) = A000590(n+7).

A059260 Triangle read by rows giving coefficient T(i,j) of x^i y^j in 1/(1-y-x*y-x^2) = 1/((1+x)(1-x-y)) for (i,j) = (0,0), (1,0), (0,1), (2,0), (1,1), (0,2), ...

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 0, 2, 2, 1, 1, 2, 4, 3, 1, 0, 3, 6, 7, 4, 1, 1, 3, 9, 13, 11, 5, 1, 0, 4, 12, 22, 24, 16, 6, 1, 1, 4, 16, 34, 46, 40, 22, 7, 1, 0, 5, 20, 50, 80, 86, 62, 29, 8, 1, 1, 5, 25, 70, 130, 166, 148, 91, 37, 9, 1, 0, 6, 30, 95, 200, 296, 314, 239, 128, 46, 10, 1
Offset: 0

Author

N. J. A. Sloane, Jan 23 2001

Keywords

Comments

Coefficients of the (left, normalized) shifted cyclotomic polynomial. Or, coefficients of the basic n-th q-series for q=-2. Indeed, let Y_n(x) = Sum_{k=0..n} x^k, having as roots all the n-th roots of unity except for 0; then coefficients in x of (-1)^n Y_n(-x-1) give exactly the n-th row of A059260 and a practical way to compute it. - Olivier Gérard, Jul 30 2002
The maximum in the (2n)-th row is T(n,n), which is A026641; also T(n,n) ~ (2/3)*binomial(2n,n). The maximum in the (2n-1)-th row is T(n-1,n), which is A014300 (but T does not have the same definition as in A026637); also T(n-1,n) ~ (1/3)*binomial(2n,n). Here is a generalization of the formula given in A026641: T(i,j) = Sum_{k=0..j} binomial(i+k-x,j-k)*binomial(j-k+x,k) for all x real (the proof is easy by induction on i+j using T(i,j) = T(i-1,j) + T(i,j-1)). - Claude Morin, May 21 2002
The second greatest term in the (2n)-th row is T(n-1,n+1), which is A014301; the second greatest term in the (2n+1)-th row is T(n+1,n) = 2*T(n-1,n+1), which is 2*A014301. - Claude Morin
Diagonal sums give A008346. - Paul Barry, Sep 23 2004
Riordan array (1/(1-x^2), x/(1-x)). As a product of Riordan arrays, factors into the product of (1/(1+x),x) and (1/(1-x),1/(1-x)) (binomial matrix). - Paul Barry, Oct 25 2004
Signed version is A239473 with relations to partial sums of sequences. - Tom Copeland, Mar 24 2014
From Robert Coquereaux, Oct 01 2014: (Start)
Columns of the triangle (cf. Example below) give alternate partial sums along nw-se diagonals of the Pascal triangle, i.e., sequences A000035, A004526, A002620 (or A087811), A002623 (or A173196), A001752, A001753, A001769, A001779, A001780, A001781, A001786, A001808, etc.
The dimension of the space of closed currents (distributional forms) of degree p on Gr(n), the Grassmann algebra with n generators, equivalently, the dimension of the space of Gr(n)-valued symmetric multilinear forms with vanishing graded divergence, is V(n,p) = 2^n T(p,n-1) - (-1)^p.
If p is odd V(n,p) is also the dimension of the cyclic cohomology group of order p of the Z2 graded algebra Gr(n).
If p is even the dimension of this cohomology group is V(n,p)+1.
Cf. A193844. (End)
From Peter Bala, Feb 07 2024: (Start)
The following remarks assume the row indexing starts at n = 1.
The sequence of row polynomials R(n,x), beginning R(1,x) = 1, R(2,x) = x, R(3,x) = 1 + x + x^2 , ..., is a strong divisibility sequence of polynomials in the ring Z[x]; that is, for all positive integers n and m, poly_gcd( R(n,x), R(m,x)) = R(gcd(n, m), x) - apply Norfleet (2005), Theorem 3. Consequently, the polynomial sequence {R(n,x): n >= 1} is a divisibility sequence; that is, if n divides m then R(n,x) divides R(m,x) in Z[x]. (End)
From Miquel A. Fiol, Oct 04 2024: (Start)
For j>=1, T(i,j) is the independence number of the (i-j)-supertoken graph FF_(i-j)(S_j) of the star graph S_j with j points.
(Given a graph G on n vertices and an integer k>=1, the k-supertoken (or reduced k-th power) FF_k(G) of G has vertices representing configurations of k indistinguishable tokens in the (not necessarily different) vertices of G, with two configurations being adjacent if one can be obtained from the other by moving one token along an edge. See an example below.)
Following the suggestion of Peter Munn, the k-supertoken graph FF_k(S_j) can also be defined as follows: Consider the Lattice graph L(k,j), whose vertices are the k^j j-vectors with elements in the set {0,..,k-1}, two being adjacent if they differ in just one coordinate by one unity. Then, FF_k(S_j) is the subgraph of L(k+1,j) induced by the vertices at distance at most k from (0,..,0). (End)

Examples

			Triangle begins
  1;
  0,  1;
  1,  1,  1;
  0,  2,  2,  1;
  1,  2,  4,  3,  1;
  0,  3,  6,  7,  4,  1;
  1,  3,  9, 13, 11,  5,  1;
  0,  4, 12, 22, 24, 16,  6,  1;
  1,  4, 16, 34, 46, 40, 22,  7,  1;
  0,  5, 20, 50, 80, 86, 62, 29,  8,  1;
Sequences obtained with _Miquel A. Fiol_'s Sep 30 2024 formula of A(n,c1,c2) for other values of (c1,c2). (In the table, rows are indexed by c1=0..6 and columns by c2=0..6):
A000007  A000012  A000027  A025747  A000292* A000332* A000389*
A059841  A008619  A087811* A002623  A001752  A001753  A001769
A193356  A008794* A005993  A005994  -------  -------  -------
-------  -------  -------  A005995  A018210  -------  A052267
-------  -------  -------  -------  A018211  A018212  -------
-------  -------  -------  -------  -------  A018213  A018214
-------  -------  -------  -------  -------  -------  A062136
*requires offset adjustment.
The 2-supertoken FF_2(S_3) of the star graph S_3 with central vertex 1 and peripheral vertices 2,3,4. (The vertex `ij' of FF_2(S_3) represents the configuration of one token in `ì' and the other token in `j'). The T(5,3)=7 independent vertices are 22, 24, 44, 23, 11, 34, and 33.
     22--12---24---14---44
          | \    / |
         23   11   34
            \  |  /
              13
               |
              33
		

Crossrefs

Cf. A059259. Row sums give A001045.
Seen as a square array read by antidiagonals this is the coefficient of x^k in expansion of 1/((1-x^2)*(1-x)^n) with rows A002620, A002623, A001752, A001753, A001769, A001779, A001780, A001781, A001786, A001808 etc. (allowing for signs). A058393 would then effectively provide the table for nonpositive n. - Henry Bottomley, Jun 25 2001

Programs

  • Maple
    read transforms; 1/(1-y-x*y-x^2); SERIES2(%,x,y,12); SERIES2TOLIST(%,x,y,12);
  • Mathematica
    t[n_, k_] := Sum[ (-1)^(n-j)*Binomial[j, k], {j, 0, n}]; Flatten[ Table[t[n, k], {n, 0, 12}, {k, 0, n}]] (* Jean-François Alcover, Oct 20 2011, after Paul Barry *)
  • PARI
    T(n, k) = sum(j=0, n, (-1)^(n - j)*binomial(j, k));
    for(n=0, 12, for(k=0, n, print1(T(n, k),", ");); print();) \\ Indranil Ghosh, Apr 11 2017
    
  • Python
    from sympy import binomial
    def T(n, k): return sum((-1)**(n - j)*binomial(j, k) for j in range(n + 1))
    for n in range(13): print([T(n, k) for k in range(n + 1)]) # Indranil Ghosh, Apr 11 2017
  • Sage
    def A059260_row(n):
        @cached_function
        def prec(n, k):
            if k==n: return 1
            if k==0: return 0
            return -prec(n-1,k-1)-sum(prec(n,k+i-1) for i in (2..n-k+1))
        return [(-1)^(n-k+1)*prec(n+1, n-k+1) for k in (1..n)]
    for n in (1..9): print(A059260_row(n)) # Peter Luschny, Mar 16 2016
    

Formula

G.f.: 1/(1-y-x*y-x^2) = 1 + y + x^2 + xy + y^2 + 2x^2y + 2xy^2 + y^3 + ...
E.g.f: (exp(-t)+(x+1)*exp((x+1)*t))/(x+2). - Tom Copeland, Mar 19 2014
O.g.f. (n-th row): ((-1)^n+(x+1)^(n+1))/(x+2). - Tom Copeland, Mar 19 2014
T(i, 0) = 1 if i is even or 0 if i is odd, T(0, i) = 1 and otherwise T(i, j) = T(i-1, j) + T(i, j-1); also T(i, j) = Sum_{m=j..i+j} (-1)^(i+j+m)*binomial(m, j). - Robert FERREOL, May 17 2002
T(i, j) ~ (i+j)/(2*i+j)*binomial(i+j, j); more precisely, abs(T(i, j)/binomial(i+j, j) - (i+j)/(2*i+j) )<=1/(4*(i+j)-2); the proof is by induction on i+j using the formula 2*T(i, j) = binomial(i+j, j)+T(i, j-1). - Claude Morin, May 21 2002
T(n, k) = Sum_{j=0..n} (-1)^(n-j)binomial(j, k). - Paul Barry, Aug 25 2004
T(n, k) = Sum_{j=0..n-k} binomial(n-j, j)*binomial(j, n-k-j). - Paul Barry, Jul 25 2005
Equals A097807 * A007318. - Gary W. Adamson, Feb 21 2007
Equals A128173 * A007318 as infinite lower triangular matrices. - Gary W. Adamson, Feb 17 2007
Equals A130595*A097805*A007318 = (inverse Pascal matrix)*(padded Pascal matrix)*(Pascal matrix) = A130595*A200139. Inverse is A097808 = A130595*(padded A130595)*A007318. - Tom Copeland, Nov 14 2016
T(i, j) = binomial(i+j, j)-T(i-1, j). - Laszlo Major, Apr 11 2017
Recurrence for row polynomials (with row indexing starting at n = 1): R(n,x) = x*R(n-1,x) + (x + 1)*R(n-2,x) with R(1,x) = 1 and R(2,x) = x. - Peter Bala, Feb 07 2024
From Miquel A. Fiol, Sep 30 2024: (Start)
The triangle can be seen as a slice of a 3-dimensional table that links it to well-known sequences as follows.
The j-th column of the triangle, T(i,j) for i >= j, equals A(n,c1,c2) = Sum_{k=0..floor(n/2)} binomial(c1+2*k-1,2*k)*binomial(c2+n-2*k-1,n-2*k) when c1=1, c2=j, and n=i-j.
This gives T(i,j) = Sum_{k=0..floor((i-j)/2)} binomial(i-2*k-1, j-1). For other values of (c1,c2), see the example below. (End)

Extensions

Formula corrected by Philippe Deléham, Jan 11 2014

A046802 T(n, k) = Sum_{j=k..n} binomial(n, j)*E1(j, j-k), where E1 are the Eulerian numbers A173018. Triangle read by rows, T(n, k) for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 7, 1, 1, 15, 33, 15, 1, 1, 31, 131, 131, 31, 1, 1, 63, 473, 883, 473, 63, 1, 1, 127, 1611, 5111, 5111, 1611, 127, 1, 1, 255, 5281, 26799, 44929, 26799, 5281, 255, 1, 1, 511, 16867, 131275, 344551, 344551, 131275, 16867, 511, 1, 1, 1023, 52905
Offset: 0

Keywords

Comments

T(n,k) is the number of positroid cells of the totally nonnegative Grassmannian G+(k,n) (cf. Postnikov/Williams). It is the triangle of the h-vectors of the stellahedra. - Tom Copeland, Oct 10 2014
See A248727 for a simple transformation of the row polynomials of this entry that produces the umbral compositional inverses of the polynomials of A074909, related to the face polynomials of the simplices. - Tom Copeland, Jan 21 2015
From Tom Copeland, Jan 24 2015: (Start)
The reciprocal of this entry's e.g.f. is [t e^(-xt) - e^(-x)] / (t-1) = 1 - (1+t) x + (1+t+t^2) x^2/2! - (1+t+t^2+t^3) x^3/3! + ... = e^(q.(0;t)x), giving the base sequence (q.(0;t))^n = q_n(0;t) = (-1)^n [1-t^(n+1)] / (1-t) for the umbral compositional inverses (q.(0;t)+z)^n = q_n(z;t) of the Appell polynomials associated with this entry, p_n(z;t) below, i.e., q_n(p.(z;t)) = z^n = p_n(q.(z;t)), in umbral notation. The relations in A133314 then apply between the two sets of base polynomials. (Inserted missing index in a formula - Mar 03 2016.)
The associated o.g.f. for the umbral inverses is Og(x) = x / (1-x q.(0:t)) = x / [(1+x)(1+tx)] = x / [1+(1+t)x+tx^2]. Applying A134264 to h(x) = x / Og(x) = 1 + (1+t) x + t x^2 leads to an o.g.f. for the Narayana polynomials A001263 as the comp. inverse Oginv(x) = [1-(1+t)x-sqrt[1-2(1+t)x+((t-1)x)^2]] / (2xt). Note that Og(x) gives the signed h-polynomials of the simplices and that Oginv(x) gives the h-polynomials of the simplicial duals of the Stasheff polynomials, or type A associahedra. Contrast this with A248727 = A046802 * A007318, which has o.g.f.s related to the corresponding f-polynomials. (End)
The Appell polynomials p_n(x;t) in the formulas below specialize to the Swiss-knife polynomials of A119879 for t = -1, so the Springer numbers A001586 are given by 2^n p_n(1/2;-1). - Tom Copeland, Oct 14 2015
The row polynomials are the h-polynomials associated to the stellahedra, whose f-polynomials are the row polynomials of A248727. Cf. page 60 of the Buchstaber and Panov link. - Tom Copeland, Nov 08 2016
The row polynomials are the h-polynomials of the stellohedra, which enumerate partial permutations according to descents. Cf. Section 10.4 of the Postnikov-Reiner-Williams reference. - Lauren Williams, Jul 05 2022
From p. 60 of the Buchstaber and Panov link, S = P * C / T where S, P, C, and T are the bivariate e.g.f.s of the h vectors of the stellahedra, permutahedra, hypercubes, and (n-1)-simplices, respectively. - Tom Copeland, Jan 09 2017
The number of Le-diagrams of type (k, n) this means the diagram uses the bounding box size k x (n-k), equivalently the number of Grassmann necklaces of type (k, n) and also the number of decorated permutations with k anti-exceedances. - Thomas Scheuerle, Dec 29 2024

Examples

			The triangle T(n, k) begins:
n\k 0   1     2      3      4      5      6     7
0:  1
1:  1   1
2:  1   3     1
3:  1   7     7      1
4:  1  15    33     15      1
5:  1  31   131    131     31      1
6:  1  63   473    883    473     63      1
7:  1 127  1611   5111   5111   1611    127     1
... Reformatted. - _Wolfdieter Lang_, Feb 14 2015
		

References

  • L. Comtet, Advanced Combinatorics, Reidel, Holland, 1974, page 245 [From Roger L. Bagula, Nov 21 2009]
  • D. Singh, The numbers L(m,n) and their relations with prepared Bernoulli and Eulerian numbers, Math. Student, 20 (1952), 66-70.

Programs

  • Maple
    T := (n, k) -> add(binomial(n, r)*combinat:-eulerian1(r, r-k), r = k .. n):
    for n from 0 to 8 do seq(T(n, k), k=0..n) od; # Peter Luschny, Jun 27 2018
  • Mathematica
    t[, 1] = 1; t[n, n_] = 1; t[n_, 2] = 2^(n-1)-1;
    t[n_, k_] = Sum[((i-k+1)^i*(k-i)^(n-i-1) - (i-k+2)^i*(k-i-1)^(n-i-1))*Binomial[n-1, i], {i, 0, k-1}];
    T[n_, k_] := t[n+1, k+1]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten
    (* Jean-François Alcover, Jan 22 2015, after Tom Copeland *)
    T[ n_, k_] := Coefficient[n! SeriesCoefficient[(1-x) Exp[t] / (1 - x Exp[(1-x) t]), {t, 0, n}] // Simplify, x, k];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] (* Michael Somos, Jan 22 2015 *)

Formula

E.g.f.: (y-1)*exp(x*y)/(y-exp((y-1)*x)). - Vladeta Jovovic, Sep 20 2003
p(t,x) = (1 - x)*exp(t)/(1 - x*exp(t*(1 - x))). - Roger L. Bagula, Nov 21 2009
With offset=0, T(n,0)=1 otherwise T(n,k) = sum_{i=0..k-1} C(n,i)((i-k)^i*(k-i+1)^(n-i) - (i-k+1)^i*(k-i)^(n-i)) (cf. Williams). - Tom Copeland, Oct 10 2014
With offset 0, T = A007318 * A123125. Second column is A000225; 3rd, appears to be A066810. - Tom Copeland, Jan 23 2015
A raising operator (with D = d/dx) associated with this entry's row polynomials is R = x + t + (1-t) / [1-t e^{(1-t)D}] = x + t + 1 + t D + (t+t^2) D^2/2! + (t+4t^2+t^3) D^3/3! + ... , containing the e.g.f. for the Eulerian polynomials of A123125. Then R^n 1 = (p.(0;t)+x)^n = p_n(x;t) are the Appell polynomials with this entry's row polynomials p_n(0;t) as the base sequence. Examples of this formalism are given in A028246 and A248727. - Tom Copeland, Jan 24 2015
With offset 0, T = A007318*(padded A090582)*(inverse of A097805) = A007318*(padded A090582)*(padded A130595) = A007318*A123125 = A007318*(padded A090582)*A007318*A097808*A130595, where padded matrices are of the form of padded A007318, which is A097805. Inverses of padded matrices are just the padded versions of inverses of the unpadded matrices. This relates the face vectors, or f-vectors, and h-vectors of the permutahedra / permutohedra to those of the stellahedra / stellohedra. - Tom Copeland, Nov 13 2016
Umbrally, the row polynomials (offset 0) are r_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A123125. - Tom Copeland, Nov 16 2016
From the previous umbral statement, OP(x,d/dy) y^n = (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = (1-x)/(1-x*exp((1-x)y)), the e.g.f. of A123125, so OP(x,d/dy) y^n evaluated at y = 1 is r_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018
Consolidating some formulas in this entry and A248727, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of this entry (A046802, the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
From Peter Luschny, Apr 30 2021: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A122045(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007047(n).
Sum_{k=0..n} T(n, n-k) = A000522(n).
Sum_{k=0..n} T(n-k, k) = Sum_{k=0..n} (n - k)^k = A026898(n-1) for n >= 1.
Sum_{k=0..n} k*T(n, k) = A036919(n) = floor(n*n!*e/2).
(End)

Extensions

More terms from Vladeta Jovovic, Sep 20 2003
First formula corrected by Wolfdieter Lang, Feb 14 2015
Offset set to 0 and edited by Peter Luschny, Apr 30 2021

A248727 A046802(x,y) --> A046802(x,y+1), transform of e.g.f. for the graded number of positroids of the totally nonnegative Grassmannians G+(k,n); enumerates faces of the stellahedra.

Original entry on oeis.org

1, 2, 1, 5, 5, 1, 16, 24, 10, 1, 65, 130, 84, 19, 1, 326, 815, 720, 265, 36, 1, 1957, 5871, 6605, 3425, 803, 69, 1, 13700, 47950, 65646, 44240, 15106, 2394, 134, 1, 109601, 438404, 707840, 589106, 267134, 63896, 7094, 263, 1
Offset: 0

Author

Tom Copeland, Oct 12 2014

Keywords

Comments

This is a transform of A046802 treating it as an array of h-vectors, so y is replaced by (y+1) in the e.g.f. for A046802.
An e.g.f. for the reversed row polynomials with signs is given by exp(a.(0;t)x) = [e^{(1+t)x} [1+t(1-e^(-x))]]^(-1) = 1 - (1+2t)x + (1+5t+5t^2)x^2/2! + ... . The reciprocal is an e.g.f. for the reversed face polynomials of the simplices A074909, i.e., exp(b.(0;t)x) = e^{(1+t)x} [1+t(1-e^(-x))] = 1 + (1+2t)x +(1+3t+3t^2) x^2/2! + ... , so the relations of A133314 apply between the two sets of polynomials. In particular, umbrally [a.(0;t)+b.(0;t)]^n vanishes except for n=0 for which it's unity, implying the two sets of Appell polynomials formed from the two bases, a_n(z;t) = (a.(0;t)+z)^n and b_n(z;t) = (b.(0;t) + z)^n, are an umbral compositional inverse pair, i.e., b_n(a.(x;t);t)= x^n = a_n(b.(x;t);t). Raising operators for these Appell polynomials are related to the polynomials of A028246, whose reverse polynomials are given by A123125 * A007318. Compare: A248727 = A007318 * A123125 * A007318 and A046802 = A007318 * A123125. See A074909 for definitions and related links. - Tom Copeland, Jan 21 2015
The o.g.f. for the umbral inverses is Og(x) = x / (1 - x b.(0;t)) = x / [(1-tx)(1-(1+t)x)] = x + (1+2t) x^2 + (1+3t+3t^2) x^3 + ... . Its compositional inverse is an o.g.f for signed A033282, the reverse f-polynomials for the simplicial duals of the Stasheff polytopes, or associahedra of type A, Oginv(x) =[1+(1+2t)x-sqrt[1+2(1+2t)x+x^2]] / (2t(1+t)x) = x - (1+2t) x^2 + (1+5t+5t^2) x^3 + ... . Contrast this with the o.g.f.s related to the corresponding h-polynomials in A046802. - Tom Copeland, Jan 24 2015
Face vectors, or coefficients of the face polynomials, of the stellahedra, or stellohedra. See p. 59 of Buchstaber and Panov. - Tom Copeland, Nov 08 2016
See A008279 for a relation between the e.g.f.s enumerating the faces of permutahedra and stellahedra. - Tom Copeland, Nov 14 2016

Examples

			The triangle T(n, k) starts:
n\k    0     1     2     3     4    5   6  7 ...
1:     1
2:     2     1
3:     5     5     1
4:    16    24    10     1
5:    65   130    84    19     1
6:   326   815   720   265    36    1
7:  1957  5871  6605  3425   803   69   1
8: 13700 47950 65646 44240 15106 2394 134  1
... reformatted, _Wolfdieter Lang_, Mar 27 2015
		

Programs

  • Mathematica
    (* t = A046802 *) t[, 1] = 1; t[n, n_] = 1; t[n_, 2] = 2^(n - 1) - 1; t[n_, k_] = Sum[((i - k + 1)^i*(k - i)^(n - i - 1) - (i - k + 2)^i*(k - i - 1)^(n - i - 1))*Binomial[n - 1, i], {i, 0, k - 1}]; T[n_, j_] := Sum[Binomial[k, j]*t[n + 1, k + 1], {k, j, n}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 23 2015, after Tom Copeland *)

Formula

Let M(n,k)= sum{i=0,..,k-1, C(n,i)[(i-k)^i*(k-i+1)^(n-i)- (i-k+1)^i*(k-i)^(n-i)]} with M(n,0)=1. Then M(n,k)= A046802(n,k), and T(n,j)= sum(k=j,..,n, C(k,j)*M(n,k)) for j>0 with T(n,0)= 1 + sum(k=1,..,n, M(n,k)) for n>0 and T(0,0)=1.
E.g.f: y * exp[x*(y+1)]/[y+1-exp(x*y)].
Row sums are A007047. Row polynomials evaluated at -1 are unity. Row polynomials evaluated at -2 are A122045.
First column is A000522. Second column appears to be A036918/2 = (A001339-1)/2 = n*A000522(n)/2.
Second diagonal is A052944. (Changed from conjecture to fact on Nov 08 2016.)
The raising operator for the reverse row polynomials with row signs is R = x - (1+t) - t e^(-D) / [1 + t(1-e^(-D))] evaluated at x = 0, with D = d/dx. Also R = x - d/dD log[exp(a.(0;t)D], or R = - d/dz log[e^(-xz) exp(a.(0;t)z)] = - d/dz log[exp(a.(-x;t)z)] with the e.g.f. defined in the comments and z replaced by D. Note that t e ^(-D) / [1+t(1-e^(-D))] = t - (t+t^2) D + (t+3t^2+2t^3) D^2/2! - ... is an e.g.f. for the signed reverse row polynomials of A028246. - Tom Copeland, Jan 23 2015
Equals A007318*(padded A090582)*A007318*A097808 = A007318*(padded (A008292*A007318))*A007318*A097808 = A007318*A130850 = A007318*(mirror of A028246). Padded means in the same way that A097805 is padded A007318. - Tom Copeland, Nov 14 2016
Umbrally, the row polynomials are p_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A130850. - Tom Copeland, Nov 16 2016
From the previous umbral statement, OP(x,d/dy) y^n = (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = x/((1+x)*exp(-x*y) - 1), the e.g.f. of A130850, so OP(x,d/dy) y^n evaluated at y = 1 is p_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018
Consolidating some formulas in this entry and A046082, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of this entry (A248727, the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020

Extensions

Title expanded by Tom Copeland, Nov 08 2016

A156644 Mirror image of triangle A080233.

Original entry on oeis.org

1, 0, 1, -1, 1, 1, -2, 0, 2, 1, -3, -2, 2, 3, 1, -4, -5, 0, 5, 4, 1, -5, -9, -5, 5, 9, 5, 1, -6, -14, -14, 0, 14, 14, 6, 1, -7, -20, -28, -14, 14, 28, 20, 7, 1, -8, -27, -48, -42, 0, 42, 48, 27, 8, 1, -9, -35, -75, -90, -42, 42, 90, 75, 35, 9, 1
Offset: 0

Author

Philippe Deléham, Feb 12 2009

Keywords

Comments

Inverse of A239473. Equals A007318*A167374. - Tom Copeland, Nov 14 2016

Examples

			Triangle begins as:
   1;
   0,  1;
  -1,  1, 1;
  -2,  0, 2, 1;
  -3, -2, 2, 3, 1;
  -4, -5, 0, 5, 4, 1; ...
		

Crossrefs

Programs

  • Magma
    A156644:= func< n,k | ((2*k-n+1)/(k+1))*Binomial(n,k) >;
    [A156644(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Feb 28 2021
  • Mathematica
    Table[Binomial[n, k] -Binomial[n, k+1], {n,0,10}, {k,0,n}]//Flatten (* Michael De Vlieger, Nov 24 2016 *)
  • Sage
    def A156644(n,k): return ((2*k-n+1)/(k+1))*binomial(n,k)
    flatten([[A156644(n,k) for k in (0..n)] for n in (0..12)]) # G. C. Greubel, Feb 28 2021
    

Formula

T(n,k) = A080233(n,n-k) = (-1)^(n-k)*A097808(n,k).
T(n,k) = ((2*k-n+1)/(k+1))*binomial(n,k).
T(n,k) = T(n-1,k-1) + T(n-1,k), k>0, with T(n,0) = 1-n = A024000(n), T(n,n) = 1.
T(n,k) = binomial(n,k) - binomial(n,k+1) = Sum_{i=-k-1..k+1} (-1)^(i+1) * binomial(n,k+1+i) * binomial(n+2,k+1-i). - Mircea Merca, Apr 28 2012
Sum_{k=0..n} T(n, k) = A000012(n) = 1^n. - G. C. Greubel, Feb 28 2021

A112466 Riordan array ((1+2*x)/(1+x), x/(1+x)).

Original entry on oeis.org

1, 1, 1, -1, 0, 1, 1, -1, -1, 1, -1, 2, 0, -2, 1, 1, -3, 2, 2, -3, 1, -1, 4, -5, 0, 5, -4, 1, 1, -5, 9, -5, -5, 9, -5, 1, -1, 6, -14, 14, 0, -14, 14, -6, 1, 1, -7, 20, -28, 14, 14, -28, 20, -7, 1, -1, 8, -27, 48, -42, 0, 42, -48, 27, -8, 1, 1, -9, 35, -75, 90, -42, -42, 90, -75, 35, -9, 1, -1, 10, -44, 110, -165, 132, 0, -132, 165, -110, 44, -10, 1
Offset: 0

Author

Paul Barry, Sep 06 2005

Keywords

Comments

Inverse is A112465.
Triangle T(n,k), 0 <= k <= n, read by rows, given by [1, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Aug 07 2006; corrected by Philippe Deléham, Dec 11 2008
Equals A097808 when the first column is removed. - Georg Fischer, Jul 26 2023

Examples

			Triangle starts
   1;
   1,  1;
  -1,  0,  1;
   1, -1, -1,  1;
  -1,  2,  0, -2,  1;
   1, -3,  2,  2, -3,  1;
  -1,  4, -5,  0,  5, -4,  1;
From _Paul Barry_, Apr 08 2011: (Start)
Production matrix begins
   1,  1;
  -2, -1,  1;
   2,  0, -1,  1;
  -2,  0,  0, -1,  1;
   2,  0,  0,  0, -1,  1;
  -2,  0,  0,  0,  0, -1,  1;
   2,  0,  0,  0,  0,  0, -1,  1; (End)
		

Crossrefs

Columns: A248157(n+2) (k=1), (-1)^n*A080956(n-2) (k=2), (-1)^(n-1)*A254749(n-2) (k=3).

Programs

  • Magma
    A112466:= func< n,k | n eq 0 select 1 else (-1)^(n+k)*(Binomial(n,k) - 2*Binomial(n-1,k)) >;
    [A112466(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Apr 30 2025
    
  • Maple
    seq(seq( (-1)^(n-k)*(2*binomial(n-1, k-1)-binomial(n, k)), k=0..n), n=0..10); # G. C. Greubel, Feb 19 2020
  • Mathematica
    {1}~Join~Table[(Binomial[n, n - k] - 2 Binomial[n - 1, n - k - 1])*(-1)^(n - k), {n, 12}, {k, 0, n}] // Flatten (* Michael De Vlieger, Feb 18 2020 *)
  • PARI
    T(n,k) = (-1)^(n-k)*(binomial(n, n-k) - 2*binomial(n-1, n-k-1)); \\ Michel Marcus, Feb 19 2020
    
  • SageMath
    def A112466(n,k): return 1 if (n==0) else (-1)^(n+k)*(binomial(n,k) - 2*binomial(n-1,k))
    print(flatten([[A112466(n,k) for k in range(n+1)] for n in range(13)])) # G. C. Greubel, Apr 30 2025

Formula

Number triangle: T(n,k) = (-1)^(n-k)*(C(n, n-k) - 2*C(n-1, n-k-1)), with T(0,0) = 1.
T(2*n, n) = 0 (main diagonal).
Sum_{k=0..n} T(n, k) = 0 + [n=0] + 2*[n=1] (row sums).
Sum_{k=0..floor(n/2)} T(n-k, k) = (-1)^(n+1)*Fibonacci(n-2) (diagonal sums).
Sum_{k=0..n} T(n,k)*x^k = (x+1)*(x-1)^(n-1), for n >= 1. - Philippe Deléham, Oct 03 2005
T(0,0) = T(1,0) = T(1,1) = 1, T(n,k) = 0 if n < 0 or if n < k, T(n,k) = T(n-1,k-1) - T(n-1,k) for n > 1. - Philippe Deléham, Nov 26 2006
G.f.: (1+2*x)/(1+x-x*y). - R. J. Mathar, Aug 11 2015
From G. C. Greubel, Apr 30 2025: (Start)
T(2*n+1, 2*n+1-k) = T(2*n+1, k) (symmetric odd n rows).
T(2*n, 2*n-k) = (-1)*T(2*n, k) (antisymmetric even n rows).
Sum_{k=0..n} (-1)^k*T(n, k) = A000007(n) (signed row sums).
Sum_{k=0..floor(n/2)} (-1)^k*T(n-k, k) = (-1)^n*A057079(n+2) (signed diagonal sums). (End)

A131085 Triangle T(n,k) (n>=0, 0<=k<=n-1) read by rows, A007318 * A129686.

Original entry on oeis.org

1, 1, 1, 0, 2, 1, -2, 2, 3, 1, -5, 0, 5, 4, 1, -9, -5, 5, 9, 5, 1, -14, -14, 0, 14, 14, 6, 1, -20, -28, -14, 14, 28, 20, 7, 1, -27, -48, -42, 0, 42, 48, 27, 8, 1, -35, -75, -90, -42, 42, 90, 75, 35, 9, 1, -44, -110, -165, -132, 0, 132, 165, 110, 44, 10, 1
Offset: 0

Author

Gary W. Adamson, Jun 14 2007

Keywords

Comments

Row sums = n.
A131085 * A000012 = A074909 starting (1, 2, 1, 3, 3, ...) instead of (1, 1, 2, 1, 3, 3, ...).

Examples

			First few rows of the triangle are:
   1;
   1,  1;
   0,  2, 1;
  -2,  2, 3, 1;
  -5,  0, 5, 4, 1;
  -9, -5, 5, 9, 5, 1;
-14, -14, 0, 14, 14, 6, 1;
-20, -28, -14, 14, 28, 20, 7, 1;
-27, -48, -42, 0, 42, 48, 27, 8, 1;
-35, -75, -90, -42, 42, 90, 75, 35, 9, 1;
   ...
		

Programs

  • PARI
    tabl(nn) = {t007318 = matrix(nn, nn, n, k, binomial(n-1, k-1)); t129686 = matrix(nn, nn, n, k, (k<=n)*((-1)^((n-k)\2)*((k==n) || (-1)*(k==(n-2))))); t131085 = t007318*t129686; for (n = 1, nn, for (k = 1, n, print1(t131085[n, k], ", ");););} \\ Michel Marcus, Feb 12 2014

Formula

Binomial transform of A129686 signed with (1, 1, 1, ...) in the main diagonal and (-1, -1, -1, ...) in the subsubdiagonal.
T(n,m) = T(n-1,m-1) + T(n-1,m). - Yuchun Ji, Dec 17 2018
T(2*k,k-1) = 0 for k > 0. - Yuchun Ji, Dec 20 2018
Comparing this triangle with the Catalan triangle A009766 one can see many similarities. For example, T(k+j,k) = A009766(k+1,j) for j < k+2. - Yuchun Ji, Dec 23 2018 [Edited by N. J. A. Sloane, Feb 11 2019]

Extensions

Missing comma corrected by Naruto Canada, Aug 26 2007
More terms from Michel Marcus, Feb 12 2014
Offset changed by N. J. A. Sloane, Feb 11 2019
Showing 1-8 of 8 results.