A065059 Primes p such that prime(p) + pi(p) is a prime.

3, 13, 53, 71, 113, 181, 397, 463, 479, 557, 593, 647, 683, 701, 719, 743, 769, 863, 941, 983, 1069, 1091, 1109, 1181, 1193, 1231, 1277, 1291, 1307, 1451, 1733, 1811, 1931, 2347, 2393, 2411, 2473, 2531, 2551, 2593, 2887, 3041, 3221, 3251, 3347, 3361

Offset: 1

Author

Robert G. Wilson v, Nov 06 2001

Keywords

• easy
• nonn

Links

• Harry J. Smith, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A065042.

Programs

Mathematica

Select[Prime[Range[500]],PrimeQ[Prime[#]+PrimePi[#]]&] (* Harvey P. Dale, Oct 12 2014 *)

PARI

{ n=0; default(primelimit, 4294965247); for (m=1, 10^9, p=prime(m); if (isprime(prime(p) + primepi(p)), write("b065059.txt", n++, " ", p); if (n==1000, return)) ) } \\ Harry J. Smith, Oct 05 2009