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,1

Comments

All terms are == {7, 13} mod 30.

Subsequence of A247275.

Links

Table of n, a(n) for n=1..4.

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 *)

Prog

(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

Crossrefs

Cf. A049492, A092120, A246842, A247269, A145741, A247273, A247275.

Sequence in context: A145741 A247273 A247275 * A127882 A240255 A045006

Adjacent sequences: A247273 A247274 A247275 * A247277 A247278 A247279

Keyword

nonn,hard,more

Author

Zak Seidov, Sep 11 2014

Status

approved