OFFSET
1,2
MATHEMATICA
usigma[1] = 1; usigma[n_] := Times @@ (1 + Power @@@ FactorInteger[n]); nusigma[n_] := DivisorSigma[1, n] - usigma[n]; rm = -1; s = {}; Do[r = nusigma[n]/n; If[r > rm, rm = r; AppendTo[s, n]], {n, 1, 10000}]; s
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Nov 23 2019
STATUS
approved