OFFSET
1,2
COMMENTS
A prime index of n is a number m such that prime(m) divides n.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
The sequence of terms together with their prime indices begins:
1: {}
2: {1}
3: {2}
4: {1,1}
5: {3}
7: {4}
8: {1,1,1}
9: {2,2}
10: {1,3}
11: {5}
13: {6}
14: {1,4}
16: {1,1,1,1}
17: {7}
19: {8}
20: {1,1,3}
21: {2,4}
22: {1,5}
23: {9}
25: {3,3}
MAPLE
filter:= proc(n) andmap(t -> not ((n/numtheory:-pi(t))::integer), numtheory:-factorset(n) minus {2}) end proc:
select(filter, [$1..200]); # Robert Israel, Mar 20 2019
MATHEMATICA
Select[Range[100], !Or@@Cases[If[#==1, {}, FactorInteger[#]], {p_, _}:>If[p==2, False, Divisible[#, PrimePi[p]]]]&]
PROG
(PARI) is(n) = my(f=factor(n)[, 1]~, idc=[]); for(k=1, #f, idc=concat(idc, [primepi(f[k])])); for(t=1, #idc, if(idc[t]==1, next); if(n%idc[t]==0, return(0))); 1 \\ Felix Fröhlich, Mar 21 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Mar 18 2019
STATUS
approved