OFFSET
1,1
COMMENTS
All terms are == {3, 7} mod 10. Subsequence of A246519.
LINKS
Zak Seidov, Table of n, a(n) for n = 1..60
MATHEMATICA
Select[Prime[Range[193*10^5]], AllTrue[#^{1, 2, 3, 5, 7}+4, PrimeQ]&] (* Harvey P. Dale, Sep 07 2024 *)
PROG
(PARI)
forprime(p=1, 10^9, if(ispseudoprime(4+p) && ispseudoprime(4+p^2) && ispseudoprime(4+p^3) && ispseudoprime(4+p^5) && ispseudoprime(4+p^7), print1(p, ", "))) \\ Derek Orr, Aug 30 2014
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Aug 29 2014
STATUS
approved