cp's OEIS Frontend

A124187 Numbers n such that 1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + ... + n^33 + n^35 is prime.

%I A124187 #17 Sep 08 2022 08:45:28
%S A124187 1,7,10,17,52,69,108,161,173,231,306,330,338,352,416,582,584,593,635,
%T A124187 767,834,855,868,892,927,944,950,1044,1060,1203,1242,1299,1302,1509,
%U A124187 1520,1551,1637,1972,2067,2078,2135,2303,2310,2366,2416,2511,2514,2556,2581
%N A124187 Numbers n such that 1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + ... + n^33 + n^35 is prime.
%H A124187 Vincenzo Librandi, <a href="/A124187/b124187.txt">Table of n, a(n) for n = 1..258</a>
%p A124187 a:= proc(n) option remember; local k;
%p A124187       for k from 1+ a(n-1) while
%p A124187         not isprime(1+(k^37-k)/(k^2-1)) do od; k
%p A124187     end: a(1):=1:
%p A124187 seq(a(n), n=1..30);  # _Alois P. Heinz_, Jun 26 2014
%t A124187 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],Print[n]],{n,1,2400}]
%t A124187 Select[Range[5000], PrimeQ[Total[#^Range[1, 35, 2]] + 1] &] (* _Vincenzo Librandi_, Jun 27 2014 *)
%o A124187 (PARI) for(n=1,10^4,if(ispseudoprime(sum(i=0,17,n^(2*i+1))+1),print1(n,", "))) \\ _Derek Orr_, Jun 24 2014
%o A124187 (Magma) [n: n in [0..4000] | IsPrime(1+n*(1+n^2)*(1+n^4+n^8)*(1+n^12+n^24))]; // _Vincenzo Librandi_, Jun 27 2014
%Y A124187 Cf. A049407.
%K A124187 nonn
%O A124187 1,2
%A A124187 _Artur Jasinski_, Dec 13 2006
%E A124187 a(45) and beyond from _Derek Orr_, Jun 24 2014