Anonymous user
Factors of a Mersenne number: Difference between revisions
→{{header|Crystal}}
Line 883:
kp2 = p2 = 2.to_big_i * p
while (kp2 - 1) ** 2 < mers_num
q = kp2 + 1
if prime?(q) && [1, 7].includes?(q % 8) && (powmod(2, p, q) == 1)
# q is a factor of 2**p-1
|