OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..755 (terms below 10^11)
MATHEMATICA
q[n_] := Module[{s = DivisorSigma[1, n/2^IntegerExponent[n, 2]] - If[OddQ[n], n, 0]}, Divisible[n, s] && OddQ[n/s]]; Select[Range[550000], # == 4 || (CompositeQ[#] && q[#]) &]
PROG
(PARI) isok(k) = if(k == 1 || isprime(k), 0, if(k == 4, 1, my(s = sigma(k >> valuation(k, 2)) - if(k%2, k)); !(k % s) && (k/s) % 2));
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Apr 27 2025
STATUS
approved
