OFFSET
1,1
COMMENTS
REFERENCES
Hagis, Peter Jr.; Some Results Concerning Exponential Divisors, Internat. J. Math. & Math. Sci., Vol. 11, No. 2, (1988), pp. 343-350.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
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
MATHEMATICA
ExponentialDivisors[1]={1}; ExponentialDivisors[n_]:=Module[{}, {pr, pows}=Transpose@FactorInteger[n]; divpowers=Distribute[Divisors[pows], List]; Sort[Times@@(pr^Transpose[divpowers])]]; se[n_]:=Plus@@ExponentialDivisors[n]-n; g[n_] := If[n > 0, se[n], 0]; eTrajectory[n_] := Most[NestWhileList[g, n, UnsameQ, All]]; ExponentialAmicableNumberQ[k_]:=If[Nest[se, k, 2]==k && !se[k]==k, True, False]; Select[Range[5 10^6], ExponentialAmicableNumberQ[ # ] &]
fun[p_, e_] := DivisorSum[e, p^# &]; esigma[1] = 1; esigma[n_] := Times @@ fun @@@ FactorInteger[n]; s = {}; Do[m = esigma[n] - n; If[m != n && esigma[m] - m == n, AppendTo[s, n]], {n, 1, 10^6}]; s (* Amiram Eldar, May 09 2019 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Ant King, Jan 25 2007
EXTENSIONS
Link corrected by Andrew Lelechenko, Dec 04 2011
More terms from Amiram Eldar, May 09 2019
STATUS
approved