讲师中心 微信公众号
AI工具推荐 视频效率加速

C++实现快速开平方根算法 _ 牛顿迭代法底层逻辑与代码【源码】

冬磊姑娘_8723

冬磊姑娘_8723

发布时间:2026-04-27 08:45:18

|

949人浏览过

|

来源于php中文网

原创

牛顿迭代法整数开方需防溢出、用long long、整数终止条件及边界特判;std::sqrt转int出错因double精度不足、浮点误差偏移及环境限制。

c++实现快速开平方根算法 _ 牛顿迭代法底层逻辑与代码【源码】

牛顿迭代法不是“套公式就行”的黑盒——它快,但错一步就溢出、死循环或返回错误整数;整数开方必须用 long long 中间计算、整数终止条件、边界特判,三者缺一不可。

为什么直接用 std::sqrt 转 int 会出错

比如对 x = 9223372030926249000LL(接近 LLONG_MAX),static_cast<int>(std::sqrt(x))</int> 可能返回 3037000499,而正确向下取整结果是 3037000498。原因有三:

  • double 只有 53 位有效精度,该数共约 63 位,高位截断后 sqrt 输入已失真
  • 浮点 sqrt 返回值可能略大于真实根(如 3037000498.9999995 → 向下取整仍得 3037000498,但若误差朝上偏移就翻车)
  • 嵌入式环境或 freestanding C++ 实现中,std::sqrt 可能根本不可用

long long 版牛顿迭代:防溢出与整数收敛的关键写法

核心不是“套 (res + x / res) / 2”,而是控制每一步不越界、不停错。以下写法经实测覆盖 [0, LLONG_MAX] 全范围:

  • 初始值用 res = x > 1 ? x >> 1 : x,避免 x == 0 或 x == 1 时除零或无效迭代
  • 所有中间计算用 long long,尤其 res * res 和 x / res 必须在 long long 下算,否则 int 溢出直接 UB
  • 停止条件不用 res != last 或浮点差值,改用:if (res > x / res) { res--; } 然后检查 res * res x —— 这才是整数平方根的定义
  • 迭代中若 res == 0,立即 break,防止后续 x / res 崩溃

浮点牛顿迭代的精度陷阱与安全退出条件

当目标是高精度浮点结果(如 double)时,常见错误是用 fabs(x1 - x0) 判定收敛。问题在于:

C++ Code Review Master
C++ Code Review Master

组合式C++代码评审方案,融合静态分析、AI推理、多轮迭代评审和C++专项检查,适用于PR审查、增量代码审查、全项目评审和代码质量评分,触发词包括review cpp、cpp代码评审、C++review、代码审查。

下载

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

  • 早期迭代差值大,后期可能卡在 eps 之外反复震荡(尤其 x 接近 0 或极大时)
  • 更稳的方式是限定最大迭代次数(如 10 次),同时辅以相对误差:fabs(x1 - x0) / fmax(1.0, fabs(x1))
  • 初始值别硬写 1.0:对 x >= 1,用 x * 0.5;对 x ,用 <code>x + 0.5,能减少 1–2 次迭代
  • 务必加 if (x ::quiet_NaN();,不处理负输入是多数手写 sqrt 的 crash 根源

整数开方 vs 浮点开方:选算法前先看需求场景

二者数学目标不同,强行混用必踩坑:

  • 要的是“≤ √x 的最大整数”(如二分查找索引、内存页对齐、算法竞赛),必须用整数牛顿或位移搜索,double 转换不可信
  • 要的是“满足 |y² − x|
  • 嵌入式无 FPU?放弃浮点牛顿,直接上位移搜索法(纯 >> 和加减),O(log n) 且零分支预测失败风险
  • 性能敏感场景(如高频数值积分内层),整数牛顿平均 6 次迭代,比二分(~63 次)快一个数量级,但要注意编译器是否把 / 优化成倒数乘法

真正难的不是写出能跑的版本,而是让 x = 0、x = 1、x = LLONG_MAX、x = LLONG_MAX - 1 全部通过,且不依赖未定义行为——这些边界才是检验是否吃透牛顿迭代的试金石。

热门AI工具

更多
蛙蛙写作

一款AI论文写作工具,主要用于超级AI智能写作助手,适合需要提升相关任务效率的用户。

立刻MV
立刻MV Hot

立刻MV是一款AI文本写作工具,AI 音乐视频(MV)创作工具。

AionClaw
AionClaw Hot

AionClaw是一款面向办公、创作和编程任务的AI桌面智能体。

咔片AIPPT

一款在线AI演示文稿制作工具,可根据主题和内容需求辅助生成PPT结构与页面,提高演示材料制作效率。

LibLibAI
LibLibAI Hot

一款AI视频创作工具,主要用于国内领先的AI创意平台,以海量模型、低门槛操作与“创作-分享-商业化”生态,让小白与专业创作者都能高效实现图文乃至视频创意表达,适合需要提升相关任务效率的用户。

WorkBuddy

一款AI办公效率工具,主要用于腾讯云推出的AI原生桌面智能体工作台,适合需要提升相关任务效率的用户。

DeepSeek

DeepSeek是一款面向对话、写作、编程和推理场景的AI大模型工具。

豆包大模型

豆包大模型是一款由字节跳动推出的企业级大语言模型服务平台。

PixTV
PixTV Hot

PixTV是一款面向AIGC内容创作的AI视频生成工具。

相关专题

更多
c++和c语言的区别有哪些
c++和c语言的区别有哪些

c++和c语言的区别:1、面向对象编程(OOP)支持不同;2、新增特性不同;3、标准库不同;4、编译方式不同;5、命名空间不同等等。想了解更多c++的相关内容,可以阅读本专题下面的文章。

2208

2024.03.14

c++和python学习顺序推荐
c++和python学习顺序推荐

一般建议先学习C++,再学习Python,因为这样可以逐步从较为底层的编程语言向更高级的语言过渡。想了解更多python的相关内容,可以阅读本专题下面的文章。

999

2024.03.14

python和c++学习性价比分析
python和c++学习性价比分析

Python易于学习,广泛应用于Web开发、数据科学和人工智能等领域,但性能较低。C语言性能高,适用于对性能要求较高的场景,如游戏开发和系统编程,但学习曲线陡峭,错误处理复杂。想了解更多python的相关内容,可以阅读本专题下面的文章。

407

2024.03.14

c语言和c++一样吗
c语言和c++一样吗

c语言和c++是两种不同的编程语言,虽然有相似之处,但存在显著差异。c语言专注于过程式编程和系统级开发,以简洁、高效著称。c++作为c语言的超集,引入了面向对象编程,增强了代码组织和管理能力,但学习曲线也更陡峭。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

307

2024.03.14

c语言和c++先学哪个好
c语言和c++先学哪个好

初学者选择学习c语言还是c++语言,需要根据个人学习目标、背景以及编程兴趣和预期应用方向来决定。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

386

2024.03.14

c语言和c++的区别和联系
c语言和c++的区别和联系

c语言和c++是计算机科学领域应用广泛的编程语言。虽然它们有着相似的基础,但它们在语言类型、语法功能和内存管理方面存在着显著差异。本专题为大家提供相关的文章、下载、课程内容,供大家免费下载体验。

580

2024.03.14

c++软件中文更改教程
c++软件中文更改教程

对于 ide,可通过打开设置,找到语言设置,选择中文,并保存更改。对于非 ide 应用程序,可查找设置或选项,选择语言设置,更改为中文,并保存更改。想了解更多c++的相关内容,可以阅读本专题下面的文章。

1389

2024.03.21

python和java和c++学习性价比分析
python和java和c++学习性价比分析

Python以其易学性、丰富的库和活跃的社区而著称,适合数据科学、人工智能和Web开发。Java以其跨平台性、企业级应用开发和Android应用开发而闻名。C++以其底层控制能力、高效性能和游戏开发而著称。选择哪种语言取决于个人兴趣、职业方向和特定需求。想了解更多python和java和c++的相关内容,可以阅读本专题下面的文章。

1197

2024.03.22

LLVM自定义Pass怎么写
LLVM自定义Pass怎么写

本专题聚焦LLVM自定义Pass开发,整理Pass类结构、run()方法、PreservedAnalyses、CMake构建、插件注册、-load-pass-plugin加载和测试用例编写流程。

100

2026.09.30

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Conan 2 Essentials 免费课程
Conan 2 Essentials 免费课程

共0课时 | 0人学习

CMake 与 Conan 集成实践
CMake 与 Conan 集成实践

共0课时 | 0人学习

Conan 2 高级依赖模型介绍
Conan 2 高级依赖模型介绍

共0课时 | 0人学习

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

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn