OFFSET
1,4
COMMENTS
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.
LINKS
EXAMPLE
The a(47) = 11 factorizations of 60 together with the corresponding multiset partitions of {1,1,2,3}:
(2*2*3*5) {{1},{1},{2},{3}}
(2*2*15) {{1},{1},{2,3}}
(2*3*10) {{1},{2},{1,3}}
(2*5*6) {{1},{3},{1,2}}
(2*30) {{1},{1,2,3}}
(3*4*5) {{2},{1,1},{3}}
(3*20) {{2},{1,1,3}}
(4*15) {{1,1},{2,3}}
(5*12) {{3},{1,1,2}}
(6*10) {{1,2},{1,3}}
(60) {{1,1,2,3}}
MATHEMATICA
nn=100;
facsusing[s_, n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facsusing[Select[s, Divisible[n/d, #]&], n/d], Min@@#>=d&]], {d, Select[s, Divisible[n, #]&]}]];
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
y=Select[Range[nn], PrimeQ[#]||CoprimeQ@@primeMS[#]&];
Table[Length[facsusing[Rest[y], n]], {n, y}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 20 2019
STATUS
approved