OFFSET
1,1
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
MAPLE
a:=proc(n)local p: p:=ithprime(n): if isprime(p+2)=true and isprime(p*(p+2)+18)=true and isprime(p*(p+2)+20)=true then p else end if end proc: seq(a(n), n= 1..40000); # Emeric Deutsch, Jul 28 2007
MATHEMATICA
Select[Prime[Range[60000]], PrimeQ[#+2] && PrimeQ[#*(#+2)+18] && PrimeQ[#*(#+2)+20] &] (* G. C. Greubel, Mar 03 2019 *)
PROG
(PARI) {isok(n) = isprime(n) && isprime(n+2) && isprime(n*(n+2)+18) && isprime(n*(n+2)+20)};
forprime(n=1, 500000, if(isok(n), print1(n", "))) \\ G. C. Greubel, Mar 03 2019
(Magma) [n: n in [1..500000] | IsPrime(n) and IsPrime(n+2) and IsPrime(n*(n+2)+18) and IsPrime(n*(n+2)+20)]; // G. C. Greubel, Mar 03 2019
(Sage) [n for n in (1..500000) if is_prime(n) and is_prime(n+2) and is_prime(n*(n+2)+18) and is_prime(n*(n+2)+20)] # G. C. Greubel, Mar 03 2019
(GAP) Filtered([1..500000], k-> IsPrime(k) and IsPrime(k+2) and IsPrime(k*(k+2)+18) and IsPrime(k*(k+2)+20)) # G. C. Greubel, Mar 03 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Ray G. Opao, Jul 06 2007
EXTENSIONS
More terms from Emeric Deutsch, Jul 28 2007
STATUS
approved