%I #14 Sep 08 2022 08:45:40
%S 17,197,30977,52901,80657,215297,287297,414737,462401,512657,739601,
%T 902501,972197,1008017,1196837,1401857,1623077,2496401,2611457,
%U 3422501,3763601,4343057,4494401,4726277,4884101,5456897,5541317,6330257
%N Primes of the form 81*n^2 - 72*n + 17.
%C Primes in A154277.
%H Vincenzo Librandi, <a href="/A154276/b154276.txt">Table of n, a(n) for n = 1..1000</a>
%p a := proc (n) if isprime(81*n^2-72*n+17) = true then 81*n^2-72*n+17 else end if end proc: seq(a(n), n = 0 .. 300); # _Emeric Deutsch_, Jan 20 2009
%t lst={};Do[If[PrimeQ[p=81*n^2-72*n+17],AppendTo[lst,p]],{n,0,7!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Jan 11 2009 *)
%t Select[Table[81*n^2-72*n+17,{n,0,2000}],PrimeQ] (* _Vincenzo Librandi_, Jul 16 2012
%o (Magma) [a: n in [0..350] | IsPrime(a) where a is 81*n^2-72*n+17]; // _Vincenzo Librandi_, Jul 16 2012
%Y Cf. A154277.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Jan 06 2009
%E Extended by _Emeric Deutsch_, Jan 20 2009
|