埃拉托斯特尼筛法

简介:埃拉托斯特尼筛法(Sieve of Eratosthenes)是一种高效计算素数的算法,它能够找出小于等于某个整数 nnn 的所有素数。该算法的基本思想是通过不断筛选合数(非素数)来获取素数,从而避免了对每个数逐一检查是否为素数的过程。

例:

#include <iostream>
using namespace std;

// 计算 1 到 a 之间的素数个数
int Yanhongyan(int a)
{
    bool b[a + 1]; // 布尔数组,用来标记是否为素数
    for (int i = 0; i <= a; i++)
    {
        b[i] = true; // 假设所有数是素数
    }
    b[0] = b[1] = false; // 0 和 1 不是素数

    // 埃拉托斯特尼筛法
    for (int i = 2; i * i <= a; i++)
    {
        if (b[i])
        {
            for (int j = i * i; j <= a; j += i)
            {
                b[j] = false; // 标记 i 的倍数为非素数
            }
        }
    }

    // 统计素数的个数
    int sum = 0;
    for (int i = 2; i <= a; i++)
    {
        if (b[i])
        {
            sum++; // 如果是素数,计数加一
        }
    }

    return sum; // 返回素数个数
}

int main()
{
    int t;
    cin >> t; // 输入测试用例的个数

    // 对每个输入的值输出其对应的素数个数
    while (t--)
    {
        int n;
        cin >> n;                      // 输入每个测试用例的值
        cout << Yanhongyan(n) << endl; // 输出对应的素数个数
    }

    return 0;
}

 

输入为:4
              5
              10
              20
              100

运行结果:

算法步骤:

  1. 初始化:创建一个布尔数组 isPrime,其中 isPrime[i] 表示数字 i 是否为素数。将数组中所有元素初始化为 true,并将 isPrime[0]isPrime[1] 设置为 false(因为 0 和 1 不是素数)。
  2. 筛选过程
    • 2 开始,检查每个数字。如果该数字是素数,则将它的所有倍数标记为非素数(false)。
    • 2 开始标记所有倍数为 false,继续到平方根范围的数字。因为对于任何一个大于 sqrt(n) 的素数,它的倍数都会在之前被处理过。
  3. 输出结果
    • 最终,isPrime[i]true 的位置表示 i 是素数。
### 埃拉托斯特尼筛法实现素数查找 埃拉托斯特尼筛法是一种高效的算法,用于生成指定范围内的素数列表。其核心思想是通过标记合数的方式逐步筛选出素数[^1]。以下是该算法的详细描述及其实现示例。 #### 算法描述 埃拉托斯特尼筛法的基本步骤如下: 1. 创建一个布尔数组 `is_prime`,长度为 `n+1`,初始值全部设为 `true`。 2. 将数组中的索引 `0` 和 `1` 的值设置为 `false`,因为它们不是素数。 3. 从第一个素数 `2` 开始,将所有 `2` 的倍数标记为 `false`(即非素数)。 4. 找到下一个未被标记为 `false` 的数字,将其视为素数,并将其所有倍数标记为 `false`。 5. 重复上述过程,直到处理完不大于根号 `n` 的所有数字。 6. 最终,数组中值为 `true` 的索引即为素数。 #### Python 实现示例 以下是一个使用 Python 实现埃拉托斯特尼筛法的代码示例: ```python def sieve_of_eratosthenes(n): is_prime = [True] * (n + 1) is_prime[0] = is_prime[1] = False # 0 和 1 不是素数 p = 2 while p * p <= n: if is_prime[p]: for i in range(p * p, n + 1, p): is_prime[i] = False p += 1 primes = [i for i in range(n + 1) if is_prime[i]] return primes # 示例:查找 50 以内的所有素数 print(sieve_of_eratosthenes(50)) ``` #### C# 实现示例 以下是一个使用 C# 实现埃拉托斯特尼筛法的代码示例[^2]: ```csharp using System; class Program { static void Main() { int max = 50; bool[] isPrime = new bool[max + 1]; for (int i = 0; i <= max; i++) isPrime[i] = true; isPrime[0] = isPrime[1] = false; for (int p = 2; p * p <= max; p++) { if (isPrime[p]) { for (int multiple = p * p; multiple <= max; multiple += p) { isPrime[multiple] = false; } } } Console.WriteLine("Primes up to " + max + ":"); for (int i = 2; i <= max; i++) { if (isPrime[i]) Console.Write(i + " "); } } } ``` #### Java 实现示例 以下是一个使用 Java 实现埃拉托斯特尼筛法的代码示例[^4]: ```java public class Eratosthenes { public static void main(String[] args) { int max = 50; boolean[] isPrime = new boolean[max + 1]; for (int i = 0; i <= max; i++) isPrime[i] = true; isPrime[0] = isPrime[1] = false; for (int p = 2; p * p <= max; p++) { if (isPrime[p]) { for (int multiple = p * p; multiple <= max; multiple += p) { isPrime[multiple] = false; } } } System.out.println("Primes up to " + max + ":"); for (int i = 2; i <= max; i++) { if (isPrime[i]) System.out.print(i + " "); } } } ``` ### 时间复杂度与空间复杂度 埃拉托斯特尼筛法的时间复杂度为 \(O(n \log(\log n))\),其中 \(n\) 是给定范围的上限[^1]。空间复杂度为 \(O(n)\),主要用于存储布尔数组 `is_prime`。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值