欧拉筛法求素数

C代码

---

#include <stdio.h>
#include <stdbool.h>

bool flag[10001] = {0};
int p[10001] = {0};

int main()
{
    int i, n, cnt = 0, j;
    scanf("%d", &n);

    for (i = 2; i <= n; ++i){
        if (flag[i] == 0)
            p[cnt++] = i;
        for (j = 0; j < cnt && i * p[j] <= n; ++j){
            flag[i * p[j]] = 1;
            if (i % p[j] == 0) break;
        }
    }

    for (i = 0; i < cnt; ++i)
        i == cnt - 1 ? printf("%d\n", p[i]) : printf("%d ", p[i]);
    return 0;
}

说明

---

• flag标记下标是否为合数
• p按顺序存放素数
• flag[i * p[j]] = 1筛掉i这一素数的素数倍数。
• 若i能整除p[j]则跳出循环，等到p[j]==i时再筛去i*p[j]这个数。即仅在合数的最大因子作为乘数时筛去这一合数。
• 时间复杂度O(n)

C++版本（实现上有所改动）

---

#include <bits/stdc++.h>
using namespace std;

const int maxn = 1e7+5;
bool np[maxn]{true,true};
vector<int> prime;

int main()
{
    int n, m, x;
    cin >> n >> m;
    for (int i = 2; i <= n; ++i)
    {
        if (!np[i]) prime.push_back(i);
        for (int j = 0; j < prime.size() && i*prime[j] <= n; ++j)
        {
            np[i*prime[j]] = true;
            if (i % prime[j] == 0) break;
        }
    }
    for (int i = 1; i <= m; ++i)
    {
        scanf("%d", &x);
        printf("%s\n", np[x] ? "No" : "Yes");
    }
    return 0;
}