OFFSET
1,1
COMMENTS
The infinitary 3-abundant numbers (analogous to A068403) are the terms that are not infinitary 3-perfect (120, 73440, 44553600, 252927360, 8402296320, 575409744000, ..., which are terms in A007358). - Amiram Eldar, Sep 16 2025
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
f[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_] := Times @@ f @@@ FactorInteger[n]; Select[Range[15000], isigma[#] >= 3*# &] (* Amiram Eldar, Sep 16 2025 *)
PROG
(PARI) isigma(n) = {my(f = factor(n), b); prod(i = 1, #f~, b = binary(f[i, 2]); prod(k = 1, #b, if(b[k], 1+f[i, 1]^(2^(#b-k)), 1))); }
isok(k) = isigma(k) >= 3*k; \\ Amiram Eldar, Sep 16 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 10 2018
EXTENSIONS
Name corrected by Amiram Eldar, Sep 16 2025
STATUS
approved
