OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := DivisorSum[e + 1, p^(# - 1) &]; mesigma[1] = 1; mesigma[n_] := Times @@ f @@@ FactorInteger[n]; meAbQ[n_] := mesigma[n] > 2*n; Select[Range[1, 10^5, 2], meAbQ]
PROG
(PARI) is(k) = if(!(k%2), 0, my(f=factor(k)); prod(i=1, #f~, sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))) > 2*k);
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Dec 14 2024
STATUS
approved