OFFSET
1,2
COMMENTS
The exponent of 2 in the expression, 521, is a Mersenne exponent.
EXAMPLE
n = 1 is a term since 2^521-1 is prime (13th Mersenne prime).
MATHEMATICA
Select[Range@ 12050, PrimeQ[# 2^521 - 1] &] (* Michael De Vlieger, Dec 09 2015 *)
PROG
(MATLAB)
if isprime(n*2^521-1)
disp(n)
end
(Magma) [n: n in [1..2*10^4] |IsPrime(n*2^521-1)]; // Vincenzo Librandi, Dec 10 2015
(PARI) is(n) = ispseudoprime(n*2^521 - 1); \\ Altug Alkan, Dec 10 2015
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vardan Semerjyan, Dec 09 2015
STATUS
approved