Factors of a Mersenne number: Difference between revisions

Line 884:
while (k * 2 * p - 1) ** 2 < mp_num
q = k * 2 * p + 1
if prime?(q) && (q % 8 == [1, || 7].includes?(q % 8 == 7) && (powmod(2, p, q) == 1)
# q is a factor of 2**p-1
return q
Anonymous user