OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..750
Eric Weisstein's World of Mathematics, Unitary Divisor.
EXAMPLE
MATHEMATICA
UnitaryDivisors[n_Integer?Positive]:=Select[Divisors[n], GCD[ #, n/# ]==1&]; sstar[n_]:=Plus@@UnitaryDivisors[n]-n; RunningMaxima[l_]:=Rest[FoldList[Max, -Infinity, l]] HighWaterMarks[l_]:=Module[{s=Split[RunningMaxima[l]]}, {First/@s, Most[FoldList[Plus, 1, Length/@s]]} ]; data1=Select[Range[20000], sstar[ # ]-#>0 &]; data2=sstar[ # ]-# &/@data1; pos=Last[HighWaterMarks[data2]]; champs=data1[[ # ]] &/@pos
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Ant King, Apr 20 2007
STATUS
approved