OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= n -> isprime(n) and not isprime(6*n-1) and not isprime(6*n+1):
select(filter, [seq(i, i=3..1000, 2)]); # Robert Israel, Jan 05 2020
MATHEMATICA
Select[Prime@Range@150, ! PrimeQ[6 # - 1] && ! PrimeQ[6 # + 1] &] (* Robert G. Wilson v, Apr 27 2017 *)
Select[Prime[Range[150]], NoneTrue[6#+{1, -1}, PrimeQ]&] (* Harvey P. Dale, Jun 10 2022 *)
PROG
(PARI) {forprime(n=3, 1000, if(!isprime(6*n-1)&&!isprime(6*n+1), print1(n", ")))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Dimitris Valianatos, Apr 26 2017
STATUS
approved