login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A124187 Numbers n such that 1 + n + n^3 + n^5 + n^7 + n^9 + n^11 + ... + n^33 + n^35 is prime. 5
1, 7, 10, 17, 52, 69, 108, 161, 173, 231, 306, 330, 338, 352, 416, 582, 584, 593, 635, 767, 834, 855, 868, 892, 927, 944, 950, 1044, 1060, 1203, 1242, 1299, 1302, 1509, 1520, 1551, 1637, 1972, 2067, 2078, 2135, 2303, 2310, 2366, 2416, 2511, 2514, 2556, 2581 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..258

MAPLE

a:= proc(n) option remember; local k;

      for k from 1+ a(n-1) while

        not isprime(1+(k^37-k)/(k^2-1)) do od; k

    end: a(1):=1:

seq(a(n), n=1..30);  # Alois P. Heinz, Jun 26 2014

MATHEMATICA

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}]

Select[Range[5000], PrimeQ[Total[#^Range[1, 35, 2]] + 1] &] (* Vincenzo Librandi, Jun 27 2014 *)

PROG

(PARI) for(n=1, 10^4, if(ispseudoprime(sum(i=0, 17, n^(2*i+1))+1), print1(n, ", "))) \\ Derek Orr, Jun 24 2014

(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

CROSSREFS

Cf. A049407.

Sequence in context: A192273 A098748 A020692 * A155027 A118420 A196939

Adjacent sequences:  A124184 A124185 A124186 * A124188 A124189 A124190

KEYWORD

nonn

AUTHOR

Artur Jasinski, Dec 13 2006

EXTENSIONS

a(45) and beyond from Derek Orr, Jun 24 2014

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy .

Last modified June 28 20:31 EDT 2017. Contains 288840 sequences.