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 35 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

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)

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

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

A128908 Riordan array (1, x/(1-x)^2).

Original entry on oeis.org

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

Author

Philippe Deléham, Apr 22 2007

Keywords

Comments

Triangle T(n,k), 0 <= k <= n, read by rows given by [0,2,-1/2,1/2,0,0,0,0,0,...] DELTA [1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.
Row sums give A088305. - Philippe Deléham, Nov 21 2007
Column k is C(n,2k-1) for k > 0. - Philippe Deléham, Jan 20 2012
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
T is the convolution triangle of the positive integers (see A357368). - Peter Luschny, Oct 19 2022

Examples

			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    2    1
   3:   0    3    4    1
   4:   0    4   10    6    1
   5:   0    5   20   21    8    1
   6:   0    6   35   56   36   10    1
   7:   0    7   56  126  120   55   12    1
   8:   0    8   84  252  330  220   78   14    1
   9:   0    9  120  462  792  715  364  105   16    1
  10:   0   10  165  792 1716 2002 1365  560  136   18    1
  ... reformatted by _Wolfdieter Lang_, Jul 31 2017
From _Peter Luschny_, Mar 06 2022: (Start)
The sequence can also be seen as a square array read by upwards antidiagonals.
   1, 1,   1,    1,    1,     1,     1,      1,      1, ...  A000012
   0, 2,   4,    6,    8,    10,    12,     14,     16, ...  A005843
   0, 3,  10,   21,   36,    55,    78,    105,    136, ...  A014105
   0, 4,  20,   56,  120,   220,   364,    560,    816, ...  A002492
   0, 5,  35,  126,  330,   715,  1365,   2380,   3876, ... (A053126)
   0, 6,  56,  252,  792,  2002,  4368,   8568,  15504, ... (A053127)
   0, 7,  84,  462, 1716,  5005, 12376,  27132,  54264, ... (A053128)
   0, 8, 120,  792, 3432, 11440, 31824,  77520, 170544, ... (A053129)
   0, 9, 165, 1287, 6435, 24310, 75582, 203490, 490314, ... (A053130)
    A27,A292, A389, A580,  A582, A1288, A10966, A10968, A165817       (End)
		

Crossrefs

Cf. A165817 (the main diagonal of the array).

Programs

  • Maple
    # Computing the rows of the array representation:
    S := proc(n,k) option remember;
    if n = k then 1 elif k < 0 or k > n then 0 else
    S(n-1, k-1) + 2*S(n-1, k) - S(n-2, k) fi end:
    Arow := (n, len) -> seq(S(n+k-1, k-1), k = 0..len-1):
    for n from 0 to 8 do Arow(n, 9) od; # Peter Luschny, Mar 06 2022
    # Uses function PMatrix from A357368.
    PMatrix(10, n -> n); # Peter Luschny, Oct 19 2022
  • Mathematica
    With[{nmax = 10}, CoefficientList[CoefficientList[Series[(1 - x)^2/(1 - (2 + y)*x + x^2), {x, 0, nmax}, {y, 0, nmax}], x], y]] // Flatten (* G. C. Greubel, Nov 22 2017 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(if(n==0 && k==0, 1, if(k==0, 0, binomial(n+k-1,2*k-1))), ", "))) \\ G. C. Greubel, Nov 22 2017
    
  • Python
    from functools import cache
    @cache
    def A128908(n, k):
        if n == k: return 1
        if (k <= 0 or k > n): return 0
        return A128908(n-1, k-1) + 2*A128908(n-1, k) - A128908(n-2, k)
    for n in range(10):
        print([A128908(n, k) for k in range(n+1)]) # Peter Luschny, Mar 07 2022
  • 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))
    A128908 = lambda n,k: T(k,n)
    for n in (0..10): print([A128908(n,k) for k in (0..n)]) # Peter Luschny, Mar 12 2016
    

Formula

T(n,0) = 0^n, T(n,k) = binomial(n+k-1, 2k-1) for k >= 1.
Sum_{k=0..n} T(n,k)*2^(n-k) = A002450(n) = (4^n-1)/3 for n>=1. - Philippe Deléham, Oct 19 2008
G.f.: (1-x)^2/(1-(2+y)*x+x^2). - Philippe Deléham, Jan 20 2012
Sum_{k=0..n} T(n,k)*x^k = (-1)^n*A001352(n), (-1)^(n+1)*A054888(n+1), (-1)^n*A008574(n), (-1)^n*A084103(n), (-1)^n*A084099(n), A163810(n), A000007(n), A088305(n) for x = -6, -5, -4, -3, -2, -1, 0, 1 respectively. - Philippe Deléham, Jan 20 2012
Riordan array (1, x/(1-x)^2). - Philippe Deléham, Jan 20 2012

A221857 Number A(n,k) of shapes of balanced k-ary trees with n nodes, where a tree is balanced if the total number of nodes in subtrees corresponding to the branches of any node differ by at most one; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 2, 1, 0, 1, 1, 3, 1, 1, 0, 1, 1, 4, 3, 4, 1, 0, 1, 1, 5, 6, 1, 4, 1, 0, 1, 1, 6, 10, 4, 9, 4, 1, 0, 1, 1, 7, 15, 10, 1, 27, 1, 1, 0, 1, 1, 8, 21, 20, 5, 16, 27, 8, 1, 0, 1, 1, 9, 28, 35, 15, 1, 96, 81, 16, 1, 0, 1, 1, 10, 36, 56, 35, 6, 25, 256, 81, 32, 1, 0
Offset: 0

Author

Alois P. Heinz, Apr 10 2013

Keywords

Examples

			: A(2,2) = 2  : A(2,3) = 3      : A(3,3) = 3          :
:   o     o   :   o    o    o   :   o      o      o   :
:  / \   / \  :  /|\  /|\  /|\  :  /|\    /|\    /|\  :
: o         o : o      o      o : o o    o   o    o o :
:.............:.................:.....................:
: A(3,4) = 6                                          :
:    o        o        o        o       o        o    :
:  /( )\    /( )\    /( )\    /( )\   /( )\    /( )\  :
: o o      o   o    o     o    o o     o   o      o o :
Square array A(n,k) begins:
  1, 1, 1,  1,   1,   1,  1,  1,  1,   1,   1, ...
  1, 1, 1,  1,   1,   1,  1,  1,  1,   1,   1, ...
  0, 1, 2,  3,   4,   5,  6,  7,  8,   9,  10, ...
  0, 1, 1,  3,   6,  10, 15, 21, 28,  36,  45, ...
  0, 1, 4,  1,   4,  10, 20, 35, 56,  84, 120, ...
  0, 1, 4,  9,   1,   5, 15, 35, 70, 126, 210, ...
  0, 1, 4, 27,  16,   1,  6, 21, 56, 126, 252, ...
  0, 1, 1, 27,  96,  25,  1,  7, 28,  84, 210, ...
  0, 1, 8, 81, 256, 250, 36,  1,  8,  36, 120, ...
		

Crossrefs

Rows n=0+1, 2-3, give: A000012, A001477, A179865.
Diagonal and upper diagonals give: A028310, A000217, A000292, A000332, A000389, A000579, A000580, A000581, A000582, A001287, A001288.
Lower diagonals give: A000012, A000290, A092364(n) for n>1.

Programs

  • Maple
    A:= proc(n, k) option remember; local m, r; if n<2 or k=1 then 1
          elif k=0 then 0 else r:= iquo(n-1, k, 'm');
          binomial(k, m)*A(r+1, k)^m*A(r, k)^(k-m) fi
        end:
    seq(seq(A(n, d-n), n=0..d), d=0..12);
  • Mathematica
    a[n_, k_] := a[n, k] = Module[{m, r}, If[n < 2 || k == 1, 1, If[k == 0, 0, {r, m} = QuotientRemainder[n-1, k]; Binomial[k, m]*a[r+1, k]^m*a[r, k]^(k-m)]]]; Table[a[n, d-n], {d, 0, 12}, {n, 0, d}] // Flatten (* Jean-François Alcover, Apr 17 2013, translated from Maple *)

A238801 Triangle T(n,k), read by rows, given by T(n,k) = C(n+1, k+1)*(1-(k mod 2)).

Original entry on oeis.org

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

Author

Philippe Deléham, Mar 05 2014

Keywords

Comments

Row sums are powers of 2.

Examples

			Triangle begins:
1;
2, 0;
3, 0, 1;
4, 0, 4, 0;
5, 0, 10, 0, 1;
6, 0, 20, 0, 6, 0;
7, 0, 35, 0, 21, 0, 1;
8, 0, 56, 0, 56, 0, 8, 0;
9, 0, 84, 0, 126, 0, 36, 0, 1;
10, 0, 120, 0, 252, 0, 120, 0, 10, 0; etc.
		

Programs

  • Mathematica
    Table[Binomial[n + 1, k + 1]*(1 - Mod[k , 2]), {n, 0, 10}, {k, 0, n}] // Flatten (* G. C. Greubel, Nov 22 2017 *)
  • PARI
    T(n,k) = binomial(n+1, k+1)*(1-(k % 2));
    tabl(nn) = for (n=0, nn, for (k=0, n, print1(T(n,k), ", ")); print); \\ Michel Marcus, Nov 23 2017

Formula

G.f.: 1/((1+(y-1)*x)*(1-(y+1)*x)).
T(n,k) = 2*T(n-1,k) + T(n-2,k-2) - T(n-2,k), T(0,0) = 1, T(1,0) = 2, T(1,1) = 0, T(n,k) = 0 if k<0 or if k>n.
Sum_{k=0..n} T(n,k)*x^k = A000027(n+1), A000079(n), A015518(n+1), A003683(n+1), A079773(n+1), A051958(n+1), A080920(n+1), A053455(n), A160958(n+1) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8 respectively.

A035927 One less than number of n-multisets chosen from a 10-set.

Original entry on oeis.org

0, 9, 54, 219, 714, 2001, 5004, 11439, 24309, 48619, 92377, 167959, 293929, 497419, 817189, 1307503, 2042974, 3124549, 4686824, 6906899, 10015004, 14307149, 20160074, 28048799, 38567099, 52451255, 70607459, 94143279
Offset: 0

Keywords

Comments

Number of distinct n-digit numbers up to permutations of digits. - Michael Somos, Jul 11 2002
Equivalently, for n > 0, a(n) = number of n-digit decimal numbers d_1 d_2 ... d_n with d_1 > 0 and d_1 >= d_2 >= ... >= d_n >= 0.. - N. J. A. Sloane, Jul 13 2023

Crossrefs

Equals A000582 - 1. Cf. A014553, A179239.

Programs

  • Maple
    binomial(10+n-1,n)-1;
  • Mathematica
    Table[Binomial[9 + n, n] - 1, {n, 0, 27}] (* Michael De Vlieger, Jul 14 2015 *)
    CoefficientList[Series[1/(1-x)^10-1/(1-x),{x,0,30}],x] (* or *) LinearRecurrence[{10,-45,120,-210,252,-210,120,-45,10,-1},{0,9,54,219,714,2001,5004,11439,24309,48619},30] (* Harvey P. Dale, Jul 11 2023 *)
  • PARI
    a(n)=if(n<0,0,binomial(n+9,9)-1)

Formula

G.f.: 1/(1-x)^10-1/(1-x). - Michael Somos, Jul 11 2002

A039948 A triangle related to A000045 (Fibonacci numbers).

Original entry on oeis.org

1, 1, 1, 4, 2, 1, 18, 12, 3, 1, 120, 72, 24, 4, 1, 960, 600, 180, 40, 5, 1, 9360, 5760, 1800, 360, 60, 6, 1, 105840, 65520, 20160, 4200, 630, 84, 7, 1, 1370880, 846720, 262080, 53760, 8400, 1008, 112, 8, 1, 19958400, 12337920, 3810240, 786240, 120960, 15120, 1512, 144, 9, 1
Offset: 0

Keywords

Examples

			Triangle begins :
    1;
    1,   1;
    4,   2,   1;
   18,  12,   3,  1;
  120,  72,  24,  4, 1;
  960, 600, 180, 40, 5, 1;
... - _Philippe Deléham_, Nov 08 2011
		

Programs

  • Magma
    [(Factorial(n)/Factorial(k))*Fibonacci(n-k+1): k in [0..n], n in [0..12]]; // G. C. Greubel, Nov 20 2022
    
  • Mathematica
    T[n_,k_]:= (n!/k!)*Fibonacci[n-k+1];
    Table[T[n,k], {n,0,12}, {k,0,n}]//Flatten (* G. C. Greubel, Nov 20 2022 *)
  • SageMath
    def A039948(n, k): return factorial(n-k)*binomial(n,k)*fibonacci(n-k+1)
    flatten([[A039948(n,k) for k in range(n+1)] for n in range(13)]) # G. C. Greubel, Nov 20 2022

Formula

T(n, m) = n!*Fibonacci(n-m+1)/m!, n >= m >= 0.
T(n, 0) = A005442(n).
T(n, 1) = A005443(n).
E.g.f. for column m: x^m/(m!*(1-x-x^2)), m >= 0.
From G. C. Greubel, Nov 20 2022: (Start)
T(n, n-1) = A000027(n).
T(n, n-2) = 4*A000217(n-1), n >= 2.
T(n, n-3) = 18*A000292(n-2), n >= 3.
T(n, n-4) = 5! * A000332(n), n >= 4.
T(n, n-5) = 8 * 5! * A000389(n), n >= 5.
T(n, n-6) = 13 * 6! * A000579(n), n >= 6.
T(n, n-7) = 21 * 7! * A000580(n), n >= 7.
T(n, n-8) = 34 * 8! * A000581(n), n >= 8.
T(n, n-9) = 55 * 9! * A000582(n), n >= 9.
T(n, n-10) = 89 * 10! * A001287(n), n >= 10.
T(n, n-11) = 12 * 12! * A001288(n), n >= 11.
T(n, n-12) = 233 * 12! * A010965(n), n >= 12.
T(n, n-13) = 89 * 13! * A010966(n), n >= 13.
Sum_{k=0..n} T(n, k) = A110313(n). (End)
Previous Showing 11-20 of 35 results. Next