OFFSET
1,1
REFERENCES
Riele, H. J. J. te; Unitary Aliquot Sequences. MR 139/72, Mathematisch Centrum, 1972, Amsterdam.
Riele, H. J. J. te; Further Results On Unitary Aliquot Sequences. NW 2/73, Mathematisch Centrum, 1973, Amsterdam.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1000
J. O. M. Pedersen, Tables of Aliquot Cycles [Broken link]
J. O. M. Pedersen, Tables of Aliquot Cycles [Via Internet Archive Wayback-Machine]
J. O. M. Pedersen, Tables of Aliquot Cycles [Cached copy, pdf file only]
EXAMPLE
a(5)=812 because the fifth non-unitary amicable number whose unitary aliquot sequence ends in a unitary amicable pair is 812.
MATHEMATICA
UnitaryDivisors[n_Integer?Positive]:=Select[Divisors[n], GCD[ #, n/# ]==1&]; sstar[n_]:=Plus@@UnitaryDivisors[n]-n; g[n_] := If[n > 0, sstar[n], 0]; UnitaryTrajectory[n_] := Most[NestWhileList[g, n, UnsameQ, All]]; UnitaryAmicableNumberQ[k_]:=If[Nest[sstar, k, 2]?k && !sstar[k]?k, True, False]; Select[Range[2500], !UnitaryAmicableNumberQ[ # ] && UnitaryAmicableNumberQ[Last[UnitaryTrajectory[ # ]]] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Ant King, Jan 25 2007
EXTENSIONS
More terms from Amiram Eldar, Apr 06 2019
STATUS
approved