OFFSET
1,2
COMMENTS
The least term k with nesigma(k)/k > m for m = 2, 3, 4, ... is 480, 332640, 1396755360, ...
MATHEMATICA
esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s[1] = 1 ; s[n_] := DivisorSigma[1, n] - esigma[n]; seq = {}; rm = -1; Do[r1 = s[n]/n; If[r1 > rm, rm = r1; AppendTo[seq, n]], {n, 1, 10^6}]; seq
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Oct 26 2021
STATUS
approved