login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A247949 Primes p such that the polynomial k^4 + k^3 + k^2 + k + p yields only primes for k = 0...5. 3

%I #34 Dec 24 2023 16:07:50

%S 7,43,79,457,877,967,1093,2437,2683,3187,5077,5923,7933,8233,11923,

%T 12889,15787,17389,19993,31543,41113,41617,42457,71359,77863,80683,

%U 91393,101719,102643,105967,107347,120163,129733,137593,151783,170263,175723,197569,210127

%N Primes p such that the polynomial k^4 + k^3 + k^2 + k + p yields only primes for k = 0...5.

%C All terms == 1 mod 6. - _Robert Israel_, Jan 11 2015

%H K. D. Bajpai, <a href="/A247949/b247949.txt">Table of n, a(n) for n = 1..6730</a>

%e a(1) = 7:

%e 0^4 + 0^3 + 0^2 + 0 + 7 = 7;

%e 1^4 + 1^3 + 1^2 + 1 + 7 = 11;

%e 2^4 + 2^3 + 2^2 + 2 + 7 = 37;

%e 3^4 + 3^3 + 3^2 + 3 + 7 = 127;

%e 4^4 + 4^3 + 4^2 + 4 + 7 = 347;

%e 5^4 + 5^3 + 5^2 + 5 + 7 = 787;

%e all six are primes.

%p select(p -> andmap(isprime, [p, p+4, p+30, p+120, p+340, p+780]), [seq(6*i+1, i=1..10^5)]); # _Robert Israel_, Jan 11 2015

%t Select[f=k^4 + k^3 + k^2 + k; k = {0, 1, 2, 3, 4, 5}; Prime[Range[2000000]], And @@ PrimeQ[#+f] &]

%t Select[Prime[Range[20000]],AllTrue[#+{4,30,120,340,780},PrimeQ]&] (* _Harvey P. Dale_, Dec 24 2023 *)

%o (PARI) forprime(p=1, 500000, if( isprime(p+0)& isprime(p+4)& isprime(p+30)& isprime(p+120)& isprime(p+340)& isprime(p+780), print1(p,", ")))

%Y Cf. A144051, A187057, A187058, A187060, A190800, A191456, A191457, A191458.

%K nonn

%O 1,1

%A _K. D. Bajpai_, Jan 11 2015

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 08:08 EDT 2024. Contains 371265 sequences. (Running on oeis4.)