%I #10 Feb 09 2017 08:27:04
%S 2,3,7,19,79,101,113,163,173,199,223,239,251,263,311,349,443,463,577,
%T 593,641,659,743,839,881,1033,1097,1109,1151,1373,1399,1429,1439,1459,
%U 1627,1693,1831,1871,2029,2069,2137,2237,2287,2423,2473,2617,2687,2713
%N Primes p such that prime(p) - pi(p) is a prime.
%H Harry J. Smith, <a href="/A065060/b065060.txt">Table of n, a(n) for n = 1..1000</a>
%p with(numtheory): A065060:=n->`if`(isprime(n) and isprime(ithprime(n)-pi(n)), n, NULL): seq(A065060(n), n=1..10^4); # _Wesley Ivan Hurt_, Feb 09 2017
%t Prime[ Select[ Range[500], PrimeQ[ Prime[ Prime[ # ]] - PrimePi[ Prime[ # ]]] & ]]
%o (PARI) { n=0; default(primelimit, 4294965247); for (m=1, 10^9, p=prime(m); if (isprime(prime(p) - primepi(p)), write("b065060.txt", n++, " ", p); if (n==1000, return)) ) } \\ _Harry J. Smith_, Oct 05 2009
%Y Cf. A000040 (primes), A000720 (pi), A065046 (prime(n) - pi(n) is a prime).
%K easy,nonn
%O 1,1
%A _Robert G. Wilson v_, Nov 06 2001