素数——详解

最新推荐文章于 2025-04-30 16:36:16 发布

枫退

最新推荐文章于 2025-04-30 16:36:16 发布

阅读量727

点赞数 10

文章标签：算法

本文链接：https://blog.csdn.net/w_w_wwwww/article/details/145272265

版权

素数与合数的定义，见数论基础。

素数计数函数：小于或等于的素数的个数，用表示。随着的增大，有这样的近似结果：。

素性测试

素性测试（Primality test）可以用于判定所给自然数是否为素数。

素性测试有两种：

确定性测试：绝对确定一个数是否为素数。常见例子包括试除法、Lucas–Lehmer 测试和椭圆曲线素性证明。
概率性测试：通常比确定性测试快很多，但有可能（尽管概率很小）错误地将合数识别为质数（尽管反之则不会）。因此，通过概率素性测试的数字被称为 可能素数，直到它们的素数可以被确定性地证明。而通过测试但实际上是合数的数字则被称为 伪素数。有许多特定类型的伪素数，最常见的是费马伪素数，它们是满足费马小定理的合数。概率性测试的常见例子包括 Miller–Rabin 测试。

试除法

暴力做法自然可以枚举从小到大的每个数看是否能整除。

参考实现

C++Python

bool isPrime(int a) {
  if (a < 2) return false;
  for (int i = 2; i < a; ++i)
    if (a % i == 0) return false;
  return true;
}

这样做是十分稳妥了，但是真的有必要每个数都去判断吗？

很容易发现这样一个事实：如果是的约数，那么也是的约数。

这个结论告诉我们，对于每一对，只需要检验其中的一个就好了。为了方便起见，我们之考察每一对里面小的那个数。不难发现，所有这些较小数就是这个区间里的数。

由于肯定是约数，所以不检验它。

参考实现

C++Python

bool isPrime(int a) {
  if (a < 2) return 0;
  for (int i = 2; (long long)i * i <= a; ++i)  // 防溢出
    if (a % i == 0) return 0;
  return 1;
}

Fermat 素性测试

Fermat 素性检验 是最简单的概率性素性检验。

我们可以根据费马小定理得出一种检验素数的思路：

基本思想是不断地选取在中的基，并检验是否每次都有。

参考实现

C++Python

bool millerRabin(int n) {
  if (n < 3) return n == 2;
  // test_time 为测试次数,建议设为不小于 8
  // 的整数以保证正确率,但也不宜过大,否则会影响效率
  for (int i = 1; i <= test_time; ++i) {
    int a = rand() % (n - 2) + 2;
    if (quickPow(a, n - 1, n) != 1) return false;
  }
  return true;
}

如果但不是素数，则被称为以为底的 伪素数。我们在实践中观察到，如果，那么通常是素数。但这里也有个反例：如果且，即使是合数，有。事实上，是最小的伪素数基数。

很遗憾，费马小定理的逆定理并不成立，换言之，满足了，也不一定是素数。甚至有些合数满足对任意满足的整数均有，这样的数称为 Carmichael 数。

Carmichael 函数

对正整数，定义 Carmichael 函数（卡迈克尔函数）为对任意满足的整数，使

恒成立的最小正整数 .

即：

Carmichael 函数有如下性质：

（Carmichael 定理）对任意素数和任意正整数，
证明

进而有：
1. 对任意正整数，有
2. 对任意正整数，，有
令的唯一分解式为，则

由中国剩余定理和 Carmichael 定理易证。

进而有：
1. 对任意正整数，，有

Carmichael 数

对于合数，如果对于所有正整数，和互素，都有同余式成立，则合数为 Carmichael 数（卡迈克尔数，OEIS:A002997）。

比如就是一个 Carmichael 数，同时也是最小的 Carmichael 数。

我们可以用如下方法判断合数是否为 Carmichael 数：

Korselt 判别法1

合数是 Carmichael 数当且仅当无平方因子且对的任意质因子均有 .

上述判别法可简化为：

Carmichael 数判别法

合数是 Carmichael 数当且仅当，其中为 Carmichael 函数。

Carmichael 数有如下性质：

Carmichael 数无平方因子且至少有个不同的质因子。
设为小于的 Carmichael 数个数，则：
1. （Alford, Granville, Pomerance. 19942）。
  
  由此可知 Carmichael 数有无限多个。
2. （Erdős. 19563），其中为常数。
  
  由此可知 Carmichael 数的分布十分稀疏。实际上，4。

注意

「若为 Carmichael 数，则也为 Carmichael 数」是错误的。

如为 Carmichael 数，考虑。

注意到，由 Korselt 判别法知，若是 Carmichael 数，则和均为的因子。

而，故，因此不是 Carmichael 数。

Miller–Rabin 素性测试

Miller–Rabin 素性测试（Miller–Rabin primality test）是更好的素数判定方法。它是由 Miller 和 Rabin 二人根据 Fermat 素性测试优化得到的。和其它概率性素数测试一样，它也只能检测出伪素数。要确保是素数，需要用慢得多的确定性算法。然而，实际上没有已知的数字通过了 Miller–Rabin 测试等高级概率性测试但实际上却是合数，因此我们可以放心使用。

在不考虑乘法的复杂度时，对数进行轮测试的时间复杂度是。Miller-Rabbin 素性测试常用于对高精度数进行测试，此时时间复杂度是，利用 FFT 等技术可以优化到

。

二次探测定理

如果是奇素数，则的解为或者。

要证明该定理，只需将上面的方程移项，再使用平方差公式，得到，即可得出上面的结论。

实现

根据 Carmichael 数的性质，可知其一定不是。

不妨将费马小定理和二次探测定理结合起来使用：

将中的指数分解为，在每轮测试中对随机出来的先求出，之后对这个值执行最多次平方操作，若发现非平凡平方根时即可判断出其不是素数，否则再使用 Fermat 素性测试判断。

还有一些实现上的小细节：

对于一轮测试，如果某一时刻，则之后的平方操作全都会得到，则可以直接通过本轮测试。
如果找出了一个非平凡平方根，则之后的平方操作全都会得到。可以选择直接返回 false，也可以放到次平方操作后再返回 false。

这样得到了较正确的 Miller Rabin：（来自 fjzzq2002）

参考实现

C++Python

bool millerRabin(int n) {
  if (n < 3 || n % 2 == 0) return n == 2;
  if (n % 3 == 0) return n == 3;
  int u = n - 1, t = 0;
  while (u % 2 == 0) u /= 2, ++t;
  // test_time 为测试次数，建议设为不小于 8
  // 的整数以保证正确率，但也不宜过大，否则会影响效率
  for (int i = 0; i < test_time; ++i) {
    // 0, 1, n-1 可以直接通过测试, a 取值范围 [2, n-2]
    int a = rand() % (n - 3) + 2, v = quickPow(a, u, n);
    if (v == 1) continue;
    int s;
    for (s = 0; s < t; ++s) {
      if (v == n - 1) break;  // 得到平凡平方根 n-1，通过此轮测试
      v = (long long)v * v % n;
    }
    // 如果找到了非平凡平方根，则会由于无法提前 break; 而运行到 s == t
    // 如果 Fermat 素性测试无法通过，则一直运行到 s == t 前 v 都不会等于 -1
    if (s == t) return 0;
  }
  return 1;
}

另外，假设广义 Riemann 猜想（generalized Riemann hypothesis, GRH）成立，则对数最多只需要测试中的全部整数即可确定数的素性，证明参见注释 7。

而在 OI 范围内，通常都是对范围内的数进行素性检验。对于范围内的数，选取三个数作为基底进行 Miller–Rabin 素性检验就可以确定素性；对于范围内的数，选取七个数作为基底进行 Miller–Rabin 素性检验就可以确定素性。参见注释 8。

也可以选取（即前个素数）检验范围内的素数。

注意如果要使用上面的数列中的数作为基底判断的素性：

所有的数都要取一遍，不能只选小于的；
把换成；
如果或，则直接通过该轮测试。

反素数

顾名思义，素数就是因子只有两个的数，那么反素数，就是因子最多的数（并且因子个数相同的时候值最小），所以反素数是相对于一个集合来说的。

一种符合直觉的反素数定义是：在一个正整数集合中，因子最多并且值最小的数，就是反素数。

反素数

对于某个正整数，如果任何小于的正数的约数个数都小于的约数个数，则称为是 反素数（anti-prime, a.k.a., highly compositive numbers）。

注意

注意区分 emirp，它表示的是逐位反转后是不同素数的素数（如 149 和 941 均为 emirp，101 不是 emirp）。

过程

那么，如何来求解反素数呢？

首先，既然要求因子数，首先要做的就是素因子分解。把分解成的形式，其中是素数，为他的指数。这样的话总因子个数就是。

但是显然质因子分解的复杂度是很高的，并且前一个数的结果不能被后面利用。所以要换个方法。

我们来观察一下反素数的特点。

反素数肯定是从开始的连续素数的幂次形式的乘积。
数值小的素数的幂次大于等于数值大的素数，即中，有。

解释：

如果不是从开始的连续素数，那么如果幂次不变，把素数变成数值更小的素数，那么此时因子个数不变，但是的数值变小了。交换到从开始的连续素数的时候值最小。
如果数值小的素数的幂次小于数值大的素数的幂，那么如果把这两个素数交换位置（幂次不变），那么所得的因子数量不变，但是的值变小。

另外还有两个问题，

对于给定的，要枚举到哪一个素数呢？

最极端的情况大不了就是，所以只要连续素数连乘到刚好小于等于就可以的呢。再大了，连全都一次幂，都用不了，当然就是用不到的啦！
我们要枚举到多少次幂呢？

我们考虑一个极端情况，当我们最小的素数的某个幂次已经比所给的（的最大值）大的话，那么展开成其他的形式，最大幂次一定小于这个幂次。unsigned long long 的最大值是，所以可以枚举到。

细节有了，那么我们具体如何具体实现呢？

我们可以把当前走到每一个素数前面的时候列举成一棵树的根节点，然后一层层的去找。找到什么时候停止呢？