15 lines
251 B
Text
15 lines
251 B
Text
print "Mersenne Primes :"
|
|
for p = 2 to 18
|
|
if lucasLehmer(p) then print "M"; p
|
|
next p
|
|
end
|
|
|
|
function lucasLehmer (p)
|
|
mp = (2 ^ p) - 1
|
|
sn = 4
|
|
for i = 2 to p-1
|
|
sn = (sn ^ 2) - 2
|
|
sn = sn - (mp * floor(sn / mp))
|
|
next
|
|
return sn = 0
|
|
end function
|