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
Links
- Jens Kruse Andersen, Table of n, a(n) for n = 1..10000
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
Comments