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.

Previous Showing 11-20 of 75 results. Next

A078812 Triangle read by rows: T(n, k) = binomial(n+k-1, 2*k-1).

Original entry on oeis.org

1, 2, 1, 3, 4, 1, 4, 10, 6, 1, 5, 20, 21, 8, 1, 6, 35, 56, 36, 10, 1, 7, 56, 126, 120, 55, 12, 1, 8, 84, 252, 330, 220, 78, 14, 1, 9, 120, 462, 792, 715, 364, 105, 16, 1, 10, 165, 792, 1716, 2002, 1365, 560, 136, 18, 1, 11, 220, 1287, 3432, 5005, 4368, 2380, 816, 171, 20, 1
Offset: 0

Views

Author

Michael Somos, Dec 05 2002

Keywords

Comments

Warning: formulas and programs sometimes refer to offset 0 and sometimes to offset 1.
Apart from signs, identical to A053122.
Coefficient array for Morgan-Voyce polynomial B(n,x); see A085478 for references. - Philippe Deléham, Feb 16 2004
T(n,k) is the number of compositions of n having k parts when there are q kinds of part q (q=1,2,...). Example: T(4,2) = 10 because we have (1,3),(1,3'),(1,3"), (3,1),(3',1),(3",1),(2,2),(2,2'),(2',2) and (2',2'). - Emeric Deutsch, Apr 09 2005
T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k (height(alpha) = |Im(alpha)|). - Abdullahi Umar, Oct 02 2008
This sequence is jointly generated with A085478 as a triangular array of coefficients of polynomials v(n,x): initially, u(1,x) = v(1,x) = 1; for n > 1, u(n,x) = u(n-1,x) + x*v(n-1)x and v(n,x) = u(n-1,x) + (x+1)*v(n-1,x). See the Mathematica section. - Clark Kimberling, Feb 25 2012
Concerning Kimberling's recursion relations, see A102426. - Tom Copeland, Jan 19 2016
Subtriangle of the triangle T(n,k), 0 <= k <= n, read by rows, given by (0, 2, -1/2, 1/2, 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, Mar 27 2012
From Wolfdieter Lang, Aug 30 2012: (Start)
With offset [0,0] the triangle with entries R(n,k) = T(n+1,k+1):= binomial(n+k+1, 2*k+1), n >= k >= 0, and zero otherwise, becomes the Riordan lower triangular convolution matrix R = (G(x)/x, G(x)) with G(x):=x/(1-x)^2 (o.g.f. of A000027). This means that the o.g.f. of column number k of R is (G(x)^(k+1))/x. This matrix R is the inverse of the signed Riordan lower triangular matrix A039598, called in a comment there S.
The Riordan matrix with entries R(n,k), just defined, provides the transition matrix between the sequence entry F(4*m*(n+1))/L(2*l), with m >= 0, for n=0,1,... and the sequence entries 5^k*F(2*m)^(2*k+1) for k = 0,1,...,n, with F=A000045 (Fibonacci) and L=A000032 (Lucas). Proof: from the inverse of the signed triangle Riordan matrix S used in a comment on A039598.
For the transition matrix R (T with offset [0,0]) defined above, row n=2: F(12*m) /L(2*m) = 3*5^0*F(2*m)^1 + 4*5^1*F(2*m)^3 + 1*5^2*F(2*m)^5, m >= 0. (End)
From R. Bagula's comment in A053122 (cf. Damianou link p. 10), this array gives the coefficients (mod sign) of the characteristic polynomials for the Cartan matrix of the root system A_n. - Tom Copeland, Oct 11 2014
For 1 <= k <= n, T(n,k) equals the number of (n-1)-length ternary words containing k-1 letters equal 2 and avoiding 01. - Milan Janjic, Dec 20 2016
The infinite sum (Sum_{i >= 0} (T(s+i,1+i) / 2^(s+2*i)) * zeta(s+1+2*i)) = 1 allows any zeta(s+1) to be expressed as a sum of rational multiples of zeta(s+1+2*i) having higher arguments. For example, zeta(3) can be expressed as a sum involving zeta(5), zeta(7), etc. The summation for each s >= 1 uses the s-th diagonal of the triangle. - Robert B Fowler, Feb 23 2022
The convolution triangle of the nonnegative integers. - Peter Luschny, Oct 07 2022

Examples

			Triangle begins, 1 <= k <= n:
                          1
                        2   1
                      3   4   1
                    4  10   6   1
                  5  20  21   8   1
                6  35  56  36  10   1
              7  56 126 120  55  12   1
            8  84 252 330 220  78  14   1
From _Peter Bala_, Feb 11 2025: (Start)
The array factorizes as an infinite product of lower triangular arrays:
  / 1               \    / 1              \ / 1              \ / 1             \
  | 2    1           |   | 2   1          | | 0  1           | | 0  1          |
  | 3    4   1       | = | 3   2   1      | | 0  2   1       | | 0  0  1       | ...
  | 4   10   6   1   |   | 4   3   2  1   | | 0  3   2  1    | | 0  0  2  1    |
  | 5   20  21   8  1|   | 5   4   3  2  1| | 0  4   3  2  1 | | 0  0  3  2  1 |
  |...               |   |...             | |...             | |...            |
Cf. A092276. (End)
		

Crossrefs

This triangle is formed from odd-numbered rows of triangle A011973 read in reverse order.
Row sums give A001906. With signs: A053122.
The column sequences are A000027, A000292, A000389, A000580, A000582, A001288 for k=1..6, resp. For k=7..24 they are A010966..(+2)..A011000 and for k=25..50 they are A017713..(+2)..A017763.

Programs

  • GAP
    Flat(List([0..12], n-> List([0..n], k-> Binomial(n+k+1, 2*k+1) ))); # G. C. Greubel, Aug 01 2019
  • Haskell
    a078812 n k = a078812_tabl !! n !! k
    a078812_row n = a078812_tabl !! n
    a078812_tabl = [1] : [2, 1] : f [1] [2, 1] where
       f us vs = ws : f vs ws where
         ws = zipWith (-) (zipWith (+) ([0] ++ vs) (map (* 2) vs ++ [0]))
                          (us ++ [0, 0])
    -- Reinhard Zumkeller, Dec 16 2013
    
  • Magma
    /* As triangle */ [[Binomial(n+k-1, 2*k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Jun 01 2018
    
  • Maple
    for n from 1 to 11 do seq(binomial(n+k-1,2*k-1),k=1..n) od; # yields sequence in triangular form; Emeric Deutsch, Apr 09 2005
    # Uses function PMatrix from A357368. Adds a row and column above and to the left.
    PMatrix(10, n -> n); # Peter Luschny, Oct 07 2022
  • Mathematica
    (* First program *)
    u[1, x_]:= 1; v[1, x_]:= 1; z = 13;
    u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
    v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x];
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%] (* A085478 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%] (* A078812 *) (* Clark Kimberling, Feb 25 2012 *)
    (* Second program *)
    Table[Binomial[n+k+1, 2*k+1], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Aug 01 2019 *)
  • Maxima
    T(n,m):=sum(binomial(2*k,n-m)*binomial(m+k,k)*(-1)^(n-m+k)*binomial(n+1,m+k+1),k,0,n-m); /* Vladimir Kruchinin, Apr 13 2016 */
    
  • PARI
    {T(n, k) = if( n<0, 0, binomial(n+k-1, 2*k-1))};
    
  • PARI
    {T(n, k) = polcoeff( polcoeff( x*y / (1 - (2 + y) * x + x^2) + x * O(x^n), n), k)};
    
  • Sage
    @cached_function
    def T(k,n):
        if k==n: return 1
        if k==0: return 0
        return sum(i*T(k-1,n-i) for i in (1..n-k+1))
    A078812 = lambda n,k: T(k,n)
    [[A078812(n,k) for k in (1..n)] for n in (1..8)] # Peter Luschny, Mar 12 2016
    
  • Sage
    [[binomial(n+k+1, 2*k+1) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019
    

Formula

G.f.: x*y / (1 - (2 + y)*x + x^2). To get row n, expand this in powers of x then expand the coefficient of x^n in increasing powers of y.
From Philippe Deléham, Feb 16 2004: (Start)
If indexing begins at 0 we have
T(n,k) = (n+k+1)!/((n-k)!*(2k+1))!.
T(n,k) = Sum_{j>=0} T(n-1-j, k-1)*(j+1) with T(n, 0) = n+1, T(n, k) = 0 if n < k.
T(n,k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j) with T(n,k) = 0 if k < 0, T(0, 0)=1 and T(0, k) = 0 for k > 0.
G.f. for the column k: Sum_{n>=0} T(n, k)*x^n = (x^k)/(1-x)^(2k+2).
Row sums: Sum_{k>=0} T(n, k) = A001906(n+1). (End)
Antidiagonal sums are A000079(n) = Sum_{k=0..floor(n/2)} binomial(n+k+1, n-k). - Paul Barry, Jun 21 2004
Riordan array (1/(1-x)^2, x/(1-x)^2). - Paul Barry, Oct 22 2006
T(0,0) = 1, T(n,k) = 0 if k < 0 or if k > n, T(n,k) = T(n-1,k-1) + 2*T(n-1,k) - T(n-2,k). - Philippe Deléham, Jan 26 2010
For another version see A128908. - Philippe Deléham, Mar 27 2012
T(n,m) = Sum_{k=0..n-m} (binomial(2*k,n-m)*binomial(m+k,k)*(-1)^(n-m+k)* binomial(n+1,m+k+1)). - Vladimir Kruchinin, Apr 13 2016
T(n, k) = T(n-1, k) + (T(n-1, k-1) + T(n-2, k-1) + T(n-3, k-1) + ...) for k >= 2 with T(n, 1) = n. - Peter Bala, Feb 11 2025
From Peter Bala, May 04 2025: (Start)
With the column offset starting at 0, the n-th row polynomial B(n, x) = 1/sqrt(x + 4) * Chebyshev_U(2*n+1, (1/2)*sqrt(x + 4)) = (-1)^n * Chebyshev_U(n, -(1/2)*(x + 2)).
B(n, x) / Product_{k = 1..2*n} (1 + 1/B(k, x)) = b(n, x), the n-th row polynomial of A085478. (End)

Extensions

Edited by N. J. A. Sloane, Apr 28 2008

A087127 This table shows the coefficients of combinatorial formulas needed for generating the sequential sums of p-th powers of triangular numbers. The p-th row (p>=1) contains a(i,p) for i=1 to 2*p-1, where a(i,p) satisfies Sum_{i=1..n} C(i+1,2)^p = 3 * C(n+2,3) * Sum_{i=1..2*p-1} a(i,p) * C(n-1,i-1)/(i+2).

Original entry on oeis.org

1, 1, 2, 1, 1, 8, 19, 18, 6, 1, 26, 163, 432, 564, 360, 90, 1, 80, 1135, 6354, 18078, 28800, 26100, 12600, 2520, 1, 242, 7291, 77400, 405060, 1210680, 2211570, 2520000, 1751400, 680400, 113400, 1, 728, 45199, 862218, 7667646, 38350080, 118848420
Offset: 1

Views

Author

André F. Labossière, Aug 11 2003

Keywords

Comments

From Peter Bala, Mar 08 2018: (Start)
The table entries T(n,k) are the coefficients when expressing the polynomial C(x+2,2)^p of degree 2*p in terms of falling factorials: C(x+2,2)^p = Sum_{k = 0..2*p} T(p,k)*C(x,k). It follows that Sum_{i = 0..n-1} C(i+2,2)^p = Sum_{k = 0..2*p} T(p,k)*C(n,k+1).
The sum of the p-th powers of the triangular numbers is also given by Sum_{i = 0..n-1} C(i+2,2)^p = Sum_{k = 2..2*p} A122193(p,k)*C(n+2,k+1) for p >= 1. (End)

Examples

			Row 3 contains 1,8,19,18,6, so Sum_{i=1..n} C(i+1,2)^3 = (n+2) * C(n+1,2) * [ a(1,3)/3 + a(2,3)*C(n-1,1)/4 + a(3,3)*C(n-1,2)/5 + a(4,3)*C(n-1,3)/6 + a(5,3)*C(n-1,4)/7 ] = [ (n+2)*(n+1)*n/2 ] * [ 1/3 + (8/4)*C(n-1,1) + (19/5)*C(n-1,2) + (18/6)*C(n-1,3) + (6/7)*C(n-1,4). Cf. A085438 for more details.
From _Peter Bala_, Mar 08 2018: (Start)
Table begins
n=0 |1
n=1 |1   2     1
n=2 |1   8    19    18      6
n=3 |1  26   163   432    564    360     90
n=4 |1  80  1135  6354  18078  28800  26100  12600  2520
...
Row 2: C(i+2,2)^2 = C(i,0) + 8*C(i,1) + 19*C(i,2) + 18*C(i,3) + 6*C(i,4). Hence, Sum_{i = 0..n-1} C(i+2,2)^2 =  C(n,1) + 8*C(n,2) + 19*C(n,3) + 18*C(n,4) + 6*C(n,5). (End)
		

Crossrefs

Programs

  • GAP
    Flat(List([0..6],n->List([0..2*n],k->Sum([0..k],i->(-1)^(k-i)*Binomial(k,i)*Binomial(i+2,2)^n)))); # Muniru A Asiru, Mar 22 2018
  • Maple
    seq(seq(add( (-1)^(k-i)*binomial(k,i)*binomial(i+2,2)^n, i = 0..k), k = 0..2*n), n = 0..8); # Peter Bala, Mar 08 2018
  • Mathematica
    a[i_, p_] := Sum[Binomial[i - 1, 2*k - 2]*Binomial[i - 2*k + 3, i - 2*k + 1]^(p - 1) - Binomial[i - 1, 2*k - 1]*Binomial[i - 2*k + 2, i - 2*k]^(p - 1), {k, 1, (2*i + 1 + (-1)^(i - 1))/4}]; Table[If[p == 1, 1, a[i, p]], {p, 1, 10}, {i, 1, 2*p - 1}]//Flatten (* G. C. Greubel, Nov 23 2017 *)
    a[i_,p_]:=(-1)^i HypergeometricPFQ[ConstantArray[3,p]~Join~{2-i},ConstantArray[1,p],1];Table[a[i,p],{p,0,10},{i,2,2 p+2}]//Flatten (* Jonathan Burns, Mar 20 2018 *)
  • PARI
    {a(i, p) = sum(k=1, (2*i + 1 + (-1)^(i - 1))/4, binomial(i - 1, 2*k - 2)*binomial(i - 2*k + 3, i - 2*k + 1)^(p - 1) - binomial(i - 1, 2*k - 1)*binomial(i - 2*k + 2, i - 2*k)^(p - 1))}; for(p=1,8, for(i=1, 2*p-1, print1(if(p==1,1,a(i,p)), ", "))) \\ G. C. Greubel, Nov 23 2017
    

Formula

a(1, p) = 1, a(2, p) = 3^(p-1)-1, a(3, p) = 3^(p-1)*[2^(p-1)-2]+1, ..., a(2*p-3, p) = [ (6*p^4-20*p^3+21*p^2-7*p)*(2*p-4)! ]/[3*2^(p-1)], a(2*p-2, p) = [ (p^2-p)*(2*p-3)! ]/2^(p-2), a(2*p-1, p) = [ (p-1)*(2*p-3)! ]/2^(p-2).
a(i, p) = Sum_{k=1..[2*i+1+(-1)^(i-1)]/4} [ C(i-1, 2*k-2)*C(i-2*k+3, i-2*k+1)^(p-1) -C(i-1, 2*k-1)*C(i-2*k+2, i-2*k)^(p-1) ]
From Peter Bala, Mar 08 2018: (Start)
The following remarks assume row and column indices start at 0.
T(n,k) = Sum_{i = 0..k} (-1)^(k-i)*C(k,i)*C(i+2,2)^n. Equivalently, let v_n denote the sequence (1, 3^n, 6^n, 10^n, ...) regarded as an infinite column vector, where 1, 3, 6, 10, ... is the sequence of triangular numbers A000217. Then the n-th row of this table is determined by the matrix product P^(-1)*v_n, where P denotes Pascal's triangle A007318. Cf. A122193.
T(n+1,k) = C(k+2,2)*T(n,k) + 2*C(k+1,2)*T(n,k-1) + C(k,2)*T(n,k-2), with boundary conditions T(n,0) = 1 for all n and T(n,k) = 0 for k > 2*n.
Let R(n,x) denote the n-th row polynomial.
R(n+1,x) = 1/2!*(1 + x)^2*(d/dx)^2 (x^2*R(n,x)).
R(n,x) = Sum_{i >= 0} binomial(i+2,2)^n*x^i/(1 + x)^(i+1).
R(n,x) = (1 + x)^2 o (1 + x)^2 o ... o (1 + x)^2 (n factors), where o denotes the black diamond product of power series defined in Dukes and White. Note the polynomial x^2 o ... o x^2 (n factors) is the n-th row polynomial of A122193.
x^2*R(n,x) = (1 + x)^2 * the n-th row polynomial of A122193 (End)

Extensions

Edited by Dean Hickerson, Aug 16 2003

A053122 Triangle of coefficients of Chebyshev's S(n,x-2) = U(n,x/2-1) polynomials (exponents of x in increasing order).

Original entry on oeis.org

1, -2, 1, 3, -4, 1, -4, 10, -6, 1, 5, -20, 21, -8, 1, -6, 35, -56, 36, -10, 1, 7, -56, 126, -120, 55, -12, 1, -8, 84, -252, 330, -220, 78, -14, 1, 9, -120, 462, -792, 715, -364, 105, -16, 1, -10, 165, -792, 1716, -2002, 1365, -560, 136, -18, 1, 11, -220, 1287, -3432, 5005, -4368, 2380, -816, 171, -20
Offset: 0

Views

Author

Keywords

Comments

Apart from signs, identical to A078812.
Another version with row-leading 0's and differing signs is given by A285072.
G.f. for row polynomials S(n,x-2) (signed triangle): 1/(1+(2-x)*z+z^2). Unsigned triangle |a(n,m)| has g.f. 1/(1-(2+x)*z+z^2) for row polynomials.
Row sums (signed triangle) A049347(n) (periodic(1,-1,0)). Row sums (unsigned triangle) A001906(n+1)=F(2*(n+1)) (even-indexed Fibonacci).
In the language of Shapiro et al. (see A053121 for the reference) such a lower triangular (ordinary) convolution array, considered as a matrix, belongs to the Bell-subgroup of the Riordan-group.
The (unsigned) column sequences are A000027, A000292, A000389, A000580, A000582, A001288 for m=0..5, resp. For m=6..23 they are A010966..(+2)..A011000 and for m=24..49 they are A017713..(+2)..A017763.
Riordan array (1/(1+x)^2,x/(1+x)^2). Inverse array is A039598. Diagonal sums have g.f. 1/(1+x^2). - Paul Barry, Mar 17 2005. Corrected by Wolfdieter Lang, Nov 13 2012.
Unsigned version is in A078812. - Philippe Deléham, Nov 05 2006
Also row n gives (except for an overall sign) coefficients of characteristic polynomial of the Cartan matrix for the root system A_n. - Roger L. Bagula, May 23 2007
From Wolfdieter Lang, Nov 13 2012: (Start)
The A-sequence for this Riordan triangle is A115141, and the Z-sequence is A115141(n+1), n>=0. For A- and Z-sequences for Riordan matrices see the W. Lang link under A006232 with details and references.
S(n,x^2-2) = sum(r(j,x^2),j=0..n) with Chebyshev's S-polynomials and r(j,x^2) := R(2*j+1,x)/x, where R(n,x) are the monic integer Chebyshv T-polynomials with coefficients given in A127672. Proof from comparing the o.g.f. of the partial sum of the r(j,x^2) polynomials (see a comment on the signed Riordan triangle A111125) with the present Riordan type o.g.f. for the row polynomials with x -> x^2. (End)
S(n,x^2-2) = S(2*n+1,x)/x, n >= 0, from the odd part of the bisection of the o.g.f. - Wolfdieter Lang, Dec 17 2012
For a relation to a generator for the Narayana numbers A001263, see A119900, whose columns are unsigned shifted rows (or antidiagonals) of this array, referring to the tables in the example sections. - Tom Copeland, Oct 29 2014
The unsigned rows of this array are alternating rows of a mirrored A011973 and alternating shifted rows of A030528 for the Fibonacci polynomials. - Tom Copeland, Nov 04 2014
Boas-Buck type recurrence for column k >= 0 (see Aug 10 2017 comment in A046521 with references): a(n, m) = (2*(m + 1)/(n - m))*Sum_{k = m..n-1} (-1)^(n-k)*a(k, m), with input a(n, n) = 1, and a(n,k) = 0 for n < k. - Wolfdieter Lang, Jun 03 2020
Row n gives the characteristic polynomial of the (n X n)-matrix M where M[i,j] = 2 if i = j, -1 if |i-j| = 1 and 0 otherwise. The matrix M is positive definite and has 2-condition number (cot(Pi/(2*n+2)))^2. - Jianing Song, Jun 21 2022
Also the convolution triangle of (-1)^(n+1)*n. - Peter Luschny, Oct 07 2022

Examples

			The triangle a(n,m) begins:
n\m   0    1    2     3     4     5     6    7    8  9
0:    1
1:   -2    1
2:    3   -4    1
3:   -4   10   -6     1
4:    5  -20   21    -8     1
5:   -6   35  -56    36   -10     1
6:    7  -56  126  -120    55   -12     1
7:   -8   84 -252   330  -220    78   -14    1
8:    9 -120  462  -792   715  -364   105  -16    1
9:  -10  165 -792  1716 -2002  1365  -560  136  -18  1
... Reformatted and extended by _Wolfdieter Lang_, Nov 13 2012
E.g., fourth row (n=3) {-4,10,-6,1} corresponds to the polynomial S(3,x-2) = -4+10*x-6*x^2+x^3.
From _Wolfdieter Lang_, Nov 13 2012: (Start)
Recurrence: a(5,1) = 35 = 1*5 + (-2)*(-20) -1*(10).
Recurrence from Z-sequence [-2,-1,-2,-5,...]: a(5,0) = -6 = (-2)*5 + (-1)*(-20) + (-2)*21 + (-5)*(-8) + (-14)*1.
Recurrence from A-sequence [1,-2,-1,-2,-5,...]: a(5,1) = 35 = 1*5  + (-2)*(-20) + (-1)*21 + (-2)*(-8) + (-5)*1.
(End)
E.g., the fourth row (n=3) {-4,10,-6,1} corresponds also to the polynomial S(7,x)/x = -4 + 10*x^2 - 6*x^4 + x^6. - _Wolfdieter Lang_, Dec 17 2012
Boas-Buck type recurrence: -56 = a(5, 2) = 2*(-1*1 + 1*(-6) - 1*21) = -2*28 = -56. - _Wolfdieter Lang_, Jun 03 2020
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
  • Theodore J. Rivlin, Chebyshev polynomials: from approximation theory to algebra and number theory, 2. ed., Wiley, New York, 1990.
  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 62.
  • Sigurdur Helgasson, Differential Geometry, Lie Groups and Symmetric Spaces, Graduate Studies in Mathematics, volume 34. A. M. S.: ISBN 0-8218-2848-7, 1978, p. 463.

Crossrefs

Cf. A285072 (version with row-leading 0's and differing signs). - Eric W. Weisstein, Apr 09 2017

Programs

  • Maple
    seq(seq((-1)^(n+m)*binomial(n+m+1,2*m+1),m=0..n),n=0..10); # Robert Israel, Oct 15 2014
    # Uses function PMatrix from A357368. Adds a row above and a column to the left.
    PMatrix(10, n -> -(-1)^n*n); # Peter Luschny, Oct 07 2022
  • Mathematica
    T[n_, m_, d_] := If[ n == m, 2, If[n == m - 1 || n == m + 1, -1, 0]]; M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}]; a = Join[M[1], Table[CoefficientList[Det[M[d] - x*IdentityMatrix[d]], x], {d, 1, 10}]]; Flatten[a] (* Roger L. Bagula, May 23 2007 *)
    (* Alternative code for the matrices from MathWorld: *)
    sln[n_] := 2IdentityMatrix[n] - PadLeft[PadRight[IdentityMatrix[n - 1], {n, n - 1}], {n, n}] - PadLeft[PadRight[IdentityMatrix[n - 1], {n - 1, n}], {n, n}] (* Roger L. Bagula, May 23 2007 *)
  • Sage
    @CachedFunction
    def A053122(n,k):
        if n< 0: return 0
        if n==0: return 1 if k == 0 else 0
        return A053122(n-1,k-1)-A053122(n-2,k)-2*A053122(n-1,k)
    for n in (0..9): [A053122(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012

Formula

a(n, m) := 0 if n
a(n, m) = -2*a(n-1, m) + a(n-1, m-1) - a(n-2, m), a(n, -1) := 0 =: a(-1, m), a(0, 0)=1, a(n, m) := 0 if n
O.g.f. for m-th column (signed triangle): ((x/(1+x)^2)^m)/(1+x)^2.
From Jianing Song, Jun 21 2022: (Start)
T(n,k) = [x^k]f_n(x), where f_{-1}(x) = 0, f_0(x) = 1, f_n(x) = (x-2)*f_{n-1}(x) - f_{n-2}(x) for n >= 2.
f_n(x) = (((x-2+sqrt(x^2-4*x))/2)^(n+1) - ((x-2-sqrt(x^2-4*x))/2)^(n+1))/sqrt(x^2-4x).
The roots of f_n(x) are 2 + 2*cos(k*Pi/(n+1)) = 4*(cos(k*Pi/(2*n+2)))^2 for 1 <= k <= n. (End)

A060540 Square array read by antidiagonals downwards: T(n,k) = (n*k)!/(k!^n*n!), (n>=1, k>=1), the number of ways of dividing nk labeled items into n unlabeled boxes with k items in each box.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 10, 15, 1, 1, 35, 280, 105, 1, 1, 126, 5775, 15400, 945, 1, 1, 462, 126126, 2627625, 1401400, 10395, 1, 1, 1716, 2858856, 488864376, 2546168625, 190590400, 135135, 1, 1, 6435, 66512160, 96197645544, 5194672859376, 4509264634875, 36212176000, 2027025, 1
Offset: 1

Author

Henry Bottomley, Apr 02 2001

Keywords

Comments

The Copeland link gives the associations of this entry with the operator calculus of Appell Sheffer polynomials, the combinatorics of simple set partitions encoded in the Faa di Bruno formula for composition of analytic functions (formal Taylor series), the Pascal matrix, and the geometry of the n-dimensional simplices (hypertriangles, or hypertetrahedra). These, in turn, are related to simple instances of the application of the exponential formula / principle / schema giving the number of not-necessarily-connected objects composed from an ensemble of connected objects. - Tom Copeland, Jun 09 2021

Examples

			Array begins:
  1,   1,       1,          1,             1,                 1, ...
  1,   3,      10,         35,           126,               462, ...
  1,  15,     280,       5775,        126126,           2858856, ...
  1, 105,   15400,    2627625,     488864376,       96197645544, ...
  1, 945, 1401400, 2546168625, 5194672859376, 11423951396577720, ...
  ...
		

Crossrefs

Main diagonal is A057599.
Related to A057599, see also A096126 and A246048.
Cf. A060358, A361948 (includes row/col 0).
Cf. A000217, A000292, A000332, A000389, A000579, A000580, A007318, A036040, A099174, A133314, A132440, A135278 (associations in Copeland link).

Programs

  • Mathematica
    T[n_, k_] := (n*k)!/(k!^n*n!);
    Table[T[n-k+1, k], {n, 1, 10}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Jun 29 2018 *)
  • PARI
    { i=0; for (m=1, 20, for (n=1, m, k=m - n + 1; write("b060540.txt", i++, " ", (n*k)!/(k!^n*n!))); ) } \\ Harry J. Smith, Jul 06 2009

Formula

T(n,k) = (n*k)!/(k!^n*n!) = T(n-1,k)*A060543(n,k) = A060538(n,k)/k!.
T(n,k) = Product_{j=2..n} binomial(j*k-1,k-1). - M. F. Hasler, Aug 22 2014

Extensions

Definition reworded by M. F. Hasler, Aug 23 2014

A086020 a(n) = Sum_(i=1..n) binomial(i+2,3)^2 [ Sequential sums of the tetragonal numbers or "tetras" (pyramidal, square) raised to power 2 (drawn from the 4th diagonal - left or right - of Pascal's Triangle) ].

Original entry on oeis.org

1, 17, 117, 517, 1742, 4878, 11934, 26334, 53559, 101959, 183755, 316251, 523276, 836876, 1299276, 1965132, 2904093, 4203693, 5972593, 8344193, 11480634, 15577210, 20867210, 27627210, 36182835, 46915011, 60266727, 76750327
Offset: 1

Author

André F. Labossière, Jul 17 2003

Keywords

Comments

Kekulé numbers for certain benzenoids (see the Cyvin-Gutman reference, p. 243; expression in (13.26) yields same sequence with offset 0). - Emeric Deutsch, Aug 02 2005
Partial sums of A001249. - R. J. Mathar, Aug 19 2008

Examples

			a(8) = Sum_{i=1..8} binomial(i+2,3)^2 = (20*(8^7) + 210*(8^6) + 854*(8^5) + 1680*(8^4) + 1610*(8^3) + 630*(8^2) + 36*8)/7! = 26334.
		

Programs

  • Magma
    [n*(n+1)*(n+2)*(n+3)*(2*n+3)*(5*n^2+15*n+1)/2520: n in [1..30]]; // G. C. Greubel, Nov 22 2017
  • Maple
    a:=n->n*(n+1)*(n+2)*(n+3)*(2*n+3)*(5*n^2+15*n+1)/2520: seq(a(n),n=1..31); # Emeric Deutsch
  • Mathematica
    Accumulate[Binomial[Range[30]+2,3]^2]  (* Harvey P. Dale, Mar 24 2011 *)
    LinearRecurrence[{8,-28,56,-70,56,-28,8,-1},{1,17,117,517,1742,4878, 11934, 26334},30] (* Harvey P. Dale, Aug 17 2014 *)
  • PARI
    a(n)=n*(n+1)*(n+2)*(n+3)*(2*n+3)*(5*n^2+15*n+1)/2520 \\ Charles R Greathouse IV, May 18 2015
    

Formula

a(n) = Sum_(i=1..n) binomial(i+2, 3)^2.
a(n) = ( C(n+3, 4)/35 )*( 35 + 84*C(n-1, 1) + 70*C(n-1, 2) + 20*C(n-1, 3) ).
a(n) = n*(n+1)*(n+2)*(n+3)*(2*n+3)(5*n^2 + 15*n + 1)/2520. - Emeric Deutsch, Aug 02 2005
O.g.f: x*(1+x)*(1 + 8*x + x^2)/(1-x)^8. - R. J. Mathar, Aug 19 2008

A085438 a(n) = Sum_{i=1..n} binomial(i+1,2)^3.

Original entry on oeis.org

1, 28, 244, 1244, 4619, 13880, 35832, 82488, 173613, 339988, 627484, 1102036, 1855607, 3013232, 4741232, 7256688, 10838265, 15838476, 22697476, 31958476, 44284867, 60479144, 81503720, 108503720, 142831845, 186075396, 240085548, 307008964, 389321839
Offset: 1

Author

André F. Labossière, Jun 30 2003

Keywords

Examples

			a(10) = (90*(10^7)+630*(10^6)+1638*(10^5)+1890*(10^4)+840*(10^3)-48*(10))/5040 = 339988.
		

References

  • Elisabeth Busser and Gilles Cohen, Neuro-Logies - "Chercher, jouer, trouver", La Recherche, April 1999, No. 319, page 97.

Programs

  • Magma
    [(90*n^7 +630*n^6 +1638*n^5 +1890*n^4+ 840*n^3 -48*n)/ Factorial(7): n in [1..30]]; // G. C. Greubel, Nov 22 2017
  • Mathematica
    Table[(90*n^7 + 630*n^6 + 1638*n^5 + 1890*n^4 + 840*n^3 - 48*n)/7!, {n, 1, 50}] (* G. C. Greubel, Nov 22 2017 *)
  • PARI
    Vec(x*(x^4+20*x^3+48*x^2+20*x+1)/(x-1)^8 + O(x^100)) \\ Colin Barker, May 02 2014
    
  • PARI
    a(n) = sum(i=1, n, binomial(i+1, 2)^3); \\ Michel Marcus, Nov 22 2017
    

Formula

a(n) = (90*n^7 +630*n^6 +1638*n^5 +1890*n^4+ 840*n^3 -48*n)/7!.
a(n) = (C(n+2, 3)/35)*(35 +210*C(n-1, 1) +399*C(n-1, 2) +315*C(n-1, 3) +90*C(n-1, 4)).
G.f.: x*(x^4+20*x^3+48*x^2+20*x+1) / (x-1)^8. - Colin Barker, May 02 2014

Extensions

More terms from Colin Barker, May 02 2014
Formula and example edited by Colin Barker, May 02 2014

A034867 Triangle of odd-numbered terms in rows of Pascal's triangle.

Original entry on oeis.org

1, 2, 3, 1, 4, 4, 5, 10, 1, 6, 20, 6, 7, 35, 21, 1, 8, 56, 56, 8, 9, 84, 126, 36, 1, 10, 120, 252, 120, 10, 11, 165, 462, 330, 55, 1, 12, 220, 792, 792, 220, 12, 13, 286, 1287, 1716, 715, 78, 1, 14, 364, 2002, 3432, 2002, 364, 14, 15, 455, 3003, 6435, 5005, 1365, 105, 1
Offset: 0

Keywords

Comments

Also triangle of numbers of n-sequences of 0,1 with k subsequences of consecutive 01 because this number is C(n+1,2*k+1). - Roger Cuculiere (cuculier(AT)imaginet.fr), Nov 16 2002
From Gary W. Adamson, Oct 17 2008: (Start)
Received from Herb Conn:
Let T = tan x, then
tan x = T
tan 2x = 2T / (1 - T^2)
tan 3x = (3T - T^3) / (1 - 3T^2)
tan 4x = (4T - 4T^3) / (1 - 6T^2 + T^4)
tan 5x = (5T - 10T^3 + T^5) / (1 - 10T^2 + 5T^4)
tan 6x = (6T - 20T^3 + 6T^5) / (1 - 15T^2 + 15T^4 - T^6)
tan 7x = (7T - 35T^3 + 21T^5 - T^7) / (1 - 21T^2 + 35T^4 - 7T^6)
tan 8x = (8T - 56T^3 + 56T^5 - 8T^7) / (1 - 28T^2 + 70T^4 - 28T^6 + T^8)
tan 9x = (9T - 84T^3 + 126T^5 - 36T^7 + T^9) / (1 - 36 T^2 + 126T^4 - 84T^6 + 9T^8)
... To get the next one in the series, (tan 10x), for the numerator add:
9....84....126....36....1 previous numerator +
1....36....126....84....9 previous denominator =
10..120....252...120...10 = new numerator
For the denominator add:
......9.....84...126...36...1 = previous numerator +
1....36....126....84....9.... = previous denominator =
1....45....210...210...45...1 = new denominator
...where numerators = A034867, denominators = A034839
(End)
Column k is the sum of columns 2k and 2k+1 of A007318. - Philippe Deléham, Nov 12 2008
Triangle, with zeros omitted, given by (2, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (0, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Dec 12 2011
The row polynomials N(n,x) = Sum_{k=0..floor((n-1)/2)} T(n-1,k)*x^k, and D(n,x) = Sum_{k=0..floor(n/2)} A034839(n,k)*x^k, n >= 1, satisfy the recurrences N(n,x) = D(n-1,x) + N(n-1,x), D(n,x) = D(n-1,x) + x*N(n-1,x), with inputs N(1,x) = 1 = D(1,x). This is due to the Pascal triangle A007318 recurrence. Q(n,x) := tan(n*x)/tan(x) satisfies the recurrence Q(n,x) = (1 + Q(n-1,x))/(1 - v(x)*Q(n-1,x)) with input Q(1,x) = 1 and v = v(x) := (tan(x))^2. This recurrence is obtained from the addition theorem for tan(n*x) using n = 1 + (n-1). Therefore Q(n,x) = N(n,-v(x))/D(n,-v(x)). This proves the Gary W. Adamson contribution from above. See also A220673. This calculation was motivated by an e-mail of Thomas Olsen. The Oliver/Prodinger and Ma references resort to HAKEM Al Memo 239, Item 16, for the tan(n*x) formula in terms of tan(x). - Wolfdieter Lang, Jan 17 2013
The infinitesimal generator (infinigen) for the Narayana polynomials A090181/A001263 can be formed from the row polynomials P(n,y) of this entry. The resulting matrix is an instance of a matrix representation of the analytic infinigens presented in A145271 for general sets of binomial Sheffer polynomials and in A001263 and A119900 specifically for the Narayana polynomials. Given the column vector of row polynomials V = (1, P(1,x) = 2x, P(2,y) = 3x + x^2, P(3,y) = 4x + 4x^2, ...), form the lower triangular matrix M(n,k) = V(n-k,n-k), i.e., diagonally multiply the matrix with all ones on the diagonal and below by the components of V. Form the matrix MD by multiplying A132440^Transpose = A218272 = D (representing derivation of o.g.f.s) by M, i.e., MD = M*D. The non-vanishing component of the first row of (MD)^n * V / (n+1)! is the n-th Narayana polynomial. - Tom Copeland, Dec 09 2015
The diagonals of this entry are A078812 (also shifted A128908 and unsigned A053122, which are embedded in A030528, A102426, A098925, A109466, A092865). Equivalently, the antidiagonals of A078812 are the rows of A034867. - Tom Copeland, Dec 12 2015
Binomial(n,2k+1) is also the number of permutations avoiding both 132 and 213 with k peaks, i.e., positions with w[i]w[i+2]. - Lara Pudwell, Dec 19 2018
Binomial(n,2k+1) is also the number of permutations avoiding both 123 and 132 with k peaks, i.e., positions with w[i]w[i+2]. - Lara Pudwell, Dec 19 2018
The row polynomial P(n, x) = Sum_{0..floor(n/2)} T(n, k)*x^k appears as numerator polynomial of the diagonal sequence m of triangle A104698 as follows. G(m, x) = P(m, x^2)/(1 - x)^(m+1), for m >= 0. - Wolfdieter Lang, May 14 2025
Number of acyclic orientations of the path graph on n+1 vertices, with k-1 sinks. - Per W. Alexandersson, Aug 15 2025

Examples

			Triangle T starts:
  n\k   0   1   2   3   4  5 ...   ----------------------------------------
0:    1
1:    2
2:    3   1
3:    4   4
4:    5  10   1
5:    6  20   6
6:    7  35  21   1
7:    8  56  56   8
8:    9  84 126  36   1
9:   10 120 252 120  10
 10:   11 165 462 330  55  1
 11:   12 220 792 792 220 12
... ... reformatted and extended by - _Wolfdieter Lang_, May 14 2025
		

References

  • A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 136.

Crossrefs

From Wolfdieter Lang, May 14 2025:(Start)
Row length A008619. Row sums A000079. Alternating row sums A009545(n+1).
Column sequences (with certain offsets): A000027, A000292, A000389, A000580, A000582, A001288, ... (End)

Programs

  • Magma
    /* as a triangle */ [[Binomial(n+1,2*k+1): k in [0..Floor(n/2)]]: n in [0..20]]; // G. C. Greubel, Mar 06 2018
  • Maple
    seq(seq(binomial(n+1,2*k+1), k=0..floor(n/2)), n=0..14); # Emeric Deutsch, Apr 01 2005
  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 12;
    u[n_, x_] := u[n - 1, x] + x*v[n - 1, x]
    v[n_, x_] := u[n - 1, x] + v[n - 1, x]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]  (* A034839 as a triangle *)
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]  (* A034867 as a triangle *)
    (* Clark Kimberling, Feb 18 2012 *)
    Table[Binomial[n+1, 2*k+1], {n,0,20}, {k,0,Floor[n/2]}]//Flatten (* G. C. Greubel, Mar 06 2018 *)
  • PARI
    for(n=0,20, for(k=0,floor(n/2), print1(binomial(n+1,2*k+1), ", "))) \\ G. C. Greubel, Mar 06 2018
    

Formula

T(n,k) = C(n+1,2k+1) = Sum_{i=k..n-k} C(i,k) * C(n-i,k).
E.g.f.: 1+(exp(x)*sinh(x*sqrt(y)))/sqrt(y). - Vladeta Jovovic, Mar 20 2005
G.f.: 1/((1-z)^2-t*z^2). - Emeric Deutsch, Apr 01 2005
T(n,k) = Sum_{j = 0..n} A034839(j,k). - Philippe Deléham, May 18 2005
Pell(n+1) = A000129(n+1) = Sum_{k=0..n} T(n,k) * 2^k = (1/n!) Sum_{k=0..n} A131980(n,k) * 2^k. - Tom Copeland, Nov 30 2007
T(n,k) = A007318(n,2k) + A007318(n,2k+1). - Philippe Deléham, Nov 12 2008
O.g.f for column k, k>=0: (1/(1-x)^2)*(x/(1-x))^(2*k). See the G.f. of this array given above by Emeric Deutsch. - Wolfdieter Lang, Jan 18 2013
T(n,k) = (x^(2*k+1))*((1+x)^n-(1-x)^n)/2. - L. Edson Jeffery, Jan 15 2014

Extensions

More terms from Emeric Deutsch, Apr 01 2005

A085442 a(n) = Sum_{i=1..n} binomial(i+1,2)^7.

Original entry on oeis.org

1, 2188, 282124, 10282124, 181141499, 1982230040, 15475158552, 93839322648, 467508775773, 1989944010148, 7445104711204, 25010673566116, 76686775501847, 217396817767472, 575714897767472, 1436257466526768, 3398894618986905, 7674255436599996, 16612972826599996
Offset: 1

Author

André F. Labossière, Jul 07 2003

Keywords

Programs

  • Magma
    [(1/823680) *n*(n+1)*(n+2)*(429*n^12 +5148*n^11 +24123*n^10 +52470*n^9 +43047*n^8 -8856*n^7 +4109*n^6 +50430*n^5 -18796*n^4 -44472*n^3 +26864*n^2 +8352*n -5568): n in [1..30]]; // G. C. Greubel, Nov 22 2017
  • Mathematica
    Table[Sum[Binomial[k+1,2]^7, {k,1,n}], {n,1,30}] (* G. C. Greubel, Nov 22 2017 *)
    LinearRecurrence[{16,-120,560,-1820,4368,-8008,11440,-12870,11440,-8008,4368,-1820,560,-120,16,-1},{1,2188,282124,10282124,181141499,1982230040,15475158552,93839322648,467508775773,1989944010148,7445104711204,25010673566116,76686775501847,217396817767472,575714897767472,1436257466526768},20] (* Harvey P. Dale, May 11 2022 *)
  • PARI
    for(n=1,30, print1(sum(k=1,n, binomial(k+1,2)^7), ", ")) \\ G. C. Greubel, Nov 22 2017
    

Formula

a(n) = (1/823680) *n*(n+1)*(n+2)*(429*n^12 +5148*n^11 +24123*n^10 +52470*n^9 +43047*n^8 -8856*n^7 +4109*n^6 +50430*n^5 -18796*n^4 -44472*n^3 +26864*n^2 +8352*n -5568). - Vladeta Jovovic, Jul 07 2003
G.f.: x*(x^12 +2172*x^11 +247236*x^10 +6030140*x^9 +49258935*x^8 +163809288*x^7 +242384856*x^6 +163809288*x^5 +49258935*x^4 +6030140*x^3 +247236*x^2 +2172*x+ 1) / (x -1)^16. - Colin Barker, May 02 2014

A254640 Third partial sums of sixth powers (A001014).

Original entry on oeis.org

1, 67, 927, 6677, 32942, 126378, 404634, 1129854, 2833479, 6515509, 13947505, 28115451, 53846156, 98669156, 173975076, 296541132, 490504893, 789878583, 1241708083, 1909993393, 2880500634, 4266609710, 6216356510, 8920844010, 12624212835, 17635378761
Offset: 1

Author

Luciano Ancora, Feb 04 2015

Keywords

Comments

This is one of a sequence of arrays that are the convolutions of the zero-padded sequences binomial(2n-1+k,k) with the Eulerian polynomials E(n,x) of A008292, represented by E(n,x) (1-x)^(-2n), which generate increasing partial sums of powers of integers:
n= 2) (1 + 4*x + x^2)/(1-x)^4 is the o.g.f. of A000578, the convolution of (1,4,1) with A000292, giving the powers of m^3.
n= 3) (1 + 11*x + 11*x^2 + x^3)/(1-x)^6 is the o.g.f. of A000538, convolution of (1,11,11,1) with A000389, giving the partial sums of m^4.
n= 4) (1 + 26*x + 66*x^2 + 26*x^3 + x^4)/(1-x)^8, the o.g.f. of A101092, convolution of (1,26,66,26,1) with A000580, the second partial sums of m^5
n= 5) (1 + 57*x + 302*x^2 + 302*x^3 + 57*x^4 + x^5)/(1-x)^10, the o.g.f. of A254460, convolution of (1,57,302,302,57,1) with A000582, giving the third partial sums of m^6. - Tom Copeland, Dec 07 2015

Programs

  • GAP
    List([1..30], n-> Binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210); # G. C. Greubel, Aug 28 2019
    
  • Magma
    [n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2-30*n+35*n^2+30*n^3+ 5*n^4)/5040: n in [1..30]]; // Vincenzo Librandi, Feb 05 2015
    
  • Maple
    seq(binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210, n=1..30); # G. C. Greubel, Aug 28 2019
  • Mathematica
    Table[n(1+n)(2+n)(3+n)(3+2n)(2 -30n +35n^2 +30n^3 +5n^4)/5040, {n, 30}] (* or *) CoefficientList[Series[(x+1)(x^4 +56x^3 +246x^2 +56x +1)/(x - 1)^10, {x, 0, 30}], x] (* Vincenzo Librandi, Feb 05 2015 *)
  • PARI
    vector(30, n, n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2-30*n+35*n^2+30*n^3+5*n^4)/5040) \\ Colin Barker, Feb 04 2015
    
  • Python
    def A254640(n): return n*(n*(n*(n*(n*(n*(n*(n*(10*n + 135) + 720) + 1890) + 2394) + 945) - 640) - 450) + 36)//5040 # Chai Wah Wu, Dec 07 2021
  • Sage
    [binomial(n+3,4)*(2*n+3)*(5*n^4 +30*n^3 +35*n^2 -30*n +2)/210 for n in (1..30)] # G. C. Greubel, Aug 28 2019
    

Formula

a(n) = n*(1+n)*(2+n)*(3+n)*(3+2*n)*(2 -30*n +35*n^2 +30*n^3 +5*n^4)/5040.
G.f.: x*(1+x)*(1 +56*x +246*x^2 +56*x^3 +x^4)/(1-x)^10. - Colin Barker, Feb 04 2015

A110555 Triangle of partial sums of alternating binomial coefficients: T(n, k) = Sum_{j=0..k} binomial(n, j)*(-1)^j, for n >= 0, 0 <= k <= n.

Original entry on oeis.org

1, 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
Offset: 0

Author

Reinhard Zumkeller, Jul 27 2005

Keywords

Examples

			Triangle T(n, k) starts:
  [0] 1;
  [1] 1,  0;
  [2] 1, -1,  0;
  [3] 1, -2,  1,   0;
  [4] 1, -3,  3,  -1,  0;
  [5] 1, -4,  6,  -4,  1,   0;
  [6] 1, -5, 10, -10,  5,  -1,  0;
  [7] 1, -6, 15, -20, 15,  -6,  1,  0;
  [8] 1, -7, 21, -35, 35, -21,  7, -1,  0.
		

Crossrefs

T(n,1) = -n + 1 for n>0;
T(n,2) = A000217(n-2) for n > 1;
T(n,3) = -A000292(n-4) for n > 2;
T(n,4) = A000332(n-1) for n > 3;
T(n,5) = -A000389(n-1) for n > 5;
T(n,6) = A000579(n-1) for n > 6;
T(n,7) = -A000580(n-1) for n > 7;
T(n,8) = A000581(n-1) for n > 8;
T(n,9) = -A000582(n-1) for n > 9;
T(n,10) = A001287(n-1) for n > 10;
T(n,11) = -A001288(n-1) for n > 11;
T(n,12) = A010965(n-1) for n > 12;
T(n,13) = -A010966(n-1) for n > 13;
T(n,14) = A010967(n-1) for n > 14;
T(n,15) = -A010968(n-1) for n > 15;
T(n,16) = A010969(n-1) for n > 16.
Cf. A071919 (variant), A000007 (row sums), A110556 (central terms).

Programs

  • Maple
    T := (n, k) -> (-1)^k * binomial(n-1, k):
    seq(print(seq(T(n, k), k = 0..n)), n = 0..7); # Peter Luschny, Apr 13 2023
  • Mathematica
    T[0, 0] := 1;  T[n_, n_] := 0; T[n_, k_] := (-1)^k*Binomial[n - 1, k]; Table[T[n, k], {n, 0, 20}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 31 2017 *)
  • PARI
    concat(1, for(n=1,10, for(k=0,n, print1(if(k != n, (-1)^k*binomial(n-1,k), 0), ", ")))) \\ G. C. Greubel, Aug 31 2017

Formula

T(n, 0) = 1, T(n, n) = 0^n, T(n, k) = -T(n-1, k-1) + T(n-1, k), for 0 < k < n.
T(n, k) = binomial(n-1, k)*(-1)^k, 0 <= k < n, T(n, n) = 0^n.
T(n, n-k-1) = -T(n, k), for 0 < k < n.
T(n, k) = A071919(n, k)*(-1)^k and A071919(n, k) = abs(T(n, k)).
Triangle T(n,k), 0 <= k <= n, read by rows, given by [1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [0, -1, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Sep 05 2005
G.f.: (1 + x*y) / (1 + x*y - x). - R. J. Mathar, Aug 11 2015

Extensions

Typo in name corrected by Andrey Zabolotskiy, Feb 22 2022
Offset corrected by Peter Luschny, Apr 13 2023
Previous Showing 11-20 of 75 results. Next