OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[p_, e_] := p^(2^(-1 + Position[Reverse@IntegerDigits[e, 2], _?(# == 1 &)])); q[1] = False; q[n_] := Times @@ (1 + 1/(2*Flatten@ (f @@@ FactorInteger[n]))) > 2; Select[Range[300000], q]
PROG
(PARI) is(n) = {my(f = factor(n), b); prod(i = 1, #f~, b = binary(f[i, 2]); prod(k=1, #b, if(b[k], 1 + 1/(2*f[i, 1]^(2^(#b-k))), 1))) > 2; }
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jul 20 2024
STATUS
approved