OFFSET
1,1
COMMENTS
All terms == 1 (mod 10). - Robert Israel, Mar 18 2020
LINKS
Robert Israel, Table of n, a(n) for n = 1..1000
MAPLE
filter:= t -> andmap(isprime, [t, t^3+6, t^3+12, t^3+18]):
select(filter, [seq(i, i=1..10^7, 10); # Robert Israel, Mar 18 2020
MATHEMATICA
Select[Prime[Range[300000]], AllTrue[#^3+{6, 12, 18}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, May 28 2019 *)
PROG
(Magma) [p: p in PrimesUpTo(1000000)|IsPrime(p^3+6) and IsPrime(p^3+12) and IsPrime(p^3+18)] // Vincenzo Librandi, Dec 13 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Nov 09 2010
STATUS
approved