|
%I #7 Sep 08 2022 08:45:28
%S 1,2,34,43,64,76,80,160,194,276,416,620,625,832,996,1022,1251,1341,
%T 1345,1356,1435,1461,1475,1672,1694,1700,1733,1769,1777,1782,1796,
%U 1808,1997,2144,2239,2245,2260,2270,2273,2411,2573,2659,2679,2829,3396,3497,3516
%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 + n^45 + n^47 + n^49 + n^51 + n^53 + n^55 + n^57 + n^59 is prime.
%H Vincenzo Librandi, <a href="/A124208/b124208.txt">Table of n, a(n) for n = 1..470</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 + n^45 + n^47 + n^49 + n^51 + n^53 + n^55 + n^57 + n^59],Print[n]],{n,1,2400}]
%t Select[Range[5000], PrimeQ[Total[#^Range[1, 59, 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..59 by 2]]]; // _Vincenzo Librandi_, Jun 28 2014
%Y Cf. A049407, similar sequences listed in A244376.
%K nonn,easy
%O 1,2
%A _Artur Jasinski_, Dec 13 2006
|
