OFFSET
1,1
COMMENTS
Primes prime(n) such that A031215(n) - n is prime.
EXAMPLE
3 is in this sequence because 3 is prime(2) and prime(2*2) - 2 = prime(4)- 2 = 7 - 2 = 5 is prime.
MAPLE
seq(`if`(isprime(ithprime(2*n)-n), ithprime(n), NULL), n=1..1000); # Robert Israel, Apr 13 2015
MATHEMATICA
Prime[Select[Range[555], PrimeQ[Prime[2*#]-#]&]] (* Ivan N. Ianakiev, Apr 14 2015 *)
PROG
(Magma) [NthPrime(n): n in [1..550] | IsPrime(NthPrime(2*n)-n)];
(PARI) for(n=1, 10^3, if(isprime(prime(2*n)-n), print1(prime(n), ", "))) \\ Derek Orr, Apr 14 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Juri-Stepan Gerasimov, Apr 13 2015
STATUS
approved