A155028 Primes p such that p^3 + 2*p^2 + p + 1 is also prime.

2, 5, 7, 13, 23, 47, 79, 89, 103, 109, 137, 139, 193, 223, 229, 257, 277, 293, 317, 347, 379, 389, 443, 457, 499, 593, 607, 653, 677, 709, 733, 823, 853, 883, 907, 953, 1013, 1033, 1087, 1103, 1117, 1283, 1297, 1307, 1327, 1367, 1409, 1489, 1523, 1549, 1597

Offset: 1

Author

Vincenzo Librandi, Jan 19 2009

Keywords

• nonn
• easy

Links

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

Crossrefs

Cf. A154686.

Programs

Magma

[p: p in PrimesUpTo(1600)| IsPrime(p^3 + 2*p^2 + p + 1)]; // Vincenzo Librandi, Oct 30 2012

Mathematica

Select[Prime[Range[1600]], PrimeQ[(#^3 + 2*#^2 + # + 1)]&] (* Vincenzo Librandi, Oct 30 2012 *)