OFFSET
1,1
MATHEMATICA
lst={}; Do[If[PrimeQ[p=81*n^2-90*n+26], If[PrimeQ[2*p+1], AppendTo[lst, p]]], {n, 8!}]; lst
Select[Table[81 n^2-90 n+26, {n, 6000}], AllTrue[{#, 2#+1}, PrimeQ]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Jun 09 2021 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 11 2009
STATUS
approved