A000258 -id:A000258 - cp's OEIS frontend

A318762 Number of permutations of a multiset whose multiplicities are the prime indices of n.

1, 1, 1, 2, 1, 3, 1, 6, 6, 4, 1, 12, 1, 5, 10, 24, 1, 30, 1, 20, 15, 6, 1, 60, 20, 7, 90, 30, 1, 60, 1, 120, 21, 8, 35, 180, 1, 9, 28, 120, 1, 105, 1, 42, 210, 10, 1, 360, 70, 140, 36, 56, 1, 630, 56, 210, 45, 11, 1, 420, 1, 12, 420, 720, 84, 168, 1, 72, 55

Offset: 1

Gus Wiseman, Sep 03 2018

This multiset is generally not the same as the multiset of prime indices of n. For example, the prime indices of 12 are {1,1,2}, while a multiset whose multiplicities are {1,1,2} is {1,1,2,3}.

			The a(12) = 12 permutations are (1123), (1132), (1213), (1231), (1312), (1321), (2113), (2131), (2311), (3112), (3121), (3211).

Alois P. Heinz, Table of n, a(n) for n = 1..20000

Cf. A000041, A000110, A000258, A000670, A005651, A008277, A008480, A056239, A112624, A112798, A124794, A215366, A300335, A305936.

a:= n-> (l-> add(i, i=l)!/mul(i!, i=l))(map(i->
       numtheory[pi](i[1])$i[2], ifactors(n)[2])):
seq(a(n), n=1..100);  # Alois P. Heinz, Sep 03 2018

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Total[primeMS[n]]!/Times@@Factorial/@primeMS[n],{n,100}]

sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
a(n)={if(n==1, 1, my(s=sig(n)); vecsum(s)!/prod(i=1, #s, s[i]!))}  \\ Andrew Howroyd, Dec 17 2018

If n = Product prime(x_i)^y_i is the prime factorization of n, then a(n) = (Sum x_i * y_i)! / Product (x_i!)^y_i.
a(n) = A008480(A181821(n)).
a(n) = A112624(n) * A124794(n). - Max Alekseyev, Oct 15 2023
Sum_{m in row n of A215366} a(m) = A005651(n).
Sum_{m in row n of A215366} a(m) * A008480(m) = A000670(n).
Sum_{m in row n of A215366} a(m) * A008480(m) / A001222(m)! = A000110(n).

A144150 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where the e.g.f. of column k is 1+g^(k+1)(x) with g = x-> exp(x)-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 3, 5, 1, 1, 1, 4, 12, 15, 1, 1, 1, 5, 22, 60, 52, 1, 1, 1, 6, 35, 154, 358, 203, 1, 1, 1, 7, 51, 315, 1304, 2471, 877, 1, 1, 1, 8, 70, 561, 3455, 12915, 19302, 4140, 1, 1, 1, 9, 92, 910, 7556, 44590, 146115, 167894, 21147, 1, 1, 1, 10, 117

Offset: 0

Alois P. Heinz, Sep 11 2008

A(n,k) is also the number of (k+1)-level labeled rooted trees with n leaves.
Number of ways to start with set {1,2,...,n} and then repeat k times: partition each set into subsets. - Alois P. Heinz, Aug 14 2015
Equivalently, A(n,k) is the number of length k+1 multichains from bottom to top in the set partition lattice of an n-set. - Geoffrey Critzer, Dec 05 2020

			Square array begins:
  1,  1,   1,    1,    1,    1,  ...
  1,  1,   1,    1,    1,    1,  ...
  1,  2,   3,    4,    5,    6,  ...
  1,  5,  12,   22,   35,   51,  ...
  1, 15,  60,  154,  315,  561,  ...
  1, 52, 358, 1304, 3455, 7556,  ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened
E. T. Bell, The Iterated Exponential Integers, Annals of Mathematics, 39(3) (1938), 539-557.
Pierpaolo Natalini and Paolo Emilio Ricci, Higher order Bell polynomials and the relevant integer sequences, in Appl. Anal. Discrete Math. 11 (2017), 327-339.
Pierpaolo Natalini and Paolo E. Ricci, Integer Sequences Connected with Extensions of the Bell Polynomials, Journal of Integer Sequences, 2017, Vol. 20, #17.10.2.
Ivar Henning Skau and Kai Forsberg Kristensen, An asymptotic Formula for the iterated exponential Bell Numbers, arXiv:1903.07979 [math.CO], 2019.
Ivar Henning Skau and Kai Forsberg Kristensen, Sets of iterated Partitions and the Bell iterated Exponential Integers, arXiv:1903.08379 [math.CO], 2019.
Index entries for sequences related to rooted trees

Columns k=0-10 give: A000012, A000110, A000258, A000307, A000357, A000405, A001669, A081624, A081629, A081697, A081740.
Rows n=0+1, 2-5 give: A000012, A000027, A000326, A005945, A005946.
First lower diagonal gives A139383.
First upper diagonal gives A346802.
Main diagonal gives A261280.
Cf. A000142, A111672, A290353.

g:= proc(p) local b; b:= proc(n) option remember; if n=0 then 1
      else (n-1)! *add(p(k)*b(n-k)/(k-1)!/(n-k)!, k=1..n) fi
    end end:
A:= (n,k)-> (g@@k)(1)(n):
seq(seq(A(n, d-n), n=0..d), d=0..12);
# second Maple program:
A:= proc(n, k) option remember; `if`(n=0 or k=0, 1,
      add(binomial(n-1, j-1)*A(j, k-1)*A(n-j, k), j=1..n))
    end:
seq(seq(A(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, Aug 14 2015
# third Maple program:
b:= proc(n, t, m) option remember; `if`(t=0, 1, `if`(n=0,
      b(m, t-1, 0), m*b(n-1, t, m)+b(n-1, t, m+1)))
    end:
A:= (n, k)-> b(n, k, 0):
seq(seq(A(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, Aug 04 2021

g[k_] := g[k] = Nest[Function[x, E^x - 1], x, k]; a[n_, k_] := SeriesCoefficient[1 + g[k + 1], {x, 0, n}]*n!; Table[a[n - k, k], {n, 0, 12}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 06 2013 *)

from sympy.core.cache import cacheit
from sympy import binomial
@cacheit
def A(n, k): return 1 if n==0 or k==0 else sum([binomial(n - 1, j - 1)*A(j, k - 1)*A(n - j, k) for j in range(1, n + 1)])
for n in range(51): print([A(k, n - k) for k in range(n + 1)]) # Indranil Ghosh, Aug 07 2017

E.g.f. of column k: 1 + g^(k+1)(x) with g = x-> exp(x)-1.

Column k+1 is Stirling transform of column k.

A122848 Exponential Riordan array (1, x(1+x/2)).

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 3, 1, 0, 0, 3, 6, 1, 0, 0, 0, 15, 10, 1, 0, 0, 0, 15, 45, 15, 1, 0, 0, 0, 0, 105, 105, 21, 1, 0, 0, 0, 0, 105, 420, 210, 28, 1, 0, 0, 0, 0, 0, 945, 1260, 378, 36, 1, 0, 0, 0, 0, 0, 945, 4725, 3150, 630, 45, 1, 0, 0, 0, 0, 0, 0, 10395, 17325, 6930, 990, 55, 1, 0, 0

Offset: 0

Paul Barry, Sep 14 2006

Entries are Bessel polynomial coefficients. Row sums are A000085. Diagonal sums are A122849. Inverse is A122850. Product of A007318 and A122848 gives A100862.
T(n,k) is the number of self-inverse permutations of {1,2,...,n} having exactly k cycles. - Geoffrey Critzer, May 08 2012
Bessel numbers of the second kind. For relations to the Hermite polynomials and the Catalan (A033184 and A009766) and Fibonacci (A011973, A098925, and A092865) matrices, see Yang and Qiao. - Tom Copeland, Dec 18 2013.
Also the inverse Bell transform of the double factorial of odd numbers Product_{k= 0..n-1} (2*k+1) (A001147). For the definition of the Bell transform see A264428 and for cross-references A265604. - Peter Luschny, Dec 31 2015

			Triangle begins:
    1
    0    1
    0    1    1
    0    0    3    1
    0    0    3    6    1
    0    0    0   15   10    1
    0    0    0   15   45   15    1
    0    0    0    0  105  105   21    1
    0    0    0    0  105  420  210   28    1
    0    0    0    0    0  945 1260  378   36    1
From _Gus Wiseman_, Jan 12 2021: (Start)
As noted above, a(n) is the number of set partitions of {1..n} into k singletons or pairs. This is also the number of set partitions of subsets of {1..n} into n - k pairs. In the first case, row n = 5 counts the following set partitions:
  {{1},{2,3},{4,5}}  {{1},{2},{3},{4,5}}  {{1},{2},{3},{4},{5}}
  {{1,2},{3},{4,5}}  {{1},{2},{3,4},{5}}
  {{1,2},{3,4},{5}}  {{1},{2,3},{4},{5}}
  {{1,2},{3,5},{4}}  {{1,2},{3},{4},{5}}
  {{1},{2,4},{3,5}}  {{1},{2},{3,5},{4}}
  {{1},{2,5},{3,4}}  {{1},{2,4},{3},{5}}
  {{1,3},{2},{4,5}}  {{1},{2,5},{3},{4}}
  {{1,3},{2,4},{5}}  {{1,3},{2},{4},{5}}
  {{1,3},{2,5},{4}}  {{1,4},{2},{3},{5}}
  {{1,4},{2},{3,5}}  {{1,5},{2},{3},{4}}
  {{1,4},{2,3},{5}}
  {{1,4},{2,5},{3}}
  {{1,5},{2},{3,4}}
  {{1,5},{2,3},{4}}
  {{1,5},{2,4},{3}}
In the second case, we have:
  {{1,2},{3,4}}  {{1,2}}  {}
  {{1,2},{3,5}}  {{1,3}}
  {{1,2},{4,5}}  {{1,4}}
  {{1,3},{2,4}}  {{1,5}}
  {{1,3},{2,5}}  {{2,3}}
  {{1,3},{4,5}}  {{2,4}}
  {{1,4},{2,3}}  {{2,5}}
  {{1,4},{2,5}}  {{3,4}}
  {{1,4},{3,5}}  {{3,5}}
  {{1,5},{2,3}}  {{4,5}}
  {{1,5},{2,4}}
  {{1,5},{3,4}}
  {{2,3},{4,5}}
  {{2,4},{3,5}}
  {{2,5},{3,4}}
(End)

G. C. Greubel, Table of n, a(n) for the first 50 rows, flattened
Peter Bala, Generalized Dobinski formulas
Richell O. Celeste, Roberto B. Corcino, and Ken Joffaniel M. Gonzales, Two Approaches to Normal Order Coefficients, Journal of Integer Sequences, Vol. 20 (2017), Article 17.3.5.
Tom Copeland, Infinitesimal Generators, the Pascal Pyramid, and the Witt and Virasoro Algebras
H. Han and S. Seo, Combinatorial proofs of inverse relations and log-concavity for Bessel numbers, Eur. J. Combinat. 29 (7) (2008) 1544-1554. [From _R. J. Mathar_, Mar 20 2009]
Robert S. Maier, Boson Operator Ordering Identities from Generalized Stirling and Eulerian Numbers, arXiv:2308.10332 [math.CO], 2023. See p. 18.
S. Yang and Z. Qiao, The Bessel Numbers and Bessel Matrices, Journal of Mathematical Research & Exposition, July, 2011, Vol. 31, No. 4, pp. 627-636. [From _Tom Copeland_, Dec 18 2013]

Cf. A001497, A008275, A039755, A096713, A104556, A113278, A130757.
Row sums are A000085.
Column sums are A001515.
Same as A049403 but with a first column k = 0.
The same set partitions counted by number of pairs are A100861.
Reversing rows gives A111924 (without column k = 0).
A047884 counts standard Young tableaux by size and greatest row length.
A238123 counts standard Young tableaux by size and least row length.
A320663/A339888 count unlabeled multiset partitions into singletons/pairs.
A322661 counts labeled covering half-loop-graphs.
A339742 counts factorizations into distinct primes or squarefree semiprimes.
Cf. A000110, A000258, A320732, A321729, A339741, A339887.

# The function BellMatrix is defined in A264428.
BellMatrix(n -> `if`(n<2,1,0), 9); # Peter Luschny, Jan 27 2016

t[n_, k_] := k!*Binomial[n, k]/((2 k - n)!*2^(n - k)); Table[ t[n, k], {n, 0, 11}, {k, 0, n}] // Flatten
(* Second program: *)
rows = 12;
t = Join[{1, 1}, Table[0, rows]];
T[n_, k_] := BellY[n, k, t];
Table[T[n, k], {n, 0, rows}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 23 2018,after Peter Luschny *)
sbs[{}]:={{}};sbs[set:{i_,_}]:=Join@@Function[s,(Prepend[#1,s]&)/@sbs[Complement[set,s]]]/@Cases[Subsets[set],{i}|{i,_}];
Table[Length[Select[sbs[Range[n]],Length[#]==k&]],{n,0,6},{k,0,n}] (* Gus Wiseman, Jan 12 2021 *)

{T(n,k)=if(2*kn, 0, n!/(2*k-n)!/(n-k)!*2^(k-n))} /* Michael Somos, Oct 03 2006 */

# uses[inverse_bell_transform from A265605]
multifact_2_1 = lambda n: prod(2*k + 1 for k in (0..n-1))
inverse_bell_matrix(multifact_2_1, 9) # Peter Luschny, Dec 31 2015

Number triangle T(n,k) = k!*C(n,k)/((2k-n)!*2^(n-k)).
T(n,k) = A001498(k,n-k). - Michael Somos, Oct 03 2006
E.g.f.: exp(y(x+x^2/2)). - Geoffrey Critzer, May 08 2012
Triangle equals the matrix product A008275*A039755. Equivalently, the n-th row polynomial R(n,x) is given by the Type B Dobinski formula R(n,x) = exp(-x/2)*Sum_{k>=0} P(n,2*k+1)*(x/2)^k/k!, where P(n,x) = x*(x-1)*...*(x-n+1) denotes the falling factorial polynomial. Cf. A113278. - Peter Bala, Jun 23 2014
From Daniel Checa, Aug 28 2022: (Start)
E.g.f. for the m-th column: (x^2/2+x)^m/m!.
T(n,k) = T(n-1,k-1) + (n-1)*T(n-2,k-1) for n>1 and k=1..n, T(0,0) = 1. (End)

A318393 Regular tetrangle where T(n,k,i) is the number of pairs of set partitions of {1,...,n} with meet of length k and join of length i.

Original entry on oeis.org

1, 1, 2, 1, 1, 6, 3, 8, 6, 1, 1, 14, 7, 48, 36, 6, 56, 44, 12, 1, 1, 30, 15, 200, 150, 25, 560, 440, 120, 10, 552, 440, 140, 20, 1, 1, 62, 31, 720, 540, 90, 3640, 2860, 780, 65, 8280, 6600, 2100, 300, 15, 7202, 5632, 1920, 340, 30, 1, 1, 126, 63, 2408, 1806

Offset: 1

Gus Wiseman, Aug 25 2018

			The T(3,3,1) = 8 pairs of set partitions:
  {{1},{2},{3}}  {{1,2,3}}
   {{1},{2,3}}  {{1,2},{3}}
   {{1},{2,3}}  {{1,3},{2}}
   {{1,2},{3}}  {{1},{2,3}}
   {{1,2},{3}}  {{1,3},{2}}
   {{1,3},{2}}  {{1},{2,3}}
   {{1,3},{2}}  {{1,2},{3}}
    {{1,2,3}}  {{1},{2},{3}}
Tetrangle begins:
   1   1     1       1            1
       2 1   6 3     14 7         30  15
             8 6 1   48 36 6      200 150 25
                     56 44 12 1   560 440 120 10
                                  552 440 140 20  1

Cf. A000110, A000258, A001247, A008277, A048994, A059849, A060639, A181939, A318389, A318390, A318391, A318392.

csm[s_]:=With[{c=Select[Tuples[Range[Length[s]],2],And[OrderedQ[#],UnsameQ@@#,Length[Intersection@@s[[#]]]>0]&]},If[c=={},s,csm[Union[Append[Delete[s,List/@c[[1]]],Union@@s[[c[[1]]]]]]]]];
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
spmeet[a_,b_]:=DeleteCases[Union@@Outer[Intersection,a,b,1],{}];spmeet[a_,b_,c__]:=spmeet[spmeet[a,b],c];
Table[Length[Select[Tuples[sps[Range[n]],2],And[Length[spmeet@@#]==k,Length[csm[Union@@#]]==j]&]],{n,6},{k,n},{j,k}]

A000307 Number of 4-level labeled rooted trees with n leaves.

Original entry on oeis.org

1, 1, 4, 22, 154, 1304, 12915, 146115, 1855570, 26097835, 402215465, 6734414075, 121629173423, 2355470737637, 48664218965021, 1067895971109199, 24795678053493443, 607144847919796830, 15630954703539323090, 421990078975569031642, 11918095123121138408128

Offset: 0

N. J. A. Sloane

J. de la Cal, J. Carcamo, Set partitions and moments of random variables, J. Math. Anal. Applic. 378 (2011) 16 doi:10.1016/j.jmaa.2011.01.002 Remark 5
J. Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353.
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).
R. P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Example 5.2.4.

Alois P. Heinz, Table of n, a(n) for n = 0..440
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Jekuthiel Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353. [Annotated scanned copy]
Gottfried Helms, Bell Numbers, 2008.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346. (Annotated scanned copy)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 293
K. A. Penson, P. Blasiak, G. Duchamp, A. Horzela and A. I. Solomon, Hierarchical Dobinski-type relations via substitution and the moment problem, arXiv:quant-ph/0312202, 2003.
K. A. Penson, P. Blasiak, G. Duchamp, A. Horzela and A. I. Solomon, Hierarchical Dobinski-type relations via substitution and the moment problem, J. Phys. A: Math.Gen 37 (2004) 3475-3487.
John Riordan, Letter, Apr 28 1976.
Kruchinin Vladimir Victorovich, Composition of ordinary generating functions, arXiv:1009.2565 [math.CO], 2010.
Index entries for sequences related to rooted trees

a(n)=|A039812(n,1)| (first column of triangle).
Cf. A000110, A000258, A000357, A000405, A001669.
Column k=3 of A144150.

g:= proc(p) local b; b:= proc(n) option remember; `if`(n=0, 1, (n-1)! *add(p(k)*b(n-k)/ (k-1)!/ (n-k)!, k=1..n)) end end: a:= g(g(g(1))): seq(a(n), n=0..30);  # Alois P. Heinz, Sep 11 2008

nn = 18; a = Exp[Exp[x] - 1]; b = Exp[a - 1];
Range[0, nn]! CoefficientList[Series[Exp[b - 1], {x, 0, nn}], x]  (*Geoffrey Critzer, Dec 28 2011*)

E.g.f.: exp(exp(exp(exp(x)-1)-1)-1).

a(n) = sum(sum(sum(stirling2(n,k) *stirling2(k,m) *stirling2(m,r), k=m..n), m=r..n), r=1..n), n>0. - Vladimir Kruchinin, Sep 08 2010

Extended with new definition by Christian G. Bower, Aug 15 1998

A000357 Number of 5-level labeled rooted trees with n leaves.

Original entry on oeis.org

1, 1, 5, 35, 315, 3455, 44590, 660665, 11035095, 204904830, 4183174520, 93055783320, 2238954627848, 57903797748386, 1601122732128779, 47120734323344439, 1470076408565099152, 48449426629560437576, 1681560512531504058350, 61293054886119796799892

Offset: 0

N. J. A. Sloane

J. de la Cal, J. Carcamo, Set partitions and moments of random variables, J. Math. Anal. Applic. 378 (2011) 16 doi:10.1016/j.jmaa.2011.01.002 Remark 5
J. Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346.
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).

Alois P. Heinz, Table of n, a(n) for n = 0..400
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Jekuthiel Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353. [Annotated scanned copy]
Gottfried Helms, Bell Numbers, 2008.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346. (Annotated scanned copy)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 294
Index entries for sequences related to rooted trees

a(n)=|A039813(n,1)| (first column of triangle).
Cf. A000110, A000258, A000307, A000405, A001669.
Column k=4 of A144150.

g:= proc(p) local b; b:=proc(n) option remember; if n=0 then 1 else (n-1)! *add(p(k)*b(n-k)/ (k-1)!/ (n-k)!, k=1..n) fi end end: a:= g(g(g(g(1)))): seq(a(n), n=0..30);  # Alois P. Heinz, Sep 11 2008

max = 17; Join[{1}, MatrixPower[Array[StirlingS2, {max, max}], 5][[All, 1]]] (* Jean-François Alcover, Mar 03 2014 *)

E.g.f.: exp(exp(exp(exp(exp(x)-1)-1)-1)-1).

Extended with new description by Christian G. Bower, Aug 15 1998

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

Christian G. Bower, Feb 15 1999

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

			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)

Tilman Piesk, First 100 rows, flattened
A. Aboud, J.-P. Bultel, A. Chouria, J.-G. Luque, and O. Mallet, Bell polynomials in combinatorial Hopf algebras, arXiv preprint arXiv:1402.2960 [math.CO], 2014-2015.
T. Mansour, A. Munagi, and M. Shattuck, Recurrence Relations and Two-Dimensional Set Partitions , J. Int. Seq. 14 (2011) # 11.4.1.

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.

# 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

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 *)

T(n, k) = sum(j=0, n, stirling(n, j, 2)*stirling(j, k, 2)); \\ Seiichi Manyama, Feb 13 2022

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)

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

A000405 Number of 6-level labeled rooted trees with n leaves.

Original entry on oeis.org

1, 1, 6, 51, 561, 7556, 120196, 2201856, 45592666, 1051951026, 26740775306, 742069051906, 22310563733864, 722108667742546, 25024187820786357, 924161461265888370, 36223781285638309482, 1501552062016443881514

Offset: 0

N. J. A. Sloane

J. de la Cal, J. Carcamo, Set partitions and moments of random variables, J. Math. Anal. Applic. 378 (2011) 16 doi:10.1016/j.jmaa.2011.01.002 Remark 5
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).

Alois P. Heinz, Table of n, a(n) for n = 0..150
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Jekuthiel Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353. [Annotated scanned copy]
Gottfried Helms, Bell Numbers, 2008.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346. (Annotated scanned copy)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 295
Index entries for sequences related to rooted trees

Cf. A000110, A000258, A000307, A000357, A001669.

Column k=5 of A144150.

g:= proc(p) local b; b:=proc(n) option remember; if n=0 then 1 else (n-1)! *add(p(k)*b(n-k)/ (k-1)!/ (n-k)!, k=1..n) fi end end: a:= g(g(g(g(g(1))))): seq(a(n), n=0..30); # Alois P. Heinz, Sep 11 2008

g[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, (n-1)!*Sum[p[k]*b[n-k]/(k-1)!/(n-k)!, {k, 1, n}]]; b]; a = Nest[g, 1&, 5]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 12 2014, after Alois P. Heinz *)

E.g.f.: exp(exp(exp(exp(exp(exp(x)-1)-1)-1)-1)-1).

Extended with new definition by Christian G. Bower, Aug 15 1998

A001669 Number of 7-level labeled rooted trees with n leaves.

Original entry on oeis.org

1, 1, 7, 70, 910, 14532, 274778, 5995892, 148154860, 4085619622, 124304629050, 4133867297490, 149114120602860, 5796433459664946, 241482353893283349, 10730629952953517859, 506500241174366575122, 25302666611855946733140

Offset: 0

N. J. A. Sloane

J. Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353.
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346.
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).

Alois P. Heinz, Table of n, a(n) for n = 0..150
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
Jekuthiel Ginsburg, Iterated exponentials, Scripta Math., 11 (1945), 340-353. [Annotated scanned copy]
T. Hogg and B. A. Huberman, Attractors on finite sets: the dissipative dynamics of computing structures, Phys. Review A 32 (1985), 2338-2346. (Annotated scanned copy)
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 296
Index entries for sequences related to rooted trees

Cf. A000110, A000258, A000307, A000357, A000405.

Column k=6 of A144150.

g:= proc(p) local b; b:= proc(n) option remember; if n=0 then 1 else (n-1)! *add(p(k)*b(n-k)/ (k-1)!/ (n-k)!, k=1..n) fi end end: a:= g(g(g(g(g(g(1)))))): seq(a(n), n=0..30); # Alois P. Heinz, Sep 11 2008

g[p_] := Module[{b}, b[n_] := b[n] = If[n == 0, 1, (n-1)!*Sum[p[k]*b[n-k]/(k-1)!/(n-k)!, {k, 1, n}]]; b]; a = Nest[g, 1&, 6]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
With[{nn=20},Join[{1},Rest[CoefficientList[Series[Nest[Exp[#]-1&,Exp[x]-1,6],{x,0,nn}],x] Range[0,nn]!]]] (* Harvey P. Dale, Mar 02 2015 *)

E.g.f.: exp(exp(exp(exp(exp(exp(exp(x)-1)-1)-1)-1)-1)-1).

Extended with new definition by Christian G. Bower, Aug 15 1998

A050345 Number of ways to factor n into distinct factors with one level of parentheses.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 3, 1, 3, 1, 6, 1, 3, 3, 4, 1, 6, 1, 6, 3, 3, 1, 13, 1, 3, 3, 6, 1, 12, 1, 7, 3, 3, 3, 15, 1, 3, 3, 13, 1, 12, 1, 6, 6, 3, 1, 25, 1, 6, 3, 6, 1, 13, 3, 13, 3, 3, 1, 31, 1, 3, 6, 12, 3, 12, 1, 6, 3, 12, 1, 37, 1, 3, 6, 6, 3, 12, 1, 25, 4, 3, 1, 31, 3, 3, 3, 13, 1, 31, 3, 6, 3, 3

Offset: 1

Christian G. Bower, Oct 15 1999

nonn

First differs from A296120 at a(36) = 15, A296120(36) = 14. - Gus Wiseman, Apr 27 2025
Each "part" in parentheses is distinct from all others at the same level. Thus (3*2)*(2) is allowed but (3)*(2*2) and (3*2*2) are not.
a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).

			12 = (12) = (6*2) = (6)*(2) = (4*3) = (4)*(3) = (3*2)*(2).
From _Gus Wiseman_, Apr 26 2025: (Start)
This is the number of ways to partition a factorization of n (counted by A001055) into a set of sets. For example, the a(12) = 6 choices are:
  {{2},{2,3}}
  {{2},{6}}
  {{3},{4}}
  {{2,6}}
  {{3,4}}
  {{12}}
(End)

R. J. Mathar, Table of n, a(n) for n = 1..2519

Cf. A050346-A050350. a(A002110)=A000258.
For multisets of multisets we have A050336.
For integer partitions we have a(p^k) = A050342(k), see A001970, A089259, A261049.
For normal multiset partitions see A116539, A292432, A292444, A381996, A382214, A382216.
The case of a unique choice (positions of 1) is A166684.
Twice-partitions of this type are counted by A358914, see A270995, A281113, A294788.
For sets of multisets we have A383310 (distinct products A296118).
For multisets of sets we have we have A383311, see A296119.
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers, distinct A050326.
A302494 gives MM-numbers of sets of sets.
A382077 counts partitions that can be partitioned into a sets of sets, ranks A382200.
A382078 counts partitions that cannot be partitioned into a sets of sets, ranks A293243.
Cf. A000009, A002865, A005117, A045782, A279785, A293511, A296120, A316439, A381441, A382201.

facs[n_]:=If[n<=1,{{}}, Join@@Table[Map[Prepend[#,d]&, Select[facs[n/d],Min@@#>=d&]],{d, Rest[Divisors[n]]}]];
sps[{}]:={{}};sps[set:{i_,_}] := Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]] /@ Cases[Subsets[set],{i,_}];
mps[set_]:=Union[Sort[Sort /@ (#/.x_Integer:>set[[x]])]& /@ sps[Range[Length[set]]]];
Table[Sum[Length[Select[mps[y], UnsameQ@@#&&And@@UnsameQ@@@#&]], {y,facs[n]}],{n,30}] (* Gus Wiseman, Apr 26 2025 *)

Dirichlet g.f.: Product_{n>=2}(1+1/n^s)^A045778(n).
a(n) = A050346(A025487^(-1)(A046523(n))), where A025487^(-1) is the inverse with A025487^(-1)(A025487(n))=n. - R. J. Mathar, May 25 2017
a(n) = A050346(A101296(n)). - Antti Karttunen, May 25 2017

cp's OEIS Frontend

A318762 Number of permutations of a multiset whose multiplicities are the prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A144150 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where the e.g.f. of column k is 1+g^(k+1)(x) with g = x-> exp(x)-1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Python

Formula

A122848 Exponential Riordan array (1, x(1+x/2)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Sage

Formula

A318393 Regular tetrangle where T(n,k,i) is the number of pairs of set partitions of {1,...,n} with meet of length k and join of length i.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

A000307 Number of 4-level labeled rooted trees with n leaves.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A000357 Number of 5-level labeled rooted trees with n leaves.

Views

Author

Keywords

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links