|
%I #11 Sep 08 2022 08:45:28
%S 2,4,8,20,22,32,36,50,54,62,64,72,78,84,86,90,92,94,96,98,112,124,134,
%T 144,146,216,224,238,240,246,250,256,262,276,294,296,298,300,314,334,
%U 370,378,382,392,400,402,406,420,430,450,472,480,482,494,510,512,526
%N Numbers n such that 1 + n + n^3 + n^5 is prime.
%C All numbers n have to be even, because sum of 3 odd + 1 is even and can't be prime >3.
%H Vincenzo Librandi, <a href="/A124154/b124154.txt">Table of n, a(n) for n = 1..1000</a>
%t Do[If[PrimeQ[1 + n + n^3 + n^5], Print[n]], {n, 1, 300}]
%t Select[Range[1000], PrimeQ[Total[#^Range[1, 5, 2]] + 1] &] (* _Vincenzo Librandi_, Jun 27 2014 *)
%o (PARI) forstep(n=2,1000,2,if(isprime(1 + n + n^3 + n^5), print1(n","))) \\ _Franklin T. Adams-Watters_, Apr 09 2009
%o (Magma) [n: n in [0..600] | IsPrime(s) where s is 1+&+[n^i: i in [1..5 by 2]]]; // _Vincenzo Librandi_, Jun 27 2014
%Y Cf. A049407, similar sequences listed in A244376.
%K nonn,easy
%O 1,1
%A _Artur Jasinski_, Dec 13 2006
%E More terms from _Franklin T. Adams-Watters_, Apr 09 2009
|
