OFFSET
1,1
COMMENTS
LINKS
Iain Fox, Table of n, a(n) for n = 1..10000 (first 4356 terms from Robert Israel)
MAPLE
select(p -> isprime(3*p+8) and isprime(3*p+2) and isprime(p+2) and isprime(p), [3, seq(i, i=5..10^6, 6)]); # Robert Israel, Nov 23 2016
MATHEMATICA
Select[Prime[Range[10000]], Union[PrimeQ/@{# + 2, 3 # + 2, 3 # + 8}] == {True}&]
PROG
(MATLAB)
P = primes(10^6);
P1 = intersect(P, P-2);
P1 = intersect(P1, (P-2)/3);
P1 = intersect(P1, (P-8)/3) % Robert Israel, Nov 23 2016
(Magma) [p: p in PrimesUpTo(100000) | IsPrime(p+2) and IsPrime(3*p+2) and IsPrime(3*p+8)]; // Vincenzo Librandi, Nov 23 2016
(PARI) isok(p) = isprime(p) && isprime(p+2) && isprime(3*p+2) && isprime(3*p+8); \\ Michel Marcus, Dec 17 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Ivan N. Ianakiev, Nov 23 2016
STATUS
approved