%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
|