A049417 -id:A049417 - cp's OEIS frontend

A126171 Number of infinitary amicable pairs (i,j) with i

0, 0, 2, 6, 22, 62, 189, 444, 1116, 2594, 6051, 14141

Offset: 1

Ant King, Dec 22 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(6)=62 because there are 62 infinitary amicable pairs (m,n) with m<n and m<=10^6

Pedersen J. M., Known amicable pairs.

Cf. A126169, A049417, A126168, A037445, A126170.

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; InfinitaryAmicableNumberQ[k_] := If[Nest[properinfinitarydivisorsum, k, 2] == k && ! properinfinitarydivisorsum[k] == k, True, False]; data1 = Select[ Range[10^6], InfinitaryAmicableNumberQ[ # ] &]; data2 = properinfinitarydivisorsum[ # ] & /@ data1; data3 = Table[{data1[[k]], data2[[k]]}, {k, 1, Length[data1]}]; data4 = Select[data3, First[ # ] < Last[ # ] &]; Table[Length[Select[data4, First[ # ] < 10^k &]], {k, 1, 6}]

Infinitary amicable pairs (m,n) satisfy isigma(m)=isigma(n)=m+n, with m

A126173 Larger element of a reduced infinitary amicable pair.

Original entry on oeis.org

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

A126174 Smaller member of an augmented infinitary amicable pair.

Original entry on oeis.org

1252216, 1754536, 2166136, 2362360, 6224890, 7626136, 7851256, 9581320, 12480160, 12494856, 13324311, 15218560, 15422536, 19028296, 29180466, 36716680, 37542190, 40682824, 45131416, 45495352, 56523810, 67195305, 71570296, 80524665, 89740456, 93182440, 101304490

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)=2166136 because 2166136 is the smaller element 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, 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; 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[First[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, n]], {n, 2, 10^9}]; s (* Amiram Eldar, Jan 20 2019 *)

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

a(9)-a(27) from Amiram Eldar, Jan 20 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

A306983 Infinitary pseudoperfect numbers: numbers n equal to the sum of a subset of their proper infinitary divisors.

Original entry on oeis.org

6, 24, 30, 40, 42, 54, 56, 60, 66, 72, 78, 88, 90, 96, 102, 104, 114, 120, 138, 150, 168, 174, 186, 210, 216, 222, 246, 258, 264, 270, 280, 282, 294, 312, 318, 330, 354, 360, 366, 378, 384, 390, 402, 408, 420, 426, 438, 440, 456, 462, 474, 480, 486, 498, 504

Offset: 1

Amiram Eldar, Mar 18 2019

nonn

Subsequence of A005835.

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

Cf. A005835, A007357, A049417, A126168, A129656, A293188, A292985, A318100, A306984.

idivs[x_] := If[x == 1, 1, Sort@Flatten@Outer[Times, Sequence @@ (FactorInteger[x] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]]; s = {}; Do[d = Most[idivs[n]]; c = SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n]; If[c > 0, AppendTo[s, n]], {n, 2, 1000}]; s

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

A323344 Numbers k whose infinitary divisors have an even sum which is larger than 2k, but they cannot be partitioned into two disjoint parts whose sums are equal.

Original entry on oeis.org

2394, 7544, 10184, 1452330, 2154584, 5021912, 5747994, 5771934, 5786298, 5800662, 5834178, 5843754, 5858118, 5886846, 5905998, 5920362, 5929938, 5992182, 6035274, 6059214, 6078366, 6087942, 6102306, 6107094, 6121458, 6174126, 6202854, 6207642, 6245946, 6265098

Offset: 1

Amiram Eldar, Jan 11 2019

nonn

The infinitary version of A171641.

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

Cf. A049417, A077609, A129656, A171641, A323341, A323342, A323343.

infdivs[x_] := If[x == 1, 1, Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[x] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]] ; 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}]]; isigma[n_] := If[n == 1, 1, Times @@ (fun @@@ FactorInteger[n])]; seq={}; Do[s=isigma[n]; If[OddQ[s] || s<=2n, Continue[]]; div = infdivs[n]; If[Coefficient[Times @@ (1 + x^div) // Expand, x, s/2] == 0, AppendTo[seq, n]], {n, 1, 100000}]; seq (* after Michael De Vlieger at A077609 *)

A324277 Infinitary untouchable numbers: numbers that are not the sum of aliquot infinitary divisors of any number.

Original entry on oeis.org

2, 3, 4, 5, 38, 68, 80, 128, 158, 164, 188, 192, 206, 212, 224, 278, 290, 308, 326, 368, 380, 398, 416, 432, 458, 518, 530, 536, 542, 548, 578, 584, 600, 626, 632, 692, 702, 710, 752, 758, 770, 782, 788, 818, 822, 836, 852, 872, 896, 902, 926, 938, 968, 998

Offset: 1

Amiram Eldar, Feb 20 2019

nonn

Amiram Eldar, Table of n, a(n) for n = 1..2256 (terms below 30000)

Cf. A049417, A126168, A005114, A063948 (unitary), A324276 (bi-unitary), A324278 (exponential).

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}]]; isigma[1] = 1; isigma[n_] := isigma[n] = Times @@ (fun @@@ FactorInteger[n]); untouchableQ[n_] := Catch[ Do[ If[n == isigma[k]-k, Throw[True]], {k, 0, (n-1)^2}]] === Null; Reap[ Table[ If[ untouchableQ[n], Sow[n]], {n, 2, 1000}]][[2, 1]] (* after Jean-François Alcover at A005114 *)

A324706 The sum of the tri-unitary divisors of n.

Original entry on oeis.org

1, 3, 4, 5, 6, 12, 8, 15, 10, 18, 12, 20, 14, 24, 24, 17, 18, 30, 20, 30, 32, 36, 24, 60, 26, 42, 40, 40, 30, 72, 32, 33, 48, 54, 48, 50, 38, 60, 56, 90, 42, 96, 44, 60, 60, 72, 48, 68, 50, 78, 72, 70, 54, 120, 72, 120, 80, 90, 60, 120, 62, 96, 80, 85, 84, 144

Offset: 1

Amiram Eldar, Mar 11 2019

A divisor d of n is tri-unitary if the greatest common bi-unitary divisor of d and n/d is 1.

Antti Karttunen, Table of n, a(n) for n = 1..20000
Graeme L. Cohen, On an integer's infinitary divisors, Mathematics of Computation, Vol. 54, No. 189 (1990), pp. 395-411.
Pentti Haukkanen, On the k-ary convolution of arithmetical functions, The Fibonacci Quarterly, Vol. 38, No. 5 (2000) pp. 440-445.

Cf. A000203, A034448, A049417, A072691, A188999, A322485, A324707, A324708, A324709.

f[p_, e_] := If[e == 3, (p^4-1)/(p-1), If[e==6, (p^8-1)/(p^2-1), p^e+1]]; a[1]=1; a[n_]:= Times @@ f @@@ FactorInteger[n]; Array[a, 100]

A324706(n) = { my(f = factor(n)); prod(i=1, #f~, if(3==f[i,2], sigma(f[i,1]^f[i,2]), if(6==f[i,2], ((f[i,1]^8)-1)/((f[i,1]^2)-1), 1+(f[i,1]^f[i,2])))); }; \\ Antti Karttunen, Mar 12 2019

Multiplicative with a(p^3) = 1 + p + p^2 + p^3, a(p^6) = 1 + p^2 + p^4 + p^6, and a(p^e) = 1 + p^e otherwise.

Sum_{k=1..n} a(k) ~ c * n^2, where c = (Pi^2/12) * Product_{p prime} (1 - 1/p^3 + 1/p^4 - 2/p^6 + 2/p^8 - 1/p^9 - 1/p^12 + 1/p^13) = 0.72189237802... . - Amiram Eldar, Nov 24 2022

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

Previous Showing 41-50 of 89 results. Next

cp's OEIS Frontend

A126171 Number of infinitary amicable pairs (i,j) with i

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A126173 Larger element of a reduced infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A126174 Smaller member of an augmented 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

A306983 Infinitary pseudoperfect numbers: numbers n equal to the sum of a subset of their proper infinitary divisors.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A126172 Smaller element of a reduced infinitary amicable pair.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A323344 Numbers k whose infinitary divisors have an even sum which is larger than 2k, but they cannot be partitioned into two disjoint parts whose sums are equal.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A324277 Infinitary untouchable numbers: numbers that are not the sum of aliquot infinitary divisors of any number.

Views