OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..1500
MAPLE
S:= [$1..10^5]: nS:= 10^5:
for k from 2 while k <= nS do
R:= [seq(i*k, i=1..floor(nS/k))];
S:= subsop(seq(r=NULL, r=R), S);
nS:= nS-floor(nS/k);
od:
select(isprime, S); # Robert Israel, Mar 01 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Aminos Abijad, Feb 01 2026
STATUS
approved
