高精度计算入门与immense数据模型精讲（七）

跳转链接高精度计算入门与immense数据模型精讲一高精度计算入门与immense数据模型精讲二高精度计算入门与immense数据模型精讲三高精度计算入门与immense数据模型精讲四高精度计算入门与immense数据模型精讲五高精度计算入门与immense数据模型精讲六高精度计算入门与immense数据模型精讲八高精度计算入门与immense数据模型精讲九高精度计算入门与immense数据模型精讲十高精度计算入门与immense数据模型精讲收官immense大实数数据库前言四则运算之上的整数-实数桥梁上一篇我们完整攻克了高精度除法完成了加减乘除四则运算的收官。但仅有四则运算还不够——很多高级算法都需要在「整数」和「实数」之间切换快速幂需要分离整数指数与小数指数三角函数需要对角度做周期约化级数迭代需要判断余项是否达到精度阈值。这些操作的底层都依赖取整函数与取模运算。取整看起来很简单不就是丢掉小数部分吗但真正严谨的取整有四种标准语义截断取整、地板取整、天花板取整、四舍五入。它们在正数下行为相似到了负数场景却截然不同很多初学者都在这里踩过坑。而取模运算更是和取整深度绑定——商的取整方式直接决定了余数的符号和大小。immense对取整的实现非常有巧思它没有针对每一位小数写复杂的截断进位逻辑而是借助Product_10小数点平移函数用「移位 → 取整 → 复位」的三步通用框架统一实现了四种取整还支持任意指定位数包括十位、百位等整数高位的取整。本篇依然遵循「通用原理先行代码实现对照」的教学节奏先把四种标准取整的数学定义、区别与适用场景讲透再讲高精度取整的通用移位法、取模运算的数学定义最后对照immense的源码逐函数拆解取整与取模的完整实现。一、通用原理取整与取模的核心方法论1.1 为什么需要标准化取整很多人对取整的理解停留在「去掉小数部分」但这只是其中一种取整方式。在数学和工程领域取整有四种严格定义的标准语义分别对应不同的应用场景取整方式英文数学定义核心方向截断取整trunc向零取整直接丢弃小数部分靠近0地板取整floor不大于原数的最大整数向负无穷天花板取整ceil不小于原数的最小整数向正无穷四舍五入round距离最近的整数0.5向绝对值大的方向入就近取整正数场景下trunc、floor、ceil 的区别还不明显一旦到了负数四种取整的结果会截然不同。我们用一组经典对比案例建立直观认知原值trunc向零floor向负无穷ceil向正无穷round四舍五入1.711221.21121-1.2-1-2-1-1-1.7-1-2-1-2教学要点最容易混淆的是 trunc 和 floor。正数下两者完全一致负数下完全相反。记住一句话floor 永远往下走trunc 永远往零走。比如 -1.2 往下走是 -2往零走是 -1。1.2 四种取整的典型应用场景每种取整都有各自的适用领域用错了会导致逻辑错误trunc最符合直觉的「砍掉小数」常用于整数除法、取模运算、小数位数截断floor数学分析、周期函数约化、向下取整计数是很多数学公式的默认取整ceil向上取整计数比如分页计算、资源分配、容量预估round日常数值展示、统计、财务近似是最符合人类习惯的取整方式1.3 高精度取整的通用方法移位取整法如果要实现「保留小数点后 Point 位」的取整直接操作小数部分的第 Point 位会比较麻烦还要处理 Point 为负数保留到十位、百位的情况。工程上有一种非常优雅的通用方案叫小数点移位法核心思想是把取整目标位移到整数部分的最低位这样取整就变成了简单的「丢掉所有小数」移位还原后就是最终结果。通用三步框架对于任意指定位数 Point 的取整移位把小数点向右移动 Point 位让目标保留位变成整数部分的最低位Point0保留小数位小数点右移小数变整数Point0保留整数高位小数点左移高位变低位取整对移位后的数执行对应语义的取整截断、地板、天花板、四舍五入复位把小数点向左移动 Point 位还原到原来的数量级举个直观例子对123.456保留 2 位小数做 floor小数点右移2位 →12345.6地板取整 →12345小数点左移2位 →123.45再比如对12345.6保留到百位Point-2做 ceil小数点左移2位 →123.456天花板取整 →124小数点右移2位 →12400这个方法的最大优势是统一框架无论保留小数位还是整数高位无论哪种取整语义外层移位逻辑完全相同中间只需要替换一步取整核心即可。几乎所有成熟的高精度实数库都采用这种实现方式。1.4 取模运算的数学定义取模求余和取整是深度绑定的。数学上取模运算的通用公式是a mod b a - b × q其中 q a / b 的商q 的取整方式决定了取模的定义不同的取整方式对应不同的取模规则截断取模trunc 模q 取 trunc(a/b)余数符号与被除数相同C/C 等绝大多数语言的整数取模都采用这种地板取模floor 模q 取 floor(a/b)余数符号与除数相同Python 等语言采用这种对于正数两种模的结果完全一致对于负数结果会有差异。工程上实现高精度取模通常会和语言原生行为保持一致采用 trunc 模。1.5 高精度取模的通用实现步骤基于 trunc 取整的高精度取模通用流程非常简洁计算商 q a / b对商做截断取整q trunc(q)余数 a - b × q这个公式对整数和实数都成立。如果是整数取模只需要先把两个操作数都截断为整数再按上述公式计算即可。1.6 常见坑点盘点负数下 floor 和 trunc 混淆现象对负数取整结果不符合预期根源默认正数思维忽略了向零和向负无穷的区别四舍五入的方向歧义现象0.5到底往哪入银行家舍入还是四舍五入常规实现0.5向绝对值大的方向入即 1.5→2、-1.5→-2取模的符号误解现象负数取模结果符号搞反根源没搞清楚商的取整方式trunc 模余数和被除数同号实数取模的精度误差现象理论上应该为0的余数实际出现极小的非零值根源除法本身的精度截断导致属于正常的数值误差二、immense 实现对照取整取模逐函数精讲immense的取整体系严格遵循移位取整法的通用框架floor、ceil、trunc三个函数共享几乎完全相同的移位外壳只有核心取整逻辑一行之差round基于floor实现取模则基于trunc实现。2.1 公共框架移位-取整-复位在拆解单个函数之前我们先提炼出三个取整函数共用的三段式骨架后面只需要关注差异部分即可。通用骨架immense 取整函数(constimmense T,intPoint){immense tT;// 1. Point 边界裁剪if(Point0Pointt.LTGn)Pointt.LTGn;if(Point0(-Point)t.DCMn)Point-t.DCMn;intL0;// 2. 移位把目标位移到整数最低位for(inti1;iPoint;i)t.Product_10(false),L;for(intiPoint;i-1;i)t.Product_10(true),L--;// 3. 核心取整逻辑三个函数唯一不同的地方...t.DCMp,t.DCMn0;// 统一清空小数// 4. 复位移回原位置for(inti1;iL;i)t.Product_10(true);for(intiL;i-1;i)t.Product_10(false);returnt;}我们先拆解公共部分第一步Point 边界裁剪if(Point0Pointt.LTGn)Pointt.LTGn;if(Point0(-Point)t.DCMn)Point-t.DCMn;对 Point 做边界保护Point 为正保留小数位不能超过整数位数避免移位过度Point 为负保留整数高位不能超过小数位数避免移位过度超出范围就裁剪到边界保证移位操作合法第二步移位与计数intL0;for(inti1;iPoint;i)t.Product_10(false),L;for(intiPoint;i-1;i)t.Product_10(true),L--;Point0小数点左移 Point 位Product_10(false)即除以10L 记录左移次数Point0小数点右移 |Point| 位L 记录右移次数负值移位完成后目标保留位刚好落在整数部分的最低位所有需要丢弃的位都变成了小数部分。第三步复位还原for(inti1;iL;i)t.Product_10(true);for(intiL;i-1;i)t.Product_10(false);用 L 记录的移位次数反向操作把小数点移回原来的位置得到最终取整结果。整个框架完全对应通用原理的「移位→取整→复位」三步法结构极其规整。2.2 trunc截断取整trunc是最简单的取整也是三个函数里唯一不需要加减修正的。完整源码immensetrunc(constimmense T,intPoint){immense tT;if(Point0Pointt.LTGn)Pointt.LTGn;if(Point0(-Point)t.DCMn)Point-t.DCMn;intL0;for(inti1;iPoint;i)t.Product_10(false),L;for(intiPoint;i-1;i)t.Product_10(true),L--;t.DCMp,t.DCMn0;for(inti1;iL;i)t.Product_10(true);for(intiL;i-1;i)t.Product_10(false);returnt;}核心取整只有一行t.DCMp,t.DCMn0;直接清空小数部分相当于向零截断。因为移位后需要保留的位都在整数里小数部分全部丢掉即可。截断取整不关心符号正数负数都是直接丢小数所以不需要额外的符号判断。2.3 floor地板取整floor向负无穷取整正数下和 trunc 一致负数下如果有小数部分要多减1。完整源码immensefloor(constimmense T,intPoint){immense tT;if(Point0Pointt.LTGn)Pointt.LTGn;if(Point0(-Point)t.DCMn)Point-t.DCMn;intL0;for(inti1;iPoint;i)t.Product_10(false),L;for(intiPoint;i-1;i)t.Product_10(true),L--;if(!t.symbolt.DCMp!)tt-immense(1);t.DCMp,t.DCMn0;for(inti1;iL;i)t.Product_10(true);for(intiL;i-1;i)t.Product_10(false);returnt;}和 trunc 唯一的区别就是多了一行核心判断if(!t.symbolt.DCMp!)tt-immense(1);!t.symbol是负数t.DCMp!有小数部分不是整数同时满足时整数部分再减1对应地板取整的定义负数有小数时要往更小的方向负无穷再走一步。比如 -1.2 移位后是 -1.xxx减1变成 -2.xxx截断小数后就是 -2符合 floor(-1.2) -2。正数情况下这行判断不触发行为和 trunc 完全一致。2.4 ceil天花板取整ceil向正无穷取整和 floor 刚好对称负数下和 trunc 一致正数下如果有小数部分要多加1。完整源码immenseceil(constimmense T,intPoint){immense tT;if(Point0Pointt.LTGn)Pointt.LTGn;if(Point0(-Point)t.DCMn)Point-t.DCMn;intL0;for(inti1;iPoint;i)t.Product_10(false),L;for(intiPoint;i-1;i)t.Product_10(true),L--;if(t.symbolt.DCMp!)ttimmense(1);t.DCMp,t.DCMn0;for(inti1;iL;i)t.Product_10(true);for(intiL;i-1;i)t.Product_10(false);returnt;}核心判断行if(t.symbolt.DCMp!)ttimmense(1);正数且有小数部分 → 加1负数不触发行为和 trunc 一致比如 1.2 移位后是 1.xxx加1变成 2.xxx截断小数后就是 2符合 ceil(1.2) 2。设计点评floor 和 ceil 完美对称一个对负数减1一个对正数加1。两个函数除了这一行判断其余代码完全相同结构极其规整维护成本极低。2.5 round四舍五入round没有重复实现移位框架而是直接复用floor采用经典的「加半取整法」。完整源码inlineimmenseround(constimmense T,constintPoint){immense abs_Tabs(T);immensehalf(0.5);immense resfloor(abs_Thalf,Point);if(!T.symbol)res-res;returnres;}实现逻辑先取绝对值统一按正数处理加上 0.5 个对应位的单位这里是加 0.5移位后对应最低位加0.5调用 floor 取整恢复原符号为什么加0.5再floor就是四舍五入小数部分 0.5加0.5后也不到下一个整数floor后不变相当于舍小数部分 ≥ 0.5加0.5后达到或超过下一个整数floor后进一位相当于入这是四舍五入的标准实现方式好处是完全复用已有的 floor 函数不需要再写一遍移位框架。注意这个实现的舍入方向是「0.5向绝对值大的方向入」即 1.5→2、-1.5→-2符合大多数编程语言的 round 行为。2.6 取模运算 operator%immense的取模运算基于 trunc 取整和 C 原生整数取模行为保持一致。完整源码inlineimmense immense::operator%(immense other)const{immense xtrunc(*this),ytrunc(other);returnx-(trunc(x/y)*y);}实现完全对应通用公式先把两个操作数都截断为整数trunc只做整数取模计算商x / y商截断取整trunc(x/y)余数 被除数 - 除数 × 取整后的商因为商是 trunc 取整所以余数的符号和被除数保持一致和 C 原生%运算符行为完全相同学习成本为零。拓展说明当前实现只支持整数取模实数取模需要去掉开头的两个 trunc。对于大多数应用场景整数取模已经足够比如快速幂的指数取模、周期约化等。三、拓展取整取模的应用与边界3.1 Point 参数的灵活用法immense的取整函数都支持 Point 参数默认是 0取整到个位但它可以正负灵活取值Point 2保留 2 位小数比如 1.234 → 1.23Point 0取整到个位默认值Point -2保留到百位比如 12345 → 12300这个设计让取整函数的适用范围大大扩展不仅能处理小数精度还能处理大数的高位约化。3.2 设计权衡点评移位法统一框架优点代码复用率极高三个取整函数共享99%的代码只改一行核心逻辑支持任意指定位数扩展性强缺点需要多次调用Product_10有少量函数调用开销综合来看是可读性、可维护性远大于性能损耗的优秀设计round 基于 floor 复用优点零冗余代码逻辑完全依赖已验证的 floor 函数不易出错缺点多一次加法和一次函数调用开销可忽略不计取模限定整数优点逻辑简单行为和原生C一致符合用户直觉缺点不支持实数取模需要的话用户可自行基于 trunc 扩展3.3 自我验证练习用以下案例推演代码逻辑检验理解程度floor(immense(-1.2), 0)→ 结果应为 -2ceil(immense(1.2), 0)→ 结果应为 2trunc(immense(-1.7), 0)→ 结果应为 -1round(immense(123.45), 1)→ 结果应为 123.5immense(10) % immense(3)→ 结果应为 1immense(-10) % immense(3)→ 结果应为 -1四、本篇小结与下篇预告4.1 核心知识点回顾通用原理部分四种标准取整trunc向零、floor向负无穷、ceil向正无穷、round四舍五入负数下行为差异显著移位取整法移位→取整→复位统一框架处理任意指定位数的取整是高精度库的通用实现方案取模运算公式a mod b a - b × trunc(a/b)trunc模下余数与被除数同号实现对照部分floor、ceil、trunc 共享移位-复位框架仅核心判断一行之差结构极其规整floor 负数有小数减1ceil 正数有小数加1两者完美对称round 采用「加半取整法」完全复用 floor 函数代码零冗余取模基于 trunc 实现行为与 C 原生整数取模一致仅支持整数取模4.2 下篇预告取整取模是高级运算的基石接下来我们就要进入真正的高级运算——幂与根。从整数次幂的快速幂算法到实数次幂的二项式展开再到平方根、任意次开方幂运算是高精度计算中非常核心的模块。immense采用了「整数指数快速幂 小数指数二项级数」的拆分方案把任意实数次幂拆解成两部分分别计算再合并结果。下一篇我们就从快速幂原理讲起一步步推导实数次幂的实现思路再对照pow_I、pow_F、pow、sqrt等函数完整拆解幂与根的实现逻辑。