login
Primes p such that 13*p + 2 is also prime.
1

%I #16 Apr 03 2023 10:36:11

%S 3,5,17,29,47,53,59,113,137,173,227,239,257,269,347,353,389,419,449,

%T 467,479,509,557,563,599,677,719,743,827,839,953,977,983,1013,1019,

%U 1103,1109,1163,1193,1217,1223,1259,1277,1289,1307,1373,1427,1487

%N Primes p such that 13*p + 2 is also prime.

%C Except for 3, all terms == 5 (mod 6). - _Robert Israel_, Sep 06 2018

%H Robert Israel, <a href="/A113169/b113169.txt">Table of n, a(n) for n = 1..10000</a>

%H Chris Caldwell <a href="https://t5k.org/lists/small/1000.txt">The First 1000 Primes</a>.

%e If p=173 then 13*p + 2 = 2251 (prime).

%p select(t -> isprime(t) and isprime(13*t+2), [3,seq(i,i=5..10000,6)]); # _Robert Israel_, Sep 06 2018

%o (Magma) [p: p in PrimesUpTo(3000)| IsPrime(13*p+2)] // _Vincenzo Librandi_, Jan 29 2011

%o (PARI) forprime(p=3, 1e4, if(isprime(13*p+2), print1(p, ", "))); \\ _Altug Alkan_, Sep 07 2018

%Y Cf. A113151, A113115.

%K nonn

%O 1,1

%A _Parthasarathy Nambi_, Jan 06 2006