筛质数(2种筛法)

最新推荐文章于 2024-10-05 22:19:52 发布

AuroraO.o

最新推荐文章于 2024-10-05 22:19:52 发布

阅读量120

点赞数 1

文章标签：算法 c++

本文链接：https://blog.csdn.net/m0_62387435/article/details/124786866

版权

给定一个正整数 n，请你求出 1∼n 中质数的个数。

输入格式
共一行，包含整数 n。

输出格式
共一行，包含一个整数，表示 1∼n 中质数的个数。

数据范围
1≤n≤106
输入样例：
8
输出样例：
4
时间复杂度
埃氏筛法 O(N * log(log(N)) )
线性筛法 O(N)

埃氏筛法

int primes[1000010], cnt;
bool st[1000010];
void get_primes(int n)
{
    int i,j,v=1,cnt=0;
   for(i=2;i<=n;i++)
   {
       if(!st[i])primes[cnt++]=i;
       else continue;
       for(j=i;j<=n/i;j++)
       {
           st[j*i]=true;
       }
   }
 }

线性筛法

const int N = 100010;
int v[N], primes[N], cnt;

void get_primes(int n)
{
    for (int i = 2; i <= n; i++)                 
    {
        if (v[i] == 0){
            v[i] = i;            // i质数 
            primes[cnt++] = i;
        }
        for (int j = 0; j < cnt; j++){               // i 的 最小质因子是v[i]  primes[j] > v[i];
            if (primes[j] > v[i] || primes[j] > n / i) 
            break;  //i有比primes[j]更小的质因子 或者超过了n的范围
            v[i * primes[j]] = primes[j];    // v[i * m] = min(v[i], v[m])  
             //i : v[i]  primes[j] : primes[j];
        }
    }
}

const int N = 100010;
int primes[N], cnt;
bool st[N];

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;
        }
    }
}