数学--质数

质数

1只包含$1$和本身两个约数

1. 判断质数

试除法，只需枚举$2\sim \sqrt{n}$即可

1. 代码：

bool is_prime(int n)
{
    if (n < 2) return false;
    for (int i = 2; i <= n / i; i ++ )
        if (x % i == 0)
            return false;
    return true;
}

2. 细节

• 不推荐写$\sqrt{n}$，因为要调用函数，太慢。
• 不推荐写$i\times i<=n$，当$n$太大时，$i\times i$有可能溢出。

2. 分解质因数

$n$最多只含一个$>\sqrt{n}$的质因子，仍使用试除法，遇到约数就约掉。如果约掉了所有$<\sqrt{n}$的质因子之后，剩余的数仍大于$1$，则剩下的数必是质因子，单独处理即可。

代码：

void divide(int n)
{
    for (int i = 2; i <= n / i; i ++ )
        if (n % i == 0)
        {
            // 约掉所有因子，易证从小到大约掉的因子必是质因子
            int s = 0;
            while (x % i == 0) x /= i, s ++ ;
            cout << i << ' ' << s << endl; // 输出质因子及其个数
        }
    if (x > 1) cout << x << ' ' << 1 << endl; // 单独处理最后一个因子
    cout << endl;
}

3. 筛质数

给定一个范围$1\sim n$，筛掉范围内所有的质数。

(1) 朴素筛法 $\to O(n\times \log n)$

1. 思路：从$2\sim n$枚举，依次筛掉$2$的倍数、$3$的倍数、$4$的倍数，以此类推。

2. 例题: 洛谷P5736 【深基7.例2】质数筛

(2) 埃氏筛法 $\to O(n\times \log (\log n))$

1. 优化：

只筛掉质数的倍数，例如在(1)中枚举到$4$，由于$4$已经在$2$的时候就被筛掉了，就跳过$4$不再筛$4$的倍数了。显然，如果枚举到某个数$i$，且$i$之前没有被筛掉，那么$i$必是质数。

2. 缺点：仍旧还是有些数被重复筛了，例如$6$会被$2$和$3$同时筛掉，效率仍待提高。

(3) 线性筛法 $\to O(n)$

1. 优化：每个数只会被自己最小的质因子筛掉，例如上面$6$只会被$2$筛掉。

2. 思路：

• $i\mathrm{mod}{p}_j=0\to p_j$一定是$i$的最小质因子，$p_j$也一定是$p_j\times i$的最小质因子。
• $i\mathrm{mod}{p}_j\ne 0\to p_j$一定小于$i$的所有质因子，$p_j$也一定是$p_j\times i$的最小质因子。
• 按照此逻辑来决定一个数一个数地筛。

3. 代码：

int primes[N], cnt;      // primes[]存储所有素数
bool st[N];           // st[x]存储x是否被筛掉

void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

4. 例题 洛谷P3383 【模板】线性筛素数