OFFSET
1,1
EXAMPLE
130001, 130003, 130007 and 130009 are all prime. Thus 13 is a member of this sequence.
PROG
(Python)
import sympy
from sympy import isprime
from sympy import prime
{print(prime(n)) for n in range(1, 10**5) if isprime(1000*prime(n)+1) and isprime(1000*prime(n)+3) and isprime(1000*prime(n)+7) and isprime(1000*prime(n)+9)}
(PARI) for(n=1, 10^5, s=prime(n); if(ispseudoprime(1000*s+1) && ispseudoprime(1000*s+3) && ispseudoprime(1000*s+7) && ispseudoprime(1000*s+9), print(s)));
CROSSREFS
KEYWORD
nonn
AUTHOR
Derek Orr, May 17 2014
STATUS
approved