OFFSET
1,1
COMMENTS
All terms are multiples of 6. - Jon E. Schoenfield, Apr 13 2024
LINKS
Zak Seidov, Table of n, a(n) for n = 1..10000 (first 1000 terms from Harry J. Smith)
EXAMPLE
12 is in the sequence because 144 +/- 37 = 181 and 107 which are both primes.
k=30 is a term: 30^2 = 900, prime(30) = 113, 900+113 = 1013 and 900-113 = 787, both primes.
MATHEMATICA
Select[ Range[10^4], PrimeQ[ #^2 + Prime[ # ]] && PrimeQ[ #^2 - Prime[ # ]] &]
PROG
(PARI) for(n=1, 20000, if(isprime(n^2+prime(n)) && isprime(n^2-prime(n)), print1(n, " ")))
(PARI) { n=0; default(primelimit, 6100000); for (m=1, 10^9, if (isprime(m^2 + prime(m)) && isprime(m^2 - prime(m)), write("b064483.txt", n++, " ", m); if (n==1000, break)) ) } \\ Harry J. Smith, Sep 16 2009
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Robert G. Wilson v and Jason Earls, Oct 05 2001
STATUS
approved