A126174 -id:A126174 - cp's OEIS frontend

A126173 Larger element of a reduced infinitary amicable pair.

2295, 75495, 817479, 1902215, 1341495, 1348935, 2226014, 2421704, 3123735, 3010215, 5644415, 4282215, 7509159, 10106504, 12900734, 24519159, 31356314, 41950359, 43321095, 80870615, 42125144, 85141719, 87689415, 87802407, 86477895, 105993657, 168669879, 129081735

Offset: 1

Ant King, Dec 23 2006

nonn

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(3)=817479 because 817479 is the largest member of the third reduced infinitary amicable pair, (573560,817479)

Amiram Eldar, Table of n, a(n) for n = 1..278
Jan Munch Pedersen, Tables of Aliquot Cycles.

Cf. A126169, A049417, A126168, A037445, A126170, A126171, A126172, 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[Last[data1[[k]]], {k, 1, Length[data1]}]
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] - 1; If[k > n && infs[k] == n + 1, AppendTo[s, k]], {n, 2, 10^5}]; s (* Amiram Eldar, Jan 22 2019 *)

The values of n for which isigma(m)=isigma(n)=m+n+1, where n>m and isigma(n) is given by A049417(n).

a(15)-a(28) from Amiram Eldar, Jan 22 2019

A126175 Larger member of an augmented infinitary amicable pair.

Original entry on oeis.org

1483785, 2479065, 2580105, 4895241, 7336455, 9100905, 10350345, 16367481, 17307105, 24829945, 15706090, 27866241, 15439545, 23872185, 53763535, 63075321, 41337555, 60923577, 51394665, 56802249, 110691295, 73809496, 89870985, 82771336, 92586585, 150672921, 108212055

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(3)=2580105 because 2580105 is the larger member of the third augmented infinitary amicable pair, (2166136,2580105).

Amiram Eldar, Table of n, a(n) for n = 1..276
Jan Munch Pedersen, Tables of Aliquot Cycles.

Cf. A126169, A049417, A126168, A037445, A126170, A126171, A126173, A126174, 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; AugmentedInfinitaryAmicableNumberQ[n_] := If[properinfinitarydivisorsum[properinfinitarydivisorsum[ n] + 1] == n - 1 && ! properinfinitarydivisorsum[n] + 1 == n, True, False]; AugmentedInfinitaryAmicablePairList[k_] := (anlist = Select[Range[k], AugmentedInfinitaryAmicableNumberQ[ # ] &]; prlist = Table[ Sort[{anlist[[n]], properinfinitarydivisorsum[anlist[[n]]] + 1}], {n, 1, Length[anlist]}]; amprlist = Union[prlist, prlist]); data = AugmentedInfinitaryAmicablePairList[10^7]; Table[Last[data[[k]]], {k, 1, Length[data]}]
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] + 1; If[k > n && infs[k] == n - 1, AppendTo[s, k]], {n, 2, 10^9}]; s (* Amiram Eldar, Jan 20 2019 *)

The values of n for which isigma(m)=isigma(n)=m+n-1, where n>m and isigma(n) is given by A049417(n).

a(9)-a(27) from Amiram Eldar, Jan 20 2019

A126172 Smaller element of a reduced infinitary amicable pair.

Original entry on oeis.org

2024, 62744, 573560, 1000824, 1173704, 1208504, 1921185, 2140215, 2198504, 2312024, 2580864, 4012184, 5416280, 9247095, 12500865, 13496840, 23939685, 26409320, 34093304, 37324584, 40818855, 52026920, 66275384, 76011992, 79842104, 101366342, 101589320, 106004024

Offset: 1

Ant King, Dec 23 2006

nonn

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(3)=573560 because 573560 is the smaller element of the third reduced infinitary amicable pair, (573560, 817479)

Amiram Eldar, Table of n, a(n) for n = 1..278
Jan Munch Pedersen, Tables of Aliquot Cycles.

Cf. A126169, A049417, A126168, A037445, A126170, A126171, 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[First[data1[[k]]], {k, 1, Length[data1]}]
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] - 1; If[k > n && infs[k] == n + 1, AppendTo[s, n]], {n, 2, 10^5}]; s (* Amiram Eldar, Jan 22 2019 *)

The values of m for which isigma(m)=isigma(n)=m+n+1, where mA049417(n).

a(15)-a(28) from Amiram Eldar, Jan 22 2019

A126176 Number of augmented infinitary amicable pairs (i,j) with i

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 8, 26, 48, 104, 227

Offset: 1

Ant King, Dec 24 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(9)=48 because there are 48 augmented infinitary amicable pairs (m,n) with m<n and m<=10^9

Pedersen J. M., Known amicable pairs.

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

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

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

A306867 Lesser of augmented bi-unitary amicable pair.

Original entry on oeis.org

434784, 1100176, 1252216, 1754536, 2166136, 2362360, 3064096, 6224890, 7626136, 7851256, 7950096, 9026235, 9581320, 12494856, 13324311, 14192080, 15218560, 15243424, 15422536, 19028296, 19466560, 19555360, 29180466, 36716680, 37542190, 40682824, 44044000, 44588896

Offset: 1

Amiram Eldar, Mar 14 2019

nonn

A pair m < n is an augmented bi-unitary amicable pair if bsigma(m) = bsigma(n) = m + n - 1, where bsigma(n) is the sum of bi-unitary divisors of n (A188999).

The larger members are in A306868.

			434784 is in the sequence since it is the lesser of the amicable pair (434784, 871585): bsigma(434784) = bsigma(871585) = 1306368 = 434784 + 871585 - 1.

Amiram Eldar, Table of n, a(n) for n = 1..509 (terms below 2*10^11, from David Moews's site).
David Moews, Perfect, amicable and sociable numbers.

Cf. A007992, A126174, A188999, A292980, A306868.

fun[p_, e_]:=If[OddQ[e], (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1)-p^(e/2)]; bsigma[1] = 1; bsigma[n_] := Times @@ (fun @@@ FactorInteger[n]); f[n_] := bsigma[n] - n + 1; s={}; Do[m = f[n]; If[m > n && f[m] == n, AppendTo[s, n]], {n, 1, 10^7}]; s

A306872 Lesser of augmented unitary amicable pair.

Original entry on oeis.org

6224890, 37542190, 56523810, 101304490, 135657795, 233990890, 5304907426, 8473747030, 8483430670, 9220653310, 11033448910, 12139959910, 13108452735, 13849730895, 16697472870, 19644687195, 20321234206, 23076788295, 40575765615, 55636542346, 89094853155, 101786530846

Offset: 1

Amiram Eldar, Mar 14 2019

nonn

A pair m < n is an augmented unitary amicable pair if usigma(m) = usigma(n) = m + n - 1, where usigma(n) is the sum of unitary divisors of n (A034460).

The larger members are in A306873.

			6224890 is in the sequence since it is the lesser of the amicable pair (6224890, 7336455): usigma(6224890) = usigma(7336455) = 13561344 = 6224890 + 7336455 - 1.

Amiram Eldar, Table of n, a(n) for n = 1..27 (terms below 10^12, from David Moews's site).
David Moews, Perfect, amicable and sociable numbers
J. O. M. Pedersen, Tables of Aliquot Cycles

Cf. A002952, A007992, A034460, A126174, A306873.

us[n_] := Times @@ (1 + Power @@@ FactorInteger[n]) - n;  s={}; Do[m = us[n] + 1; If[m > n && us[m] == n - 1, AppendTo[s, n]], {n, 1, 10^9}]; s

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

Showing 1-7 of 7 results.

cp's OEIS Frontend

A126173 Larger element of a reduced infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A126175 Larger member of an augmented infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A126172 Smaller element of a reduced infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A126176 Number of augmented infinitary amicable pairs (i,j) with i

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A306867 Lesser of augmented bi-unitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A306872 Lesser of augmented unitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula