OFFSET
1,1
FORMULA
MATHEMATICA
Select[Prime[Range[150]], Mod[#, 6] == 5 && \[Not] PrimeQ[2 # + 1] &]
PROG
(PARI) forprime(p=1, 900, if(Mod(p, 6)==5 && !ispseudoprime(2*p+1), print1(p, ", "))) \\ Felix Fröhlich, Aug 08 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
Hilko Koning, Aug 08 2016
STATUS
approved