A201107 Primes of the form 2k^3+1.
3, 17, 251, 433, 1459, 2663, 3457, 16001, 18523, 35153, 39367, 48779, 54001, 65537, 85751, 170369, 370387, 410759, 432001, 715823, 746497, 913067, 1272113, 1557377, 1714751, 1769473, 1940599, 2450087, 2735263, 3456001, 4812209, 5488001, 6615899, 6750001
Offset: 1
Links
- Michael S. Branicky, Table of n, a(n) for n = 1..10000 (terms 1..3000 from Vincenzo Librandi)
Programs
-
Magma
[a: n in [0..200] | IsPrime(a) where a is 2*n^3+1];
-
Mathematica
Select[Table[2n^3+1,{n,0,4000}],PrimeQ]
-
Python
from sympy import isprime print(list(filter(isprime, (2*k**3+1 for k in range(152))))) # Michael S. Branicky, Jun 17 2021
Formula
a(n) = 2*A168550(n)^3+1. - R. J. Mathar, Aug 20 2019