OFFSET
1,6
COMMENTS
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}.
LINKS
Robert Price, Table of n, a(n) for n = 1..1000
EXAMPLE
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.
MATHEMATICA
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}]
CROSSREFS
Before taking sums we had A045778.
Other multiset partitions of prime indices:
A000040 lists the primes.
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
