OFFSET
1,1
COMMENTS
Primes p such that 4*p+3, 16*p+15 and 64*p+63 are also primes. - Vincenzo Librandi, Aug 04 2010
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
MAPLE
f:=proc(x) options operator, arrow: 4*x+3 end proc: a:=proc(n) if isprime(n)= true and isprime(f(n))=true and isprime(f(f(n)))=true and isprime(f(f(f(n)))) =true then n else end if end proc: seq(a(n), n=1..20000); # Emeric Deutsch, Jan 01 2008
MATHEMATICA
Select[Prime@ Range@ 2100, Times @@ Boole@ PrimeQ@ Rest@ NestList[4 # + 3 &, #, 3] > 0 &] (* Michael De Vlieger, Sep 19 2016 *)
PROG
(Magma) [n: n in [1..150000] | IsPrime(n) and IsPrime(4*n+3) and IsPrime(16*n+15) and IsPrime(64*n+63)] // Vincenzo Librandi, Aug 04 2010
(PARI) is(n)=isprime(n) && isprime(4*n+3) && isprime(16*n+15) && isprime(64*n+63) \\ Charles R Greathouse IV, Sep 20 2016
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved