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-10 of 13 results. Next

A028575 Row sums of triangle A011801.

Original entry on oeis.org

1, 5, 49, 721, 14177, 349141, 10334689, 357361985, 14137664833, 629779342213, 31195027543505, 1700812505769169, 101218448336028193, 6528869281965115541, 453720852957751220353, 33796334125623555379969, 2686138908337714715560577, 226908450494953996837748869
Offset: 1

Views

Author

Keywords

Crossrefs

Sequences with e.g.f. exp(1-(1-m*x)^(1/m)) - 1: A000012 (m=1), A001515 (m=2), A015735 (m=3), A016036 (m=4), this sequence (m=5), A028844 (m=6).
Cf. A011801.

Programs

  • Magma
    R:=PowerSeriesRing(Rationals(), 30); Coefficients(R!(Laplace( Exp(1-(1-5*x)^(1/5)) - 1 ))); // G. C. Greubel, Oct 02 2023
    
  • Mathematica
    Rest[With[{nn=20},CoefficientList[Series[Exp[1-(1-5x)^(1/5)]-1, {x,0,nn}], x] Range[0,nn]!]] (* Harvey P. Dale, Aug 02 2016 *)
  • SageMath
    def A028575_list(prec):
        P. = PowerSeriesRing(QQ, prec)
        return P( exp(1-(1-5*x)^(1/5)) -1 ).egf_to_ogf().list()
    a=A028575_list(40); a[1:] # G. C. Greubel, Oct 02 2023

Formula

E.g.f.: exp(1 - (1-5*x)^(1/5)) - 1.
a(n) = D^n(exp(x)) evaluated at x = 0, where D is the operator 1/(1-x)^4*d/dx. Cf. A001515, A015735 and A016036. - Peter Bala, Nov 25 2011
D-finite with recurrence: a(n) -20*(n-3)*a(n-1) +30*(5*n^2-35*n +62)*a(n-2) -100*(n-4)*(5*n^2-40*n+81)*a(n-3) +(5*n-22)*(5*n-21)*(5*n-24)*(5*n-23)*a(n-4) -a(n-5) = 0. - R. J. Mathar, Jan 28 2020
From Seiichi Manyama, Jan 20 2025: (Start)
a(n) = Sum_{k=0..n} (-1)^k * 5^(n-k) * |Stirling1(n,k)| * A000587(k).
a(n) = e * (-5)^n * n! * Sum_{k>=0} (-1)^k * binomial(k/5,n)/k!. (End)

A144267 Partition number array, called M32(-4), related to A011801(n,m)= |S2(-4;n,m)| ( generalized Stirling triangle).

Original entry on oeis.org

1, 4, 1, 36, 12, 1, 504, 144, 48, 24, 1, 9576, 2520, 1440, 360, 240, 40, 1, 229824, 57456, 30240, 12960, 7560, 8640, 960, 720, 720, 60, 1, 6664896, 1608768, 804384, 635040, 201096, 211680, 90720, 60480, 17640, 30240, 6720, 1260, 1680, 84, 1, 226606464, 53319168
Offset: 1

Views

Author

Wolfdieter Lang, Oct 09 2008

Keywords

Comments

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k)=:M32(-4;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42, ...].
a(n,k) enumerates special unordered forests related to the k-th partition of n in the A-St order. The k-th partition of n is given by the exponents enk =(e(n,k,1),...,e(n,k,n)) of 1,2,...n. The number of parts is m = sum(e(n,k,j),j=1..n). The special (enk)-forest is composed of m rooted increasing (r+3)-ary trees if the outdegree is r >= 0.
If M32(-4;n,k) is summed over those k with fixed number of parts m one obtains triangle A011801(n,m)= |S2(-4;n,m)|, a generalization of Stirling numbers of the second kind. For S2(K;n,m), K from the integers, see the reference under A035342.

Examples

			a(4,3)=48. The relevant partition of 4 is (2^2). The 48 unordered (0,2,0,0)-forests are composed of the following 2 rooted increasing trees 1--2,3--4; 1--3,2--4 and 1--4,2--3. The trees are quaternary because r=1 vertices are quaternary (4-ary) and for the leaves (r=0) the arity does not matter. Each of the three differently labeled forests comes therefore in 4^2=16 versions due to the two quaternary root vertices.
		

Crossrefs

Cf. A143173 (M32(-3) array), A144268 (M32(-5) array).

Formula

a(n,k) = (n!/product(e(n,k,j)!*j!^(e(n,k,j),j=1..n))*product(|S2(-4,j,1)|^e(n,k,j),j=1..n) = M3(n,k)*product(|S2(-4,j,1)|^e(n,k,j),j=1..n), with |S2(-4,n,1)|= A008546(n-1) = (5*n-6)(!^5) (5-factorials) for n>=2 and 1 if n=1 and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n. Exponents 0 can be omitted due to 0!=1. M3(n,k):= A036040(n,k), k=1..p(n), p(n):= A000041(n).

A144284 Partition number array, called M32hat(-4)= 'M32(-4)/M3'= 'A144267/A036040', related to A011801(n,m)= |S2(-4;n,m)| (generalized Stirling triangle).

Original entry on oeis.org

1, 4, 1, 36, 4, 1, 504, 36, 16, 4, 1, 9576, 504, 144, 36, 16, 4, 1, 229824, 9576, 2016, 1296, 504, 144, 64, 36, 16, 4, 1, 6664896, 229824, 38304, 18144, 9576, 2016, 1296, 576, 504, 144, 64, 36, 16, 4, 1, 226606464, 6664896, 919296, 344736, 254016, 229824, 38304, 18144
Offset: 1

Views

Author

Wolfdieter Lang Oct 09 2008

Keywords

Comments

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k) =: M32hat(-4;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42,...].
If M32hat(-4;n,k) is summed over those k with fixed number of parts m one obtains triangle S2hat(-4):= A144285(n,m).

Examples

			a(4,3)= 16 = |S2(-4,2,1)|^2. The relevant partition of 4 is (2^2).
		

Crossrefs

A144279 (M32hat(-3) array). A144341 (M32hat(-5) array)

Formula

a(n,k)= product(|S2(-4,j,1)|^e(n,k,j),j=1..n) with |S2(-4,n,1)|= A008546(n-1) = (5*n-6)(!^5) (5-factorials) for n>=2 and 1 if n=1 and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n.
Formally a(n,k)= 'M32(-4)/M3' = 'A144267/A036040' (elementwise division of arrays).

A144341 Partition number array, called M32hat(-5)= 'M32(-5)/M3'= 'A144268/A036040', related to A011801(n,m)= |S2(-4;n,m)| (generalized Stirling triangle).

Original entry on oeis.org

1, 5, 1, 55, 5, 1, 935, 55, 25, 5, 1, 21505, 935, 275, 55, 25, 5, 1, 623645, 21505, 4675, 3025, 935, 275, 125, 55, 25, 5, 1, 21827575, 623645, 107525, 51425, 21505, 4675, 3025, 1375, 935, 275, 125, 55, 25, 5, 1, 894930575, 21827575, 3118225, 1182775, 874225, 623645
Offset: 1

Views

Author

Wolfdieter Lang Oct 09 2008

Keywords

Comments

Each partition of n, ordered as in Abramowitz-Stegun (A-St order; for the reference see A134278), is mapped to a nonnegative integer a(n,k) =: M32hat(-5;n,k) with the k-th partition of n in A-St order.
The sequence of row lengths is A000041 (partition numbers) [1, 2, 3, 5, 7, 11, 15, 22, 30, 42,...].
If M32hat(-5;n,k) is summed over those k with fixed number of parts m one obtains triangle S2hat(-5):= A144342(n,m).

Examples

			a(4,3)= 25 = |S2(-5,2,1)|^2. The relevant partition of 4 is (2^2).
		

Crossrefs

A144284 (M32hat(-4) array).

Formula

a(n,k)= product(|S2(-5,j,1)|^e(n,k,j),j=1..n) with |S2(-5,n,1)|= A008543(n-1) = (6*n-7)(!^6) (6-factorials) for n>=2 and 1 if n=1 and the exponent e(n,k,j) of j in the k-th partition of n in the A-St ordering of the partitions of n.
Formally a(n,k)= 'M32(-5)/M3' = 'A144268/A036040' (elementwise division of arrays).

A144347 Second column (m=2) of triangle S2p(-4) = A011801.

Original entry on oeis.org

1, 12, 192, 3960, 100656, 3048192, 107255232, 4302305280, 193836779136, 9693022090752, 532784148728832, 31930395433896960, 2072320885985366016, 144803002560595771392, 10838696766561262190592, 865256088684973495910400, 73383436891415208719056896
Offset: 0

Views

Author

Wolfdieter Lang, Oct 09 2008

Keywords

Crossrefs

Cf. A011801, A008546(n-1) (m=1 column), A144348 (m=3 column).

Formula

a(n) = A011801(n+2,2), n>=0.

A144348 Third column (m=3) of triangle S2p(-4) = A011801.

Original entry on oeis.org

1, 24, 600, 17160, 563976, 21095424, 887785920, 41589313920, 2148534533376, 121416817826304, 7453542764828160, 494050046853242880, 35173674025152638976, 2677307532371585777664, 216991376759173367070720, 18658548982937189595709440, 1696677198406950703545778176
Offset: 0

Views

Author

Wolfdieter Lang, Oct 09 2008

Keywords

Crossrefs

Cf. A011801, A144347 (m=2 column).

Formula

a(n) = A011801(n+3,3), n>=0.

A105278 Triangle read by rows: T(n,k) = binomial(n,k)*(n-1)!/(k-1)!.

Original entry on oeis.org

1, 2, 1, 6, 6, 1, 24, 36, 12, 1, 120, 240, 120, 20, 1, 720, 1800, 1200, 300, 30, 1, 5040, 15120, 12600, 4200, 630, 42, 1, 40320, 141120, 141120, 58800, 11760, 1176, 56, 1, 362880, 1451520, 1693440, 846720, 211680, 28224, 2016, 72, 1, 3628800, 16329600
Offset: 1

Views

Author

Miklos Kristof, Apr 25 2005

Keywords

Comments

T(n,k) is the number of partially ordered sets (posets) on n elements that consist entirely of k chains. For example, T(4, 3)=12 since there are exactly 12 posets on {a,b,c,d} that consist entirely of 3 chains. Letting ab denote a<=b and using a slash "/" to separate chains, the 12 posets can be given by a/b/cd, a/b/dc, a/c/bd, a/c/db, a/d/bc, a/d/cb, b/c/ad, b/c/da, b/d/ac, b/d/ca, c/d/ab and c/d/ba, where the listing of the chains is arbitrary (e.g., a/b/cd = a/cd/b =...cd/b/a). - Dennis P. Walsh, Feb 22 2007
Also the matrix product |S1|.S2 of Stirling numbers of both kinds.
This Lah triangle is a lower triangular matrix of the Jabotinsky type. See the column e.g.f. and the D. E. Knuth reference given in A008297. - Wolfdieter Lang, Jun 29 2007
The infinitesimal matrix generator of this matrix is given in A132710. See A111596 for an interpretation in terms of circular binary words and generalized factorials. - Tom Copeland, Nov 22 2007
Three combinatorial interpretations: T(n,k) is (1) the number of ways to split [n] = {1,...,n} into a collection of k nonempty lists ("partitions into sets of lists"), (2) the number of ways to split [n] into an ordered collection of n+1-k nonempty sets that are noncrossing ("partitions into lists of noncrossing sets"), (3) the number of Dyck n-paths with n+1-k peaks labeled 1,2,...,n+1-k in some order. - David Callan, Jul 25 2008
Given matrices A and B with A(n,k) = T(n,k)*a(n-k) and B(n,k) = T(n,k)*b(n-k), then A*B = D where D(n,k) = T(n,k)*[a(.)+b(.)]^(n-k), umbrally. - Tom Copeland, Aug 21 2008
An e.g.f. for the row polynomials of A(n,k) = T(n,k)*a(n-k) is exp[a(.)* D_x * x^2] exp(x*t) = exp(x*t) exp[(.)!*Lag(.,-x*t,1)*a(.)*x], umbrally, where [(.)! Lag(.,x,1)]^n = n! Lag(n,x,1) is a normalized Laguerre polynomial of order 1. - Tom Copeland, Aug 29 2008
Triangle of coefficients from the Bell polynomial of the second kind for f = 1/(1-x). B(n,k){x1,x2,x3,...} = B(n,k){1/(1-x)^2,...,(j-1)!/(1-x)^j,...} = T(n,k)/(1-x)^(n+k). - Vladimir Kruchinin, Mar 04 2011
The triangle, with the row and column offset taken as 0, is the generalized Riordan array (exp(x), x) with respect to the sequence n!*(n+1)! as defined by Wang and Wang (the generalized Riordan array (exp(x), x) with respect to the sequence n! is Pascal's triangle A007318, and with respect to the sequence n!^2 is A021009 unsigned). - Peter Bala, Aug 15 2013
For a relation to loop integrals in QCD, see p. 33 of Gopakumar and Gross and Blaizot and Nowak. - Tom Copeland, Jan 18 2016
Also the Bell transform of (n+1)!. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016
Also the number of k-dimensional flats of the n-dimensional Shi arrangement. - Shuhei Tsujie, Apr 26 2019
The numbers T(n,k) appear as coefficients when expanding the rising factorials (x)^k = x(x+1)...(x+k-1) in the basis of falling factorials (x)k = x(x-1)...(x-k+1). Specifically, (x)^n = Sum{k=1..n} T(n,k) (x)k. - _Jeremy L. Martin, Apr 21 2021

Examples

			T(1,1) = C(1,1)*0!/0! = 1,
T(2,1) = C(2,1)*1!/0! = 2,
T(2,2) = C(2,2)*1!/1! = 1,
T(3,1) = C(3,1)*2!/0! = 6,
T(3,2) = C(3,2)*2!/1! = 6,
T(3,3) = C(3,3)*2!/2! = 1,
Sheffer a-sequence recurrence: T(6,2)= 1800 = (6/3)*120 + 6*240.
B(n,k) =
   1/(1-x)^2;
   2/(1-x)^3,  1/(1-x)^4;
   6/(1-x)^4,  6/(1-x)^5,  1/(1-x)^6;
  24/(1-x)^5, 36/(1-x)^6, 12/(1-x)^7, 1/(1-x)^8;
The triangle T(n,k) begins:
  n\k      1       2       3      4      5     6    7  8  9 ...
  1:       1
  2:       2       1
  3:       6       6       1
  4:      24      36      12      1
  5:     120     240     120     20      1
  6:     720    1800    1200    300     30     1
  7:    5040   15120   12600   4200    630    42    1
  8:   40320  141120  141120  58800  11760  1176   56  1
  9:  362880 1451520 1693440 846720 211680 28224 2016 72  1
  ...
Row n=10: [3628800, 16329600, 21772800, 12700800, 3810240, 635040, 60480, 3240, 90, 1]. - _Wolfdieter Lang_, Feb 01 2013
From _Peter Bala_, Feb 24 2025: (Start)
The array factorizes as an infinite product (read from right to left):
  /  1                \        /1             \^m /1           \^m /1           \^m
  |  2    1            |      | 0   1          |  |0  1         |  |1  1         |
  |  6    6   1        | = ...| 0   0   1      |  |0  1  1      |  |0  2  1      |
  | 24   36  12   1    |      | 0   0   1  1   |  |0  0  2  1   |  |0  0  3  1   |
  |120  240 120  20   1|      | 0   0   0  2  1|  |0  0  0  3  1|  |0  0  0  4  1|
  |...                 |      |...             |  |...          |  |...          |
where m = 2. Cf. A008277 (m = 1), A035342 (m = 3), A035469 (m = 4), A049029 (m = 5) A049385 (m = 6), A092082 (m = 7), A132056 (m = 8), A223511 - A223522 (m = 9 through 20), A001497 (m = -1), A004747 (m = -2), A000369 (m = -3), A011801 (m = -4), A013988 (m = -5). (End)
		

Crossrefs

Triangle of Lah numbers (A008297) unsigned.
Cf. A111596 (signed triangle with extra n=0 row and m=0 column).
Cf. A130561 (for a natural refinement).
Cf. A094638 (for differential operator representation).
Cf. A248045 (central terms), A002868 (row maxima).
Cf, A059110.
Cf. A089231 (triangle with mirrored rows).
Cf. A271703 (triangle with extra n=0 row and m=0 column).

Programs

  • GAP
    Flat(List([1..10],n->List([1..n],k->Binomial(n,k)*Factorial(n-1)/Factorial(k-1)))); # Muniru A Asiru, Jul 25 2018
  • Haskell
    a105278 n k = a105278_tabl !! (n-1) !! (k-1)
    a105278_row n = a105278_tabl !! (n-1)
    a105278_tabl = [1] : f [1] 2 where
       f xs i = ys : f ys (i + 1) where
         ys = zipWith (+) ([0] ++ xs) (zipWith (*) [i, i + 1 ..] (xs ++ [0]))
    -- Reinhard Zumkeller, Sep 30 2014, Mar 18 2013
    
  • Magma
    /* As triangle */ [[Binomial(n,k)*Factorial(n-1)/Factorial(k-1): k in [1..n]]: n in [1.. 15]]; // Vincenzo Librandi, Oct 31 2014
    
  • Maple
    The triangle: for n from 1 to 13 do seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n) od;
    the sequence: seq(seq(binomial(n,k)*(n-1)!/(k-1)!,k=1..n),n=1..13);
    # The function BellMatrix is defined in A264428.
    # Adds (1, 0, 0, 0, ...) as column 0.
    BellMatrix(n -> (n+1)!, 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    nn = 9; a = x/(1 - x); f[list_] := Select[list, # > 0 &]; Flatten[Map[f, Drop[Range[0, nn]! CoefficientList[Series[Exp[y a], {x, 0, nn}], {x, y}], 1]]] (* Geoffrey Critzer, Dec 11 2011 *)
    nn = 9; Flatten[Table[(j - k)! Binomial[j, k] Binomial[j - 1, k - 1], {j, nn}, {k, j}]] (* Jan Mangaldan, Mar 15 2013 *)
    rows = 10;
    t = Range[rows]!;
    T[n_, k_] := BellY[n, k, t];
    Table[T[n, k], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
    T[n_, n_] := 1; T[n_, k_] /;0Oliver Seipel, Dec 06 2024 *)
  • Perl
    use ntheory ":all"; say join ", ", map { my $n=$; map { stirling($n,$,3) } 1..$n; } 1..9; # Dana Jacobsen, Mar 16 2017
    

Formula

T(n,k) = Sum_{m=n..k} |S1(n,m)|*S2(m,k), k>=n>=1, with Stirling triangles S2(n,m):=A048993 and S1(n,m):=A048994.
T(n,k) = C(n,k)*(n-1)!/(k-1)!.
Sum_{k=1..n} T(n,k) = A000262(n).
n*Sum_{k=1..n} T(n,k) = A103194(n) = Sum_{k=1..n} T(n,k)*k^2.
E.g.f. column k: (x^(k-1)/(1-x)^(k+1))/(k-1)!, k>=1.
Recurrence from Sheffer (here Jabotinsky) a-sequence [1,1,0,...] (see the W. Lang link under A006232): T(n,k)=(n/k)*T(n-1,m-1) + n*T(n-1,m). - Wolfdieter Lang, Jun 29 2007
The e.g.f. is, umbrally, exp[(.)!* L(.,-t,1)*x] = exp[t*x/(1-x)]/(1-x)^2 where L(n,t,1) = Sum_{k=0..n} T(n+1,k+1)*(-t)^k = Sum_{k=0..n} binomial(n+1,k+1)* (-t)^k / k! is the associated Laguerre polynomial of order 1. - Tom Copeland, Nov 17 2007
For this Lah triangle, the n-th row polynomial is given umbrally by
n! C(B.(x)+1+n,n) = (-1)^n C(-B.(x)-2,n), where C(x,n)=x!/(n!(x-n)!),
the binomial coefficient, and B_n(x)= exp(-x)(xd/dx)^n exp(x), the n-th Bell / Touchard / exponential polynomial (cf. A008277). E.g.,
2! C(-B.(-x)-2,2) = (-B.(x)-2)(-B.(x)-3) = B_2(x) + 5*B_1(x) + 6 = 6 + 6x + x^2.
n! C(B.(x)+1+n,n) = n! e^(-x) Sum_{j>=0} C(j+1+n,n)x^j/j! is a corresponding Dobinski relation. See the Copeland link for the relation to inverse Mellin transform. - Tom Copeland, Nov 21 2011
The row polynomials are given by D^n(exp(x*t)) evaluated at x = 0, where D is the operator (1+x)^2*d/dx. Cf. A008277 (D = (1+x)*d/dx), A035342 (D = (1+x)^3*d/dx), A035469 (D = (1+x)^4*d/dx) and A049029 (D = (1+x)^5*d/dx). - Peter Bala, Nov 25 2011
T(n,k) = Sum_{i=k..n} A130534(n-1,i-1)*A008277(i,k). - Reinhard Zumkeller, Mar 18 2013
Let E(x) = Sum_{n >= 0} x^n/(n!*(n+1)!). Then a generating function is exp(t)*E(x*t) = 1 + (2 + x)*t + (6 + 6*x + x^2)*t^2/(2!*3!) + (24 + 36*x + 12*x^2 + x^3)*t^3/(3!*4!) + ... . - Peter Bala, Aug 15 2013
P_n(x) = L_n(1+x) = n!*Lag_n(-(1+x);1), where P_n(x) are the row polynomials of A059110; L_n(x), the Lah polynomials of A105278; and Lag_n(x;1), the Laguerre polynomials of order 1. These relations follow from the relation between the iterated operator (x^2 D)^n and ((1+x)^2 D)^n with D = d/dx. - Tom Copeland, Jul 23 2018
Dividing each n-th diagonal by n!, where the main diagonal is n=1, generates the Narayana matrix A001263. - Tom Copeland, Sep 23 2020
T(n,k) = A089231(n,n-k). - Ron L.J. van den Burg, Dec 12 2021
T(n,k) = T(n-1,k-1) + (n+k-1)*T(n-1,k). - Bérénice Delcroix-Oger, Jun 25 2025

Extensions

Stirling comments and e.g.f.s from Wolfdieter Lang, Apr 11 2007

A008546 Quintuple factorial numbers: Product_{k = 0..n-1} (5*k + 4).

Original entry on oeis.org

1, 4, 36, 504, 9576, 229824, 6664896, 226606464, 8837652096, 388856692224, 19053977918976, 1028914807624704, 60705973649857536, 3885182313590882304, 268077579637770878976, 19837740893195045044224, 1567181530562408558493696, 131643248567242318913470464
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Crossrefs

Programs

  • GAP
    List([0..20], n-> Product([0..n-1], k-> 5*k+4 )); # G. C. Greubel, Aug 20 2019
  • Magma
    [1] cat [(&*[5*k+4: k in [0..n-1]]): n in [1..20]]; // G. C. Greubel, Aug 20 2019
    
  • Maple
    f:= n-> product(5*k+4, k=0..n-1);
  • Mathematica
    FoldList[Times, 1, 5Range[0, 20] + 4] (* Vincenzo Librandi, Jun 10 2013 *)
    CoefficientList[Series[(1 - 5x)^(-4/5), {x, 0, 20}], x] Range[0, 20]! (* Vaclav Kotesovec, Jan 28 2015 *)
    Table[5^n Pochhammer[4/5, n], {n, 0, 20}] (* G. C. Greubel, Aug 20 2019 *)
  • PARI
    vector(20, n, n--; prod(j=0,n-1, 5*j+4) ) \\ G. C. Greubel, Aug 20 2019
    
  • Sage
    [5^n*rising_factorial(4/5, n) for n in (0..20)] # G. C. Greubel, Aug 20 2019
    

Formula

a(n) = 4*A034301(n) = (5*n - 1)(!^5), n >= 1, with a(0) = 1.
a(n) = A011801(n + 1, 1) (first column of triangle).
a(n) ~ (sqrt(2*Pi)/Gamma(4/5))*n^(n + 3/10)*(5/e)^n*(1 + 1/(300*n) + ...). - Joe Keane (jgk(AT)jgk.org), Nov 24 2001
G.f.: 1/(1 - 4*x/(1 - 5*x/(1 - 9*x/(1 - 10*x/(1 - 14*x/(1 - 15*x/(1 - 19*x/(1 - 20*x/(1 - 24*x/(1 - ... (continued fraction). - Philippe Deléham, Jan 08 2012
a(n) = (-1)^n*Sum_{k = 0..n} 5^k*s(n + 1, n + 1 - k), where s(n, k) are the Stirling numbers of the first kind, A048994. - Mircea Merca, May 03 2012
G.f.: ( 1 - 1/Q(0) )/x where Q(k) = 1 - x*(5*k - 1)/(1 - x*(5*k + 5)/Q(k + 1) ); (continued fraction); e.g.f. (1 - 5*x)^(-4/5). - Sergei N. Gladkovskii, Mar 20 2013
G.f.: 1/x - G(0)/(2*x), where G(k) = 1 + 1/(1 - x*(5*k - 1)/(x*(5*k - 1) + 1/G(k + 1))); (continued fraction). - Sergei N. Gladkovskii, May 27 2013
a(n) = 5^n * Gamma(n + 4/5) / Gamma(4/5). - Vaclav Kotesovec, Jan 28 2015
a(n) + (-5*n + 1)*a(n - 1) = 0. - R. J. Mathar, Sep 04 2016
G.f.: 1/(1 - 4*x - 20*x^2/(1 - 14*x - 90*x^2/(1 - 24*x - 210*x^2/(1 - 34*x - 380*x^2/(1 - 44*x - 600*x^2/(1 - 54*x - 870*x^2/(1 - ...))))))) (Jacobi continued fraction). - Nikolaos Pantelidis, Feb 29 2020
Sum_{n>=0} 1/a(n) = 1 + (e/5)^(1/5)*(Gamma(4/5) - Gamma(4/5, 1/5)). - Amiram Eldar, Dec 19 2022

A004747 Triangle read by rows: the Bell transform of the triple factorial numbers A008544 without column 0.

Original entry on oeis.org

1, 2, 1, 10, 6, 1, 80, 52, 12, 1, 880, 600, 160, 20, 1, 12320, 8680, 2520, 380, 30, 1, 209440, 151200, 46480, 7840, 770, 42, 1, 4188800, 3082240, 987840, 179760, 20160, 1400, 56, 1, 96342400, 71998080, 23826880, 4583040, 562800, 45360, 2352, 72, 1
Offset: 1

Views

Author

Keywords

Comments

Previous name was: Triangle of numbers related to triangle A048966; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
T(n,m) = S2p(-2; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n, m) (Stirling 2nd kind). T(n,1)= A008544(n-1).
T(n,m), n>=m>=1, enumerates unordered n-vertex m-forests composed of m plane (aka ordered) increasing (rooted) trees where vertices of out-degree r>=0 come in r+1 different types (like an (r+1)-ary vertex). Proof from the e.g.f. of the first column Y(z) = 1 - (1-3*x)^(1/3) and the F. Bergeron et al. eq. (8) Y'(z)= phi(Y(z)), Y(0) = 0, with out-degree o.g.f. phi(w)=1/(1-w)^2. - Wolfdieter Lang, Oct 12 2007
Also the Bell transform of the triple factorial numbers A008544 which adds a first column (1,0,0 ...) on the left side of the triangle. For the definition of the Bell transform see A264428. See A051141 for the triple factorial numbers A032031 and A203412 for the triple factorial numbers A007559 as well as A039683 and A132062 for the case of double factorial numbers. - Peter Luschny, Dec 21 2015

Examples

			Triangle begins:
       1;
       2,      1;
      10,      6,     1;
      80,     52,    12,    1;
     880,    600,   160,   20,   1;
   12320,   8680,  2520,  380,  30,  1;
  209440, 151200, 46480, 7840, 770, 42, 1;
Tree combinatorics for T(3,2)=6: Consider first the unordered forest of m=2 plane trees with n=3 vertices, namely one vertex with out-degree r=0 (root) and two different trees with two vertices (one root with out-degree r=1 and a leaf with r=0). The 6 increasing labelings come then from the forest with rooted (x) trees x, o-x (1,(3,2)), (2,(3,1)) and (3,(2,1)) and similarly from the second forest x, x-o (1,(2,3)), (2,(1,3)) and (3,(1,2)).
		

Crossrefs

Cf. A015735 (row sums).
Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), this sequence (m=3), A000369 (m=4), A011801 (m=5), A013988 (m=6).

Programs

  • Magma
    function T(n,k) // T = A004747
      if k eq 0 then return 0;
      elif k eq n then return 1;
      else return (3*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
      end if;
    end function;
    [T(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Oct 03 2023
  • Maple
    T := (n, m) -> 3^n/m!*(1/3*m*GAMMA(n-1/3)*hypergeom([1-1/3*m, 2/3-1/3*m, 1/3-1/3*m], [2/3, 4/3-n], 1)/GAMMA(2/3)-1/6*m*(m-1)*GAMMA(n-2/3)*hypergeom( [1-1/3*m, 2/3-1/3*m, 4/3-1/3*m], [4/3, 5/3-n], 1)/Pi*3^(1/2)*GAMMA(2/3)):
    for n from 1 to 6 do seq(simplify(T(n,k)),k=1..n) od;
    # Karol A. Penson, Feb 06 2004
    # The function BellMatrix is defined in A264428.
    # Adds (1,0,0,0, ..) as column 0.
    BellMatrix(n -> mul(3*k+2, k=(0..n-1)), 9); # Peter Luschny, Jan 29 2016
  • Mathematica
    (* First program *)
    T[1,1]= 1; T[, 0]= 0; T[0, ]= 0; T[n_, m_]:= (3*(n-1)-m)*T[n-1, m]+T[n-1, m-1];
    Flatten[Table[T[n, m], {n,12}, {m,n}] ][[1 ;; 45]] (* Jean-François Alcover, Jun 16 2011, after recurrence *)
    (* Second program *)
    f[n_, m_]:= m/n Sum[Binomial[k, n-m-k] 3^k (-1)^(n-m-k) Binomial[n+k-1, n-1], {k, 0, n-m}]; Table[n! f[n, m]/(m! 3^(n-m)), {n,12}, {m,n}]//Flatten (* Michael De Vlieger, Dec 23 2015 *)
    (* Third program *)
    rows = 12;
    T[n_, m_]:= BellY[n, m, Table[Product[3k+2, {k, 0, j-1}], {j, 0, rows}]];
    Table[T[n, m], {n,rows}, {m,n}]//Flatten (* Jean-François Alcover, Jun 22 2018 *)
  • Sage
    # uses [bell_transform from A264428]
    triplefactorial = lambda n: prod(3*k+2 for k in (0..n-1))
    def A004747_row(n):
        trifact = [triplefactorial(k) for k in (0..n)]
        return bell_transform(n, trifact)
    [A004747_row(n) for n in (0..10)] # Peter Luschny, Dec 21 2015
    

Formula

T(n, m) = n!*A048966(n, m)/(m!*3^(n-m));
T(n+1, m) = (3*n-m)*T(n, m)+ T(n, m-1), for n >= m >= 1, with T(n, m) = 0, for n
E.g.f. of m-th column: ( 1 - (1-3*x)^(1/3) )^m/m!.
Sum_{k=1..n} T(n, k) = A015735(n).
For a formula expressed as special values of hypergeometric functions 3F2 see the Maple program below. - Karol A. Penson, Feb 06 2004
T(n,1) = A008544(n-1). - Peter Luschny, Dec 23 2015

Extensions

New name from Peter Luschny, Dec 21 2015

A013988 Triangle read by rows, the inverse Bell transform of n!*binomial(5,n) (without column 0).

Original entry on oeis.org

1, 5, 1, 55, 15, 1, 935, 295, 30, 1, 21505, 7425, 925, 50, 1, 623645, 229405, 32400, 2225, 75, 1, 21827575, 8423415, 1298605, 103600, 4550, 105, 1, 894930575, 358764175, 59069010, 5235405, 271950, 8330, 140, 1, 42061737025, 17398082625, 3016869625, 289426830, 16929255, 621810, 14070, 180, 1
Offset: 1

Keywords

Comments

Previous name was: Triangle of numbers related to triangle A049224; generalization of Stirling numbers of second kind A008277, Bessel triangle A001497.
T(n, m) = S2p(-5; n,m), a member of a sequence of triangles including S2p(-1; n,m) = A001497(n-1,m-1) (Bessel triangle) and ((-1)^(n-m))*S2p(1; n,m) = A008277(n,m) (Stirling 2nd kind). T(n, 1) = A008543(n-1).
For the definition of the Bell transform see A264428 and the link. - Peter Luschny, Jan 16 2016

Examples

			Triangle begins as:
          1;
          5,         1;
         55,        15,        1;
        935,       295,       30,       1;
      21505,      7425,      925,      50,      1;
     623645,    229405,    32400,    2225,     75,     1;
   21827575,   8423415,  1298605,  103600,   4550,   105,    1;
  894930575, 358764175, 59069010, 5235405, 271950,  8330,  140,   1;
		

Crossrefs

Cf. A028844 (row sums).
Triangles with the recurrence T(n,k) = (m*(n-1)-k)*T(n-1,k) + T(n-1,k-1): A010054 (m=1), A001497 (m=2), A004747 (m=3), A000369 (m=4), A011801 (m=5), this sequence (m=6).

Programs

  • Magma
    function T(n,k) // T = A013988
      if k eq 0 then return 0;
      elif k eq n then return 1;
      else return (6*(n-1)-k)*T(n-1,k) + T(n-1,k-1);
      end if;
    end function;
    [T(n,k): k in [1..n], n in [1..12]]; // G. C. Greubel, Oct 03 2023
  • Mathematica
    (* First program *)
    rows = 10;
    b[n_, m_] := BellY[n, m, Table[k! Binomial[5, k], {k, 0, rows}]];
    A = Table[b[n, m], {n, 1, rows}, {m, 1, rows}] // Inverse // Abs;
    A013988 = Table[A[[n, m]], {n, 1, rows}, {m, 1, n}] // Flatten (* Jean-François Alcover, Jun 22 2018 *)
    (* Second program *)
    T[n_, k_]:= T[n, k]= If[k==0, 0, If[k==n, 1, (6*(n-1) -k)*T[n-1,k] +T[n-1, k-1]]];
    Table[T[n,k], {n,12}, {k,n}]//Flatten (* G. C. Greubel, Oct 03 2023 *)
  • Sage
    # uses[inverse_bell_matrix from A264428]
    # Adds 1,0,0,0, ... as column 0 at the left side of the triangle.
    inverse_bell_matrix(lambda n: factorial(n)*binomial(5, n), 8) # Peter Luschny, Jan 16 2016
    

Formula

T(n, m) = n!*A049224(n, m)/(m!*6^(n-m));
T(n+1, m) = (6*n-m)*T(n, m) + T(n, m-1), for n >= m >= 1, with T(n, m) = 0, n
E.g.f. of m-th column: ((1 - (1-6*x)^(1/6))^m)/m!.
Sum_{k=1..n} T(n, k) = A028844(n).

Extensions

New name from Peter Luschny, Jan 16 2016
Showing 1-10 of 13 results. Next