login
Integers whose infinitary aliquot sequences end in an infinitary perfect number (A007357).
5

%I #18 Mar 11 2023 05:43:00

%S 6,30,42,54,60,66,72,78,90,100,140,148,152,192,194,196,208,220,238,

%T 244,252,268,274,292,296,298,300,336,348,350,360,364,372,374,380,382,

%U 386,400,416,420,424,476,482,492,516,520,532,540,542,544,550,572,576,578,586

%N Integers whose infinitary aliquot sequences end in an infinitary perfect number (A007357).

%H Amiram Eldar, <a href="/A127662/b127662.txt">Table of n, a(n) for n = 1..75</a>

%H Graeme L. Cohen, <a href="http://dx.doi.org/10.1090/S0025-5718-1990-0993927-5">On an integer's infinitary divisors</a>, Math. Comp., 54 (1990), 395-411.

%H J. O. M. Pedersen, <a href="http://amicable.homepage.dk/tables.htm">Tables of Aliquot Cycles</a>. [Broken link]

%H J. O. M. Pedersen, <a href="http://web.archive.org/web/20140502102524/http://amicable.homepage.dk/tables.htm">Tables of Aliquot Cycles</a>. [Via Internet Archive Wayback-Machine]

%H J. O. M. Pedersen, <a href="/A063990/a063990.pdf">Tables of Aliquot Cycles</a>. [Cached copy, pdf file only]

%e a(5) = 60 because the fifth number whose infinitary aliquot sequence ends in an infinitary perfect number is 60.

%e 6 -> 6 ...

%e 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 60 -> 60 -> ..

%e 66 -> 78 -> 90 -> 90 -> ..

%e 72 -> 78 -> 90 -> 90 -> ..

%e 78 -> 90 -> 90 -> ..

%e 90 -> 90 -> ..

%e 100 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 102 -> 114 -> 126 -> 114 -> .. cycle but not in the sequence

%e 114 -> 126 -> 114 -> .. cycle but not in the sequence

%e 126 -> 114 -> 126 -> ..

%e 140 -> 100 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 148 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 152 -> 148 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 192 -> 148 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 194 -> 100 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 196 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 208 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 210 -> 366 -> 378 -> 582 -> 594 -> 846 -> 594 -> ..

%e 220 -> 140 -> 100 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 238 -> 194 -> 100 -> 30 -> 42 -> 54 -> 66 -> 78 -> 90 -> 90 -> ..

%e 244 -> 66 -> 78 -> 90 -> 90 -> ..

%e 246 -> 258 -> 270 -> 450 -> 330 -> 534 -> 546 -> 798 -> 1122 -> 1470 -> 2130 -> 3054 -> 3066 -> 4038 -> 4050 -> 2346 -> 2838 -> 3498 -> 4278 -> 4938 -> 4950 -> 4410 -> 4590 -> 8370 -> 14670 -> 14850 -> 22590 -> 22770 -> 29070 -> 35730 -> 35910 -> 79290 -> 79470 -> 79650 -> 107550 -> 79650 -> ..

%p isA007357 := proc(n)

%p A049417(n) = 2*n ;

%p simplify(%) ;

%p end proc:

%p isA127662 := proc(n)

%p local trac,x;

%p x := n ;

%p trac := [x] ;

%p while true do

%p x := A049417(x)-trac[-1] ;

%p if x = 0 then

%p return false ;

%p elif x in trac then

%p return isA007357(x) ;

%p end if;

%p trac := [op(trac),x] ;

%p end do:

%p end proc:

%p for n from 1 do

%p if isA127662(n) then

%p printf("%d,\n",n) ;

%p end if;

%p end do: # _R. J. Mathar_, Oct 05 2017

%t ExponentList[n_Integer,factors_List]:={#,IntegerExponent[n,# ]}&/@factors;InfinitaryDivisors[1]:={1}; InfinitaryDivisors[n_Integer?Positive]:=Module[ { factors=First/@FactorInteger[n], d=Divisors[n] }, d[[Flatten[Position[ Transpose[ Thread[Function[{f,g}, BitOr[f,g]==g][ #,Last[ # ]]]&/@ Transpose[Last/@ExponentList[ #,factors]&/@d]],_?(And@@#&),{1}]] ]] ] Null;properinfinitarydivisorsum[k_]:=Plus@@InfinitaryDivisors[k]-k;g[n_] := If[n > 0,properinfinitarydivisorsum[n], 0];iTrajectory[n_] := Most[NestWhileList[g, n, UnsameQ, All]];InfinitaryPerfectNumberQ[0]=False;InfinitaryPerfectNumberQ[k_Integer] :=If[properinfinitarydivisorsum[k]==k,True,False];Select[Range[500],InfinitaryPerfectNumberQ[Last[iTrajectory[ # ]]] &]

%t s[n_] := Times @@ (1 + Power @@@ FactorInteger[n]) - n; s[0] = s[1] = 0; q[n_] := Module[{v = NestWhileList[s, n, UnsameQ, All]}, v[[-1]] != n && v[[-2]] == v[[-1]] > 0]; Select[Range[3200], q] (* _Amiram Eldar_, Mar 11 2023 *)

%Y Cf. A007357, A126168, A127661 - A127667.

%K hard,nonn

%O 1,1

%A _Ant King_, Jan 26 2007

%E More terms from _Amiram Eldar_, Mar 11 2023