A236952 - OEIS

%I #16 Sep 08 2022 08:46:06

%S 13,79,3571,3739,5023,5443,7459,7621,10243,13339,14251,17359,17551,

%T 17863,17971,18061,19483,21481,27631,32611,37501,38821,48463,49711,

%U 54709,56443,57073,57751,69313,71353,72883,74293,81883,82051,84223

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

%C Intersection of A236940 and A236071.

%H Vincenzo Librandi, <a href="/A236952/b236952.txt">Table of n, a(n) for n = 1..900</a>

%e 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.

%t Select[Prime[Range[10000]], PrimeQ[#^4 - # - 1] && PrimeQ[#^4 - # + 1]&] (* _Vincenzo Librandi_, Feb 14 2014 *)

%o (Python)

%o import sympy

%o from sympy import isprime

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

%o (Magma) [p: p in PrimesUpTo(90000) | IsPrime(p^4-p-1) and IsPrime(p^4-p+1)]; // _Vincenzo Librandi_, Feb 14 2014

%Y Cf. A236940, A236071.

%K nonn

%O 1,1

%A _Derek Orr_, Feb 01 2014