OFFSET
1,1
COMMENTS
Primes in A154277.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
MAPLE
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
MATHEMATICA
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 *)
Select[Table[81*n^2-72*n+17, {n, 0, 2000}], PrimeQ] (* Vincenzo Librandi, Jul 16 2012 *)
PROG
(Magma) [a: n in [0..350] | IsPrime(a) where a is 81*n^2-72*n+17]; // Vincenzo Librandi, Jul 16 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jan 06 2009
EXTENSIONS
Extended by Emeric Deutsch, Jan 20 2009
STATUS
approved