素数的判断

• 一个数 $a$ 想成为另外一个数 $b$ 的因子，能整除另一个数，该 $a\leq 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$ 是一个质数的话，那么对于任意一个数 $a$，则 $n|a^n-a$（$|$：表示整除）。其逆否命题便是，如果 $n$ 不能被 $a^n-a$ 整除，则 $n$ 一定不是素数。例为了判断 209 是不是质数，我们随便选择一个数，比如最简单的 $2$，结果发现 $2^{209}-2\;\%\;209=70$（即使把 209 换成上百位的大数，利用计算机也很快能算出这个余数来）。于是，209肯定不是素数。
但反过来说，满足 $n|a^n-a$，并不能保证 $n$ 就是素数。我们再举一个例子。为了判断 221 是不是质数，我们随机选择 a ，比如说还是 38 吧。你会发现 38221 – 38 除以 221 正好除尽。那么， 221 是否就一定是质数了呢？麻烦就麻烦在这里：这并不能告诉我们 221 是质数，因为 Fermat 小定理毕竟只说了对一切质数都成立，但没说对其他的数成不成立。万一 221 根本就不是质数，但 a = 38 时碰巧也符合 Fermat 小定理呢？为了保险起见，我们不妨再选一个不同的 a 值。比方说，令 a = 26 ，可以算出 $26^221 – 26$ 除以 221 余 169 ，因而 221 果然并不是质数。这个例子告诉了我们，如果运气不好的话，所选的 a 值会让不是质数的数也能骗过检测，虽然这个概率其实并不大。因此，我们通常的做法便是，多选几个不同的 a ，只要有一次没通过测试，被检测的数一定不是质数，如果都通过测试了，则被检测的数很可能是质数。没错， Fermat 素性测试的效率非常高，但它是基于一定概率的，有误报的可能。如果发现某个数 n 不满足 Fermat 小定理，它一定不是质数；但如果发现某个数 n 总能通过 Fermat 小定理的检验，只能说明它有很大的几率是质数。

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算法