筛法——求n以内质数个数

最新推荐文章于 2024-01-07 20:44:29 发布

OGMWA

最新推荐文章于 2024-01-07 20:44:29 发布

阅读量3.4k

点赞数 4

文章标签：算法

本文链接：https://blog.csdn.net/qq_30603195/article/details/86677762

版权

博客介绍了不同筛法求质数的原理和复杂度，重点讲解了欧氏筛法，指出其复杂度为O(n)且避免了重复筛选的问题。还探讨了判断质数的优化算法，并分享了在实际编程中遇到的问题和解决经验，包括输出效率和数组大小的设定。

摘要由CSDN通过智能技术生成

文章目录

预备知识：
判断n是否为素数
- 朴素筛
- 优化：
续筛法部分
区间筛法

预备知识：

若一个数可以进行因数分解，则得到的两个数一定是有一个≤ $\sqrt x$ ，另一个≥ $\sqrt x$
任何一个数都可以拆分为若干个质因数的乘积，如6=2*3，7=1*7 。

朴素筛

复杂度O(n $\sqrt n$ )

#include<iostream>
using namespace std;
int plainSieve(int n)
{
   
    int p=0;
    if(n==2)
        return 1;
    if(n<2)
        return 0;
    p++;
    // i+=2优化步长
    for(int i=3;i<=n;i+=2)
    {
   
        int j;
        // j*j优化范围
        for(j=2;j*j<=i;j++)
            if(i%j==0)
                break;
        if(j*j>i)
            p++;
    }
    return p;
}
int main()
{
   
    int n;
    while(cin>>n)
        cout<<plainSieve(n)<<endl;
    return 0;
}

埃氏筛

复杂度 O(nlglgn)

解释：当prime[i]为0时，说明 i 没有被筛到过，即 i 为质数。质数总数 p++
之后要将 i 的n以内的倍数 j 全部置为合数，prime[j]标记为1。
此后若筛到 j ，prime为1，j 为合数，直接跳过，进行下个数的判断

#define Max 1000010
int prime[Max]={
   0};
int EraSieve(int n)
{
   
    if(n<2)
        return 0;
    if(n==2)
        return 1;
    int p=0;
    prime[0]=prime[1]=1;
    for(int i=2;i<=n;i++)
    {
   
        if(prime[i]==0)
        {
   
            p++;
            for(int j=i*2;j<=n;j+=i)
                prime[j]=1;
        }
    }
    return p;
}