筛质数(暴力法、埃氏筛、欧拉筛)

筛质数(暴力法、埃氏筛、欧拉筛)

暴力法

思路分析:

直接双for循环来求解质数

如果不设置标记只是简单地执行了break会导致内部循环(由j控制)而不是立即打印i或者跳过它。如果打印语句写到内部循环中,也会导致每个

非素数也被打印出来多次(在找到其因子之前)。

code:

cpp 复制代码
#include <iostream>
using namespace std;
int n;
int main() {
    cin >> n;
    for (int i = 2; i <= n; i++) { // 从2开始,因为1不是素数
        bool isprime = true; // 假设i是素数
        for (int j = 2; j * j <= i; j++) { // 只需要检查到sqrt(i),因为如果i有一个大于sqrt(i)的因子,那么它必然还有一个小于或等于sqrt(i)的因子
            if (i % j == 0) {
                isprime = false; // 如果找到一个因子,则i不是素数
                break;
            }
        }
        if (isprime) { // 如果i是素数,则打印它
            cout << i << " ";
        }
    }
    return 0;
}

埃氏筛

思路分析:

如果一个数是素数,那么它的倍数就不是素数

code:

cpp 复制代码
#include <iostream>
using namespace std;
int main(){
    int p[100010]={1,1};//f[i]=1代表不是素数
    int n;
    cin>>n;
    for(int i=2;i<=n;i++){
        if(p[i]==1) continue;
        for(int j=2;i*j<=n;j++){
            p[i*j]=1;
        }
    }
    for(int i=1;i<=n;i++){
        if(p[i]==0){
            cout<<i<<' ';
        }
    }
}

欧拉筛

思路分析:

基本思想:确保每个数只被其除1以外的最小质因数标记。

实现:<需要一个数组把它存储起来,

逆向思考:从最大因数开始找到最小的质因数,

在处理一个数时,需要遍历所有已经找到的质数,并当该质数是当前数的因数时退出循环,这是因为对于更大的数,当前数已不再是其最大的因数,无需继续处理。>

code:

cpp 复制代码
#include <iostream>
using namespace std;
int main(){
    int d=0;
    int p[10010]={0};
    int f[10010]={1,1};
    int n;
    cin>>n;
    for(int i=2;i<=n;i++){
        if(f[i]==0){//如果没有被标记过,那么i是质数
            p[d++]=i;//先运行后加
        }
        for(int j=0;j<d;j++){
            if(p[j]*i<=n){//标记以i为最大因数的数不是素数
                f[p[j]*i]=1;
            }
            else{
                break;
            }
            if(i%p[j]==0){//如果p[j]是i的因数,那么后面的数都不是以i为最大因数的
                break;
            }
        }
    }
    for(int i=0;i<d;i++){
        cout<<p[i]<<' ';
    }
}

如果 i 能够被 pj 整除(即 i % pj == 0),那么说明 pj 是 i 的一个因数。由于我们已经用 pj 标记了所有小于或等于n的 pj 的倍数,所以没有必要再用更大的质数去标记 i 的倍数。因此,可以跳出内层循环。

相关推荐
luj_176815 分钟前
残熵算法实时化三大瓶颈突破
c语言·开发语言·网络·经验分享·算法
mCell1 小时前
你以为短链接只是 Hash + 301/302?
后端·算法·架构
远离UE42 小时前
UE5 compute shader 原子加
开发语言·c++·ue5
C+-C资深大佬2 小时前
C++ 显式类型转换详解:static_cast、dynamic_cast、const_cast、reinterpret_cast
开发语言·c++
KaMeidebaby2 小时前
卡梅德生物技术快报|抗体亲和力成熟工业化调控新机制:差异性浆细胞增殖工艺优化思路
java·开发语言·人工智能·算法·机器学习·架构·spark
luj_17683 小时前
心形曲线轨迹控制三大关键技术
c语言·开发语言·c++·经验分享·算法
取地址符3 小时前
C++学习笔记(基于learn-cxx)(1)
c++·经验分享·笔记·学习
tmlx3I0813 小时前
高光谱拼接算法(六)RANSAC 误匹配剔除
人工智能·算法·机器学习
不相心 -w-4 小时前
基础-滑动窗口-(定长 | 不定长)
算法
思麟呀4 小时前
C++17(三)if constexpr+折叠表达式
开发语言·c++