OFFSET
1,1
COMMENTS
All terms are odd.
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.
EXAMPLE
The prime indices of 87 are {2,10}, with divisors {{1,2},{1,2,5,10}}, with intersection {1,2}, so 87 is in the sequence.
The prime indices of 91 are {4,6}, with divisors {{1,2,4},{1,2,3,6}}, with intersection {1,2}, so 91 is in the sequence.
The terms together with their prime indices begin:
3: {2}
5: {3}
9: {2,2}
11: {5}
17: {7}
21: {2,4}
25: {3,3}
27: {2,2,2}
31: {11}
39: {2,6}
41: {13}
57: {2,8}
59: {17}
63: {2,2,4}
65: {3,6}
67: {19}
81: {2,2,2,2}
MATHEMATICA
prix[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
Select[Range[2, 100], Length[Intersection@@Divisors/@prix[#]]==2&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Aug 21 2025
STATUS
approved
