|
%I #20 Sep 08 2022 08:45:28
%S 1,7,10,171,304,322,357,385,418,436,448,592,606,628,642,858,1065,1186,
%T 1438,1543,1639,1843,1848,1879,1899,2110,2241,2517,2544,2578,2668,
%U 3195,3207,3306,3307,3657,4137,4177,4329,4617,4837,4945,4956,5026,5281,5509
%N Numbers n such that 1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + n^13 + n^15 + n^17 + n^19 + n^21 + n^23 + n^25 + n^27 + n^29 + n^31 + n^33 + n^35 + n^37 + n^39 + n^41 + n^43 is prime.
%H Vincenzo Librandi, <a href="/A124200/b124200.txt">Table of n, a(n) for n = 1..620</a>
%t Do[If[PrimeQ[1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + n^13 + n^15 + n^17 + n^19 + n^21 + n^23 + n^25 + n^27 + n^29 + n^31 + n^33 + n^35 + n^37 + n^39 + n^41 + n^43], Print[n]], {n, 1, 2400}]
%t Select[Range[6000], PrimeQ[Total[#^Range[1, 43, 2]] + 1] &] (* _Vincenzo Librandi_, Jun 28 2014 *)
%o (Magma) [n: n in [0..5000] | IsPrime(s) where s is 1+&+[n^i: i in [1..43 by 2]]]; // _Vincenzo Librandi_, Jun 28 2014
%o (PARI) is(n)=isprime(1+n+n^3+n^5+n^7+n^9+n^11+n^13+n^15+n^17+n^19+n^21+n^23 + n^25+n^27+n^29+n^31+n^33+n^35+n^37+n^39+n^41+n^43) \\ _Charles R Greathouse IV_, Feb 20 2017
%Y Cf. A049407, similar sequences listed in A244376.
%K nonn,easy
%O 1,2
%A _Artur Jasinski_, Dec 13 2006
|
