OFFSET
1,1
COMMENTS
LINKS
Harvey P. Dale and K. D. Bajpai, Table of n, a(n) for n = 1..10000 (first 1000 terms from Harvey P. Dale)
EXAMPLE
For p=2, 2*2^3 -+ 3 = (13, 19), both prime, so 2 is a term of the sequence.
For p=13, 2*13^3 -+ 3 = (4391, 4397), both prime, so 13 is a term of the sequence.
MAPLE
select(p -> andmap(isprime, [p, 2*p^3+3, 2*p^3-3]), [seq(p, p=1.. 10^5)]); # K. D. Bajpai, Nov 28 2019
MATHEMATICA
Select[Prime[Range[5000]], And@@PrimeQ[2 #^3+{3, -3}]&] (* Harvey P. Dale, Jan 25 2013 *)
PROG
(Magma) [p: p in PrimesUpTo(100000)|IsPrime(2*p^3-3) and IsPrime(2*p^3+3)]
(PARI) forprime(p=1, 55000, if(ispseudoprime(2*p^3-3) && ispseudoprime(2*p^3+3), print1(p, ", "))) \\ Felix Fröhlich, Nov 27 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Vincenzo Librandi, Mar 17 2010
STATUS
approved