OFFSET
1,4
COMMENTS
Also the number of factorizations of n into prime powers > 1 with equal sums of prime indices.
FORMULA
EXAMPLE
The prime indices of 144 are {1,1,1,1,2,2}, with the following 2 multiset partitions into constant blocks with a common sum:
{{2,2},{1,1,1,1}}
{{2},{2},{1,1},{1,1}}
so a(144) = 2.
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[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]& /@ sps[Range[Length[set]]]];
Table[Length[Select[mps[prix[n]], SameQ@@Total/@#&&And@@SameQ@@@#&]], {n, 100}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 19 2025
STATUS
approved
