OFFSET
1,1
COMMENTS
For the corresponding exponents n see A032356.
LINKS
Ray Ballinger, Wilfried Keller, List of primes k*2^n + 1 for k < 300.
FORMULA
MAPLE
a:=(n, k)->`if`(isprime(k*2^n+1), k*2^n+1, NULL):
seq(a(n, 13), n=1..316);
MATHEMATICA
Select[Table[13 2^n + 1, {n, 400}], PrimeQ] (* Vincenzo Librandi, Mar 06 2018 *)
PROG
(Magma) [a: n in [1..400] | IsPrime(a) where a is 13*2^n + 1]; // Vincenzo Librandi, Mar 06 2018
(GAP) Filtered(List([1..500], n->13*2^n + 1), IsPrime); # Muniru A Asiru, Mar 06 2018
(PARI) lista(nn) = {for(k=1, nn, if(ispseudoprime(p=13*2^k+1), print1(p, ", "))); } \\ Altug Alkan, Mar 29 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Martin Renner, Mar 05 2018
STATUS
approved