OFFSET
1,1
COMMENTS
All these primes are congruent to 1 mod 12.
LINKS
Edward Jiang, Table of n, a(n) for n = 1..10000
MAPLE
select(t -> isprime(t) and isprime((13+t)/2), [seq(12*k+1, k=1..100)]); # Robert Israel, Aug 05 2014
MATHEMATICA
aa = {}; k = 13; Do[If[PrimeQ[(k + Prime[n])/2], AppendTo[aa, Prime[n]]], {n, 1, 500}]; aa
PROG
(PARI)
forprime(p=3, 10^4, if(isprime((13+p)/2), print1(p, ", "))) \\ Derek Orr, Aug 05 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 11 2008
STATUS
approved