0

0

C++怎么生成素数 C++求质数的高效算法实现【精析】

裘德小鎮的故事

裘德小鎮的故事

发布时间:2026-03-15 13:18:13

|

912人浏览过

|

来源于php中文网

原创

埃氏筛配合std::vector预处理10⁶内素数最实用,兼顾简洁性、可读性与性能;单个大数素性判断用Miller-Rabin(配7个确定性底数);查素数表优先用std::lower_bound。

c++怎么生成素数 c++求质数的高效算法实现【精析】

std::vector + 埃氏筛预处理 10⁶ 以内素数最实用

绝大多数 C++ 实际场景(比如 OJ 题、算法预处理、小范围质数枚举)不需要实时判断单个大数是否为素数,而是要快速获取一段连续范围内的所有素数。这时候埃拉托斯特尼筛法(埃氏筛)配合 std::vector<bool></bool> 是平衡简洁性、可读性和性能的首选。

常见错误是手写 is_prime(int n) 对每个数试除——10⁵ 个查询就可能超时;或者用 vector<int></int> 存筛出的素数但忘了初始化大小,导致反复 push_back 引发内存重分配开销。

实操建议:

  • 筛前先 vector<bool> is_prime(n + 1, true)</bool>,注意索引 0 和 1 要设为 false
  • 外层循环只需到 sqrt(n),且从 i * i 开始标记(i * 2 可能已被更小的素数筛过)
  • 如果只要素数列表,最后遍历一次 is_prime 收集下标即可,别边筛边存——缓存不友好
  • vector<bool></bool> 是特化类型,空间省但访问略慢;若对性能极端敏感且内存充足,可用 vector<char></char>
int n = 1000000;
vector<bool> is_prime(n + 1, true);
is_prime[0] = is_prime[1] = false;
for (int i = 2; i * i <= n; ++i) {
    if (is_prime[i]) {
        for (int j = i * i; j <= n; j += i) {
            is_prime[j] = false;
        }
    }
}

判断单个大数是否为素数:用 Miller-Rabin 比试除快得多

当输入是单个 64 位整数(比如 long long),且值可能接近 10¹⁸ 时,试除法最多要跑 √n ≈ 10⁹ 次,肯定超时。Miller-Rabin 是概率性算法,但对 64 位整数,用固定底数集(如 {2, 325, 9375, 28178, 450775, 9780504, 1795265022})可做到确定性正确。

立即学习C++免费学习笔记(深入)”;

容易踩的坑:

  • 快速幂中乘法溢出:必须用 __int128(GCC)或自写模乘,不能直接 (a * b) % mod
  • n == 2n < 2 的边界漏掉,导致对 1、2、偶数处理出错
  • 误以为“通过一轮测试 = 是素数”,其实需对多个底数测试;不过对 uint64_t,上述 7 个底数已覆盖全部

使用场景:密码学参数生成、大数分解前置判断、交互题中动态质数验证。

ChatDOC
ChatDOC

ChatDOC是一款基于chatgpt的文件阅读助手,可以快速从pdf中提取、定位和总结信息

下载

std::lower_bound 查素数表比线性扫快得多

如果你已经用埃氏筛生成了 vector<int> primes</int>(升序),后续要频繁查“≤x 的最大素数”或“第 k 小素数”,千万别用 for 循环遍历——O(n) 太慢。

实操建议:

  • 查第一个 ≥ x 的素数位置:用 lower_bound(primes.begin(), primes.end(), x)
  • 查最后一个 ≤ x 的素数:用 upper_bound 然后减一,注意判越界
  • 确保 primes 是全局或静态构造,避免重复筛;C++17 可用 inline constexpr vector(需编译器支持)
  • 若只查几十次,预处理成本可能得不偿失;但查上百次,O(log π(n)) ≈ O(log n) 就明显占优

Windows 下 clang++ 编译 Miller-Rabin 可能缺 __int128

Clang for Windows(尤其 MinGW-w64 工具链)默认不启用 __int128,而 Miller-Rabin 的模乘又极度依赖它。此时会编译报错:error: use of undeclared identifier '__int128',或运行时乘法溢出导致误判。

解决办法很具体:

  • 加编译选项 -D__SIZEOF_INT128__ 强制定义(仅限 GCC/Clang 兼容环境)
  • 改用三模乘(mulmod(a, b, m) 分三段模拟),但代码变长、常数变大
  • 换用 MSVC?不行,它原生不支持 __int128,且无对应内置函数
  • 最稳方案:Linux/macOS 下开发,或用 WSL;生产环境部署也建议避开 Windows 原生 Clang 做大数素性检测

这个点非常容易被忽略——算法逻辑全对,只因平台差异,is_prime(1000000000000000000LL) 返回 false。

热门AI工具

更多
DeepSeek
DeepSeek

幻方量化公司旗下的开源大模型平台

豆包大模型
豆包大模型

字节跳动自主研发的一系列大型语言模型

WorkBuddy
WorkBuddy

腾讯云推出的AI原生桌面智能体工作台

腾讯元宝
腾讯元宝

腾讯混元平台推出的AI助手

文心一言
文心一言

文心一言是百度开发的AI聊天机器人,通过对话可以生成各种形式的内容。

讯飞写作
讯飞写作

基于讯飞星火大模型的AI写作工具,可以快速生成新闻稿件、品宣文案、工作总结、心得体会等各种文文稿

即梦AI
即梦AI

一站式AI创作平台,免费AI图片和视频生成。

ChatGPT
ChatGPT

最最强大的AI聊天机器人程序,ChatGPT不单是聊天机器人,还能进行撰写邮件、视频脚本、文案、翻译、代码等任务。

相关专题

更多
scripterror怎么解决
scripterror怎么解决

scripterror的解决办法有检查语法、文件路径、检查网络连接、浏览器兼容性、使用try-catch语句、使用开发者工具进行调试、更新浏览器和JavaScript库或寻求专业帮助等。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

493

2023.10.18

500error怎么解决
500error怎么解决

500error的解决办法有检查服务器日志、检查代码、检查服务器配置、更新软件版本、重新启动服务、调试代码和寻求帮助等。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

382

2023.10.25

string转int
string转int

在编程中,我们经常会遇到需要将字符串(str)转换为整数(int)的情况。这可能是因为我们需要对字符串进行数值计算,或者需要将用户输入的字符串转换为整数进行处理。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

1051

2023.08.02

int占多少字节
int占多少字节

int占4个字节,意味着一个int变量可以存储范围在-2,147,483,648到2,147,483,647之间的整数值,在某些情况下也可能是2个字节或8个字节,int是一种常用的数据类型,用于表示整数,需要根据具体情况选择合适的数据类型,以确保程序的正确性和性能。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

616

2024.08.29

c++怎么把double转成int
c++怎么把double转成int

本专题整合了 c++ double相关教程,阅读专题下面的文章了解更多详细内容。

335

2025.08.29

C++中int的含义
C++中int的含义

本专题整合了C++中int相关内容,阅读专题下面的文章了解更多详细内容。

235

2025.08.29

windows查看端口占用情况
windows查看端口占用情况

Windows端口可以认为是计算机与外界通讯交流的出入口。逻辑意义上的端口一般是指TCP/IP协议中的端口,端口号的范围从0到65535,比如用于浏览网页服务的80端口,用于FTP服务的21端口等等。怎么查看windows端口占用情况呢?php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

1518

2023.07.26

查看端口占用情况windows
查看端口占用情况windows

端口占用是指与端口关联的软件占用端口而使得其他应用程序无法使用这些端口,端口占用问题是计算机系统编程领域的一个常见问题,端口占用的根本原因可能是操作系统的一些错误,服务器也可能会出现端口占用问题。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

1172

2023.07.27

TypeScript类型系统进阶与大型前端项目实践
TypeScript类型系统进阶与大型前端项目实践

本专题围绕 TypeScript 在大型前端项目中的应用展开,深入讲解类型系统设计与工程化开发方法。内容包括泛型与高级类型、类型推断机制、声明文件编写、模块化结构设计以及代码规范管理。通过真实项目案例分析,帮助开发者构建类型安全、结构清晰、易维护的前端工程体系,提高团队协作效率与代码质量。

69

2026.03.13

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
C# 教程
C# 教程

共94课时 | 11.3万人学习

C 教程
C 教程

共75课时 | 5.4万人学习

C++教程
C++教程

共115课时 | 21.9万人学习

关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号 技术交流群
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn | 湘ICP备2023035733号