A248368 Primes p such that 52*p + 1 is prime.

3, 13, 31, 43, 73, 109, 139, 151, 181, 193, 211, 223, 229, 283, 349, 379, 409, 421, 463, 523, 601, 619, 691, 769, 823, 853, 1021, 1033, 1069, 1153, 1231, 1279, 1303, 1453, 1459, 1471, 1531, 1663, 1693, 1723, 1741, 1783, 1831, 1873, 1933, 2029, 2131, 2251, 2269, 2293, 2593, 2671, 2749, 2791

Offset: 1

Author

Zak Seidov, Oct 05 2014

Keywords

• nonn
• easy

Comments

Or, primes in A248221. Subsequence of A248221. Note that a(1..6) coincide with A171517(1..6).

Links

• Jens Kruse Andersen, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A033210, A142508, A171517, A248221.

Programs

Maple

A248368:=n->`if`(isprime(52*n+1) and isprime(n), n, NULL): seq(A248368(n), n=1..4000); # Wesley Ivan Hurt, Oct 05 2014

Mathematica

s = {}; Do[If[PrimeQ[1 + 52*(p = Prime[n])], AppendTo[s, p]], {n, 500}]; s
Select[Prime[Range[500]],PrimeQ[52#+1]&] (* Harvey P. Dale, Aug 15 2017 *)

PARI

forprime(p=1,10^4,if(isprime(52*p+1),print1(p,", "))) \\ Derek Orr, Oct 05 2014