%I #10 Sep 08 2022 08:46:06
%S 2,4,11,14,15,21,31,35,41,45,111,130,136,140,155,176,189,221,230,239,
%T 274,316,406,414,441,465,466,504,521,561,570,580,584,591,686,689,696,
%U 759,834,836,860,869,904,960,1026,1159,1379,1539,1614,1625,1660
%N Numbers n such that n^3 - n +/- 1 are twin primes.
%H Vincenzo Librandi, <a href="/A236524/b236524.txt">Table of n, a(n) for n = 1..2600</a>
%e 561^3 - 561 + 1 (176557921) and 561^3 - 561 - 1 (176557919) are twin primes. Thus, 561 is a member of this sequence.
%t Select[Range[2000], PrimeQ[#^3 - # - 1] && PrimeQ[#^3 - # + 1] &] (* _Vincenzo Librandi_, Jan 30 2018 *)
%o (Python)
%o import sympy
%o from sympy import isprime
%o {print(n) for n in range(10**4) if isprime(n**3-n-1) and isprime(n**3-n+1)}
%o (Magma) [n: n in [1..2000] | IsPrime(n^3-n-1) and IsPrime(n^3-n+1)]; // _Vincenzo Librandi_, Jan 30 2018
%Y Intersection of A126421 and A236477.
%K nonn
%O 1,1
%A _Derek Orr_, Jan 27 2014