OFFSET
1,1
COMMENTS
All terms are == 0 (mod 6).
LINKS
Zak Seidov, Table of n, a(n) for n = 1..1278
EXAMPLE
n = 948, p = 898691 = prime(71194), q = 898717 = prime(71195);
n = 1134, p = 1285943 = prime(99033), q = 1285969 = prime(99034).
MAPLE
with(numtheory): A248785:=n->`if`(isprime(n^2-13) and isprime(n^2+13) and pi(n^2+13) = pi(n^2-13)+1, n, NULL): seq(A248785(n), n=1..2*10^4); # Wesley Ivan Hurt, Oct 13 2014
MATHEMATICA
Select[Range[17000], PrimeQ[#^2-13]&&NextPrime[#^2-13]==#^2+13&] (* Harvey P. Dale, Aug 14 2020 *)
PROG
(PARI) isok(n) = isprime(p=n^2-13) && isprime(q=n^2+13) && (q==nextprime(p+1)); \\ Michel Marcus, Oct 14 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Oct 13 2014
EXTENSIONS
More terms from Michel Marcus, Oct 14 2014
STATUS
approved