A126170 -id:A126170 - cp's OEIS frontend

A357496 Greater 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.

1136, 11696, 22256, 25472, 43424, 73664, 131355, 304336, 267968, 492608, 612704, 674920, 640305, 788697, 691292, 705344, 723392, 813728, 809776, 1117395, 1258335, 1559696, 1518570, 1598368, 1821376, 2218250, 2058944, 2678752, 2744288, 2765024, 2848864, 2610656, 3134224

Offset: 1

Amiram Eldar, Oct 01 2022

nonn

Analogous to amicable numbers (A002025 and A002046) with nonsquarefree divisors.
The terms are ordered according to their lesser counterparts (A357495).
Both members of each pair are necessarily nonsquarefree numbers.

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

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

Cf. A162296, A325314, A322609, A357493, A357494, A357495.
Subsequence of A013929.
Similar sequences: A002046, A002953, A126166, A126170, A259039, A292981, A322542, A324709, A348344.

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, m]], {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

A333930 Larger 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

284, 378, 2924, 4584, 5564, 16632, 16728, 28752, 30912, 53692, 76084, 69552, 87633, 124155, 139815, 179118, 168730, 225096, 202444, 256338, 245904, 266568, 365084, 389924, 320016, 430402, 391656, 353616, 387720, 393528, 486178, 525915, 555216, 642720, 814698, 682896

Offset: 1

Amiram Eldar, Apr 10 2020

nonn

The terms are ordered according to their lesser counterparts (A333929).

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

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

Cf. A333926, A333927, A333929.

Analogous sequences: A002046, A002953 (unitary), A126166 (exponential), A126170 (infinitary), A292981 (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, m]], {n, 1, 10^5}]; seq

A371420 Greater 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

14, 62, 124, 189, 254, 508, 2032, 16382, 32764, 131056, 262142, 524284, 524224, 1048574, 2097148, 2097136, 8388592, 8388544, 33554368, 536866816, 2147479552, 4294967294, 8589934588, 34359738352, 34359672832, 137438953408

Offset: 1

Amiram Eldar, Mar 23 2024

The terms are ordered according to their lesser counterparts (A371419).

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

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. A371418, A371419.

Similar sequences: A002046, A002953, A126166, A126170, A259039, A292981, A322542, A324709, A348344.

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, m]], {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(m, ", ")));}

A348603 Larger 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

204, 19332, 168730, 1099390, 1292570, 1598470, 2062570, 2429030, 3077354, 3903012, 4488910, 6135962, 5504110, 5812130, 7158710, 8221598, 9627915, 10893230, 10043690, 11049730, 10273670, 18087818, 19150222, 17578785, 23030090, 32174506, 35997346, 40117714, 39944086

Offset: 1

Amiram Eldar, Oct 25 2021

nonn

The terms are ordered according to their smaller counterparts (A348602).

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

Cf. A160135, A348601, A348602.

Similar sequences: A002046, A002953, A126166, A126170, A259039, A292981, A322542, A324709, A348344.

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, m]], {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-17 of 17 results.

cp's OEIS Frontend

A357496 Greater 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

A333930 Larger 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

A371420 Greater 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

A348603 Larger 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

Crossrefs

Programs

Mathematica

PARI