OFFSET
1,6
COMMENTS
First differs from A050326 at a(210) = 13, A050326(210) = 15. This comes from the set systems {{3},{1,2,4}} and {{1,2},{3,4}}, and from {{4},{1,2,3}} and {{1,3},{2,4}}.
Also the number of multisets that can be obtained by taking the sums of prime indices of each factor in a strict factorization of n into squarefree numbers > 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 ranked 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 sets are generally not transitive. For example, we have arrows: {{1},{1,2}}: {1,1,2} -> {1,3} and {{1,3}}: {1,3} -> {4}, but there is no set of sets {1,1,2} -> {4}.
LINKS
Robert Price, Table of n, a(n) for n = 1..1000
EXAMPLE
The prime indices of 60 are {1,1,2,3}, with partitions into sets of sets:
{{1},{1,2,3}}
{{1,2},{1,3}}
{{1},{2},{1,3}}
{{1},{3},{1,2}}
with block-sums: {1,6}, {3,4}, {1,2,4}, {1,3,3}, which are all different, so a(60) = 4.
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
facs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Table[Length[Union[Sort[Total/@prix/@#]&/@Select[facs[n], UnsameQ@@#&&And@@SquareFreeQ/@#&]]], {n, 100}]
CROSSREFS
Positions of 0 are A293243.
Positions of 1 are A293511.
Other multiset partitions of prime indices:
A003963 gives product of prime indices.
A122111 represents conjugation in terms of Heinz numbers.
A265947 counts refinement-ordered pairs of integer partitions.
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 06 2025
STATUS
approved
