OFFSET
1,1
COMMENTS
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]; seq[kmax_] := Module[{v = {}}, Do[If[s[k] > 2, AppendTo[v, k]], {k, 1, kmax, 2}]; v]; seq[2*10^10]
PROG
(PARI) is(k) = if(!(k % 2), 0, my(f = factor(k)); 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)) > 2);
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 19 2024
STATUS
approved