|
|
A061067
|
|
m-th prime prime(m) is included iff prime(m) + m is also prime.
|
|
15
|
|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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: A273814 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
|
|
|
|