OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..1000
MAPLE
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
MATHEMATICA
Prime[ Select[ Range[500], PrimeQ[ Prime[ Prime[ # ]] - PrimePi[ Prime[ # ]]] & ]]
PROG
(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
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Robert G. Wilson v, Nov 06 2001
STATUS
approved