A005651 -id:A005651 - 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).

A226873 Number A(n,k) of n-length words w over a k-ary alphabet {a1,a2,...,ak} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 0, where #(w,x) counts the letters x in word w; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 3, 1, 0, 1, 1, 3, 4, 1, 0, 1, 1, 3, 10, 11, 1, 0, 1, 1, 3, 10, 23, 16, 1, 0, 1, 1, 3, 10, 47, 66, 42, 1, 0, 1, 1, 3, 10, 47, 126, 222, 64, 1, 0, 1, 1, 3, 10, 47, 246, 522, 561, 163, 1, 0, 1, 1, 3, 10, 47, 246, 882, 1821, 1647, 256, 1, 0

Offset: 0

Alois P. Heinz, Jun 21 2013

			A(4,3) = 23: aaaa, aaab, aaba, aabb, aabc, aacb, abaa, abab, abac, abba, abca, acab, acba, baaa, baab, baac, baba, baca, bbaa, bcaa, caab, caba, cbaa.
Square array A(n,k) begins:
  1, 1,  1,   1,    1,    1,    1,     1, ...
  0, 1,  1,   1,    1,    1,    1,     1, ...
  0, 1,  3,   3,    3,    3,    3,     3, ...
  0, 1,  4,  10,   10,   10,   10,    10, ...
  0, 1, 11,  23,   47,   47,   47,    47, ...
  0, 1, 16,  66,  126,  246,  246,   246, ...
  0, 1, 42, 222,  522,  882, 1602,  1602, ...
  0, 1, 64, 561, 1821, 3921, 6441, 11481, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0-10 give: A000007, A000012, A027306, A092255, A092429, A226875, A226876, A226877, A226878, A226879, A226880.
Main diagonal gives: A005651.
Cf. A131632, A182172.

b:= proc(n, i, t) option remember;
      `if`(t=1, 1/n!, add(b(n-j, j, t-1)/j!, j=i..n/t))
    end:
A:= (n, k)-> `if`(k=0, `if`(n=0, 1, 0), n!*b(n, 0, k)):
seq(seq(A(n, d-n), n=0..d), d=0..14);

b[n_, i_, t_] := b[n, i, t] = If[t == 1, 1/n!, Sum[b[n-j, j, t-1]/j!, {j, i, n/t}]]; a[n_, k_] := If[k == 0, If[n == 0, 1, 0], n!*b[n, 0, k]]; Table[Table[a[n, d-n], {n, 0, d}], {d, 0, 14}] // Flatten (* Jean-François Alcover, Dec 13 2013, translated from Maple *)

A(n,k) = Sum_{i=0..min(n,k)} A226874(n,i).

A321912 Tetrangle where T(n,H(u),H(v)) is the coefficient of m(v) in e(u), where u and v are integer partitions of n, H is Heinz number, m is monomial symmetric functions, and e is elementary symmetric functions.

Original entry on oeis.org

1, 0, 1, 1, 2, 0, 0, 1, 0, 1, 3, 1, 3, 6, 0, 0, 0, 0, 1, 0, 1, 0, 2, 6, 0, 0, 0, 1, 4, 0, 2, 1, 5, 12, 1, 6, 4, 12, 24, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 5, 0, 0, 0, 1, 0, 3, 10, 0, 0, 1, 5, 2, 12, 30, 0, 0, 0, 2, 1, 7, 20, 0, 1, 3, 12, 7, 27, 60, 1, 5

Offset: 1

Gus Wiseman, Nov 22 2018

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Also the coefficient of f(v) in h(u), where f is forgotten symmetric functions and h is homogeneous symmetric functions.

			Tetrangle begins (zeroes not shown):
  (1):  1
.
  (2):      1
  (11):  1  2
.
  (3):          1
  (21):      1  3
  (111):  1  3  6
.
  (4):                 1
  (22):       1     2  6
  (31):             1  4
  (211):      2  1  5 12
  (1111):  1  6  4 12 24
.
  (5):                        1
  (41):                    1  5
  (32):              1     3 10
  (221):          1  5  2 12 30
  (311):             2  1  7 20
  (2111):      1  3 12  7 27 60
  (11111):  1  5 10 30 20 60 20
For example, row 14 gives: e(32) = m(221) + 3m(2111) + 10m(11111).

Wikipedia, Symmetric polynomial

This is a regrouping of the triangle A321742.

Cf. A005651, A008480, A056239, A124794, A124795, A215366, A318284, A318360, A319191, A319193, A321854, A321738, A321913-A321935.

A247551 Decimal expansion of Product_{k>=2} 1/(1-1/k!).

Original entry on oeis.org

2, 5, 2, 9, 4, 7, 7, 4, 7, 2, 0, 7, 9, 1, 5, 2, 6, 4, 8, 1, 8, 0, 1, 1, 6, 1, 5, 4, 2, 5, 3, 9, 5, 4, 2, 4, 1, 1, 7, 8, 7, 0, 2, 3, 9, 4, 8, 4, 5, 9, 9, 7, 3, 3, 7, 5, 8, 4, 9, 3, 4, 9, 8, 2, 5, 0, 0, 2, 1, 1, 8, 7, 8, 0, 0, 8, 6, 6, 9, 9, 1, 2, 1, 9, 9, 8, 8, 3, 6, 4, 6, 2, 5, 2, 6, 1, 4, 9, 5, 5, 1, 6, 4, 3, 2

Offset: 1

Vaclav Kotesovec, Sep 19 2014

			2.5294774720791526481801161542539542411787023948459973375849349825...

Steven R. Finch, Mathematical Constants II, Encyclopedia of Mathematics and Its Applications, Cambridge University Press, Cambridge, 2018, p. 269.
A. Knopfmacher, A. M. Odlyzko, B. Pittel, L. B. Richmond, D. Stark, G. Szekeres, and N. C. Wormald, The Asymptotic Number of Set Partitions with Unequal Block Sizes, The Electronic Journal of Combinatorics, Volume 6 (1999), Research Paper #R2.

Cf. A005651, A035341, A072605, A140585, A183239, A209668, A261600, A282529, A292713, A320566, A327827, A335643.

evalf(1/product(1-1/k!,k=2..infinity),100); # Vaclav Kotesovec, Sep 19 2014

digits = 105;
RealDigits[NProduct[1/(1-1/k!), {k, 2, Infinity}, WorkingPrecision -> digits+10, NProductFactors -> digits], 10, digits][[1]] (* Jean-François Alcover, Nov 02 2020 *)

default(realprecision,150); 1/prodinf(k=2,1 - 1/k!) \\ Vaclav Kotesovec, Sep 21 2014

Product_{k>=2} 1/(1-1/k!).
Equals lim_{n -> oo} A005651(n) / n!.
Equals 1/A282529. - Amiram Eldar, Sep 15 2023

A226874 Number T(n,k) of n-length words w over a k-ary alphabet {a1, a2, ..., ak} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 1, where #(w,x) counts the letters x in word w; triangle T(n,k), n >= 0, 0 <= k <= n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 1, 3, 6, 0, 1, 10, 12, 24, 0, 1, 15, 50, 60, 120, 0, 1, 41, 180, 300, 360, 720, 0, 1, 63, 497, 1260, 2100, 2520, 5040, 0, 1, 162, 1484, 6496, 10080, 16800, 20160, 40320, 0, 1, 255, 5154, 20916, 58464, 90720, 151200, 181440, 362880

Offset: 0

Alois P. Heinz, Jun 21 2013

T(n,k) is the sum of multinomials M(n; lambda), where lambda ranges over all partitions of n into parts that form a multiset of size k.

			T(4,2) = 10: aaab, aaba, aabb, abaa, abab, abba, baaa, baab, baba, bbaa.
T(4,3) = 12: aabc, aacb, abac, abca, acab, acba, baac, baca, bcaa, caab, caba, cbaa.
T(5,2) = 15: aaaab, aaaba, aaabb, aabaa, aabab, aabba, abaaa, abaab, ababa, abbaa, baaaa, baaab, baaba, babaa, bbaaa.
Triangle T(n,k) begins:
  1;
  0,  1;
  0,  1,   2;
  0,  1,   3,    6;
  0,  1,  10,   12,   24;
  0,  1,  15,   50,   60,   120;
  0,  1,  41,  180,  300,   360,   720;
  0,  1,  63,  497, 1260,  2100,  2520,  5040;
  0,  1, 162, 1484, 6496, 10080, 16800, 20160, 40320;
  ...

Alois P. Heinz, Rows n = 0..140, flattened
Wikipedia, Iverson bracket
Wikipedia, Multinomial coefficients
Wikipedia, Partition (number theory)

Columns k=0-10 give: A000007, A057427, A226881, A226882, A226883, A226884, A226885, A226886, A226887, A226888, A226889.
Main diagonal gives: A000142.
Row sums give: A005651.
T(2n,n) gives A318796.
Cf. A131632, A285824, A292222, A327803.

b:= proc(n, i, t) option remember;
      `if`(t=1, 1/n!, add(b(n-j, j, t-1)/j!, j=i..n/t))
    end:
T:= (n, k)-> `if`(n*k=0, `if`(n=k, 1, 0), n!*b(n, 1, k)):
seq(seq(T(n, k), k=0..n), n=0..12);
# second Maple program:
b:= proc(n, i) option remember; expand(
      `if`(n=0, 1, `if`(i<1, 0, add(x^j*b(n-i*j, i-1)*
      combinat[multinomial](n, n-i*j, i$j), j=0..n/i))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n$2)):
seq(T(n), n=0..12);

b[n_, i_, t_] := b[n, i, t] = If[t == 1, 1/n!, Sum[b[n - j, j, t - 1]/j!, {j, i, n/t}]]; t[n_, k_] := If[n*k == 0, If[n == k, 1, 0], n!*b[n, 1, k]]; Table[Table[t[n, k], {k, 0, n}], {n, 0, 12}] // Flatten (* Jean-François Alcover, Dec 13 2013, translated from first Maple *)

T(n)={Vec(serlaplace(prod(k=1, n, 1/(1-y*x^k/k!) + O(x*x^n))))}
{my(t=T(10)); for(n=1, #t, for(k=0, n-1, print1(polcoeff(t[n], k), ", ")); print)} \\ Andrew Howroyd, Dec 20 2017

T(n,k) = A226873(n,k) - [k>0] * A226873(n,k-1).

A261719 Number T(n,k) of partitions of n where each part i is marked with a word of length i over a k-ary alphabet whose letters appear in alphabetical order and all k letters occur at least once in the partition; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 2, 3, 0, 3, 12, 10, 0, 5, 40, 81, 47, 0, 7, 104, 396, 544, 246, 0, 11, 279, 1751, 4232, 4350, 1602, 0, 15, 654, 6528, 25100, 44475, 36744, 11481, 0, 22, 1577, 23892, 136516, 369675, 512787, 352793, 95503, 0, 30, 3560, 80979, 666800, 2603670, 5413842, 6170486, 3641992, 871030

Offset: 0

Alois P. Heinz, Aug 29 2015

T(n,k) is defined for n,k >= 0. The triangle contains only the terms with k<=n. T(n,k) = 0 for k>n.

			A(3,2) = 12: 3aab, 3abb, 2aa1b, 2ab1a, 2ab1b, 2bb1a, 1a1a1b, 1a1b1a, 1a1b1b, 1b1a1a, 1b1a1b, 1b1b1a.
Triangle T(n,k) begins:
  1
  0,  1;
  0,  2,    3;
  0,  3,   12,    10;
  0,  5,   40,    81,     47;
  0,  7,  104,   396,    544,    246;
  0, 11,  279,  1751,   4232,   4350,   1602;
  0, 15,  654,  6528,  25100,  44475,  36744,  11481;
  0, 22, 1577, 23892, 136516, 369675, 512787, 352793, 95503;
  ...

Alois P. Heinz, Rows n = 0..140, flattened

Columns k=0-10 give: A000007, A000041 (for n>0), A293366, A293367, A293368, A293369, A293370, A293371, A293372, A293373, A293374.
Row sums give A035341.
Main diagonal gives A005651.
T(2n,n) gives A261732.
Cf. A060642, A261718, A261781 (same for compositions).

b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<1, 0,
      b(n, i-1, k)+`if`(i>n, 0, b(n-i, i, k)*binomial(i+k-1, k-1))))
    end:
T:= (n, k)-> add(b(n$2, k-i)*(-1)^i*binomial(k, i), i=0..k):
seq(seq(T(n, k), k=0..n), n=0..10);

b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0, b[n, i - 1, k] + If[i > n, 0, b[n - i, i, k]*Binomial[i + k - 1, k - 1]]]]; T[n_, k_] := Sum[b[n, n, k - i]*(-1)^i*Binomial[k, i], {i, 0, k}]; Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jan 21 2017, translated from Maple *)

T(n,k) = Sum_{i=0..k} (-1)^i * C(k,i) * A261718(n,k-i).

A319225 Number of acyclic spanning subgraphs of a cycle graph, where the sizes of the connected components are given by the prime indices of n.

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 4, 1, 2, 4, 5, 4, 6, 5, 5, 1, 7, 5, 8, 5, 6, 6, 9, 5, 3, 7, 2, 6, 10, 12, 11, 1, 7, 8, 7, 9, 12, 9, 8, 6, 13, 14, 14, 7, 7, 10, 15, 6, 4, 7, 9, 8, 16, 7, 8, 7, 10, 11, 17, 21, 18, 12, 8, 1, 9, 16, 19, 9, 11, 16, 20, 14, 21, 13, 8, 10, 9, 18

Offset: 1

Gus Wiseman, Sep 13 2018

nonn

a(1) = 1 by convention.

A prime index of n is a number m such that prime(m) divides n.

			Of the cycle ({1,2,3}, {(1,2),(2,3),(3,1)}) the spanning subgraphs where the sizes of connected components are (2,1) are: ({1,2,3}, {(1,2)}), ({1,2,3}, {(2,3)}), ({1,2,3}, {(3,1)}). Since the prime indices of 6 are (2,1), we conclude a(6) = 3.

Wikipedia, Expressing power sums in terms of elementary symmetric polynomials
Gus Wiseman, Enumeration of paths and cycles and e-coefficients of incomparability graphs, arXiv:0709.0430 [math.CO], 2007.

Different orderings with signs are A115131, A210258, A263916.

Cf. A005651, A008480, A048994, A056239, A124794, A124795, A135278, A215366, A318762, A319191, A319193, A319226.

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]]]]]]]]];
Table[Length[With[{m=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]},Select[Subsets[Partition[Range[Total[m]],2,1,1],{Total[m]-PrimeOmega[n]}],Sort[Length/@csm[Union[#,List/@Range[Total[m]]]]]==m&]]],{n,30}]

a(n) = A056239(n) * (Omega(n) - 1)! / Product c_i! where c_i is the multiplicity of prime(i) in the prime factorization of n.

A386585 Triangle read by rows where T(n,k) is the number of integer partitions y of n into k = 0..n parts such that any multiset whose multiplicities are the parts of y is separable.

Original entry on oeis.org

1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 1, 1, 0, 0, 1, 2, 1, 1, 0, 0, 1, 2, 2, 1, 1, 0, 0, 1, 3, 3, 2, 1, 1, 0, 0, 1, 3, 4, 3, 2, 1, 1, 0, 0, 1, 5, 5, 5, 3, 2, 1, 1, 0, 0, 1, 4, 7, 6, 5, 3, 2, 1, 1

Offset: 0

Gus Wiseman, Aug 02 2025

We say that such partitions are of separable type.

A multiset is separable iff it has a permutation without any adjacent equal parts.

			Row n = 8 counts the following partitions:
  .  .  44  431  4211  41111  311111  2111111  11111111
            422  3311  32111  221111
            332  3221  22211
                 2222
with the following separable multisets:
  . . 11112222 11112223 11112234 11112345 11123456 11234567 12345678
               11112233 11122234 11122345 11223456
               11122233 11122334 11223345
                        11223344
Triangle begins:
  1
  0  1
  0  0  1
  0  0  1  1
  0  0  1  1  1
  0  0  1  2  1  1
  0  0  1  2  2  1  1
  0  0  1  3  3  2  1  1
  0  0  1  3  4  3  2  1  1
  0  0  1  5  5  5  3  2  1  1
  0  0  1  4  7  6  5  3  2  1  1

This is the separable type case of A072233 or A008284.
Row sums are A336106, ranks A335127.
For separable instead of separable type we have A386583, inseparable A386584.
For inseparable instead of separable we have A386586, sums A025065, ranks A335126.
A003242 and A335452 count anti-runs, ranks A333489, patterns A005649.
A239455 counts Look-and-Say partitions, ranks A351294.
A279790 counts disjoint families on strongly normal multisets.
A325534 counts separable multisets, ranks A335433.
A325535 counts inseparable multisets, ranks A335448.
A336103 counts normal separable multisets, inseparable A336102.
A351293 counts non-Look-and-Say partitions, ranks A351295.
A386633 counts separable set partitions, row sums of A386635.
A386634 counts inseparable set partitions, row sums of A386636.
Cf. A005651, A106351, A111133, A238130, A335434, A386575, A386576, A386579, A386580, A386581, A386582.

sepQ[y_]:=Select[Permutations[y],Length[Split[#]]==Length[y]&]!={};
mst[y_]:=Join@@Table[ConstantArray[k,y[[k]]],{k,Length[y]}];
Table[Length[Select[IntegerPartitions[n,{k}],sepQ[mst[#]]&]],{n,0,5},{k,0,n}]

a(n) = A072233(n) - A386586(n).

A386586 Triangle read by rows where T(n,k) is the number of integer partitions y of n into k parts such that any multiset whose multiplicities are the parts of y is inseparable.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 1, 2, 1, 0, 0, 0, 0, 0, 1, 3, 2, 1, 0, 0, 0, 0, 0, 1, 3, 2, 1, 0, 0, 0, 0, 0, 0, 1, 4, 4, 2, 1, 0, 0, 0, 0, 0

Offset: 0

Gus Wiseman, Aug 05 2025

We say that such partitions are of inseparable type. This is different from inseparable partitions (see A386584). A multiset is separable iff it has a permutation without any adjacent equal parts.

			The partition y = (7,2,1) is the multiplicities of the multiset {1,1,1,1,1,1,1,2,2,3}, which is inseparable, so y is counted under T(10,3).
Row n = 10 counts the following partitions (A = 10):
  .  A  91  811  7111  61111  .  .  .  .  .
        82  721  6211
        73  631
        64  622
Triangle begins:
  0
  0 0
  0 1 0
  0 1 0 0
  0 1 1 0 0
  0 1 1 0 0 0
  0 1 2 1 0 0 0
  0 1 2 1 0 0 0 0
  0 1 3 2 1 0 0 0 0
  0 1 3 2 1 0 0 0 0 0
  0 1 4 4 2 1 0 0 0 0 0

This is the inseparable type case of A008284 or A072233.
Row sums shifted left once are A025065 (ranks A335126), separable version A336106 (ranks A335127).
For separable instead of inseparable type we have A386583.
For integer partitions instead of normal multisets we have A386584.
For separable type instead of inseparable type we have A386585.
A003242 and A335452 count anti-runs, ranks A333489, patterns A005649.
A239455 counts Look-and-Say partitions, ranks A351294.
A325534 counts separable multisets, ranks A335433.
A325535 counts inseparable multisets, ranks A335448.
A336103 counts normal separable multisets, inseparable A336102.
A351293 counts non-Look-and-Say partitions, ranks A351295.
Cf. A000670, A005651, A106351, A238130, A279790, A386575, A386576, A386579, A386582, A386634.

insepQ[y_]:=Select[Permutations[y],Length[Split[#]]==Length[y]&]=={};
ptm[y_]:=Join@@Table[ConstantArray[k,y[[k]]],{k,Length[y]}];
Table[Length[Select[IntegerPartitions[n,{k}],insepQ[ptm[#]]&]],{n,0,5},{k,0,n}]

a(n) = A072233(n) - A386585(n).

A178682 The number of functions f:{1,2,...,n}->{1,2,...,n} such that the number of elements that are mapped to m is divisible by m.

Original entry on oeis.org

1, 1, 2, 5, 13, 42, 150, 576, 2266, 9966, 47466, 237019, 1224703, 6429152, 35842344, 212946552, 1325810173, 8488092454, 55276544436, 362961569008, 2465240278980, 17538501945077, 130454679958312, 1002493810175093, 7838007702606372, 61789072382062638

Offset: 0

Geoffrey Critzer, Dec 25 2010

a(n) is also the number of partitions of n where each block of part i with multiplicity j is marked with a word of length i*j over an n-ary alphabet whose letters appear in alphabetical order and all n letters occur exactly once in the partition. a(3) = 5: 3abc, 2ab1c, 2ac1b, 2bc1a, 111abc. There is a simple bijection between the marked partitions and the functions f. - Alois P. Heinz, Aug 30 2015

			a(3) = 5 because there are 5 such functions: (1,1,1), (1,2,2), (2,1,2), (2,2,1), (3,3,3).
G.f. = 1 + x + 2*x^2 + 5*x^3 + 13*x^4 + 42*x^5 + 150*x^6 + 576*x^7 + ...

Alois P. Heinz, Table of n, a(n) for n = 0..680

Cf. A000110, A005651, A007837, A261774.
Main diagonal of A326500, A326616, A326617.
Row sums of A364285, A364310.

m:=25; R:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( (&*[(&+[x^(k*j)/Factorial(k*j): k in [0..m]]): j in [1..m]]) )); [Factorial(n-1)*b[n]: n in [1..m]]; // G. C. Greubel, Jan 26 2019

b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<1, 0,
      add(b(n-i*j, i-1)*binomial(n, i*j), j=0..n/i)))
    end:
a:= n-> b(n$2):
seq(a(n), n=0..30);  # Alois P. Heinz, Aug 30 2015

Range[0,20]! CoefficientList[Series[Product[Sum[x^(j i)/(j i)!,{i,0,20}],{j,1,20}],{x,0,20}],x]

m=30; my(x='x+O('x^m)); Vec(serlaplace(prod(j=1, m, sum(k=0,m, x^(k*j)/(k*j)!)))) \\ G. C. Greubel, Jan 26 2019

m = 30; T = taylor(product(sum(x^(k*j)/factorial(k*j) for k in (0..m)) for j in (1..m)), x, 0, m); [factorial(n)*T.coefficient(x, n) for n in (0..m)] # G. C. Greubel, Jan 26 2019

E.g.f.: Product_{j>=1} Sum_{i>=0} x^(j*i)/(j*i)!.

a(21)-a(25) from Alois P. Heinz, Aug 30 2015

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

A226873 Number A(n,k) of n-length words w over a k-ary alphabet {a1,a2,...,ak} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 0, where #(w,x) counts the letters x in word w; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A321912 Tetrangle where T(n,H(u),H(v)) is the coefficient of m(v) in e(u), where u and v are integer partitions of n, H is Heinz number, m is monomial symmetric functions, and e is elementary symmetric functions.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A247551 Decimal expansion of Product_{k>=2} 1/(1-1/k!).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A226874 Number T(n,k) of n-length words w over a k-ary alphabet {a1, a2, ..., ak} such that #(w,a1) >= #(w,a2) >= ... >= #(w,ak) >= 1, where #(w,x) counts the letters x in word w; triangle T(n,k), n >= 0, 0 <= k <= n, read by rows.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A261719 Number T(n,k) of partitions of n where each part i is marked with a word of length i over a k-ary alphabet whose letters appear in alphabetical order and all k letters occur at least once in the partition; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A319225 Number of acyclic spanning subgraphs of a cycle graph, where the sizes of the connected components are given by the prime indices of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica