OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
840 is in the sequence since isigma(840) = 2880 > 3 * 840.
MATHEMATICA
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}]]]]]; properinfinitarydivisorsum[k_] := Plus @@ InfinitaryDivisors[k] - k; Infinitary3AbundantNumberQ[k_] := If[properinfinitarydivisorsum[k] >= 2 k, True, False]; Select[Range[15000], Infinitary3AbundantNumberQ[#] &] (* after Ant King at A129656 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 10 2018
STATUS
approved