A201107 - OEIS

%I #17 Sep 08 2022 08:46:00

%S 3,17,251,433,1459,2663,3457,16001,18523,35153,39367,48779,54001,

%T 65537,85751,170369,370387,410759,432001,715823,746497,913067,1272113,

%U 1557377,1714751,1769473,1940599,2450087,2735263,3456001,4812209,5488001,6615899,6750001

%N Primes of the form 2k^3+1.

%H Michael S. Branicky, <a href="/A201107/b201107.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..3000 from Vincenzo Librandi)

%F a(n) = 2*A168550(n)^3+1. - _R. J. Mathar_, Aug 20 2019

%t Select[Table[2n^3+1,{n,0,4000}],PrimeQ]

%o (Magma) [a: n in [0..200] | IsPrime(a) where a is 2*n^3+1]

%o (Python)

%o from sympy import isprime

%o print(list(filter(isprime, (2*k**3+1 for k in range(152))))) # _Michael S. Branicky_, Jun 17 2021

%Y Cf. A168550, A033562.

%K nonn,easy

%O 1,1

%A _Vincenzo Librandi_, Nov 27 2011