OFFSET
1,1
COMMENTS
EXAMPLE
The sequence of terms together with their prime indices begins:
2: {1}
4: {1,1}
8: {1,1,1}
15: {2,3}
16: {1,1,1,1}
32: {1,1,1,1,1}
33: {2,5}
35: {3,4}
45: {2,2,3}
51: {2,7}
55: {3,5}
64: {1,1,1,1,1,1}
69: {2,9}
75: {2,3,3}
77: {4,5}
85: {3,7}
93: {2,11}
95: {3,8}
99: {2,2,5}
119: {4,7}
MATHEMATICA
primeMS[n_]:=If[n==1, {}, Flatten[Cases[FactorInteger[n], {p_, k_}:>Table[PrimePi[p], {k}]]]];
stableQ[u_, Q_]:=!Apply[Or, Outer[#1=!=#2&&Q[#1, #2]&, u, u, 1], {0, 1}];
Select[Range[100], GCD@@primeMS[#]==1&&stableQ[primeMS[#], Divisible]&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Oct 30 2019
STATUS
approved