OFFSET
1,2
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, sum A056239. The multiset of multisets with MM-number n is formed by taking the multiset of prime indices of each part of the multiset of prime indices of n. For example, the prime indices of 78 are {1,2,6}, so the multiset of multisets with MM-number 78 is {{},{1},{1,2}}.
LINKS
Robert Price, Table of n, a(n) for n = 1..1091
EXAMPLE
The terms together with their prime indices of prime indices begin:
1: {}
2: {{}}
3: {{1}}
4: {{},{}}
5: {{2}}
7: {{1,1}}
8: {{},{},{}}
9: {{1},{1}}
11: {{3}}
16: {{},{},{},{}}
17: {{4}}
19: {{1,1,1}}
23: {{2,2}}
25: {{2},{2}}
27: {{1},{1},{1}}
31: {{5}}
32: {{},{},{},{},{}}
35: {{2},{1,1}}
41: {{6}}
49: {{1,1},{1,1}}
53: {{1,1,1,1}}
59: {{7}}
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[100], SameQ@@Total/@prix/@prix[#] && And@@SameQ@@@prix/@prix[#]&]
PROG
(PARI) is(k) = my(f=factor(k)[, 1]~, k, p, v=vector(#f, i, primepi(f[i]))); for(i=1, #v, k=isprimepower(v[i], &p); if(k||v[i]==1, v[i]=k*primepi(p), return(0))); #Set(v)<2; \\ Jinyuan Wang, Apr 02 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 21 2025
STATUS
approved
