OFFSET
1,1
COMMENTS
Also called nonomatic primes. There is probably an infinity of them. There seems to be no prime number with a similar property using 5 or a larger factor in the polynomials.
LINKS
Vincenzo Librandi and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (first 1000 terms from Librandi)
EXAMPLE
a(2) = 29 since (3*29)^2 + 29^2 + 3^2 = 8419 and (3*29)^2 - 29^2 - 3^2 = 6719 are both primes.
MATHEMATICA
Select[Prime@Range[1, 2000], PrimeQ[9 #^2 + #^2 + 9] && PrimeQ[9 #^2 - #^2 - 9] &] (* Vincenzo Librandi, Oct 18 2012 *)
PROG
(Magma) [p: p in PrimesUpTo(15000)| IsPrime( (3*p)^2 + p^2 + 3^2) and IsPrime((3*p)^2 - p^2 - 3^2)]; // Vincenzo Librandi, Oct 18 2012
(PARI) is(n)=isprime(10*n^2+9) && isprime(8*n^2-9) && isprime(n) \\ Charles R Greathouse IV, Jun 10 2015
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Olivier Gérard, Feb 19 2003
STATUS
approved