1,1

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.

Pedersen J. M., Knowm amicable pairs.

The values of m for which isigma(m)=isigma(n)=m+n, and m<n

a(5)=5940 because the fifth infinitary amicable pair is (5940,8460), and 5940 is its smallest member

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[First[data4[[k]]], {k, 1, Length[data4]}]

Cf. A049417, A126168, A037445.

Sequence in context: A095619 A043403 A122279 this_sequence A002952 A108344 A112485

Adjacent sequences: A126166 A126167 A126168 this_sequence A126170 A126171 A126172

hard,nonn

Ant King (mathstutoring(AT)ntlworld.com), Dec 21 2006