素数（算法、C++）

        素数定义相信大家很多都知道了，简单来说，就是只能被1和它本身除的数，就叫做素数。以下是一些素数常见的题目。

        一、判断一个数是否是素数

        我们一般采用试除法来判断。（题目比较简单，就不附上注释了）

        题目

---

时间复杂度：O(sqrt(n))

代码：

#include<iostream>

using namespace std;

int n;

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

int main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        int x;
        scanf("%d",&x);
        bool result=isPrime(x);
        if(result) printf("Yes\n");
        else printf("No\n");
    }
    return 0;
}

二、分解质因数

我们也是采用试除法来解决。题目也比较简单，这里就不附上注释了。

题目

---

时间复杂度：O(sqrt(n))

代码：

#include<iostream>

using namespace std;

int n;

void divide(int n){
    for(int i=2;i<= n/i ;i++){
        if(n%i==0){
            int res=0;
            while(n%i==0){
                n /= i;
                res++;
            }
            printf("%d %d\n",i,res);
        }
    }
    if(n>1) printf("%d %d\n",n,1);
    puts("");
}
int main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++){
    int x;
    scanf("%d",&x);
    divide(x);
    }
    return 0;
}

---

三、筛素数

题目

---

筛素数的话，这里有两种筛选方法。一种是埃氏筛法，还有一种是线性筛法。

（1）埃氏筛法

时间复杂度O(n*logn*logn)

埃式筛选法是一种在内找出2 − N 所有质数的算法,基本思路是遍历2 − N ,如果当前数未被筛选那么他就是素数,之后淘汰当前数的所有倍数,当所有数字被遍历完毕后,最后剩下的就是素数
 

for(int i=2;i<=N;i++){
	//被剔除了,直接pass
	if(flag[i]) continue;
	//未被剔除,当前数i是素数,淘汰它的所有倍数
	for(int j=i+i;j<=N;j=j+i){
		flag[j]=true;
	}
}

---

有一种加速的形式,吧j=i+i改为j=ii,因为2i,3*i,…,(i-1)*i都在前面被遍历过了。

for(int i=2;i<=N;i++){
	//被剔除了,直接pass
	if(flag[i]) continue;
	//未被剔除,当前数i是素数,淘汰它的所有倍数
	for(int j=i*i;j<=N;j=j+i){
		flag[j]=true;
	}
}

附上题目代码：

#include <iostream>
#include <vector>

using namespace std;

int countPrimes(int n) {
    vector<bool> isPrime(n + 1, true);
    int count = 0;
    
    for (int i = 2; i * i <= n; ++i) {
        if (isPrime[i]) {
            for (int j = i * i; j <= n; j += i) {
                isPrime[j] = false;
            }
        }
    }
    
    for (int i = 2; i <= n; ++i) {
        if (isPrime[i]) {
            count++;
        }
    }
    
    return count;
}

int main() {
    int n;
    cin >> n;
    
    cout << countPrimes(n) << endl;
    
    return 0;
}

---

（2）线性筛法

线性筛法的核心：每个数只有一个最小质因子，且会被他的最小质因子筛掉。

//首先定义一个数组primes[]用来存放质数，cnt用来记录质数个数
//st[]用来标记这个数是否被筛掉
void get_primes(int n){
    for(int i = 2; i <= n; i++){//筛选出2~n之间的质数
        if(!st[i])    primes[cnt++] = i;//如果没有被标记，则它一定是质数，将其放入primes[]数组
        for(int j = 0; primes[j] <= n / i; j++){//从小到大枚举的所有质数
            st[primes[j]*i] = true;
            if(i % primes[j] == 0)    break;//若成立，primes[j]一定是i的最小质因子
        }
    }
}        

---

首先需要判断 i 是质数还是合数，如果是质数就将它放入primes[ ] 数组。

在第二个for循环中，将primes[ j ]的 i 倍筛掉（i 是从小到大依次遍历）

若存在一个合数X，设primes[ j ] 是它的最小质因子，当 i 枚举到X之前，一定会先枚举到X / primes[ j ]，而在那时，就已经先将X筛掉了。所以任何一个合数一定会被筛掉，是因为我们只用最小质因子来筛，但每个数只有一个最小质因子，所以每个数都只会筛一遍。例如 X = 12，2是12的最小质因子，在枚举到12之前一定会先枚举到6，并执行了st[ 2 * 6] = true，所以 i == 12 时，不会进入primes[ j ]数组  ；

最后的if语句有两种情况：

（1）若 i % primes[ j ] == 0 ，则说明 primes[ j ] 是 i 的最小质因子，那么primes[ j ] 也一定是primes[ j ] * i 的最小质因子。例如：4 % 2 == 0 ，2 是 4的最小质因子，且 2 也是2 * 4 的最小质因子。

（2）若i % primes[ j ] != 0 ，由于我们是从小到大枚举的所有的质数，并且我们没有枚举到 i 的任何一个质因子，则此时primes[ j ] 一定小于 i 的所有质因子，但是primes[ j ] 也一定是primes[ j ] * i 的最小质因子。

例如：9 % 2 ！= 0，2 一定小于 9 的最小质因子，但 2 一定是2 * 9 的最小质因子。（9的最小质因子是3，9要枚举到3才break，但9 * 2 的最小质因子是 2，当枚举到2时就可以break ）
 

附上题目代码：

#include <iostream>
#include <vector>

using namespace std;

int countPrimes(int n) {
    vector<bool> isPrime(n + 1, true);
    vector<int> primes;

    for (int i = 2; i <= n; ++i) {
        if (isPrime[i]) {
            primes.push_back(i);
        }
        for (int j = 0; j < primes.size() && i * primes[j] <= n; ++j) {
            isPrime[i * primes[j]] = false;
            if (i % primes[j] == 0) {
                break;
            }
        }
    }

    return primes.size();
}

int main() {
    int n;
    cin >> n;

    int result = countPrimes(n);
    cout << result << endl;

    return 0;
}