OFFSET
1,1
COMMENTS
Integers such that A247015(n) = 1.
EXAMPLE
sigma(8)/8 is greater than all sigma(x)/x when x < 8 except 6; so 8 is here.
MAPLE
M1:= 3/2: M2:= 1: c1:= 1:
Res:= NULL: count:= 0:
for n from 3 while count < 20 do
v:= numtheory:-sigma(n)/n;
if v > M1 then M2:= M1; M1:= v; c1:= 1
elif v = M1 then
c1:= c1+1
elif c1 = 1 and v >= M2 then
M2:= v;
Res:= Res, n: count:= count+1
fi
od:
Res; # Robert Israel, Jul 28 2020
PROG
(PARI) lista(nn) = {my(t=1, x=3/2, y); for(m=3, nn, if((g=sigma(m)/m)>x, t=1; y=x; x=g, if(g==x, t=0, if(g>=y&&t, y=g; print1(m, ", "))))); } \\ Jinyuan Wang, Jul 28 2020
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Sep 09 2014
EXTENSIONS
a(15)-a(21) from Robert Israel, Jun 08 2018
Corrected and name changed by Robert Israel, Jul 28 2020
More terms from Jinyuan Wang, Jul 28 2020
STATUS
approved