This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A301633 #30 Dec 08 2024 05:46:59
%S A301633 0,1,33,123,235,513,3441,5836,71071
%N A301633 Numbers k such that 4^(2*k) + 2*k + 1 is a prime.
%C A301633 Next term, if it exists, is greater than 10000. - _Vaclav Kotesovec_, Mar 25 2018
%p A301633 a:=n->`if`(isprime(4^(2*k)+2*k+1),k,NULL): seq(a(k),k=0..5000); # _Muniru A Asiru_, Mar 25 2018
%t A301633 Flatten[{0, Select[Range[1000], PrimeQ[16^# + 2*# + 1] &]}] (* _Vaclav Kotesovec_, Mar 25 2018 *)
%o A301633 (PARI) for(n=0, 500, if(isprime(4^(2*n)+2*n+1), print1(n", ")))
%Y A301633 A301637 gives the primes.
%Y A301633 Cf. A173053, A301632.
%K A301633 nonn,more
%O A301633 1,3
%A A301633 _Seiichi Manyama_, Mar 25 2018
%E A301633 a(7) from _Muniru A Asiru_, Mar 25 2018
%E A301633 a(8) from _Vaclav Kotesovec_, Mar 25 2018
%E A301633 a(9) from _Michael S. Branicky_, Dec 08 2024