半秒求一亿内的素数：经典筛选法&线性筛选法的介绍、改进、实现与性能分析

最新推荐文章于 2024-01-21 14:56:32 发布

orange_wolf

最新推荐文章于 2024-01-21 14:56:32 发布

阅读量2.5k

点赞数 4

分类专栏： C++ 算法

本文链接：https://blog.csdn.net/my_precious/article/details/82789144

版权

本文介绍了经典筛选法和线性筛选法的原理与C++实现，探讨了两种筛选法的改进策略，减少了重复计算和内存消耗。通过性能对比，展示在Intel Core i3 8100平台上，改进后的筛选法运行时间减少，并实现了内存占用的减半。寄存器优化进一步提升了20%的性能。

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

埃拉托斯特尼筛法，简称埃氏筛，也有人称素数筛。这是一种简单且历史悠久的筛法，用来找出一定范围内所有的素数。
本文主要内容为四种筛选法的原理介绍和具体实现（C++），并对比它们之间的性能差距。这四种筛选法分别是：经典筛选法、采用了欧拉函数的线性筛选法，以及个人对这两种筛选法的改进。

经典筛选法

原理简介

关于经典筛选法，维基百科上的动图能很好地解释其原理。先用2去筛，即把2留下，把2的倍数剔除掉；再用下一个素数，也就是3筛，把3留下，把3的倍数剔除掉；接下去用下一个素数5筛，把5留下，把5的倍数剔除掉；不断重复下去…

具体实现

#define MAX 100000000  // 求一亿范围内的素数
long prime_numbers[6000000], n = 0;  // 一亿内有500多万个素数

void eratosthenes() {
   
    static bool not_prime_flag[MAX];  // 假设所有数都为素数
    for (long num = 2; num < MAX; num++) {
   
        if (!not_prime_flag[num]) {
   
            prime_numbers[n++] = num;
            for (long times = 2; times * num < MAX; times += 1) {
   
                not_prime_flag[times * num] = true;  // 将素数的倍数剔除掉
            }
        }
    }
}<