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,1

Links

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

Formula

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

Example

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

Mathematica

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

Prog

(PARI) { n=-1; m=0; forprime (p=2, 109597, if (isprime(p + m++), write("b061067.txt", n++, " ", p)) ) } \\ Harry J. Smith, Jul 17 2009
(Magma) [NthPrime(n): n in [1..250] | IsPrime(NthPrime(n)+ n)]; // Vincenzo Librandi, Jan 19 2015

Crossrefs

Cf. A064402, A061068, A227420.

Sequence in context: A361988 A085872 A070858 * A171434 A104366 A104367

Adjacent sequences: A061064 A061065 A061066 * A061068 A061069 A061070

Keyword

nonn

Author

Labos Elemer, May 28 2001

Extensions

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

Status

approved