A061067 m-th prime prime(m) is included iff prime(m) + m is also prime.

2, 3, 7, 13, 61, 79, 89, 101, 131, 139, 181, 223, 317, 349, 359, 421, 479, 503, 521, 541, 557, 569, 577, 593, 619, 641, 701, 733, 827, 863, 881, 929, 1021, 1069, 1223, 1231, 1249, 1277, 1361, 1399, 1439, 1487, 1733, 1831, 1889, 1987, 2069, 2137, 2143, 2203

Offset: 0

Author

Labos Elemer, May 28 2001

Keywords

• nonn

Examples

			5th term here is 61 = prime(18) and 61 + 18 = 79.
		

Links

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

Crossrefs

Cf. A064402, A061068, A227420.

Programs

Magma

[NthPrime(n): n in [1..250] | IsPrime(NthPrime(n)+ n)]; // Vincenzo Librandi, Jan 19 2015

Mathematica

Prime[Select[Range[500], PrimeQ[Prime[ # ] + # ] &]] (* Stefan Steinerberger, Jul 21 2006 *)
Select[Prime[Range[400]],PrimeQ[#+PrimePi[#]]&] (* Harvey P. Dale, Oct 03 2016 *)

PARI

{ n=-1; m=0; forprime (p=2, 109597, if (isprime(p + m++), write("b061067.txt", n++, " ", p)) ) } \\ Harry J. Smith, Jul 17 2009

Formula

a(n) + A064402(n+1) = A061068(n+1). [corrected by Martin Fuller]

Extensions

Edited by N. J. A. Sloane at the suggestion of Andrew S. Plewe, Aug 23 2007