A126169 -id:A126169 - cp's OEIS frontend

A127665 Numbers whose infinitary aliquot sequences end in an infinitary amicable pair.

102, 114, 126, 210, 246, 258, 270, 318, 330, 342, 354, 366, 378, 388, 390, 408, 426, 436, 438, 450, 474, 484, 486, 498, 510, 522, 534, 536, 546, 552, 570, 582, 594, 600, 606, 618, 630, 642, 648, 654, 666, 672, 702, 726, 738, 750, 760, 762, 774, 786, 798

Offset: 1

Ant King, Jan 26 2007

Sometimes called the infinitary 2-cycle attractor set.

			a(5)=246 because 246 is the fifth number whose infinitary aliquot sequence ends in an infinitary amicable pair.

Graeme L. Cohen, On an integer's infinitary divisors, Math. Comp., 54 (1990), 395-411.
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]

Cf. A126168, A007357, A127661, A127662, A127663, A127664, A127666, A127667, A126169, A126170, A126171.

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]];InfinitaryAmicableNumberQ[k_]:=If[Nest[properinfinitarydivisorsum,k,2]==k && !properinfinitarydivisorsum[k]==k,True,False];Select[Range[820],InfinitaryAmicableNumberQ[Last[iTrajectory[ # ]]] &]

A357495 Lesser of a pair of amicable numbers k < m such that s(k) = m and s(m) = k, where s(k) = A162296(k) - k is the sum of aliquot divisors of k that have a square factor.

Original entry on oeis.org

880, 10480, 20080, 24928, 42976, 69184, 110565, 252080, 267712, 489472, 566656, 569240, 603855, 626535, 631708, 687424, 705088, 741472, 786896, 904365, 1100385, 1234480, 1280790, 1425632, 1749824, 1993750, 2012224, 2401568, 2439712, 2496736, 2542496, 2573344, 2671856

Offset: 1

Amiram Eldar, Oct 01 2022

nonn

Analogous to amicable numbers (A002025 and A002046) with nonsquarefree divisors.
The larger counterparts are in A357496.
Both members of each pair are necessarily nonsquarefree numbers.

			880 is a term since s(880) = 1136 and s(1136) = 880.

Amiram Eldar, Table of n, a(n) for n = 1..1000

Cf. A162296, A325314, A322609, A357493, A357494, A357496.
Subsequence of A013929.
Similar sequences: A002025, A002952, A126165, A126169, A259038, A292980, A322541, A324708, A348343.

s[n_] := Module[{f = FactorInteger[n], p, e}, p = f[[;; , 1]]; e = f[[;; , 2]]; Times @@ ((p^(e + 1) - 1)/(p - 1)) - Times @@ (p + 1) - n]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 2, 3*10^6}]; seq

A127664 Infinitary amicable numbers.

Original entry on oeis.org

114, 126, 594, 846, 1140, 1260, 4320, 5940, 7920, 8460, 8640, 10744, 10856, 11760, 12285, 13500, 14595, 17700, 25728, 35712, 43632, 44772, 45888, 49308, 60858, 62100, 62700, 67095, 67158, 71145, 73962, 74784, 79296, 79650, 79750, 83142, 83904, 86400, 88730

Offset: 1

Ant King, Jan 26 2007

nonn

			a(5)=1140 because 1140 is the fifth infinitary amicable number.

Graeme L. Cohen, On an integer's infinitary divisors, Math. Comp., 54 (1990), 395-411.
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]

Cf. A007357, A126168, A127661, A127662, A127663, A127665, A127666, A127667, A126169, A126170, A126171.

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]];InfinitaryAmicableNumberQ[k_]:=If[Nest[properinfinitarydivisorsum,k,2]==k && !properinfinitarydivisorsum[k]==k,True,False];Select[Range[50000],InfinitaryAmicableNumberQ[ # ] &]
fun[p_, e_] := Module[{b = IntegerDigits[e, 2]}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; infs[n_] := Times @@ (fun @@@ FactorInteger[n]) - n; s = {}; Do[k = infs[n]; If[k != n && infs[k] == n, AppendTo[s, n]], {n, 2, 10^5}]; s (* Amiram Eldar, Mar 16 2019 *)

Non-infinitary perfect numbers which satisfy A126168(A126168(n)) = n.

More terms from Amiram Eldar, Mar 16 2019

A333929 Lesser of recursive amicable numbers pair: numbers m < k such that m = s(k) and k = s(m), where s(k) = A333926(k) - k is the sum of proper recursive divisors of k.

Original entry on oeis.org

220, 366, 2620, 3864, 5020, 16104, 16536, 26448, 29760, 43524, 63020, 67344, 69615, 100485, 122265, 142290, 142310, 196248, 196724, 198990, 239856, 240312, 280540, 308620, 309264, 319550, 326424, 341904, 348840, 366792, 469028, 522405, 537744, 580320, 647190, 661776

Offset: 1

Amiram Eldar, Apr 10 2020

nonn

The larger counterparts are in A333930.

			220 is a terms since A333926(220) - 220 = 284 and A333926(284) - 284 = 220.

Amiram Eldar, Table of n, a(n) for n = 1..1000

Cf. A333926, A333927, A333930.

Analogous sequences: A002025, A002952 (unitary), A126165 (exponential), A126169 (infinitary), A292980 (bi-unitary).

recDivQ[n_, 1] = True; recDivQ[n_, d_] := recDivQ[n, d] = Divisible[n, d] && AllTrue[FactorInteger[d], recDivQ[IntegerExponent[n, First[#]], Last[#]] &]; recDivs[n_] := Select[Divisors[n], recDivQ[n, #] &]; f[p_, e_] := 1 + Total[p^recDivs[e]]; recDivSum[1] = 1; recDivSum[n_] := Times @@ (f @@@ FactorInteger[n]); s[n_] := recDivSum[n] - n; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 1, 10^5}]; seq

A371419 Lesser member of Carmichael's variant of amicable pair: numbers k < m such that s(k) = m and s(m) = k, where s(k) = A371418(k).

Original entry on oeis.org

12, 48, 112, 160, 192, 448, 1984, 12288, 28672, 126976, 196608, 458752, 520192, 786432, 1835008, 2031616, 8126464, 8323072, 33292288, 536805376, 2147221504, 3221225472, 7516192768, 33285996544, 34359476224, 136365211648

Offset: 1

Amiram Eldar, Mar 23 2024

Analogous to amicable numbers (A002025 and A002046) with the largest aliquot divisor of the sum of divisors (A371418) instead of the sum of aliquot divisors (A001065).
Carmichael (1921) proposed this function (A371418) for the purpose of studying periodic chains that are formed by repeatedly applying the mapping x -> A371418(x). The chains of cycle 2 are analogous to amicable numbers.
Carmichael noted that if q < p are two different Mersenne exponents (A000043), then 2^(p-1)*(2^q-1) and 2^(q-1)*(2^p-1) are an amicable pair. With the 51 Mersenne exponents that are currently known it is possible to calculate 51 * 50 / 2 = 1275 amicable pairs. (160, 189) is a pair that is not of this "Mersenne form". Are there any other pairs like it? There are no other such pairs with lesser member below a(26).
a(27) <= 8795019280384.
The greater counterparts are in A371420.

			12 is a term since A371418(12) = 14 > 12, and A371418(14) = 12.

Robert D. Carmichael, Empirical Results in the Theory of Numbers, The Mathematics Teacher, Vol. 14, No. 6 (1921), pp. 305-310; alternative link. See p. 309.

Cf. A000043, A001065, A002046, A337874, A337876, A371418, A371420.

Similar sequences: A002025, A002952, A126165, A126169, A259038, A292980, A322541, A324708, A348343.

r[n_] := n/FactorInteger[n][[1, 1]]; s[n_] := r[DivisorSigma[1, n]]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 1, 10^6}]; seq

f(n) = {my(s = sigma(n)); if(s == 1, 1, s/factor(s)[1, 1]);}
lista(nmax) = {my(m); for(n = 1, nmax, m = f(n); if(m > n && f(m) == n, print1(n, ", ")));}

A348602 Smaller member of a nonexponential amicable pair: numbers (k, m) such that nesigma(k) = m and nesigma(m) = k, where nesigma(k) is the sum of the nonexponential divisors of k (A160135).

Original entry on oeis.org

198, 18180, 142310, 1077890, 1156870, 1511930, 1669910, 2236570, 2728726, 3776580, 4246130, 4532710, 5123090, 5385310, 6993610, 7288930, 8619765, 8754130, 8826070, 9478910, 10254970, 14426230, 17041010, 17257695, 21448630, 30724694, 34256222, 35361326, 37784810

Offset: 1

Amiram Eldar, Oct 25 2021

nonn

The larger counterparts are in A348603.

			198 is a term since A160135(198) = 204 and A160135(204) = 198.

Cf. A160135, A348601, A348603.

Similar sequences: A002025, A002952, A126165, A126169, A259038, A292980, A322541, A324708, A348343.

esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s[n_] := DivisorSigma[1, n] - esigma[n]; seq = {}; Do[m = s[n]; If[m > n && s[m] == n, AppendTo[seq, n]], {n, 1, 1.7*10^6}]; seq

A124663 Number of reduced infinitary amicable pairs (i,j) with i

Original entry on oeis.org

0, 0, 0, 1, 2, 3, 14, 25, 51, 112, 213

Offset: 1

Ant King, Dec 23 2006

A divisor of n is called infinitary if it is a product of divisors of the form p^{y_a 2^a}, where p^y is a prime power dividing n and sum_a y_a 2^a is the binary representation of y.

			a(7)=14 because there are 14 reduced infinitary amicable pairs (m,n) with m<n and m<=10^7

Pedersen J. M., Known amicable pairs.

Cf. A126169, A049417, A126168, A037445, A126170, A126171, A126172, A126173, A126174, A126175, A126176.

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; ReducedInfinitaryAmicableNumberQ[n_] := If[properinfinitarydivisorsum[properinfinitarydivisorsum[ n] - 1] == n + 1 && n > 1, True, False]; ReducedInfinitaryAmicablePairList[k_] := (anlist = Select[Range[k], ReducedInfinitaryAmicableNumberQ[ # ] &]; prlist = Table[Sort[{anlist[[n]], properinfinitarydivisorsum[anlist[[n]]] - 1}], {n, 1, Length[anlist]}]; amprlist = Union[prlist, prlist]); data1 = ReducedInfinitaryAmicablePairList[10^7]; Table[Length[Select[data1, First[ # ] < 10^k &]], {k, 1, 7}]

Reduced infinitary amicable pairs (m,n) satisfy isigma(m)=isigma(n)=m+n+1, with m

A361811 Smallest members of infinitary sociable quadruples.

Original entry on oeis.org

1026, 10098, 10260, 41800, 45696, 100980, 241824, 685440, 4938136, 13959680, 14958944, 25581600, 28158165, 32440716, 36072320, 55204500, 74062944, 81128632, 149589440, 178327008, 192793770, 209524210, 283604220, 319848642, 498215416, 581112000, 740629440, 1236402232

Offset: 1

Amiram Eldar, Mar 25 2023

nonn

The first 8 terms were found by Cohen (1990).

			1026 is a term since the iterations of the sum of aliquot infinitary divisors function (A126168) that start with 1026 are cyclic with period 4: 1026, 1374, 1386, 1494, 1026, ..., and 1026 is the smallest member of the quadruple.
The first five quadruples are {1026, 1374, 1386, 1494}, {10098, 15822, 19458, 15102}, {10260, 13740, 13860, 14940}, {41800, 51800, 66760, 83540}, {45696, 101184, 94656, 88944}.

Graeme L. Cohen, On an integer's infinitary divisors, Mathematics of Computation, Vol. 54, No. 189 (1990), pp. 395-411.
Jan Munch Pedersen, Known Infinitary Sociable Numbers of order four.

Cf. A049417, A126168.
Cf. A007357 (period 1), A126169 and A126170 (period 2).
Subsequence of A004607 (all cycles of length > 2).
Similar sequences: A090615 (all divisors), A319902 (unitary), A319915 (bi-unitary).

f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]]>0, 1 + p^(2^(m-j)), 1], {j, 1, m}]]; infs[n_] := Times @@ f @@@ FactorInteger[n] - n;  infs[1] = 0; seq[n_] := NestList[infs, n, 4][[2;; 5]] ; q[n_] := Module[{s = seq[n]}, n == Min[s] && Count[s, n] == 1]; Select[Range[10^6], q]

infs(n) = {my(f = factor(n), b); prod(i=1, #f~, b = binary(f[i, 2]); prod(k=1, #b, if(b[k], f[i, 1]^(2^(#b-k)) + 1, 1))) - n; }
is(n) = {my(m = n); for(k = 1, 4, m = infs(m); if(k < 4 && m <= n, return(0))); m == n; }

Previous Showing 11-18 of 18 results.

cp's OEIS Frontend

A127665 Numbers whose infinitary aliquot sequences end in an infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A357495 Lesser of a pair of amicable numbers k < m such that s(k) = m and s(m) = k, where s(k) = A162296(k) - k is the sum of aliquot divisors of k that have a square factor.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A127664 Infinitary amicable numbers.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A333929 Lesser of recursive amicable numbers pair: numbers m < k such that m = s(k) and k = s(m), where s(k) = A333926(k) - k is the sum of proper recursive divisors of k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A371419 Lesser member of Carmichael's variant of amicable pair: numbers k < m such that s(k) = m and s(m) = k, where s(k) = A371418(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A348602 Smaller member of a nonexponential amicable pair: numbers (k, m) such that nesigma(k) = m and nesigma(m) = k, where nesigma(k) is the sum of the nonexponential divisors of k (A160135).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A124663 Number of reduced infinitary amicable pairs (i,j) with i

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A361811 Smallest members of infinitary sociable quadruples.

Views

Author

Keywords

Comments

Examples

Links