A247276 Primes p such that p + m^2 is prime for all m in {2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 26}.
163, 409333, 13491637509487, 19802478368863
Offset: 1
Programs
-
Mathematica
Select[Prime[Range[35000]],AllTrue[#+{4,16,36,64,100,144,196,256,324,400,484,576,676},PrimeQ]&] (* The program generates the first two terms of the sequence. To generate a(3) and a(4), increase the Range constant to 67*10^10 but the program will take a very long time to run. *) (* Harvey P. Dale, Mar 05 2025 *) -
PARI
forprime(p=1, 10^12, c=0; for(i=1, 13, if(ispseudoprime(p+(2*i)^2), c++);if(!ispseudoprime(p+(2*i)^2),break)); if(c==13, print1(p, ", "))) \\ Derek Orr, Sep 11 2014
Comments