OFFSET
1,2
COMMENTS
The least term k with A348271(k)/k > m for m = 1, 2, 3, .... is 36, 3600, 1587600, ...
MATHEMATICA
f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := DivisorSigma[1, n] - isigma[n]; seq = {}; rm = -1; Do[r1 = s[n]/n; If[r1 > rm, rm = r1; AppendTo[seq, n]], {n, 1, 10^6}]; seq
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Oct 09 2021
STATUS
approved