%I #21 May 03 2024 19:52:16
%S 577,31249,40897,64081,77617,126001,186049,208657,581041,982801,
%T 1367857,1616401,2101249,2567377,3836449,4038481,4141441,4566241,
%U 4898449,5359537,7691041,8413201,8711137,10941841,12152449,13240657,14183137
%N Primes of the form 648n^2 - 72n + 1.
%C The primes occur at n = 1, 7, 8, 10, 11, 14, 17, 18, 30, 39, 46, 50, 57,..
%H Vincenzo Librandi, <a href="/A154511/b154511.txt">Table of n, a(n) for n = 1..1000</a>
%t lst={};Do[p=648*n^2-72*n+1;If[PrimeQ[p],AppendTo[lst,p]],{n,0,5*5!}];lst (* _Vladimir Joseph Stephan Orlovsky_, Jan 27 2009 *)
%t Select[Table[648n^2-72n+1,{n,0,6001}],PrimeQ] (* _Vincenzo Librandi_, Jul 23 2012 *)
%o (Magma) [a: n in [0..200] | IsPrime(a) where a is 648*n^2-72*n+1]; // _Vincenzo Librandi_, Jul 23 2012
%o (PARI) (for (n=0, 200, if (isprime (k=648*n^2-72*n+1), print1 (k, “, “)))) ; \\ _Vincenzo Librandi_, Jul 23 2012
%Y Cf. A154510.
%K nonn,easy
%O 1,1
%A _Vincenzo Librandi_, Jan 11 2009
%E Incorrect formula replaced with a comment by _R. J. Mathar_, Jul 06 2010
|