A317776 -id:A317776 - cp's OEIS frontend

A318286 Number of strict multiset partitions of a multiset whose multiplicities are the prime indices of n.

1, 1, 1, 2, 2, 3, 2, 5, 5, 5, 3, 9, 4, 7, 9, 15, 5, 18, 6, 16, 14, 10, 8, 31, 17, 14, 40, 25, 10, 34, 12, 52, 21, 19, 27, 70, 15, 25, 31, 59, 18, 57, 22, 38, 80, 33, 27, 120, 46, 67, 44, 56, 32, 172, 42, 100, 61, 43, 38, 141, 46, 55, 143, 203, 64, 91, 54, 80

Offset: 1

Gus Wiseman, Aug 23 2018

nonn

Andrew Howroyd, Table of n, a(n) for n = 1..500

Cf. A001055, A007716, A045778, A181821, A305936, A316980, A317776.

Cf. A318283, A318284, A318285, A318287, A318360, A318361.

nrmptn[n_]:=Join@@MapIndexed[Table[#2[[1]],{#1}]&,If[n==1,{},Flatten[Cases[FactorInteger[n]//Reverse,{p_,k_}:>Table[PrimePi[p],{k}]]]]];
strfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strfacs[n/d],Min@@#>d&]],{d,Rest[Divisors[n]]}]];
Table[Length[strfacs[Times@@Prime/@nrmptn[n]]],{n,60}]

permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
sig(n)={my(f=factor(n)); concat(vector(#f~, i, vector(f[i, 2], j, primepi(f[i, 1]))))}
count(sig)={my(r=0, A=O(x*x^vecmax(sig))); for(n=1, vecsum(sig)+1, my(s=0); forpart(p=n, my(q=1/prod(i=1, #p, 1 - x^p[i] + A)); s+=prod(i=1, #sig, polcoef(q, sig[i]))*(-1)^#p*permcount(p)); r+=(-1)^n*s/n!); r/2}
a(n)={if(n==1, 1, count(sig(n)))} \\ Andrew Howroyd, Dec 18 2018

a(n) = A045778(A181821(n)).

a(prime(n)^k) = A219585(n, k). - Andrew Howroyd, Dec 17 2018

A317752 Number of multiset partitions of normal multisets of size n such that the blocks have empty intersection.

Original entry on oeis.org

0, 1, 8, 49, 305, 1984, 13686, 100124, 776885, 6386677, 55532358, 509549386, 4921352952, 49899820572, 529807799836, 5876162077537, 67928460444139, 816764249684450, 10195486840926032, 131896905499007474, 1765587483656124106, 24419774819813602870

Offset: 1

Gus Wiseman, Aug 06 2018

nonn

A multiset is normal if it spans an initial interval of positive integers.

			The a(3) = 8 multiset partitions with empty intersection:
  {{2},{1,1}}
  {{1},{2,2}}
  {{1},{2,3}}
  {{2},{1,3}}
  {{3},{1,2}}
  {{1},{1},{2}}
  {{1},{2},{2}}
  {{1},{2},{3}}

Andrew Howroyd, Table of n, a(n) for n = 1..100

Cf. A007716, A255903, A255906, A317073, A281116, A317077, A317532, A317533.

Cf. A317748, A317751, A317755, A317757, A317776.

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]]]];
allnorm[n_]:=Function[s,Array[Count[s,y_/;y<=#]+1&,n]]/@Subsets[Range[n-1]+1];
Table[Length[Join@@Table[Select[mps[m],Intersection@@#=={}&],{m,allnorm[n]}]],{n,6}]

P(n,k)={1/prod(i=1, n, (1 - x^i*y + O(x*x^n))^binomial(k+i-1, k-1))}
R(n,k)={my(p=P(n,k), q=p/(1-y+O(y*y^n))); Vec(sum(i=2, n, polcoef(p,i,y) + polcoef(q,i,y)*sum(j=1, n\i, (-1)^j*binomial(k,j)*x^(i*j))), -n)}
seq(n)={sum(k=2, n, R(n, k)*sum(r=k, n, binomial(r, k)*(-1)^(r-k)) )} \\ Andrew Howroyd, Feb 05 2021

Terms a(9) and beyond from Andrew Howroyd, Feb 05 2021

A381452 Number of multisets that can be obtained by partitioning the prime indices of n into a set of multisets and taking their sums.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Mar 06 2025

nonn

First differs from A045778 at a(24) = 4, A045778(24) = 5.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a factorization of n into distinct factors > 1.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
A multiset partition can be regarded as an arrow in the poset of integer partitions. For example, we have {{1},{1,2},{1,3},{1,2,3}}: {1,1,1,1,2,2,3,3} -> {1,3,4,6}, or (33221111) -> (6431) (depending on notation).
Sets of multisets are generally not transitive. For example, we have arrows: {{1},{2},{1,2}}: {1,1,2,2} -> {1,2,3} and {{1,2},{3}}: {1,2,3} -> {3,3}, but there is no set of multisets {1,1,2,2} -> {3,3}.

			The prime indices of 24 are {1,1,1,2}, with 5 partitions into a set of multisets:
  {{1,1,1,2}}
  {{1},{1,1,2}}
  {{2},{1,1,1}}
  {{1,1},{1,2}}
  {{1},{2},{1,1}}
with block-sums: {5}, {1,4}, {2,3}, {2,3}, {1,2,2}, of which 4 are distinct, so a(24) = 4.

Before taking sums we had A045778.
If each block is a set we have A381441, before sums A050326.
For distinct block-sums instead of blocks we have A381637, before sums A321469.
Other multiset partitions of prime indices:
- For multisets of constant multisets (A000688) see A381455 (upper), A381453 (lower).
- For multiset partitions (A001055) see A317141 (upper), A300383 (lower).
- For set multipartitions (A050320) see A381078 (upper), A381454 (lower).
- For sets of constant multisets (A050361) see A381715.
- For set systems with distinct sums (A381633) see A381634, zeros A293243.
- For sets of constant multisets with distinct sums (A381635) see A381716, A381636.
More on sets of multisets: A261049, A317776, A317775, A296118, A318286.
A000041 counts integer partitions, strict A000009.
A000040 lists the primes.
A003963 gives product of prime indices.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
Cf. A000720, A001222, A001970, A002846, A066328, A213385, A213427, A299200, A299202, A300385, A317142.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mps[mset_]:=Union[Sort[Sort/@(#/.x_Integer:>mset[[x]])]&/@sps[Range[Length[mset]]]];
Table[Length[Union[Sort[Total/@#]&/@Select[mps[prix[n]],UnsameQ@@#&]]],{n,100}]

a(A002110(n)) = A066723(n).

A327116 Number T(n,k) of colored integer partitions of n using all colors of a k-set such that all parts have different color patterns and a pattern for part i has i colors in (weakly) increasing order; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 1, 2, 0, 2, 6, 5, 0, 2, 15, 27, 15, 0, 3, 32, 102, 124, 52, 0, 4, 65, 319, 656, 600, 203, 0, 5, 124, 897, 2780, 4210, 3084, 877, 0, 6, 230, 2346, 10305, 23040, 27567, 16849, 4140, 0, 8, 414, 5818, 34864, 108135, 188284, 186095, 97640, 21147

Offset: 0

Alois P. Heinz, Sep 13 2019

			T(3,2) = 6; 3aab, 3abb, 2aa1b, 2ab1a, 2ab1b, 2bb1a.
Triangle T(n,k) begins:
  1;
  0, 1;
  0, 1,   2;
  0, 2,   6,    5;
  0, 2,  15,   27,    15;
  0, 3,  32,  102,   124,     52;
  0, 4,  65,  319,   656,    600,    203;
  0, 5, 124,  897,  2780,   4210,   3084,    877;
  0, 6, 230, 2346, 10305,  23040,  27567,  16849,  4140;
  0, 8, 414, 5818, 34864, 108135, 188284, 186095, 97640, 21147;
  ...

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

Columns k=0-2 give: A000007, A000009 (for n>0), A327598.
Main diagonal gives A000110.
Row sums give A317776.
T(2n,n) gives A327556.
Cf. A255903, A326914, A326962, A327117, A327557, A309973.

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

c = Binomial;
b[n_, i_, k_] := b[n, i, k] = If[n == 0, 1, If[i < 1, 0, Sum[b[n - i j, Min[n - i j, i - 1], k] c[c[k + i - 1, i], j], {j, 0, n/i}]]];
T[n_, k_] := Sum[b[n, n, i] (-1)^(k - i) c[k, i], {i, 0, k}];
Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, Apr 27 2020, after Alois P. Heinz *)

Sum_{k=1..n} k * T(n,k) = A327557(n).

A317775 Number of strict multiset partitions of strongly normal multisets of size n, where a multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

Original entry on oeis.org

1, 3, 10, 36, 136, 596, 2656, 13187, 68226, 381572, 2233091, 13940407, 90981030, 626911429, 4509031955, 33987610040, 266668955183, 2180991690286, 18512572760155, 163103174973092, 1487228204311039, 14027782824491946, 136585814043190619, 1371822048393658001, 14190528438090988629

Offset: 1

Gus Wiseman, Aug 06 2018

nonn

			The a(3) = 10 strict multiset partitions:
  {{1,1,1}}, {{1},{1,1}},
  {{1,1,2}}, {{1},{1,2}}, {{2},{1,1}},
  {{1,2,3}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}, {{1},{2},{3}}.

Andrew Howroyd, Table of n, a(n) for n = 1..50

Cf. A001055, A007716, A045778, A281116, A317449, A317584, A317654, A317755, A317757, A317776.

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]]]];
strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
Table[Length[Select[Join@@mps/@strnorm[n],UnsameQ@@#&]],{n,6}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1,-n)/prod(i=1, #v, i^v[i]*v[i]!)}
seq(n)={my(s); for(k=1, n, forpart(p=k, s+=(-1)^(k+#p)*D(p,n))); s[n]+=1; s/2} \\ Andrew Howroyd, Dec 30 2020

Terms a(10) and beyond from Andrew Howroyd, Dec 30 2020

A318287 Number of non-isomorphic strict multiset partitions of a multiset whose multiplicities are the prime indices of n.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 2, 3, 4, 5, 3, 7, 4, 7, 9, 5, 5, 12, 6, 12, 14, 10, 8, 13, 12, 14, 14, 18, 10, 34

Offset: 1

Gus Wiseman, Aug 23 2018

			Non-isomorphic representatives of the a(20) = 12 strict multiset partitions of {1,1,1,2,3}:
  {{1,1,1,2,3}}
  {{1},{1,1,2,3}}
  {{2},{1,1,1,3}}
  {{1,1},{1,2,3}}
  {{1,2},{1,1,3}}
  {{2,3},{1,1,1}}
  {{1},{2},{1,1,3}}
  {{1},{1,1},{2,3}}
  {{1},{1,2},{1,3}}
  {{2},{3},{1,1,1}}
  {{2},{1,1},{1,3}}
  {{1},{2},{3},{1,1}}

Cf. A001055, A007716, A045778, A181821, A305936, A316980, A317533, A317776, A317791.

Cf. A318283, A318284, A318285, A318286, A318357, A318362, A318371.

a(n) = A318357(A181821(n)).

A317449 Regular triangle where T(n,k) is the number of multiset partitions of strongly normal multisets of size n into k blocks, where a multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

Original entry on oeis.org

1, 2, 2, 3, 6, 3, 5, 21, 16, 5, 7, 52, 72, 32, 7, 11, 141, 306, 216, 65, 11, 15, 327, 1113, 1160, 512, 113, 15, 22, 791, 4033, 6052, 3737, 1154, 199, 22, 30, 1780, 13586, 28749, 24325, 10059, 2317, 323, 30, 42, 4058, 45514, 133642, 151994, 82994, 24854, 4493, 523, 42

Offset: 1

Gus Wiseman, Aug 06 2018

			The T(3,2) = 6 multiset partitions are {{1},{1,1}}, {{1},{1,2}}, {{2},{1,1}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}.
Triangle begins:
    1
    2    2
    3    6    3
    5   21   16    5
    7   52   72   32    7
   11  141  306  216   65   11
   15  327 1113 1160  512  113   15
   ...

Andrew Howroyd, Table of n, a(n) for n = 1..1275 (rows 1..50)

Row sums are A035310. First and last columns are both A000041.

Cf. A001055, A007716, A045778, A255906, A281116, A317584, A317654, A317755, A317775, A317776.

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]]]];
strnorm[n_]:=Flatten[MapIndexed[Table[#2,{#1}]&,#]]&/@IntegerPartitions[n];
Table[Length[Select[Join@@mps/@strnorm[n],Length[#]==k&]],{n,6},{k,n}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
D(p, n)={my(v=vector(n)); for(i=1, #p, v[p[i]]++); my(u=EulerT(v)); Vec(1/prod(k=1, n, 1 - u[k]*x^k + O(x*x^n))-1, -n)/prod(i=1, #v, i^v[i]*v[i]!)}
U(m, n)={my(s=0); forpart(p=m, s+=D(p,n)); s}
M(n)={Mat(vector(n,k,(U(k,n)-U(k-1,n))~))}
{ my(A=M(8)); for(n=1, #A~, print(A[n,1..n])) } \\ Andrew Howroyd, Dec 30 2020

Terms a(46) and beyond from Andrew Howroyd, Dec 30 2020

A330787 Triangle read by rows: T(n,k) is the number of strict multiset partitions of normal multisets of size n into k blocks, where a multiset is normal if it spans an initial interval of positive integers.

Original entry on oeis.org

1, 2, 1, 4, 8, 1, 8, 32, 18, 1, 16, 124, 140, 32, 1, 32, 444, 888, 432, 50, 1, 64, 1568, 5016, 4196, 1060, 72, 1, 128, 5440, 26796, 34732, 15064, 2224, 98, 1, 256, 18768, 138292, 262200, 174240, 44348, 4172, 128, 1, 512, 64432, 698864, 1870840, 1781884, 692668, 112424, 7200, 162, 1

Offset: 1

Andrew Howroyd, Dec 31 2019

			Triangle begins:
    1;
    2,    1;
    4,    8,     1;
    8,   32,    18,     1;
   16,  124,   140,    32,     1;
   32,  444,   888,   432,    50,    1;
   64, 1568,  5016,  4196,  1060,   72,  1;
  128, 5440, 26796, 34732, 15064, 2224, 98, 1;
  ...
The T(3,1) = 4 multiset partitions are {{1,1,1}}, {{1,1,2}}, {{1,2,2}}, {{1,2,3}}.
The T(3,2) = 8 multiset partitions are {{1},{1,1}}, {{1},{2,2}}, {{2},{1,2}}, {{1},{1,2}}, {{2},{1,1}}, {{1},{2,3}}, {{2},{1,3}}, {{3},{1,2}}.
The T(3,3) = 1 multiset partition is {{1},{2},{3}}.

Andrew Howroyd, Table of n, a(n) for n = 1..1275 (first 50 rows)

Row sums are A317776.
Column 1 is A000079(n-1).
Main diagonal is A000012.
Cf. A317532, A327116.

B[n_, k_] := Sum[Binomial[r, k] (-1)^(r-k), {r, k, n}];
row[n_] := Sum[B[n, j] SeriesCoefficient[ Product[(1 + x^k y)^Binomial[k + j - 1, j - 1], {k, 1, n}], {x, 0, n}], {j, 1, n}]/y + O[y]^n // CoefficientList[#, y]&;
Array[row, 10] // Flatten (* Jean-François Alcover, Dec 17 2020, after Andrew Howroyd *)

\\ here B(n, k) is A239473(n, k)
B(n,k)={sum(r=k, n, binomial(r, k)*(-1)^(r-k))}
Row(n)={Vecrev(sum(j=1, n, B(n,j)*polcoef(prod(k=1, n, (1 + x^k*y + O(x*x^n))^binomial(k+j-1,j-1)), n))/y)}
{ for(n=1, 10, print(Row(n))) }

A368099 Triangle read by rows where T(n,k) is the number of non-isomorphic k-element sets of finite nonempty multisets with cardinalities summing to n, or strict multiset partitions of weight n and length k.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 3, 4, 1, 0, 5, 12, 5, 1, 0, 7, 28, 22, 5, 1, 0, 11, 66, 83, 31, 5, 1, 0, 15, 134, 252, 147, 34, 5, 1, 0, 22, 280, 726, 620, 203, 35, 5, 1, 0, 30, 536, 1946, 2283, 1069, 235, 35, 5, 1, 0, 42, 1043, 4982, 7890, 5019, 1469, 248, 35, 5, 1

Offset: 0

Gus Wiseman, Dec 31 2023

			Triangle begins:
    1
    0    1
    0    2    1
    0    3    4    1
    0    5   12    5    1
    0    7   28   22    5    1
    0   11   66   83   31    5    1
    0   15  134  252  147   34    5    1
    0   22  280  726  620  203   35    5    1
    0   30  536 1946 2283 1069  235   35    5    1
    0   42 1043 4982 7890 5019 1469  248   35    5    1
    ...
Row n = 4 counts the following representatives:
  .  {{1,1,1,1}}  {{1},{1,1,1}}  {{1},{2},{1,1}}  {{1},{2},{3},{4}}
     {{1,1,1,2}}  {{1},{1,1,2}}  {{1},{2},{1,2}}
     {{1,1,2,2}}  {{1},{1,2,2}}  {{1},{2},{1,3}}
     {{1,1,2,3}}  {{1},{1,2,3}}  {{1},{2},{3,3}}
     {{1,2,3,4}}  {{1},{2,2,2}}  {{1},{2},{3,4}}
                  {{1},{2,2,3}}
                  {{1},{2,3,4}}
                  {{1,1},{1,2}}
                  {{1,1},{2,2}}
                  {{1,1},{2,3}}
                  {{1,2},{1,3}}
                  {{1,2},{3,4}}

Andrew Howroyd, Table of n, a(n) for n = 0..1325 (rows 0..50)

Row sums are A316980, connected case A319557.
For multiset partitions we have A317533, connected A322133.
Counting connected components instead of edges gives A321194.
For normal multiset partitions we have A330787, row sums A317776.
For set multipartitions we have A334550.
For set-systems we have A368096, row-sums A283877 (connected A300913).
A000110 counts set-partitions, non-isomorphic A000041.
A003465 counts covering set-systems, unlabeled A055621.
A007716 counts non-isomorphic multiset partitions, connected A007718.
A049311 counts non-isomorphic set multipartitions, connected A056156.
A058891 counts set-systems, unlabeled A000612, connected A323818.
Cf. A255903, A296122, A302545, A306005, A317532, A317775, A317794, A317795, A319560, A368094, A368095.

sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]&/@sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
Table[Length[Union[brute /@ Select[mpm[n],UnsameQ@@#&&Length[#]==k&]]], {n,0,5},{k,0,n}]

EulerT(v)={Vec(exp(x*Ser(dirmul(v, vector(#v, n, 1/n))))-1, -#v)}
permcount(v) = {my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
K(q, t, k)={EulerT(Vec(sum(j=1, #q, my(g=gcd(t, q[j])); g*x^(q[j]/g)) + O(x*x^k), -k))}
G(n)={my(s=0); forpart(q=n, my(p=sum(t=1, n, y^t*subst(x*Ser(K(q, t, n\t))/t, x, x^t))); s+=permcount(q)*exp(p-subst(subst(p, x, x^2), y, y^2))); s/n!}
T(n)={[Vecrev(p) | p <- Vec(G(n))]}
{ my(A=T(10)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 11 2024

A383310 Number of ways to choose a strict multiset partition of a factorization of n into factors > 1.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 5, 2, 3, 1, 8, 1, 3, 3, 9, 1, 8, 1, 8, 3, 3, 1, 20, 2, 3, 5, 8, 1, 12, 1, 19, 3, 3, 3, 24, 1, 3, 3, 20, 1, 12, 1, 8, 8, 3, 1, 46, 2, 8, 3, 8, 1, 20, 3, 20, 3, 3, 1, 38, 1, 3, 8, 37, 3, 12, 1, 8, 3, 12, 1, 67, 1, 3, 8, 8, 3, 12, 1, 46, 9, 3

Offset: 1

Gus Wiseman, Apr 26 2025

nonn

			The a(36) = 24 choices:
  {{2,2,3,3}}  {{2},{2,3,3}}  {{2},{3},{2,3}}
  {{2,2,9}}    {{3},{2,2,3}}  {{2},{3},{6}}
  {{2,3,6}}    {{2,2},{3,3}}
  {{2,18}}     {{2},{2,9}}
  {{3,3,4}}    {{9},{2,2}}
  {{3,12}}     {{2},{3,6}}
  {{4,9}}      {{3},{2,6}}
  {{6,6}}      {{6},{2,3}}
  {{36}}       {{2},{18}}
               {{3},{3,4}}
               {{4},{3,3}}
               {{3},{12}}
               {{4},{9}}

The case of a unique choice (positions of 1) is A008578.
This is the strict case of A050336.
For distinct strict blocks we have A050345.
For integer partitions we have A261049, strict case of A001970.
For strict blocks that are not necessarily distinct we have A296119.
Twice-partitions of this type are counted by A296122.
For normal multisets we have A317776, strict case of A255906.
A001055 counts factorizations, strict A045778.
A050320 counts factorizations into squarefree numbers, distinct A050326.
A281113 counts twice-factorizations, strict A296121, see A296118, A296120.
Cf. A000009, A005117, A045782, A050342, A279785, A293243, A293511, A302494, A316439, A358914, A381992, 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@@#&]],{y,facs[n]}],{n,30}]

cp's OEIS Frontend

A318286 Number of strict multiset partitions of a multiset whose multiplicities are the prime indices of n.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A317752 Number of multiset partitions of normal multisets of size n such that the blocks have empty intersection.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A381452 Number of multisets that can be obtained by partitioning the prime indices of n into a set of multisets and taking their sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A327116 Number T(n,k) of colored integer partitions of n using all colors of a k-set such that all parts have different color patterns and a pattern for part i has i colors in (weakly) increasing order; triangle T(n,k), n>=0, 0<=k<=n, read by rows.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A317775 Number of strict multiset partitions of strongly normal multisets of size n, where a multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A318287 Number of non-isomorphic strict multiset partitions of a multiset whose multiplicities are the prime indices of n.

Views

Author

Keywords

Examples

Crossrefs

Formula

A317449 Regular triangle where T(n,k) is the number of multiset partitions of strongly normal multisets of size n into k blocks, where a multiset is strongly normal if it spans an initial interval of positive integers with weakly decreasing multiplicities.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A330787 Triangle read by rows: T(n,k) is the number of strict multiset partitions of normal multisets of size n into k blocks, where a multiset is normal if it spans an initial interval of positive integers.

Views

Author

Keywords

Examples

Links

Crossrefs