OFFSET
1,2
COMMENTS
Because of the bound d(m) <= 2*sqrt(m), in order for k to be in the sequence it suffices that d(k)/k >= d(m)/m for k < m < (2*k/d(k))^2. - Robert Israel, Jan 23 2023
REFERENCES
David desJardins, Posting to Math Fun Mailing List, Jun 21 2022.
MAPLE
N:= 10^6:
Q:= [seq(numtheory:-tau(k)/k, k=1..N)]:
V:= Vector(10^6):
r:= 2/10^3:
for n from 10^6 to 1 by -1 do
r:= max(Q[n], r);
V[n]:= r;
od:
select(i -> Q[i] >= V[i+1], [$1..10^6-1]); # Robert Israel, Jan 23 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Jun 21 2022
EXTENSIONS
More terms from Robert Israel, Jan 23 2023
STATUS
approved