cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A127660 Integers whose exponential aliquot sequences end in an exponential amicable pair.

Original entry on oeis.org

90972, 100548, 454860, 502740, 937692, 968436, 1000692, 1106028, 1182636, 1307124, 1383732, 1536416, 1546524, 1709316, 2092356, 2312604, 2502528, 2638188, 2690100, 2820132, 2915892, 3116988, 3365964, 3720276, 3729852, 3907008, 3911796, 4122468, 4248552, 4275684
Offset: 1

Views

Author

Ant King, Jan 25 2007

Keywords

Comments

Sometimes called the exponential 2-cycle attractor set. The first 10 terms of this sequence are the same as the first 10 terms of A127659.

Examples

			a(11) = 1383732 because the eleventh integer whose exponential aliquot sequence ends in an exponential amicable pair is 1383732.
		

Crossrefs

Subsequences: A127659, A126165, A126166.

Programs

  • 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[Last[eTrajectory[ # ]]] &]
    f[p_, e_] := DivisorSum[e, p^# &]; s[0] = s[1] = 0; s[n_] := Times @@ f @@@ FactorInteger[n] - n; q[n_] := Module[{v = NestWhileList[s, n, UnsameQ, All]}, v[[-2]] != v[[-1]] > 0 && v[[-3]] == v[[-1]]]; Select[Range[10^6], q] (* Amiram Eldar, Mar 11 2023 *)