[蓝桥杯]质数数目

质数数目

题目描述

实现一个算法求质数数目。介绍如下：

• 质数是指在大于 1 的自然数中，除了 1 和它本身以外不再有其他因数的自然数。
• 给定数字 nn，对于从 0 至 nn 的数字，需要判断这个数字是不是质数。然后输出质数总数。

输入描述

输入一个数字 n (1≤N≤105)n (1≤N≤105)，含义见题干。

输出描述

输出一行，为 0 至 nn 之间质数的个数。

输入输出样例

示例

输入

5

输出

3

运行限制

• 最大运行时间：1s
• 最大运行内存: 256M

总通过次数: 3378  |  总提交次数: 3783  |  通过率: 89.3%

难度: 困难   标签: 数学, 质数

C++质数数目计算：埃拉托斯特尼筛法实现

📝 算法思路（埃拉托斯特尼筛法）

1. ​​初始化标记数组​​：创建大小为n+1的布尔数组isPrime，初始化所有元素为true（默认都是质数）
2. ​​标记非质数​​：
   • 0和1不是质数 → 标记为false
   • 从i=2开始遍历到√n：
     • 若isPrime[i]为true，则将所有i的倍数标记为false
     • ​​优化​​：从i*i开始标记（更小的倍数已被处理）
3. ​​统计质数​​：遍历数组，统计isPrime[i] == true的数量

graph TD
    A[开始] --> B[初始化isPrime数组]
    B --> C[标记0和1为false]
    C --> D[i从2到√n循环]
    D --> E{isPrime[i]为真?}
    E --是--> F[j从i*i到n, 步长i]
    F --> G[标记isPrime[j]=false]
    E --否--> D
    G --> D
    D --> H[统计true数量]
    H --> I[输出结果]

⚙️ C++代码实现

#include <iostream>
#include <vector>
#include <cmath>
using namespace std;

int countPrimes(int n) {
    if (n < 2) return 0;  // 0和1无质数
    
    vector<bool> isPrime(n + 1, true);  // 初始化标记数组
    isPrime[0] = isPrime[1] = false;   // 标记非质数
    
    // 埃氏筛核心算法
    for (int i = 2; i * i <= n; i++) {
        if (isPrime[i]) {
            // 从i*i开始标记i的倍数
            for (int j = i * i; j <= n; j += i) {
                isPrime[j] = false;
            }
        }
    }
    
    // 统计质数数量
    int count = 0;
    for (int i = 2; i <= n; i++) {
        if (isPrime[i]) count++;
    }
    return count;
}

int main() {
    int n;
    cin >> n;
    cout << countPrimes(n) << endl;
    return 0;
}

🔍 代码解析

1. ​​时间复杂度​​：O(n log log n)

   • 内层循环只对质数执行，且从i*i开始跳跃标记
   • 10⁵数据量下仅需约0.01秒

2. ​​空间复杂度​​：O(n)

   • 使用vector<bool>特化优化（1bit/元素）
   • 10⁵数据仅需12.5KB内存

3. ​​关键优化点​​：

   • 外层循环终止于√n（i*i <= n）
   • 跳过偶数处理（但本实现未显式优化）
   • 避免重复标记（从i*i开始）

🧪 实例验证

输入: 5 → 质数: [2, 3, 5] → 输出: 3
输入: 10 → 质数: [2, 3, 5, 7] → 输出: 4
输入: 100 → 输出: 25（验证通过）[1](@ref)

⚠️ 注意事项

1. ​​边界处理​​：

   • n < 2时直接返回0
   • 数组大小设为n+1避免越界

2. ​​特殊测试点​​：

   输入	预期输出	说明
   0	0	空范围
   1	0	无质数
   2	1	最小质数
   100000	9592	最大边界值 2

3. ​​常见错误​​：

   • 未初始化isPrime[0]和isPrime[1]
   • 外层循环终止条件错误（应i*i <= n非i <= n）
   • 内层循环从2*i开始（导致重复标记）

🔧 优化建议

1. ​​空间优化​​：

   // 奇偶压缩存储（节省50%空间）
   vector<bool> isPrime((n+1)/2, true);  
   // 映射关系：数字i对应索引i/2（i为奇数）

2. ​​时间优化​​：

   // 单独处理偶数
   if (n >= 2) count++;  
   for (int i = 3; i <= n; i += 2) {  // 仅检查奇数
       if (isPrime[i/2]) count++;
   }

3. ​​分块筛法​​：

   • 适用于n > 10⁷的超大范围
   • 将区间分成小块，逐块筛质数

性能对比（10⁵数据量）