OFFSET
1,2
EXAMPLE
a(5) = 7 is a term because 2*7 = 14 = 3+11 = 7+7, with 3*11 == 5 (mod 14) and 7*7 == 7 (mod 14), and both 5 and 7 are prime.
5 is not a term because 2*5 = 10 = 3+7 = 5+5, but 3*7 == 1 (mod 10) and 1 is not prime.
MAPLE
filter:= proc(k) local p;
p:= 1;
while p <= k do
p:= nextprime(p);
if isprime(2*k-p) and not isprime(-p^2 mod 2*k) then return false fi
od;
true
end proc:
select(filter, [$1..1000]);
MATHEMATICA
q[k_] := AllTrue[Select[Range[2, 2*k], PrimeQ], ! PrimeQ[2*k - #] || PrimeQ[Mod[#*(2*k - #), 2*k]] &]; Select[Range[100], q] (* Amiram Eldar, Dec 28 2021 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J. M. Bergot and Robert Israel, Dec 28 2021
STATUS
approved