A301633 Numbers k such that 4^(2*k) + 2*k + 1 is a prime.
0, 1, 33, 123, 235, 513, 3441, 5836, 71071
Offset: 1
Programs
-
Maple
a:=n->`if`(isprime(4^(2*k)+2*k+1),k,NULL): seq(a(k),k=0..5000); # Muniru A Asiru, Mar 25 2018
-
Mathematica
Flatten[{0, Select[Range[1000], PrimeQ[16^# + 2*# + 1] &]}] (* Vaclav Kotesovec, Mar 25 2018 *) -
PARI
for(n=0, 500, if(isprime(4^(2*n)+2*n+1), print1(n", ")))
Extensions
a(7) from Muniru A Asiru, Mar 25 2018
a(8) from Vaclav Kotesovec, Mar 25 2018
a(9) from Michael S. Branicky, Dec 08 2024
Comments