一个数 a 想成为另外一个数
b 的因子,能整除另一个数,该 a≤b/2 ,因为 b 的因子如果有的话,它最大为b/2 互质的两个数不必都是质数,比如 91 与 1000
试除法
为了判断一个数是不是素数(质数),最笨的方法就是试除法——看它能不能被 2 整除,如果不能的话再看它能不能被 3 整除,这样不断试除下去,直到除遍所有比它小的数(一种稍微启发的做法是一直判断到该数的一半,再一种启发的做法是一直判断到该数的二次方根),都还不能把它分解开来,它就是质数了。
from math import sqrt
def is_prime(n):
for i in range(2, int(sqrt(n))+1):
if n%i == 0:
return false
return True
Fermat 素性测试
Fermat 是一种比较常用的高效方法,它基于如下原理:Fermat 小定理对一切质数都是成立的。Fermat 小定理的内容是:如果
n
是一个质数的话,那么对于任意一个数
但反过来说,满足
n|an−a
,并不能保证
n
就是素数。我们再举一个例子。为了判断 221 是不是质数,我们随机选择 a ,比如说还是 38 吧。你会发现 38221 – 38 除以 221 正好除尽。那么, 221 是否就一定是质数了呢?麻烦就麻烦在这里:这并不能告诉我们 221 是质数,因为 Fermat 小定理毕竟只说了对一切质数都成立,但没说对其他的数成不成立。万一 221 根本就不是质数,但 a = 38 时碰巧也符合 Fermat 小定理呢?为了保险起见,我们不妨再选一个不同的 a 值。比方说,令 a = 26 ,可以算出
def fermet_test(n, a):
if (a**n-a)%n:
return False
Fermat 素性测试真正麻烦的地方就是,居然有这么一种极其特殊的数,它不是质数,但对于任意的 a 值,它都能通过测试。这样的数叫做 Carmichael 数,最小的一个是 561 ,接下来的几个则是 1105, 1729, 2465, 2821, 6601, 8911… 虽然不多,但很致命。因此,在实际应用时,我们通常会选用 Miller-Rabin 素性测试算法。这个算法以 Gary Miller 的研究成果为基础,由 Michael Rabin 提出,时间大约是 1975 年。它可以看作是对 Fermat 素性测试的改良。如果选用了 k 个不同的 a 值,那么 Miller-Rabin 素性测试算法出现误判的概率不会超过 1 / 4k ,足以应付很多现实需要了。
有没有什么高效率的、确定性的质数判定算法呢?有,不过这已经是很后来的事情了。 2002 年, Manindra Agrawal 、 Neeraj Kayal 和 Nitin Saxena 发表了一篇重要的论文 PRIMES is in P ,给出了第一个高效判断质数的确定性算法,并以三人名字的首字母命名,叫做 AKS 素性测试。不过,已有的质数判断算法已经做得很好了,因此对于 AKS 来说,更重要的是它的理论意义。
References
[1] 跨越千年的RSA算法