Anonymous user
Miller–Rabin primality test: Difference between revisions
Improved python code: cleaned up style, added comments, added doctests
(added perl) |
(Improved python code: cleaned up style, added comments, added doctests) |
||
Line 554:
=={{header|Python}}==
<lang python>
import random
if n<=3: return True▼
_mrpt_num_trials = 5 # number of bases to test
if n % 2 == 0: return False▼
d = n - 1▼
def is_probable_prime(n):
"""
Miller-Rabin primality test.
A return value of False means n is certainly not prime. A return value of
True means n is very likely a prime.
>>> is_probable_prime(1)
Traceback (most recent call last):
AssertionError
>>> is_probable_prime(2)
True
>>> is_probable_prime(3)
True
>>> is_probable_prime(4)
False
>>> is_probable_prime(5)
True
>>> is_probable_prime(123456789)
False
>>>
>>> len(primes_under_1000)
168
>>> primes_under_1000[-10:]
[937, 941, 947, 953, 967, 971, 977, 983, 991, 997]▼
>>> is_probable_prime(6438080068035544392301298549614926991513861075340134\
3291807343952413826484237063006136971539473913409092293733259038472039\
7133335969549256322620979036686633213903952966175107096769180017646161\
851573147596390153)
True
>>> is_probable_prime(7438080068035544392301298549614926991513861075340134\
3291807343952413826484237063006136971539473913409092293733259038472039\
7133335969549256322620979036686633213903952966175107096769180017646161\
851573147596390153)
False
"""
assert n >= 2
# special case 2
# ensure n is odd
if n % 2 == 0:
return False
# write n-1 as 2**s * d
# repeatedly try to divide n-1 by 2
s = 0
while
quotient, remainder = divmod(d, 2)
if remainder == 1:
s += 1
assert(2**s *
▲ x = pow(a, d, n)
# test the base a to see whether it is a witness for the compositeness of n
▲ if x != 1:
def
if pow(a, d,
return
for i in
if
return False
return True # n is definitely composite
for i in range(_mrpt_num_trials):
a = random.randint(2, n-1)
if try_composite(a):
return True # no base tested showed n as composite
▲>>> [ i for i in range(1,1000) if miller_rabin(i, 10)][-10:]
▲[937, 941, 947, 953, 967, 971, 977, 983, 991, 997]
▲>>></lang>
=={{header|Ruby}}==
|