素数测试

欧几里得算法

输入:m，n
输出:m和n的最大公约数
步骤：
1. 如果$n=0$,返回$m$的值作为结果，结束。否则进入2；
2. 计算m/n, 将余数赋值给r。
3. 将n的值赋值给m，将r的值赋值给n， 返回第一步。
其实这里的思想十分清晰，被除数与除数的最大公约数等于除数与余数的最大公约数。

实际上根据上面的描述我们将代码写成递归形式：

int Euclid_gcd(int m, int n){
    if (n == 0) return m;
    return Euclid_gcd(n, m % n);
}

也可以写成递推形式：

int Euclid_gcd1(int m, int n){
    while (n != 0){
        int r = m % n;
        m = n;
        n = r;
    }
    return m;
}

素数测试1

如果两个整数相乘为N，则必有一个数小于等于$\sqrt N$
所以我们只需要测试：

for i in range(2,sqrt(n)+1):
   remain = N % i
   if(remain==0):
      return not prime
return prime

素数测试2

公式：
if:

$$x(mod N)=x'(mod N) \\ y(mod N)=y'(mod N)$$
then:
$$x+y(mod N)=x'+y'(mod N) \\ xy(mod N) = x'y'(mod N)$$

假设$p$是一个素数，$1<=a < p$,可以得到$\{1,2,3,4,.......,p-1\}$mod p的值就是$\{1,2,3,4,.......,p-1\}$,而$\{a,2a,3a,4a,.......,(p-1)a\}$mod p的值也在$\{1,2,3,4,.......,p-1\}$，理由如下：如果上面(p-1)个有两个数的余数相等，设为$n$和$m$且($n> m$)，那么n和m，就需要满足$n-m = kp$,(k为大于0的整数)，同时由于$n=k_1a, m =k_2a$,所以$n-m = k_3a$,所以我们就可以得出$kp=k_3a$,从而推出$p=\frac{k_3}{k}a$,（不管怎样都可以分解成两个整数相乘）从这里可以推出p可以由两个整数相乘可以，而$p$是素数，这与条件不符，所以不存在两个余数相等的数，所以$\{a,2a,3a,4a,.......,(p-1)a\}$中每一个的余数不相同。根据两个组数同余，那么两组数的乘数也是同余。

$$(p-1)!=a^{p-1}(p-1)!mod p \\ 1=a^{p-1}mod p$$
上面就推出了费马小定理。
所以我们判定一个数是不是素数，只需要选择一个 $a$，然后看$a^{p-1}$的模是不是1，如果不是1，说明是合数，但是是1，也不能说它就是素数(少量的合数也满足上述公式)，所以为了筛出合数，我们一般会选择几个 a <script type="math/tex" id="MathJax-Element-6328">a</script>.