%I #38 Sep 08 2022 08:45:03
%S 2,3,7,13,61,79,89,101,131,139,181,223,317,349,359,421,479,503,521,
%T 541,557,569,577,593,619,641,701,733,827,863,881,929,1021,1069,1223,
%U 1231,1249,1277,1361,1399,1439,1487,1733,1831,1889,1987,2069,2137,2143,2203
%N m-th prime prime(m) is included iff prime(m) + m is also prime.
%H Harry J. Smith, <a href="/A061067/b061067.txt">Table of n, a(n) for n = 0..1000</a>
%F a(n) + A064402(n+1) = A061068(n+1). [corrected by _Martin Fuller_]
%e 5th term here is 61 = prime(18) and 61 + 18 = 79.
%t Prime[Select[Range[500], PrimeQ[Prime[ # ] + # ] &]] (* _Stefan Steinerberger_, Jul 21 2006 *)
%t Select[Prime[Range[400]],PrimeQ[#+PrimePi[#]]&] (* _Harvey P. Dale_, Oct 03 2016 *)
%o (PARI) { n=-1; m=0; forprime (p=2, 109597, if (isprime(p + m++), write("b061067.txt", n++, " ", p)) ) } \\ _Harry J. Smith_, Jul 17 2009
%o (Magma) [NthPrime(n): n in [1..250] | IsPrime(NthPrime(n)+ n)]; // _Vincenzo Librandi_, Jan 19 2015
%Y Cf. A064402, A061068, A227420.
%K nonn
%O 0,1
%A _Labos Elemer_, May 28 2001
%E Edited by _N. J. A. Sloane_ at the suggestion of _Andrew S. Plewe_, Aug 23 2007
|