OFFSET
1,2
COMMENTS
MATHEMATICA
f[p_, e_] := p^e - 1; uphi[1] = 1; uphi[n_] := Times @@ f @@@ FactorInteger[n]; r[n_] := Plus @@ DeleteDuplicates[uphi /@ Select[Divisors[n], CoprimeQ[#, n/#] &]]/n; rm = 2; seq = {}; Do[r1 = r[n]; If[r1 < rm, rm = r1; AppendTo[seq, n]], {n, 1, 2*10^5}]; seq
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Oct 04 2021
STATUS
approved