A236952 Primes p such that p^4 - p +/- 1 are twin primes.

13, 79, 3571, 3739, 5023, 5443, 7459, 7621, 10243, 13339, 14251, 17359, 17551, 17863, 17971, 18061, 19483, 21481, 27631, 32611, 37501, 38821, 48463, 49711, 54709, 56443, 57073, 57751, 69313, 71353, 72883, 74293, 81883, 82051, 84223

Offset

1,1

Comments

Intersection of A236940 and A236071.

Links

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

Example

13 is prime and 13^4-13-1 (28547) and 13^4-13+1 (28549) are twin primes. So, 13 is a member of this sequence.

Mathematica

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

Prog

(Python)
import sympy
from sympy import isprime
{print(n) for n in range(10**6) if isprime(n) and isprime(n**4-n-1) and isprime(n**4-n+1)}
(Magma) [p: p in PrimesUpTo(90000) | IsPrime(p^4-p-1) and IsPrime(p^4-p+1)]; // Vincenzo Librandi, Feb 14 2014

Crossrefs

Cf. A236940, A236071.

Sequence in context: A032625 A120782 A032652 * A136373 A071614 A045552

Adjacent sequences: A236949 A236950 A236951 * A236953 A236954 A236955

Keyword

nonn

Author

Derek Orr, Feb 01 2014

Status

approved