OFFSET
1,1
COMMENTS
LINKS
Martin Ehrenstein, Table of n, a(n) for n = 1..10000
EXAMPLE
a(4) = 223 is a term because 223, 223+4 = 227, 3*223+4 = 673 and 3*223+8 = 677 are all prime.
MAPLE
select(p -> isprime(p) and isprime(p+4) and isprime(3*p+4) and isprime(3*p+8), [3, seq(i, i=13..10^6, 30)]);
MATHEMATICA
Select[Range[200000], AllTrue[{#, # + 4, 3*# + 4, 3*# + 8}, PrimeQ] &] (* Amiram Eldar, Mar 07 2022 *)
PROG
(Python)
from sympy import sieve, isprime
for p in sieve.primerange(0, 10**6):
if(all(isprime(q) for q in [p+4, 3*p+4, 3*p+8])):
print (p, end=", ") # Martin Ehrenstein, Mar 09 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Mar 07 2022
STATUS
approved