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.

A119258 Triangle read by rows: T(n,0) = T(n,n) = 1 and for 0

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 7, 1, 1, 7, 17, 15, 1, 1, 9, 31, 49, 31, 1, 1, 11, 49, 111, 129, 63, 1, 1, 13, 71, 209, 351, 321, 127, 1, 1, 15, 97, 351, 769, 1023, 769, 255, 1, 1, 17, 127, 545, 1471, 2561, 2815, 1793, 511, 1, 1, 19, 161, 799, 2561, 5503, 7937, 7423, 4097, 1023, 1
Offset: 0

Views

Author

Reinhard Zumkeller, May 11 2006

Keywords

Comments

From Richard M. Green, Jul 26 2011: (Start)
T(n,n-k) is the (k-1)-st Betti number of the subcomplex of the n-dimensional half cube obtained by deleting the interiors of all half-cube shaped faces of dimension at least k.
T(n,n-k) is the (k-2)-nd Betti number of the complement of the k-equal real hyperplane arrangement in R^n.
T(n,n-k) gives a lower bound for the complexity of the problem of determining, given n real numbers, whether some k of them are equal.
T(n,n-k) is the number of nodes used by the Kronrod-Patterson-Smolyak cubature formula in numerical analysis. (End)

Examples

			Triangle begins as:
  1;
  1, 1;
  1, 3,  1;
  1, 5,  7,  1;
  1, 7, 17, 15,  1;
  1, 9, 31, 49, 31, 1;
		

Crossrefs

A145661, A119258 and A118801 are all essentially the same (see the Shattuck and Waldhauser paper). - Tamas Waldhauser, Jul 25 2011

Programs

  • Haskell
    a119258 n k = a119258_tabl !! n !! k
    a119258_row n = a119258_tabl !! n
    a119258_list = concat a119258_tabl
    a119258_tabl = iterate (\row -> zipWith (+)
       ([0] ++ init row ++ [0]) $ zipWith (+) ([0] ++ row) (row ++ [0])) [1]
    -- Reinhard Zumkeller, Nov 15 2011
    
  • Magma
    function T(n,k)
      if k eq 0 or k eq n then return 1;
      else return 2*T(n-1,k-1) + T(n-1,k);
      end if;
      return T;
    end function;
    [T(n,k): k in [0..n], n in [0..12]]; // G. C. Greubel, Nov 18 2019
    
  • Maple
    # Case m = 2 of the more general:
    A119258 := (n,k,m) -> (1-m)^(-n+k)-m^(k+1)*pochhammer(n-k, k+1) *hypergeom([1,n+1],[k+2],m)/(k+1)!;
    seq(seq(round(evalf(A119258(n,k,2))),k=0..n), n=0..10); # Peter Luschny, Jul 25 2014
  • Mathematica
    T[n_, k_] := Binomial[n, k] Hypergeometric2F1[-k, n-k, n-k+1, -1];
    Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Sep 10 2017 *)
  • PARI
    T(n,k) = if(k==0 || k==n, 1, 2*T(n-1, k-1) + T(n-1,k) ); \\ G. C. Greubel, Nov 18 2019
    
  • Sage
    @CachedFunction
    def T(n, k):
        if (k==0 or k==n): return 1
        else: return 2*T(n-1, k-1) + T(n-1, k)
    [[T(n, k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Nov 18 2019

Formula

T(2*n,n-1) = T(2*n-1,n) for n > 0;
central terms give A119259; row sums give A007051;
T(n,0) = T(n,n) = 1;
T(n,1) = A005408(n-1) for n > 0;
T(n,2) = A056220(n-1) for n > 1;
T(n,n-4) = A027608(n-4) for n > 3;
T(n,n-3) = A055580(n-3) for n > 2;
T(n,n-2) = A000337(n-1) for n > 1;
T(n,n-1) = A000225(n) for n > 0.
T(n,k) = [k<=n]*(-1)^k*Sum_{i=0..k} (-1)^i*C(k-n,k-i)*C(n,i). - Paul Barry, Sep 28 2007
From Richard M. Green, Jul 26 2011: (Start)
T(n,k) = [k<=n] Sum_{i=n-k..n} (-1)^(n-k-i)*2^(n-i)*C(n,i).
T(n,k) = [k<=n] Sum_{i=n-k..n} C(n,i)*C(i-1,n-k-1).
G.f. for T(n,n-k): x^k/(((1-2x)^k)*(1-x)). (End)
T(n,k) = R(n,k,2) where R(n, k, m) = (1-m)^(-n+k)-m^(k+1)*Pochhammer(n-k,k+1)* hyper2F1([1,n+1], [k+2], m)/(k+1)!. - Peter Luschny, Jul 25 2014
From Peter Bala, Mar 05 2018: (Start)
The n-th row polynomial R(n,x) equals the n-th degree Taylor polynomial of the function (1 + 2*x)^n/(1 + x) about 0. For example, for n = 4 we have (1 + 2*x)^4/(1 + x) = 1 + 7*x + 17*x^2 + 15*x^3 + x^4 + O(x^5).
Row reverse of A112857. (End)

A119259 Central terms of the triangle in A119258.

Original entry on oeis.org

1, 3, 17, 111, 769, 5503, 40193, 297727, 2228225, 16807935, 127574017, 973168639, 7454392321, 57298911231, 441739706369, 3414246490111, 26447737520129, 205272288591871, 1595964714385409, 12427568655368191, 96905907580960769, 756583504975757311, 5913649000782757889
Offset: 0

Views

Author

Reinhard Zumkeller, May 11 2006

Keywords

Comments

The Gauss congruences a(n*p^k) == a(n^p^(k-1)) (mod p^k) hold for prime p and positive integers n and k. - Peter Bala, Jan 06 2022

References

  • R. P. Stanley, Enumerative Combinatorics Volume 2, Cambridge Univ. Press, 1999, Theorem 6.33, p. 197.

Crossrefs

Programs

  • Haskell
    a119259 n = a119258 (2 * n) n  -- Reinhard Zumkeller, Aug 06 2014
    
  • Mathematica
    Table[Binomial[2k - 1, k] Hypergeometric2F1[-2k, -k, 1 - 2k, -1], {k, 0, 10}] (* Vladimir Reshetnikov, Feb 16 2011 *)
  • Python
    from itertools import count, islice
    def A119259_gen(): # generator of terms
        yield from (1,3)
        a, c = 2, 1
        for n in count(1):
            yield (a<>1
    A119259_list = list(islice(A119259_gen(),20)) # Chai Wah Wu, Apr 26 2023

Formula

a(n) = A119258(2*n,n).
a(n) = Sum_{k=0..n} C(2*n,k)*C(2*n-k-1,n-k). - Paul Barry, Sep 28 2007
a(n) = Sum_{k=0..n} C(n+k-1,k)*2^k. - Paul Barry, Sep 28 2007
2*a(n) = A064062(n)+A178792(n). - Joseph Abate, Jul 21 2010
G.f.: (4*x^2+3*sqrt(1-8*x)*x-5*x)/(sqrt(1-8*x)*(2*x^2+x-1)-8*x^2-7*x+1). - Vladimir Kruchinin, Aug 19 2013
a(n) = (-1)^n - 2^(n+1)*binomial(2*n,n-1)*hyper2F1([1,2*n+1],[n+2],2). - Peter Luschny, Jul 25 2014
a(n) = (-1)^n + 2^(n+1)*A014300(n). - Peter Luschny, Jul 25 2014
a(n) = [x^n] ( (1 + x)^2/(1 - x) )^n. Exp( Sum_{n >= 1} a(n)*x^n/n ) = 1 + 3*x + 13*x^2 + 67*x^3 + ... is essentially the o.g.f. for A064062. - Peter Bala, Oct 01 2015
The o.g.f. is the diagonal of the bivariate rational function 1/(1 - t*(1 + x)^2/(1 - x)) and hence is algebraic by Stanley 1999, Theorem 6.33, p.197. - Peter Bala, Aug 21 2016
n*(3*n-4)*a(n) +(-21*n^2+40*n-12)*a(n-1) -4*(3*n-1)*(2*n-3)*a(n-2)=0. - R. J. Mathar, Aug 09 2017
From Peter Bala, Mar 23 2020: (Start)
a(p) == 3 ( mod p^3 ) for prime p >= 5. Cf. A002003, A103885 and A156894.
More generally, we conjecture that a(n*p^k) == a(n*p^(k-1)) ( mod p^(3*k) ) for prime p >= 5 and positive integers n and k. (End)
G.f.: (8*x)/(sqrt(1-8*x)*(1+4*x)-1+8*x). - Fabian Pereyra, Jul 20 2024
a(n) = 2^(n+1)*binomial(2*n,n) - A178792(n). - Akiva Weinberger, Dec 06 2024
a(n) = Sum_{k=0..n} 2^k * (-1)^(n-k) * binomial(2*n,k). - Seiichi Manyama, Jul 31 2025

A112857 Triangle T(n,k) read by rows: number of Green's R-classes in the semigroup of order-preserving partial transformations (of an n-element chain) consisting of elements of height k (height(alpha) = |Im(alpha)|).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 5, 1, 1, 15, 17, 7, 1, 1, 31, 49, 31, 9, 1, 1, 63, 129, 111, 49, 11, 1, 1, 127, 321, 351, 209, 71, 13, 1, 1, 255, 769, 1023, 769, 351, 97, 15, 1, 1, 511, 1793, 2815, 2561, 1471, 545, 127, 17, 1, 1, 1023, 4097, 7423, 7937, 5503, 2561, 799, 161, 19, 1
Offset: 0

Views

Author

Abdullahi Umar, Aug 25 2008

Keywords

Comments

Sum of rows of T(n, k) is A007051; T(n,k) = |A118801(n,k)|.
Row-reversed variant of A119258. - R. J. Mathar, Jun 20 2011
Pairwise sums of row terms starting from the right yields triangle A038207. - Gary W. Adamson, Feb 06 2012
Riordan array (1/(1 - x), x/(1 - 2*x)). - Philippe Deléham, Jan 17 2014
Appears to coincide with the triangle T(n,m) (n >= 1, 1 <= m <= n) giving number of set partitions of [n], avoiding 1232, with m blocks [Crane, 2015]. See also A250118, A250119. - N. J. A. Sloane, Nov 25 2014
(A007318)^2 = A038207 = T*|A167374|. See A118801 for other relations to the Pascal matrix. - Tom Copeland, Nov 17 2016

Examples

			T(3,2) = 5 because in a regular semigroup of transformations the Green's R-classes coincide with convex partitions of subsets of {1,2,3} with convex classes (modulo the subsets): {1}, {2}/{1}, {3}/{2}, {3}/{1,2}, {3}/{1}, {2,3}
Triangle T(n,k) (with rows n >= 0 and columns k = 0..n) begins:
  1;
  1,    1;
  1,    3,    1;
  1,    7,    5,    1;
  1,   15,   17,    7,    1;
  1,   31,   49,   31,    9,    1;
  1,   63,  129,  111,   49,   11,    1;
  1,  127,  321,  351,  209,   71,   13,   1;
  1,  255,  769, 1023,  769,  351,   97,  15,   1;
  1,  511, 1793, 2815, 2561, 1471,  545, 127,  17,  1;
  1, 1023, 4097, 7423, 7937, 5503, 2561, 799, 161, 19, 1;
  ...
As to matrix M, top row of M^3 = (1, 7, 5, 1, 0, 0, 0, ...)
		

Crossrefs

Programs

  • Maple
    A112857 := proc(n,k) if k=0 or k=n then 1; elif k <0 or k>n then 0; else 2*procname(n-1,k)+procname(n-1,k-1) ; end if; end proc: # R. J. Mathar, Jun 20 2011
  • Mathematica
    Table[Abs[1 + (-1)^k*2^(n - k + 1)*Sum[ Binomial[n - 2 j - 2, k - 2 j - 1], {j, 0, Floor[k/2]}]] - 4 Boole[And[n == 1, k == 0]], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Nov 24 2016 *)

Formula

T(n,k) = Sum_{j = k..n} C(n,j)*C(j-1,k-1).
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) for n >= 2 and 1 <= k <= n-1 with T(n,0) = 1 = T(n,n) for n >= 0.
n-th row = top row of M^n, deleting the zeros, where M is an infinite square production matrix with (1,1,1,...) as the superdiagonal and (1,2,2,2,...) as the main diagonal. - Gary W. Adamson, Feb 06 2012
From Peter Bala, Mar 05 2018 (Start):
The following remarks are particular cases of more general results for Riordan arrays of the form (f(x), x/(1 - k*x)).
Let R(n,x) denote the n-th row polynomial of this triangle. The polynomial R(n,2*x) has the e.g.f. Sum_{k = 0..n} T(n,k)*(2*x)^k/k!. The e.g.f. for the n-th diagonal of the triangle (starting at n = 0 for the main diagonal) equals exp(x) * the e.g.f. for the polynomial R(n,2*x). For example, when n = 3 we have exp(x)*(1 + 7*(2*x) + 5*(2*x)^2/2! + (2*x)^3/3!) = 1 + 15*x + 49*x^2/2! + 111*x^3/3! + 209*x^4/4! + ....
Let P(n,x) = Sum_{k = 0..n} T(n,k)*x^(n-k) denote the n-th row polynomial in descending powers of x. Then P(n,x) is the n-th degree Taylor polynomial of the function (1 + 2*x)^n/(1 + x) about 0. For example, for n = 4 we have (1 + 2*x)^4/(1 + x) = x^4 + 15*x^3 + 17*x^2 + 7*x + 1 + O(x^5).
See A118801 for a signed version of this triangle and A145661 for a signed version of the row reversed triangle. (End)
Bivariate o.g.f.: Sum_{n,k>=0} T(n,k)*x^n*y^k = (1 - 2*x)/((1 - x)*(1 - 2*x - x*y)). - Petros Hadjicostas, Feb 14 2021
The matrix inverse of the Lucas triangle A029635 is -T(n, k)/(-2)^(n-k+1). - Peter Luschny, Dec 22 2024

A055584 Triangle of partial row sums (prs) of triangle A055252.

Original entry on oeis.org

1, 5, 1, 19, 6, 1, 63, 25, 7, 1, 192, 88, 32, 8, 1, 552, 280, 120, 40, 9, 1, 1520, 832, 400, 160, 49, 10, 1, 4048, 2352, 1232, 560, 209, 59, 11, 1, 10496, 6400, 3584, 1792, 769, 268, 70, 12, 1, 26624, 16896, 9984, 5376, 2561, 1037, 338, 82, 13, 1, 66304, 43520
Offset: 0

Views

Author

Wolfdieter Lang, May 26 2000

Keywords

Comments

In the language of the Shapiro et al. reference (given in A053121) such a lower triangular (ordinary) convolution array, considered as matrix, belongs to the Riordan-group. The G.f. for the row polynomials p(n,x) (increasing powers of x) is (((1-z)^3)/(1-2*z)^4)/(1-x*z/(1-z)).
This is the fourth member of the family of Riordan-type matrices obtained from A007318(n,m) (Pascal's triangle read as lower triangular matrix) by repeated application of the prs-procedure.
The column sequences appear as A049612(n+1), A055585, A001794, A001789(n+3), A027608, A055586 for m=0..5.

Examples

			[0] 1
[1] 5, 1
[2] 19, 6, 1
[3] 63, 25, 7, 1
[4] 192, 88, 32, 8, 1
[5] 552, 280, 120, 40, 9, 1
[6] 1520, 832, 400, 160, 49, 10, 1
[7] 4048, 2352, 1232, 560, 209, 59, 11, 1
Fourth row polynomial (n=3): p(3, x)= 63 + 25*x + 7*x^2 + x^3.
		

Crossrefs

Cf. A007318, A055248, A055249, A055252. Row sums: A049600(n+1, 4).

Programs

  • Maple
    T := (n, k) -> binomial(n, k)*hypergeom([4, k - n], [k + 1], -1):
    for n from 0 to 7 do seq(simplify(T(n, k)), k = 0..n) od; # Peter Luschny, Sep 23 2024

Formula

a(n, m)=sum(A055252(n, k), k=m..n), n >= m >= 0, a(n, m) := 0 if n
Column m recursion: a(n, m)= sum(a(j, m), j=m..n-1)+ A055252(n, m), n >= m >= 0, a(n, m) := 0 if n
G.f. for column m: (((1-x)^3)/(1-2*x)^4)*(x/(1-x))^m, m >= 0.
T(n, k) = binomial(n, k)*hypergeom([4, k - n], [k + 1], -1). - Peter Luschny, Sep 23 2024

A193844 Triangular array: the fission of ((x+1)^n) by ((x+1)^n); i.e., the self-fission of Pascal's triangle.

Original entry on oeis.org

1, 1, 3, 1, 5, 7, 1, 7, 17, 15, 1, 9, 31, 49, 31, 1, 11, 49, 111, 129, 63, 1, 13, 71, 209, 351, 321, 127, 1, 15, 97, 351, 769, 1023, 769, 255, 1, 17, 127, 545, 1471, 2561, 2815, 1793, 511, 1, 19, 161, 799, 2561, 5503, 7937, 7423, 4097, 1023
Offset: 0

Author

Clark Kimberling, Aug 07 2011

Keywords

Comments

See A193842 for the definition of fission of two sequences of polynomials or triangular arrays.
A193844 is also the fission of (p1(n,x)) by (q1(n,x)), where p1(n,x)=x^n+x^(n-1)+...+x+1 and q1(n,x)=(x+2)^n.
Essentially A119258 but without the main diagonal. - Peter Bala, Jul 16 2013
From Robert Coquereaux, Oct 02 2014: (Start)
This is also a rectangular array A(n,p) read down the antidiagonals:
1 1 1 1 1 1 1 1 1
3 5 7 9 11 13 15 17 19
7 17 31 49 71 97 127 161 199
15 49 111 209 351 545 799 1121 1519
31 129 351 769 1471 2561 4159 6401 9439
...
Calling Gr(n) the Grassmann algebra with n generators, A(n,p) is the dimension of the space of Gr(n)-valued symmetric multilinear forms with vanishing graded divergence. If p is odd A(n,p) is 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 A(n,p)+1. A(n,p) = 2^n*A059260(p,n-1)-(-1)^p.
(End)
The n-th row are also the coefficients of the polynomial P=sum_{k=0..n} (X+2)^k (in falling order, i.e., that of X^n first). - M. F. Hasler, Oct 15 2014

Examples

			First six rows:
1
1....3
1....5....7
1....7....17....15
1....9....31....49....31
1....11...49....111...129...63
		

Crossrefs

A145661 is an essentially identical triangle.

Programs

  • Maple
    A193844 := (n,k) -> 2^k*binomial(n+1,k)*hypergeom([1,-k],[-k+n+2],1/2);
    for n from 0 to 5 do seq(round(evalf(A193844(n,k))),k=0..n) od; # Peter Luschny, Jul 23 2014
    # Alternatively
    p := (n,x) -> add(x^k*(1+2*x)^(n-k), k=0..n): for n from 0 to 7 do [n], PolynomialTools:-CoefficientList(p(n,x), x) od; # Peter Luschny, Jun 18 2017
  • Mathematica
    z = 10;
    p[n_, x_] := (x + 1)^n;
    q[n_, x_] := (x + 1)^n
    p1[n_, k_] := Coefficient[p[n, x], x^k];
    p1[n_, 0] := p[n, x] /. x -> 0;
    d[n_, x_] := Sum[p1[n, k]*q[n - 1 - k, x], {k, 0, n - 1}]
    h[n_] := CoefficientList[d[n, x], {x}]
    TableForm[Table[Reverse[h[n]], {n, 0, z}]]
    Flatten[Table[Reverse[h[n]], {n, -1, z}]]  (* A193844 *)
    TableForm[Table[h[n], {n, 0, z}]]
    Flatten[Table[h[n], {n, -1, z}]]  (* A193845 *)
  • Sage
    # uses[fission from A193842]
    p = lambda n,x: (x+1)^n
    A193844_row = lambda n: fission(p, p, n)
    for n in range(7): print(A193844_row(n)) # Peter Luschny, Jul 23 2014

Formula

From Peter Bala, Jul 16 2013: (Start)
T(n,k) = sum {i = 0..k} (-1)^i*binomial(n+1,k-i)*2^(k-i).
O.g.f.: 1/( (1 - x*t)*(1 - (2*x + 1)*t) ) = 1 + (1 + 3*x)*t + (1 + 5*x + 7*x^2)*t^2 + ....
The n-th row polynomial R(n,x) = 1/(x+1)*( (2*x+1)^(n+1) - x^(n+1) ). (End)
T(n,k) = T(n-1,k) + 3*T(n-1,k-1) - T(n-2,k-1) - 2*T(n-2,k-2), T(0,0) = 1, T(1,0) = 1, T(1,1) = 3, T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Jan 17 2014
T(n,k) = 2^k*binomial(n+1,k)*hyper2F1(1,-k,-k+n+2, 1/2). - Peter Luschny, Jul 23 2014

A193845 Mirror of the triangle A193844.

Original entry on oeis.org

1, 3, 1, 7, 5, 1, 15, 17, 7, 1, 31, 49, 31, 9, 1, 63, 129, 111, 49, 11, 1, 127, 321, 351, 209, 71, 13, 1, 255, 769, 1023, 769, 351, 97, 15, 1, 511, 1793, 2815, 2561, 1471, 545, 127, 17, 1, 1023, 4097, 7423, 7937, 5503, 2561, 799, 161, 19, 1
Offset: 0

Author

Clark Kimberling, Aug 07 2011

Keywords

Comments

This triangle is obtained by reversing the rows of the triangle A193844.
From Philippe Deléham, Jan 17 2014: (Start)
Subtriangle of the triangle in A112857.
T(n,0) = A000225(n+1).
T(n,1) = A000337(n).
T(n+2,2) = A055580(n).
T(n+3,3) = A027608(n).
T(n+4,4) = A211386(n).
T(n+5,5) = A211388(n).
T(n,n) = A000012(n).
T(n+1,n) = A005408(n).
T(n+2,n) = A056220(n+2).
T(n+3,n) = A199899(n+1).
Row sums are A003462(n+1).
Diagonal sums are A048739(n).
Riordan array (1/((1-2*x)*(1-x)), x/(1-2*x)). (End)
Consider the transformation 1 + x + x^2 + x^3 + ... + x^n = A_0*(x-2)^0 + A_1*(x-2)^1 + A_2*(x-2)^2 + ... + A_n*(x-2)^n. This sequence gives A_0, ... A_n as the entries in the n-th row of this triangle, starting at n = 0. - Derek Orr, Oct 14 2014
The n-th row lists the coefficients of the polynomial sum_{k=0..n} (X+2)^k, in order of increasing powers. - M. F. Hasler, Oct 15 2014

Examples

			First six rows:
1
3....1
7....5....1
15...17...7....1
31...49...31...9...1
63...129..111..49..11..1
		

Crossrefs

Programs

  • Mathematica
    z = 10;
    p[n_, x_] := (x + 1)^n;
    q[n_, x_] := (x + 1)^n
    p1[n_, k_] := Coefficient[p[n, x], x^k];
    p1[n_, 0] := p[n, x] /. x -> 0;
    d[n_, x_] := Sum[p1[n, k]*q[n - 1 - k, x], {k, 0, n - 1}]
    h[n_] := CoefficientList[d[n, x], {x}]
    TableForm[Table[Reverse[h[n]], {n, 0, z}]]
    Flatten[Table[Reverse[h[n]], {n, -1, z}]]  (* A193844 *)
    TableForm[Table[h[n], {n, 0, z}]]
    Flatten[Table[h[n], {n, -1, z}]]  (* A193845 *)
    Table[2^k*Binomial[n + 1, k]*Hypergeometric2F1[1, -k, -k + n + 2, 1/2], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Michael De Vlieger, Nov 09 2021 *)
  • PARI
    for(n=0,20,for(k=0,n,print1(1/k!*sum(i=0,n,(2^(i-k)*prod(j=0,k-1,i-j))),", "))) \\ Derek Orr, Oct 14 2014

Formula

T(n,k) = A193844(n,n-k).
T(n,k) = 3*T(n-1,k) + T(n-1,k-1) - 2*T(n-2,k) - T(n-2,k-1), T(0,0) = 1, T(1,0) = 3, T(1,1) = 1, T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Jan 17 2014

A055586 Sixth column of triangle A055584.

Original entry on oeis.org

1, 10, 59, 268, 1037, 3598, 11535, 34832, 100369, 278546, 749587, 1966100, 5046293, 12714006, 31522839, 77070360, 186122265, 444596250, 1051721755, 2466250780, 5737807901, 13254000670, 30417092639, 69390565408
Offset: 0

Author

Wolfdieter Lang, May 26 2000

Keywords

Crossrefs

Cf. A055584.
Partial sums of A027608.

Programs

  • Mathematica
    CoefficientList[Series[1/(((1-x)^2)(1-2x)^4),{x,0,30}],x] (* or *) LinearRecurrence[{10,-41,88,-104,64,-16},{1,10,59,268,1037,3598},30] (* Harvey P. Dale, Jul 31 2025 *)
  • PARI
    Vec(1/(((1-x)^2)*(1-2*x)^4) + O(x^40)) \\ Michel Marcus, Dec 11 2015

Formula

G.f.: 1/(((1-x)^2)*(1-2*x)^4).
a(n) = A055584(n+5, 5).
a(n) = Sum_{j=0..n-1} a(j) + A055582(n) for n >= 1.
E.g.f.: exp(x)*(x + 9) + 8*exp(2*x)*(2*x^3 + 3*x^2 + 6*x - 3)/3. - Stefano Spezia, Sep 24 2024

A375550 Triangle read by rows: T(m, n, k) = binomial(n + 1, n - k)*hypergeom([m, k - n], [k + 2], -1) for m = 4.

Original entry on oeis.org

1, 6, 1, 25, 7, 1, 88, 32, 8, 1, 280, 120, 40, 9, 1, 832, 400, 160, 49, 10, 1, 2352, 1232, 560, 209, 59, 11, 1, 6400, 3584, 1792, 769, 268, 70, 12, 1, 16896, 9984, 5376, 2561, 1037, 338, 82, 13, 1, 43520, 26880, 15360, 7937, 3598, 1375, 420, 95, 14, 1
Offset: 0

Author

Peter Luschny, Sep 23 2024

Keywords

Comments

Triangle T(m,n,k) is a Riordan array of the form ((1-x)^(m-1)*(1-2x)^(-m-1), x/(1-x)), for m = 3. - Igor Victorovich Statsenko, Feb 08 2025

Examples

			Triangle starts:
  [0]     1;
  [1]     6,     1;
  [2]    25,     7,     1;
  [3]    88,    32,     8,    1;
  [4]   280,   120,    40,    9,    1;
  [5]   832,   400,   160,   49,   10,    1;
  [6]  2352,  1232,   560,  209,   59,   11,   1;
  [7]  6400,  3584,  1792,  769,  268,   70,  12,  1;
  [8] 16896,  9984,  5376, 2561, 1037,  338,  82, 13,  1;
  [9] 43520, 26880, 15360, 7937, 3598, 1375, 420, 95, 14, 1;
  ...
Seen as an array of the columns:
  [0] 1,  6, 25,  88,  280,  832,  2352,  6400,  16896, ...
  [1] 1,  7, 32, 120,  400, 1232,  3584,  9984,  26880, ...
  [2] 1,  8, 40, 160,  560, 1792,  5376, 15360,  42240, ...
  [3] 1,  9, 49, 209,  769, 2561,  7937, 23297,  65537, ...
  [4] 1, 10, 59, 268, 1037, 3598, 11535, 34832, 100369, ...
  [5] 1, 11, 70, 338, 1375, 4973, 16508, 51340, 151709, ...
  [6] 1, 12, 82, 420, 1795, 6768, 23276, 74616, 226325, ...
		

Crossrefs

Column k: A055585 (k=0), A001794 (k=1), A001789 (k=2), A027608 (k=3), A055586 (k=4).
Cf. A145018 (diagonal n-2), A375549 (row sums), A049612 (alternating row sums), A122433.

Programs

  • Maple
    T := (m, n, k) -> binomial(n + 1, n - k)*hypergeom([m, k - n], [k + 2], -1);
    for n from 0 to 9 do seq(simplify(T(4, n, k)), k = 0..n) od;
    # As a binomial sum:
    T := (m, n, k) -> add(binomial(m + j, m)*binomial(n + 1, n - (j + k)), j = 0..n-k):
    for n from 0 to 9 do [n], seq(T(3, n, k), k = 0..n) od;
    # Alternative, generating the array of the columns:
    cgf := k -> (1 - x)^(2 - k) / (1 - 2*x)^4:
    ser := (k, len) -> series(cgf(k), x, len + 2):
    Tcol := (k, len) -> seq(coeff(ser(k, len), x, j), j = 0..len):
    seq(lprint([k], Tcol(k, 8)), k = 0..6);

Formula

T(m, n, k) = Sum_{j=0..n-k} binomial(m + j, m)*binomial(n + 1, n - (j + k)) for m = 3.
G.f. of column k: (1 - x)^(2 - k) / (1 - 2*x)^4.
Showing 1-8 of 8 results.