|
%I #23 Sep 08 2022 08:46:20
%S 53,3329,13313,13631489,3489660929,62864142619960717084721153,
%T 5100145160001678120616578906356228963083163798627028041729,
%U 6779255729241169695101387251026410519979286814120235842117075415451380965612384558178346467329,1735489466685739441945955136262761093114697424414780375581971306355553527196770446893656695635969
%N Primes of the form 13*2^n + 1.
%C For the corresponding exponents n see A032356.
%H Ray Ballinger, Wilfried Keller, <a href="http://www.prothsearch.com/riesel1.html">List of primes k*2^n + 1 for k < 300</a>.
%F a(n) = A168596(A032356(n)). - _Michel Marcus_, Mar 29 2018
%p a:=(n,k)->`if`(isprime(k*2^n+1), k*2^n+1, NULL):
%p seq(a(n,13), n=1..316);
%t Select[Table[13 2^n + 1, {n, 400}], PrimeQ] (* _Vincenzo Librandi_, Mar 06 2018 *)
%o (Magma) [a: n in [1..400] | IsPrime(a) where a is 13*2^n + 1]; // _Vincenzo Librandi_, Mar 06 2018
%o (GAP) Filtered(List([1..500],n->13*2^n + 1),IsPrime); # _Muniru A Asiru_, Mar 06 2018
%o (PARI) lista(nn) = {for(k=1, nn, if(ispseudoprime(p=13*2^k+1), print1(p, ", ")));} \\ _Altug Alkan_, Mar 29 2018
%Y Cf. A019434, A039687, A032356, A050527, A050528, A050529, A173236, A195745, A050526, A300407, A300408.
%K nonn
%O 1,1
%A _Martin Renner_, Mar 05 2018
|
