A168274 Primes p such that 10*p^3 + 1 is prime.

3, 31, 37, 43, 73, 109, 163, 211, 241, 313, 337, 409, 499, 673, 739, 757, 907, 1033, 1039, 1069, 1237, 1327, 1447, 1483, 1489, 1597, 1663, 1741, 1777, 1933, 2083, 2143, 2251, 2437, 2683, 2797, 3001, 3181, 3307, 3319, 3463, 3739, 3793, 4051

Offset: 1

Author

Eva-Maria Zschorn (e-m.zschorn(AT)zaschendorf.km3.de), Nov 22 2009

Keywords

• nonn

Comments

Subsequence of A168219.

Examples

			10*3^3+1 = 271 is prime, so prime 3 is in the sequence.
10*313^3+1 = 306642971 is prime, so prime 313 is in the sequence.
		

Links

• G. C. Greubel, Table of n, a(n) for n = 1..1000

Crossrefs

Cf. A168219, A168147.

Programs

Magma

[ p: p in PrimesUpTo(5000) | IsPrime(10*p^3+1) ]; // Vincenzo Librandi, Jan 29 2011

Mathematica

Select[Prime[Range[100]], PrimeQ[10*#^3 + 1] &] (* G. C. Greubel, Jul 16 2016 *)

PARI

select(p->isprime(10*p^3+1), primes(1000)) \\ Charles R Greathouse IV, Jul 16 2016

Extensions

Edited by Charles R Greathouse IV, Apr 24 2010