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

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

暴力法

思路分析:

直接双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 的倍数。因此,可以跳出内层循环。

相关推荐
坤坤子吖13 分钟前
C++继承:友元、静态成员、多继承与组合
开发语言·c++·笔记·学习·1024程序员节
qz_Serene34 分钟前
练习算法题:双指针 -cpp
c++·算法·leetcode·1024程序员节
扶尔魔ocy1 小时前
【C++开发】从裸指针到智能指针
c++·智能指针·1024程序员节
Nuanyt1 小时前
Hot 100 刷题笔记 · 数组与哈希专题(3 题)
java·笔记·算法·leetcode·哈希算法·1024程序员节
碧海银沙音频科技研究院1 小时前
基于炬芯平台ATS308X多模型兼容语音识别算法
人工智能·算法·语音识别
爱读书的小胖2 小时前
LeetCode 274 H 指数 C 语言解题分析
c语言·算法·leetcode·1024程序员节
basketball6162 小时前
动态规划第一章:引论
算法·动态规划·1024程序员节
星夜离尘562 小时前
MATLAB 因果采样保持:用历史匹配和数据年龄防止未来信息泄漏
数据结构·算法·matlab
思麟呀3 小时前
BTree的原理与 C++ 完整实现
c语言·数据结构·c++·算法·1024程序员节
jimy13 小时前
C++ 中 xvalue 的 identity 存在哪里?和 lvalue 的 identity 有什么区别?
开发语言·c++·1024程序员节