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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A124181 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 is prime. 16
1, 3, 69, 86, 104, 110, 138, 146, 210, 238, 247, 260, 264, 269, 316, 436, 572, 600, 621, 654, 666, 715, 737, 740, 744, 754, 779, 1056, 1156, 1159, 1216, 1218, 1221, 1343, 1419, 1434, 1442, 1524, 1580, 1603, 1676, 1680, 1731, 1742, 1804, 1952, 1956, 1985 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

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

FORMULA

1 together with numbers n such that (n^25-n)/(n^2-1) + 1 is prime. - Charles R Greathouse IV, Jul 02 2013

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], Print[n]], {n, 1, 1400}]

Select[Range[3000], PrimeQ[Total[#^Range[1, 23, 2]] + 1] &] (* Vincenzo Librandi, Jun 28 2014 *)

PROG

(PARI) is(n)=n==1 || isprime((n^25-n)/(n^2-1)+1) \\ Charles R Greathouse IV, Jul 02 2013

(Sage)

i, n = var('i, n')

[n for n in (1..2000) if is_prime(1+(n^(2*i+1)).sum(i, 0, 11))] # Bruno Berselli, Jun 27 2014

(MAGMA) [n: n in [0..2000] | IsPrime(s) where s is 1+&+[n^i: i in [1..23 by 2]]]; // Vincenzo Librandi, Jun 28 2014

CROSSREFS

Cf. A049407, similar sequences listed in A244376.

Sequence in context: A073163 A279491 A264700 * A046432 A232326 A270869

Adjacent sequences:  A124178 A124179 A124180 * A124182 A124183 A124184

KEYWORD

nonn,easy

AUTHOR

Artur Jasinski, Dec 13 2006

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 April 29 11:38 EDT 2017. Contains 285604 sequences.