login
Primes p such that 4+p, 4+p^2, 4+p^3, 4+p^5, and 4+p^7 are all prime.
3

%I #9 Sep 07 2024 17:56:25

%S 7,469363,2552713,3378103,6595597,6629683,39837517,46024063,46167307,

%T 97371007,97629403,105528217,136983307,169483033,202953613,213792193,

%U 216520987,216738043,221705647,304033927,317502193,359133553

%N Primes p such that 4+p, 4+p^2, 4+p^3, 4+p^5, and 4+p^7 are all prime.

%C All terms are == {3, 7} mod 10. Subsequence of A246519.

%H Zak Seidov, <a href="/A246562/b246562.txt">Table of n, a(n) for n = 1..60</a>

%t Select[Prime[Range[193*10^5]],AllTrue[#^{1,2,3,5,7}+4,PrimeQ]&] (* _Harvey P. Dale_, Sep 07 2024 *)

%o (PARI)

%o 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

%Y Cf. A007591, A073573, A125260, A172367, A246519.

%K nonn

%O 1,1

%A _Zak Seidov_, Aug 29 2014