OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
7 is a term because 7*2^1 + 1 = 15 and 7*2^7 + 1 = 897 are composites for divisors 1 and 7 of 7.
MAPLE
filter:= proc(k) andmap(d -> not isprime(k*2^d+1), NumberTheory:-Divisors(k)) end proc:
select(filter, [$1..200]); # Robert Israel, Mar 29 2026
MATHEMATICA
q[k_] := AllTrue[Divisors[k], CompositeQ[k*2^# + 1] &]; Select[Range[161], q] (* Amiram Eldar, Feb 24 2026 *)
PROG
(Magma) [k: k in [1..200] | #[d: d in Divisors(k) | IsPrime(k*2^d+1)] eq 0];
(PARI) isok(k) = fordiv(k, d, if (isprime(k*2^d+1), return(0))); 1; \\ Michel Marcus, Feb 18 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Feb 13 2026
STATUS
approved
