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 A164041 #36 Sep 08 2022 08:45:47
%S A164041 17,31,71,127,647,1151,2887,3527,7687,12799,19207,20807,23327,34847,
%T A164041 39199,49927,53791,73727,79999,103967,117127,145799,172871,194687,
%U A164041 220447,279751,294911,323207,336199,387199,394271,419527,438047,587527,649799,724807
%N A164041 Primes of the form 2*p^2 + 4*p + 1, where p is also prime.
%C A164041 A subsequence of the primes of the form 2k^2+4k+1 = 2*(k+1)^2-1, A066436. - _R. J. Mathar_, Aug 10 2009
%H A164041 Vincenzo Librandi, <a href="/A164041/b164041.txt">Table of n, a(n) for n = 1..1000</a>
%F A164041 a(n) = 2*(A164042(n))^2 + 4*A164042(n) + 1.
%t A164041 lst={}; Do[p=Prime@n; a=2*p^2+4*p+1; If[PrimeQ@a,AppendTo[lst,a]],{n,7!}]; lst (* _Vladimir Joseph Stephan Orlovsky_, Aug 12 2009 *)
%t A164041 Select[Table[2p^2+4p+1,{p,Prime[Range[200]]}],PrimeQ] (* _Harvey P. Dale_, Aug 25 2019 *)
%o A164041 (Magma) [a: p in PrimesUpTo(700)|IsPrime(a) where a is 2*p^2+4*p+1 ] // _Vincenzo Librandi_, Sep 01 2012
%o A164041 (PARI) lista(nn) = {forprime(p=2, nn, if(isprime(q=2*p^2+4*p+1), print1(q, ", ")));} \\ _Altug Alkan_, Mar 29 2018
%K A164041 nonn,easy
%O A164041 1,1
%A A164041 _Vincenzo Librandi_, Aug 08 2009
%E A164041 a(29) corrected by _R. J. Mathar_, Aug 11 2009
%E A164041 Edited by _N. J. A. Sloane_, Aug 11 2009