login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A236940 Primes p such that p^4-p+1 is prime. 2
3, 13, 19, 43, 79, 109, 313, 379, 613, 709, 1171, 1213, 1399, 1543, 1693, 1759, 1861, 1933, 2089, 2239, 2341, 2371, 2503, 2521, 2731, 2749, 3001, 3061, 3229, 3433, 3571, 3739, 3769, 4219, 4801, 4933, 4951, 4993, 5011, 5023, 5209, 5281 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

Primes in the sequence A236761.

LINKS

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

EXAMPLE

1213 is prime and 1213^4 - 1213 + 1 = 2164926732949 is prime. Thus, 1213 is a member of this sequence.

MATHEMATICA

Select[Prime[Range[10000]], PrimeQ[#^4 - # + 1]&] (* Vincenzo Librandi, Feb 14 2014 *)

PROG

(Python)

import sympy

from sympy import isprime

{print(n) for n in range(10**4) if isprime(n) and isprime(n**4-n+1)}

(PARI)

s=[]; forprime(p=2, 6000, if(isprime(p^4-p+1), s=concat(s, p))); s \\ Colin Barker, Feb 05 2014

(MAGMA) [p: p in PrimesUpTo(6000) | IsPrime(p^4-p+1)]; // Vincenzo Librandi, Feb 14 2014

CROSSREFS

Cf. A236761.

Sequence in context: A229494 A211438 A018621 * A024469 A007697 A055202

Adjacent sequences:  A236937 A236938 A236939 * A236941 A236942 A236943

KEYWORD

nonn

AUTHOR

Derek Orr, Feb 01 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
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 18 08:40 EDT 2021. Contains 345098 sequences. (Running on oeis4.)