OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
MATHEMATICA
f[p_, e_] := ((e+1)*p^2 - (e+2)*p + p^(-e))/((e+1)*(p-1)^2); s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[4*10^7], s[#] > 3 &]
PROG
(PARI) is(n) = {my(f = factor(n)); prod(i = 1, #f~, p=f[i, 1]; e=f[i, 2]; (-2*p - e*p + p^2 + e*p^2 + p^(-e))/((e + 1)*(p - 1)^2)) > 3; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 19 2024
STATUS
approved