%I #23 Sep 08 2022 08:46:03
%S 2,3,5,17,29,41,53,71,83,131,179,191,239,263,311,389,491,509,557,569,
%T 593,653,701,719,743,797,821,863,887,953,971,977,1019,1049,1097,1109,
%U 1277,1301,1319,1373,1427,1481,1523,1559,1601,1607,1613,1667,1787,1823
%N Primes p such that p^3 + p + 1 is prime.
%H Alois P. Heinz, <a href="/A216061/b216061.txt">Table of n, a(n) for n = 1..1000</a>
%p A := {}; for n to 1000 do p := ithprime(n); if isprime(p^3+p+1) then A := `union`(A, {p}) end if end do; A := A
%t Select[Prime[Range[400]], PrimeQ[#^3 + # + 1] &] (* _Bruno Berselli_, Sep 01 2012 *)
%o (Magma) [p: p in PrimesUpTo(2000)  IsPrime(p^3+p+1)]; // _Bruno Berselli_, Sep 01 2012
%Y Cf. A053182.
%Y Subsequence of A045309.
%K nonn,easy
%O 1,1
%A _César Eliud Lozada_, Aug 31 2012
