OFFSET
1,1
COMMENTS
First differs from A383100 in lacking 108.
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, sum A056239.
Also numbers that cannot be written as a product of prime powers with equal sums of prime indices.
Partitions of this type are counted by A381993.
EXAMPLE
The terms together with their prime indices begin:
6: {1,2}
10: {1,3}
14: {1,4}
15: {2,3}
18: {1,2,2}
20: {1,1,3}
21: {2,4}
22: {1,5}
24: {1,1,1,2}
26: {1,6}
28: {1,1,4}
30: {1,2,3}
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]]]];
Select[Range[100], Select[mps[prix[#]], SameQ@@Total/@#&&And@@SameQ@@@#&]=={}&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 13 2025
STATUS
approved
