OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..666
MAPLE
R:= NULL: count:= 0: p:= 1:
for n from 1 while count < 100 do
p:= nextprime(p);
if (p - n*numtheory:-tau(n))^2 < n then
R:= R, n; count:= count+1;
fi
od:
R; # Robert Israel, Jul 23 2025
MATHEMATICA
Reap[For[n = 2, n <= 40000, n++, If[Abs[Prime[n] - n*DivisorSigma[0, n]] < Sqrt[n], Sow[n]]]][[2, 1]] (* Jean-François Alcover, Nov 07 2016 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Mar 22 2002
STATUS
approved
