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

A039810 Matrix square of Stirling2 triangle A008277: 2-levels set partitions of [n] into k first-level subsets.

Original entry on oeis.org

1, 2, 1, 5, 6, 1, 15, 32, 12, 1, 52, 175, 110, 20, 1, 203, 1012, 945, 280, 30, 1, 877, 6230, 8092, 3465, 595, 42, 1, 4140, 40819, 70756, 40992, 10010, 1120, 56, 1, 21147, 283944, 638423, 479976, 156072, 24570, 1932, 72, 1, 115975, 2090424, 5971350, 5660615, 2350950, 487704, 53550, 3120, 90, 1
Offset: 1

Views

Author

Christian G. Bower, Feb 15 1999

Keywords

Comments

This triangle groups certain generalized Stirling numbers of the second kind A000558, A000559, ... They can also be interpreted in terms of trees of height 3 with n leaves and constraints on the order of the root.
From Peter Bala, Jul 19 2014: (Start)
The (n,k)-th entry in this table gives the number of double partitions of the set [n] = {1,2,...,n} into k blocks. To form a double partition of [n] we first write [n] as a disjoint union X_1 U...U X_k of k nonempty subsets (blocks) X_i of [n]. Then each block X_i is further partitioned into sub-blocks to give a double partition. For instance, {1,2,4} U {3,5} is a partition of [5] into 2 blocks and {{1,4},{2}} U {{3},{5}} is a refinement of this partition to a double partition of [5] into 2 blocks (and 4 sub-blocks).
Compare the above interpretation for the (n,k)-th entry of this table with the interpretation of the (n,k)-th entry of A013609 (the square of Pascal's triangle but with the rows read in reverse order) as counting the pairs (X,Y) of subsets of [n] such that |Y| = k and X is contained in Y. (End)
Also the Bell transform of the shifted Bell numbers B(n+1) without column 0. For the definition of the Bell transform see A264428. - Peter Luschny, Jan 28 2016
T(n,k) is the number of partitions of an n-set into colored blocks, such that exactly k colors are used and the colors are introduced in increasing order. T(3,2) = 6: 1a|23b, 13a|2b, 12a|3b, 1a|2a|3b, 1a|2b|3a, 1a|2b|3b. - Alois P. Heinz, Aug 27 2019

Examples

			Triangle begins:
      k = 1    2    3    4    5          sum
  n
  1       1                                1
  2       2    1                           3
  3       5    6    1                     12
  4      15   32   12    1                60
  5      52  175  110   20    1          358
Matrix multiplication Stirling2 * Stirling2:
                  1  0  0  0
                  1  1  0  0
                  1  3  1  0
                  1  7  6  1
.
  1  0  0  0      1  0  0  0
  1  1  0  0      2  1  0  0
  1  3  1  0      5  6  1  0
  1  7  6  1     15 32 12  1
From _Peter Bala_, Jul 19 2014: (Start)
T(5,2) = 175: A 5-set can be partitioned into 2 blocks as either a union of a 3-set and a 2-set or as a union of a 4-set and a singleton set.
In the first case there are 10 ways of partitioning a 5-set into a 3-set and a 2-set. Each 3-set can be further partitioned into sub-blocks in Bell(3) = 5 ways and each 2-set can be further partitioned into sub-blocks in Bell(2) = 2 ways. So altogether we obtain 10*5*2 = 100 double partitions of this type.
In the second case, there are 5 ways of partitioning a 5-set into a 4-set and a 1-set. Each 4-set can be further partitioned in Bell(4) = 15 ways and each 1-set can be further partitioned in Bell(1) = 1 way. So altogether we obtain 5*15*1 = 75 double partitions of this type.
Hence, in total, T(5,2) = 100 + 75 = 175. (End)
		

Crossrefs

Cf. A039811, A039814, A039813 (other products of Stirling matrices).
T(n, 1) = A000110(n) (first column) (Bell numbers).
T(n, 2) = A000558(n) 2-levels set partitions with 2 first-level classes.
T(n, n-1) = A002378(n-1) = n*(n-1) = 2*C(n,2) = set-partitions into (n-2) singletons and one of the two possible set partitions of [2].
Sum is A000258(n), 2-levels set partitions.
Another version with offset 0: A130191.
Horizontal mirror triangle is A046817.
T(2n,n) gives A321712.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    # Adds (1,0,0,0, ..) as column 0.
    BellMatrix(n -> combinat:-bell(n+1), 10); # Peter Luschny, Jan 28 2016
  • Mathematica
    Flatten[Table[Sum[StirlingS2[n,i]*StirlingS2[i,k],{i,k,n}],{n,1,10},{k,1,n}]] (* Indranil Ghosh, Feb 22 2017 *)
    rows = 10;
    t = Table[BellB[n+1], {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 *)
  • PARI
    T(n, k) = sum(j=0, n, stirling(n, j, 2)*stirling(j, k, 2)); \\ Seiichi Manyama, Feb 13 2022

Formula

S2 = A008277 (Stirling numbers of the second kind).
T = (S2)^2.
T(n,k) = Sum_{i=k..n} S2(n,i) * S2(i,k).
E.g.f. of k-th column: (exp(exp(x)-1)-1)^k/k!. [corrected by Seiichi Manyama, Feb 12 2022]
From Peter Bala, Jul 19 2014: (Start)
T(n,k) = Sum_{disjoint unions X_1 U...U X_k = [n]} Bell(|X_1|)*...*Bell(|X_k|), where Bell(n) = A000110(n).
Recurrence equation: T(n+1,k+1) = Sum_{j = k..n} Bell(n+1-j)*binomial(n,j)* T(j,k).
Row sums [1,3,12,60,358,...] = A000258. (End)

Extensions

Definition and interpretation edited by Olivier Gérard, Jul 31 2011

A000558 Generalized Stirling numbers of second kind.

Original entry on oeis.org

1, 6, 32, 175, 1012, 6230, 40819, 283944, 2090424, 16235417, 132609666, 1135846062, 10175352709, 95108406130, 925496853980, 9357279554071, 98118527430960, 1065259283215810, 11956366813630835, 138539436100687988, 1655071323662574756, 20361556640795422729
Offset: 2

Views

Author

Keywords

Comments

From Olivier Gérard, Mar 25 2009: (Start)
a(n) is the number of hierarchical partitions of a set of n elements into two second level classes : k>1 subsets of [n] are further grouped in two classes.
a(n) is equivalently the number of trees of uniform height 3 with n labeled leaves, and a root of order two. (End)

Examples

			From _Olivier Gérard_, Mar 25 2009: (Start)
a(2) = 1, since there is only one partition of {1,2} into two classes, and only one way to partition those classes.
a(4) = 32 = 7*1 + 6*3 + 1*7 since there are 7 ways of partitioning {1,2,3,4} into two classes (which cannot be grouped further), 6 ways of partitioning a set of 4 elements into three classes and three ways to partition three classes into two super-classes, etc. (End)
		

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Column k=2 of A130191.
Cf. A001861 for the related bicolor set partitions. - Olivier Gérard, Mar 25 2009

Programs

  • Mathematica
    nn = 22; t = Range[0, nn]! CoefficientList[Series[1/2*(Exp[Exp[x] - 1] - 1)^2, {x, 0, nn}], x]; Drop[t, 2] (* T. D. Noe, Aug 10 2012 *)
    a[n_] := Sum[StirlingS2[n, k] (2^(k-1)-1), {k, 0, n}];
    a /@ Range[2, 100] (* Jean-François Alcover, Mar 30 2021 *)

Formula

E.g.f.: (1/2) * (exp(exp(x) - 1) - 1)^2. - Vladeta Jovovic, Sep 28 2003
a(n) = Sum_{k=0..n} Stirling2(n,k) * Stirling2(k,2). - Olivier Gérard, Mar 25 2009
a(n) = Sum_{k=1..n-1} binomial(n-1,k) * Bell(k) * Bell(n-k). - Ilya Gutkovskiy, Feb 15 2021

Extensions

More terms from David W. Wilson, Jan 13 2000

A130191 Square of the Stirling2 matrix A048993.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 5, 6, 1, 0, 15, 32, 12, 1, 0, 52, 175, 110, 20, 1, 0, 203, 1012, 945, 280, 30, 1, 0, 877, 6230, 8092, 3465, 595, 42, 1, 0, 4140, 40819, 70756, 40992, 10010, 1120, 56, 1, 0, 21147, 283944, 638423, 479976, 156072, 24570, 1932, 72, 1
Offset: 0

Views

Author

Wolfdieter Lang, Jun 01 2007

Keywords

Comments

Without row n=0 and column k=0 this is triangle A039810.
This is an associated Sheffer matrix with e.g.f. of the m-th column ((exp(f(x))-1)^m)/m! with f(x)=:exp(x)-1.
The triangle is also called the exponential Riordan array [1, exp(exp(x)-1)]. - Peter Luschny, Apr 19 2015
Also the Bell transform of shifted Bell numbers A000110(n+1). For the definition of the Bell transform see A264428. - Peter Luschny, Jan 27 2016

Examples

			Triangle starts:
  1;
  0,   1;
  0,   2,    1;
  0,   5,    6,    1;
  0,  15,   32,   12,    1;
  0,  52,  175,  110,   20,   1;
  0, 203, 1012,  945,  280,  30,  1;
  0, 877, 6230, 8092, 3465, 595, 42, 1;
		

Crossrefs

Columns k=0..3 give A000007, A000110 (for n > 0), A000558, A000559.
Row sums: A000258.
Alternating row sums: A130410.
T(2n,n) gives A321712.
Cf. A039810 (another version), A048993.

Programs

  • Maple
    # The function BellMatrix is defined in A264428.
    BellMatrix(n -> combinat:-bell(n+1), 9); # Peter Luschny, Jan 27 2016
  • Mathematica
    BellMatrix[f_Function, len_] := With[{t = Array[f, len, 0]}, Table[BellY[n, k, t], {n, 0, len - 1}, {k, 0, len - 1}]];
    rows = 10;
    M = BellMatrix[BellB[# + 1]&, rows];
    Table[M[[n, k]], {n, 1, rows}, {k, 1, n}] // Flatten (* Jean-François Alcover, Jun 23 2018, after Peter Luschny *)
    a[n_, m_]:= Sum[StirlingS2[n, k]*StirlingS2[k, m], {k,m,n}]; Table[a[n, m], {n, 0, 100}, {m, 0, n}]//Flatten (* G. C. Greubel, Jul 10 2018 *)
  • PARI
    for(n=0, 9, for(k=0, n, print1(sum(j=k, n, stirling(n, j, 2)*stirling(j, k, 2)), ", "))) \\ G. C. Greubel, Jul 10 2018
  • Sage
    # uses[riordan_array from A256893]
    riordan_array(1, exp(exp(x) - 1), 8, exp=true) # Peter Luschny, Apr 19 2015
    

Formula

a(n,k) = Sum_{j=k..n} S2(n,j) * S2(j,k), n>=k>=0.
E.g.f. row polynomials with argument x: exp(x*f(f(z))).
E.g.f. column k: ((exp(exp(x) - 1) - 1)^k)/k!.

A341588 E.g.f.: -log(1 + log(1 - x))^3 / 6.

Original entry on oeis.org

1, 12, 130, 1485, 18508, 253400, 3805723, 62437500, 1113510409, 21479997957, 446094038806, 9930796412082, 236037249893092, 5968192832899412, 160007282538148508, 4534905316824903144, 135500246340709682692, 4257646241716404353684, 140366073694357927723936, 4845119946789226304526392
Offset: 3

Views

Author

Ilya Gutkovskiy, Feb 15 2021

Keywords

Crossrefs

Programs

  • Mathematica
    nmax = 22; CoefficientList[Series[-Log[1 + Log[1 - x]]^3/6, {x, 0, nmax}], x] Range[0, nmax]! // Drop[#, 3] &
    Table[Sum[Abs[StirlingS1[n, k] StirlingS1[k, 3]], {k, 3, n}], {n, 3, 22}]

Formula

a(n) = Sum_{k=3..n} |Stirling1(n, k) * Stirling1(k, 3)|.
a(n) ~ (n-1)! * log(n)^2 / (2 * (1 - exp(-1))^n) * (1 + (2*gamma - 2*log(exp(1) - 1)) / log(n) + (gamma^2 - Pi^2/6 - 2*log(exp(1) - 1)*gamma + log(exp(1)-1)^2) / log(n)^2). - Vaclav Kotesovec, Jun 04 2022

A046817 Triangle of generalized Stirling numbers of 2nd kind.

Original entry on oeis.org

1, 1, 2, 1, 6, 5, 1, 12, 32, 15, 1, 20, 110, 175, 52, 1, 30, 280, 945, 1012, 203, 1, 42, 595, 3465, 8092, 6230, 877, 1, 56, 1120, 10010, 40992, 70756, 40819, 4140, 1, 72, 1932, 24570, 156072, 479976, 638423, 283944, 21147, 1, 90, 3120, 53550, 487704, 2350950, 5660615, 5971350
Offset: 0

Views

Author

Keywords

Examples

			Triangle begins:
      k = 0    1    2    3    4          sum
n
1         1                                1
2         1    2                           3
3         1    6    5                     12
4         1   12   32   15                60
5         1   20  110  175   52          358
		

Crossrefs

Diagonals give A000558, A000559, A000110, A002378, etc.
Row sums give A000258.
Horizontal mirror triangle is A039810 (matrix square of Stirling2).

Programs

  • Mathematica
    a[n_, k_] = Sum[StirlingS2[n, i]*StirlingS2[i, k], {i, k, n}]; Flatten[Table[a[n, k], {n, 1, 10}, {k, n, 1, -1}]][[1 ;; 53]]  (* Jean-François Alcover, Apr 26 2011 *)

Formula

a(n, k) = Sum_{i=k..n} S2(n, i)*S2(i, k).
E.g.f.: exp(exp(exp(x*y)-1)-1)^(1/y). - Vladeta Jovovic, Dec 14 2003

Extensions

More terms from David W. Wilson, Jan 13 2000

A383208 Expansion of e.g.f. (exp(f(x)) - 1)^3 / 6, where f(x) = (exp(2*x) - 1)/2.

Original entry on oeis.org

0, 0, 0, 1, 18, 245, 3120, 39697, 517790, 6999785, 98520060, 1445923149, 22129416210, 352932509085, 5859167661256, 101122879922313, 1811960841148774, 33662625853200337, 647550189266734452, 12881675626292023173, 264677402162135670554, 5610552395871699336453
Offset: 0

Views

Author

Seiichi Manyama, Apr 19 2025

Keywords

Crossrefs

Column k=3 of A383206.
Cf. A000559.

Programs

  • PARI
    a(n) = sum(k=3, n, 2^(n-k)*stirling(n, k, 2)*stirling(k, 3, 2));

Formula

a(n) = Sum_{k=3..n} 2^(n-k) * Stirling2(n,k) * Stirling2(k,3).
Showing 1-6 of 6 results.