OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Unitary Divisor
MATHEMATICA
mx = -1; t = {}; Do[u = DivisorSum[n, # &, GCD[#, n/#] == 1 &] - n; If[u > mx, mx = u; AppendTo[t, n]], {n, 8000}]; t
PROG
(PARI) sumud(n) = sumdiv(n, d, if (gcd(d, n/d)==1, d)) - n;
lista(nn) = {lasts = -1; for (n=1, nn, if ((news = sumud(n)) > lasts, print1(n, ", "); lasts = news); ); } \\ Michel Marcus, Aug 04 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Ilya Gutkovskiy, Aug 03 2017
STATUS
approved