OFFSET
1,1
COMMENTS
Primes p such that p*pi(p)+1 is prime (see Crossrefs).
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
EXAMPLE
13 is in the list because, being the 6th prime, 13*6+1=79 is prime.
MATHEMATICA
Select[Prime[Range[700]], PrimeQ[# PrimePi[#]+1]&] (* Harvey P. Dale, May 28 2012 *)
PROG
(Magma) [ NthPrime(n): n in [1..650] | IsPrime(NthPrime(n)*n+1) ]; // Bruno Berselli, Jul 05 2011
(PARI) k=0; forprime(p=2, 1e4, if(isprime(k++*p+1), print1(p", "))) \\ Charles R Greathouse IV, Jul 05 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrea Raffetti, Jul 05 2011
STATUS
approved