login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A238185
Primes p such that prime(prime(p^2)) + 2 is also prime.
1
2, 23, 97, 163, 463, 491, 557, 659, 677, 977, 1033, 1151, 1187, 1429, 1439, 1511, 1579, 1663, 1933, 2111, 2113, 2141, 2381, 2969, 3301, 3491, 3803, 3929, 4201, 4421, 4447, 4513, 4547, 4789, 5039, 5273, 5281, 5303, 5309, 5449, 5669, 5741, 5939, 5981, 6053
OFFSET
1,1
LINKS
EXAMPLE
23 is in the sequence because 23 is prime and prime(prime(23^2)) + 2 = 35803 is also prime.
97 is in the sequence because 97 is prime and prime(prime(97^2)) + 2 = 1269643 is also prime.
MAPLE
KD := proc() local a, b, d; a:=ithprime(n); b:=ithprime(ithprime(a^2))+2; if isprime (b) then RETURN (a); fi; end: seq(KD(), n=1..300);
MATHEMATICA
n=0; Do[If[PrimeQ[Prime[Prime[Prime[k]^2]]+2], n=n+1; Print[n, " ", Prime[k]]], {k, 1, 5000}]
Select[Prime[Range[800]], PrimeQ[Prime[Prime[#^2]]+2]&] (* Harvey P. Dale, Dec 19 2014 *)
CROSSREFS
KEYWORD
nonn,less
AUTHOR
K. D. Bajpai, Feb 19 2014
STATUS
approved