Lucas-Lehmer test: Difference between revisions

m
minor formatting changes
(Added support up to M127)
m (minor formatting changes)
Line 1,497:
is mersenne prime = p .EQ. 2
ELSE
n = 22_16**INT(p, 16) - 1
s = 4
DO 10 i = 1, p - 2
Line 1,524:
b = b / 2
IF (a .LT. m - a) THEN
a = a * 2
ELSE
a = a - m + a
357

edits