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

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

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

Prog

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

Crossrefs

Cf. A033210, A142508, A171517, A248221.

Sequence in context: A347988 A275081 A097443 * A171517 A179026 A179027

Adjacent sequences: A248365 A248366 A248367 * A248369 A248370 A248371

Keyword

nonn,easy

Author

Zak Seidov, Oct 05 2014

Status

approved