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

A256890 Triangle T(n,k) = t(n-k, k); t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.

Original entry on oeis.org

1, 2, 2, 4, 12, 4, 8, 52, 52, 8, 16, 196, 416, 196, 16, 32, 684, 2644, 2644, 684, 32, 64, 2276, 14680, 26440, 14680, 2276, 64, 128, 7340, 74652, 220280, 220280, 74652, 7340, 128, 256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172, 256, 512, 72076, 1637860, 10978444, 27227908, 27227908, 10978444, 1637860, 72076, 512
Offset: 0

Views

Author

Dale Gerdemann, Apr 12 2015

Keywords

Comments

Related triangles may be found by varying the function f(x). If f(x) is a linear function, it can be parameterized as f(x) = a*x + b. With different values for a and b, the following triangles are obtained:
a\b 1.......2.......3.......4.......5.......6
The row sums of these, and similarly constructed number triangles, are shown in the following table:
a\b 1.......2.......3.......4.......5.......6.......7.......8.......9
The formula can be further generalized to: t(n,m) = f(m+s)*t(n-1,m) + f(n-s)*t(n,m-1), where f(x) = a*x + b. The following table specifies triangles with nonzero values for s (given after the slash).
a\b 0 1 2 3
-2 A130595/1
-1
0
With the absolute value, f(x) = |x|, one obtains A038221/3, A038234/4, A038247/5, A038260/6, A038273/7, A038286/8, A038299/9 (with value for s after the slash).
If f(x) = A000045(x) (Fibonacci) and s = 1, the result is A010048 (Fibonomial).
In the notation of Carlitz and Scoville, this is the triangle of generalized Eulerian numbers A(r, s | alpha, beta) with alpha = beta = 2. Also the array A(2,1,4) in the notation of Hwang et al. (see page 31). - Peter Bala, Dec 27 2019

Examples

			Array, t(n, k), begins as:
   1,    2,      4,        8,        16,         32,          64, ...;
   2,   12,     52,      196,       684,       2276,        7340, ...;
   4,   52,    416,     2644,     14680,      74652,      357328, ...;
   8,  196,   2644,    26440,    220280,    1623964,    10978444, ...;
  16,  684,  14680,   220280,   2643360,   27227908,   251195000, ...;
  32, 2276,  74652,  1623964,  27227908,  381190712,  4677894984, ...;
  64, 7340, 357328, 10978444, 251195000, 4677894984, 74846319744, ...;
Triangle, T(n, k), begins as:
    1;
    2,     2;
    4,    12,      4;
    8,    52,     52,       8;
   16,   196,    416,     196,      16;
   32,   684,   2644,    2644,     684,      32;
   64,  2276,  14680,   26440,   14680,    2276,     64;
  128,  7340,  74652,  220280,  220280,   74652,   7340,   128;
  256, 23172, 357328, 1623964, 2643360, 1623964, 357328, 23172,   256;
		

Crossrefs

Programs

  • Magma
    A256890:= func< n,k | (&+[(-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n: j in [0..k]]) >;
    [A256890(n,k): k in [0..n], n in [0..10]]; // G. C. Greubel, Oct 18 2022
    
  • Mathematica
    Table[Sum[(-1)^(k-j)*Binomial[j+3, j] Binomial[n+4, k-j] (j+2)^n, {j,0,k}], {n,0, 9}, {k,0,n}]//Flatten (* Michael De Vlieger, Dec 27 2019 *)
  • PARI
    t(n,m) = if ((n<0) || (m<0), 0, if ((n==0) && (m==0), 1, (m+2)*t(n-1, m) + (n+2)*t(n, m-1)));
    tabl(nn) = {for (n=0, nn, for (k=0, n, print1(t(n-k, k), ", ");); print(););} \\ Michel Marcus, Apr 14 2015
    
  • SageMath
    def A256890(n,k): return sum((-1)^(k-j)*Binomial(j+3,j)*Binomial(n+4,k-j)*(j+2)^n for j in range(k+1))
    flatten([[A256890(n,k) for k in range(n+1)] for n in range(11)]) # G. C. Greubel, Oct 18 2022

Formula

T(n,k) = t(n-k, k); t(0,0) = 1, t(n,m) = 0 if n < 0 or m < 0 else t(n,m) = f(m)*t(n-1,m) + f(n)*t(n,m-1), where f(x) = x + 2.
Sum_{k=0..n} T(n, k) = A001715(n).
T(n,k) = Sum_{j = 0..k} (-1)^(k-j)*binomial(j+3,j)*binomial(n+4,k-j)*(j+2)^n. - Peter Bala, Dec 27 2019
Modified rule of Pascal: T(0,0) = 1, T(n,k) = 0 if k < 0 or k > n else T(n,k) = f(n-k) * T(n-1,k-1) + f(k) * T(n-1,k), where f(x) = x + 2. - Georg Fischer, Nov 11 2021
From G. C. Greubel, Oct 18 2022: (Start)
T(n, n-k) = T(n, k).
T(n, 0) = A000079(n). (End)

A196347 Triangle T(n, k) read by rows, T(n, k) = n!*binomial(n, k).

Original entry on oeis.org

1, 1, 1, 2, 4, 2, 6, 18, 18, 6, 24, 96, 144, 96, 24, 120, 600, 1200, 1200, 600, 120, 720, 4320, 10800, 14400, 10800, 4320, 720, 5040, 35280, 105840, 176400, 176400, 105840, 35280, 5040, 40320, 322560, 1128960, 2257920, 2822400, 2257920, 1128960, 322560, 40320
Offset: 0

Views

Author

Philippe Deléham, Oct 28 2011

Keywords

Comments

Unsigned version of A021012.
Equal to A136572*A007318.

Examples

			Triangle begins:
    1;
    1,   1;
    2,   4,    2;
    6,  18,   18,    6;
   24,  96,  144,   96,  24;
  120, 600, 1200, 1200, 600, 120;
  ...
		

Crossrefs

Programs

  • Magma
    /* As triangle */ [[Factorial(n)*Binomial(n, k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Sep 28 2015
  • Mathematica
    Table[n!*Binomial[n, j], {n, 0, 30}, {j, 0, n}] (* G. C. Greubel, Sep 27 2015 *)
  • Sage
    factorial(n)*binomial(n,k) # Danny Rorabaugh, Sep 27 2015
    

Formula

T(n,k) is given by (1,1,2,2,3,3,4,4,5,5,6,6,...) DELTA (1,1,2,2,3,3,4,4,5,5,6,6, ...) where DELTA is the operator defined in A084938.
Sum_{k>=0} T(m,k)*T(n,k) = (m+n)!.
T(2n,n) = A122747(n).
Sum_{k>=0} T(n,k)^2 = A010050(n) = (2n)!.
Sum_{k>=0} T(n,k)*x^k = A000007(n), A000142(n), A000165(n), A032031(n), A047053(n), A052562(n), A047058(n), A051188(n), A051189(n), A051232(n), A051262(n), A196258(n), A145448(n) for x = -1,0,1,2,3,4,5,6,7,8,9,10,11 respectively.
The row polynomials have the form (x + 1) o (x + 2) o ... o (x + n), where o denotes the black diamond multiplication operator of Dukes and White. See example E10 in the Bala link. - Peter Bala, Jan 18 2018

Extensions

Name exchanged with a formula by Peter Luschny, Feb 01 2015

A048176 Generalized Stirling number triangle of first kind.

Original entry on oeis.org

1, -10, 1, 200, -30, 1, -6000, 1100, -60, 1, 240000, -50000, 3500, -100, 1, -12000000, 2740000, -225000, 8500, -150, 1, 720000000, -176400000, 16240000, -735000, 17500, -210, 1, -50400000000, 13068000000, -1313200000, 67690000, -1960000, 32200, -280, 1, 4032000000000, -1095840000000
Offset: 1

Views

Author

Keywords

Comments

a(n,m)= R_n^m(a=0,b=10) in the notation of the given reference.
a(n,m) is a Jabotinsky matrix, i.e., the monic row polynomials E(n,x) := sum(a(n,m)*x^m,m=1..n) = product(x-10*j,j=0..n-1), n >= 1, E(0,x) := 1, are exponential convolution polynomials (see A039692 for the definition and a Knuth reference).
Also the Bell transform of the sequence (-1)^n*A051262(n) without column 0. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016

Examples

			{1}; {-10,1}; {200,-30,1}; {-6000,1100,-60,1}; ... E(3,x) = 200*x-30*x^2+x^3.
		

References

  • Mitrinovic, D. S.; Mitrinovic, R. S.; Tableaux d'une classe de nombres relies aux nombres de Stirling. Univ. Beograd. Pubi. Elektrotehn. Fak. Ser. Mat. Fiz. No. 77 1962, 77 pp.

Crossrefs

First (m=1) (unsigned) column sequence is: A051262(n-1). Row sums (signed triangle): A049212(n-1)*(-1)^(n-1). Row sums (unsigned triangle): A045757(n). b=8: A051187, b=9: A051231.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    # Adds (1,0,0,0, ..) as column 0.
    BellMatrix(n -> (-1)^n*n!*10^n, 9); # Peter Luschny, Jan 28 2016
  • Mathematica
    rows = 9;
    t = Table[(-1)^n*n!*10^n, {n, 0, rows}];
    T[n_, k_] := BellY[n, k, t];
    Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018, after Peter Luschny *)

Formula

a(n, m) = a(n-1, m-1) - 10*(n-1)*a(n-1, m), n >= m >= 1; a(n, m) := 0, n

A196258 a(n) = 11^n*n!.

Original entry on oeis.org

1, 11, 242, 7986, 351384, 19326120, 1275523920, 98215341840, 8642950081920, 855652058110080, 94121726392108800, 11388728893445164800, 1503312213934761753600, 214973646592670930764800, 33105941575271323337779200
Offset: 0

Author

Philippe Deléham, Oct 27 2011

Keywords

Programs

Formula

a(n) = 11^n*n!.
E.g.f.: 1/(1-11*x).
From Amiram Eldar, Jun 25 2020: (Start)
Sum_{n>=0} 1/a(n) = e^(1/11).
Sum_{n>=0} (-1)^n/a(n) = e^(-1/11). (End)
Showing 1-4 of 4 results.