OFFSET
0,1
COMMENTS
Since a(15) = 0, M(17) = 131071 is prime.
LINKS
Eric Weisstein's World of Mathematics, Lucas Lehmer Test.
Wikipedia, Lucas Lehmer Primality Test.
FORMULA
a(0) = 4, a(n) = a(n-1)^2 mod 2^p-1. Last term: a(p-2).
EXAMPLE
a(15) = 130559^2 - 2 mod 131071 = 0.
CROSSREFS
KEYWORD
fini,full,nonn
AUTHOR
Sergio Pimentel, Apr 04 2007
STATUS
approved