%I #8 Sep 08 2022 08:45:40
%S 2,5,7,13,23,47,79,89,103,109,137,139,193,223,229,257,277,293,317,347,
%T 379,389,443,457,499,593,607,653,677,709,733,823,853,883,907,953,1013,
%U 1033,1087,1103,1117,1283,1297,1307,1327,1367,1409,1489,1523,1549,1597
%N Primes p such that p^3 + 2*p^2 + p + 1 is also prime.
%H Vincenzo Librandi, <a href="/A155028/b155028.txt">Table of n, a(n) for n = 1..1000</a>
%t Select[Prime[Range[1600]], PrimeQ[(#^3 + 2*#^2 + # + 1)]&] (* _Vincenzo Librandi_, Oct 30 2012 *)
%o (Magma) [p: p in PrimesUpTo(1600)| IsPrime(p^3 + 2*p^2 + p + 1)]; // _Vincenzo Librandi_, Oct 30 2012
%Y Cf. A154686.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Jan 19 2009
|