惯性聚合 高效追踪和阅读你感兴趣的博客、新闻、科技资讯
阅读原文 在惯性聚合中打开

推荐订阅源

H
Help Net Security
Scott Helme
Scott Helme
爱范儿
爱范儿
WordPress大学
WordPress大学
博客园 - 三生石上(FineUI控件)
阮一峰的网络日志
阮一峰的网络日志
博客园 - Franky
V
V2EX
腾讯CDC
博客园_首页
博客园 - 司徒正美
酷 壳 – CoolShell
酷 壳 – CoolShell
T
Tailwind CSS Blog
freeCodeCamp Programming Tutorials: Python, JavaScript, Git & More
OSCHINA 社区最新新闻
OSCHINA 社区最新新闻
小众软件
小众软件
J
Java Code Geeks
大猫的无限游戏
大猫的无限游戏
月光博客
月光博客
Microsoft Azure Blog
Microsoft Azure Blog
B
Blog
雷峰网
雷峰网
Stack Overflow Blog
Stack Overflow Blog
IT之家
IT之家
罗磊的独立博客
Recorded Future
Recorded Future
博客园 - 聂微东
O
OpenAI News
S
Secure Thoughts
Hacker News: Ask HN
Hacker News: Ask HN
S
Schneier on Security
Hacker News - Newest:
Hacker News - Newest: "LLM"
Y
Y Combinator Blog
C
Cyber Attacks, Cyber Crime and Cyber Security
Project Zero
Project Zero
宝玉的分享
宝玉的分享
K
Kaspersky official blog
N
Netflix TechBlog - Medium
T
The Exploit Database - CXSecurity.com
Google Online Security Blog
Google Online Security Blog
cs.CL updates on arXiv.org
cs.CL updates on arXiv.org
cs.CV updates on arXiv.org
cs.CV updates on arXiv.org
Webroot Blog
Webroot Blog
云风的 BLOG
云风的 BLOG
Simon Willison's Weblog
Simon Willison's Weblog
C
Check Point Blog
D
Darknet – Hacking Tools, Hacker News & Cyber Security
L
LINUX DO - 热门话题
美团技术团队
L
Lohrmann on Cybersecurity

博客园_首页

Plist 二进制格式 Milvus 和 PGVector,哪个更好? OpenClaw 已过时?在 VS Code 中运行 Hermes Agent! 第30篇文章:一个大三计科生的自白 Manim如何在数学公式中完美显示中文? Docker 部署 RocketMQ 5 并发编程核心概念辨析 C#事务处理最佳实践:别再让“主表存了、明细丢了”的破事发生 CLI 是什么?为什么大厂突然集体卷命令行? 【从0到1构建一个ClaudeAgent】协作-自主Agent UIImageView 设置图片不生效的原因排查 最小二乘问题详解20:无先验约束下的增量式SFM自由网平差 痞子衡嵌入式:大话双核i.MXRT1180之XIP应用里借助MU实现可靠Flash IAP的方法 AI Chat 封装, SemanticKerne.AiProvider.Unified 已发布 Windows下右键编辑js文件无法打开记事本——在注册表中使用环境变量 在后台服务中使用 Scoped 服务,为什么总是报错? H200 安装驱动并使用sglang启动模型 wireshark 抓包Trap上报告警内容 我用 AI 辅助开发了一系列小工具(2):图片压缩工具 [A Primer On MC and CC] 2.1 Memory Consistency 1 - 指令重排序和 SC 模型 Oracle数据库SCN推进技术详解与实践指南 玩转控件:封装个带图片的Label控件 Claude Code 4.7 真正该升级的不是模型,而是你的工作流 前端小白一句话,AI 帮我做了个颜值拉满的桌面媒体播放器。当代码不再是门槛,一句话编程就是现实。 5. WorkBuddy: 小龙虾的灵魂三件套,让你的小龙虾不只是工具 SQLite 分片方案实战:三种分片策略的深度对比 告别简陋 UI!一款基于 Fluent Design 和基于 WinUI 的开源免费、现代化的 Avalonia UI 控件库 关于二进制排列组合枚举的总结 AI开发-python-LangGraph框架(3-27-LangGraph从零实现大模型智能决策工作流) ElasticSearch主分片和副本分片概念详解 【002】HTTPS 粗解:证书、TLS 握手与对后端配置的影响 Hermes Agent 一周暴涨五万 Star,但我劝你别急着追 明明连接的是Redis的DB0,为什么能查到DB3的数据? 【从0到1构建一个ClaudeAgent】协作-Agent团队 熟悉电子元器件之后,电子小白下一步该怎么走? MAF快速入门(23)通过C#类定义Skills .NET 高级开发 | 手写一个对象映射框架 FastAPI数据库ORM怎么选?我肝了三个Demo后,终于不再纠结了 mysqldump 参数拾遗:在遗忘与铭记之间 C# .NET 周刊|2026年3月5期 Claude code入门 - 陈彦斌 一文学习入门 ThingsBoard 开源物联网平台 GitHub 热门项目 | 2026年04月16日 如何为GIT设置全局勾子,为每次提交追加信息 Number.isFinite和isFinite与isNaN()和Number.isNaN的区别 PortSwigger SQL注入LAB2 推荐一个测试人必备的Skills,从功能到性能全搞定(附详细实操和安装下载方式) 筑基期:掌握Odoo基础核心知识点02(Odoo XML 开发方式详解) GLM模型这么火,咱们用vllm也咧一个呗! 深入理解 AbortController:从底层原理到跨语言设计哲学 字符串学习笔记 多租户系统框架的基础模块设计和分析设计 Apache SeaTunnel Zeta 为什么能做到“又快又稳”? AI开发-python-LangGraph框架(3-26-LangGraph基本概念及第一个简单样例) Vue 3 组件通信,别只会用 Props 和 Emits 了,这几个狠活儿你得看看 ElasticSearch7.X版本配置密码 用Manim实现动态交点计算--从一个动点问题说起 团结引擎+Addressable+Instant Game打包抖音小游戏 function call 实战:让 LLM 自动判断 pod 异常、调用日志工具并完成故障分析 bubseek —— 让 Agent 的足迹,变成团队的洞察 通过 C# 读取并导出 PDF 书签 如何用 GitHub Actions 实现 Steam 自动化发布 【从0到1构建一个ClaudeAgent】并发-后台任务 .NET 高级开发 | 定制 ASP.NET Core 框架 电子小白:什么是运算放大器(运放) zero2Agent:面向大厂面试的 Agent 工程教程,从概念到生产的完整学习路线 堆上的ORW HC32F460 USB CDC通信异常:非对齐访问异常排查 20260413-Hyperbridge 攻击事件:发生在默克尔山上的验证绕过 那些喊着AI 要淘汰你的人,正在靠你的焦虑赚大钱! 深度学习进阶(八)Swin Transformer 最小二乘问题详解19:带先验约束的增量式SFM优化与实现 SnapTranslate 3.0 正式发布:全局划词翻译 + 完整英语学习闭环,一站式搞定查词、记词、复习 工作的意义、工作的困难认知再思考 .NET + AI 进阶实战:基于类的技能开发 - 打造可治理的 Agent 能力模块 【从0到1构建一个ClaudeAgent】规划与协调-技能 上周热点回顾(4.6-4.12) 电子小白的工具三件套:面包板、杜邦线、万能板 单表五亿数据的查询优化 | Mysql、StarRocks 2. WorkBuddy:从“我是谁”到“帮我干活” C# 如何减少代码运行时间:7 个实战技巧 基于HelixToolkit.SharpDX 渲染3D模型 - 笺上知微 从零开始的双臂具身VLA起源及现阶段发展综述 - SkyXZ 记对 xonsh shell 的使用, 脚本编写, 迁移及调优 - pluvium27 受够了Vibe Coding的失控?换个起点,让AI事半功倍 从开始配置漏洞环境到漏洞复现流程 - 難しい 关于10年工作经验的程序员对OpenClaw的实战经验分享以及看法 - 虚无境 Any metadata 的内存布局 C# .NET 周刊|2026年3月2期 - InCerry 我帮你测过了,测试圈排名第二的 Skill 依然很牛逼 Skill Discovery | 无监督技能发现的经典工作总结 - MoonOut 上下文工程是什么?过时了么?一文讲明白! - 一枫说码 开了 TUN 模式还是直连?90% 的人都踩过这个坑 AScript扩展多种脚本语言 - rockey627 AI 学习笔记:Agent 的记忆机制 你能被装进一个文件里吗?——7 万人把同事"蒸馏"成了 AI - 我没有三颗心脏 Claude Code 通关手册(七):给 AI 装上技能包——Skills 完全指南 - 暮色之狐 在浏览器中快速编辑代码:VSCode Web 集成实践 - Newbe36524 蒸馏自己 skill?基于 Deepseek 的蒸馏器,丐版蒸馏方式,简单便捷 - To_Carpe_Diem Spring AI Aliababa和AgentScope,哪个更好? - 苏三说技术
高精度算法
myLv · 2026-06-17 · via 博客园_首页

一、为什么需要高精度运算

普通的整数类型(int、long long)都有固定的取值范围,例如64位 long long 最大只能存储约 9×10¹⁸ 的数值。当我们需要处理几十位、上百位甚至更长的大整数时(比如大数阶乘、密码学计算、超长数值运算等场景),内置类型会发生溢出,无法正确存储和计算。高精度运算的核心思路就是用字符串或数组存储大整数的每一位数字,模拟人类手算的方式完成加减乘除,从而突破数值范围的限制。

二、基本概念

高精度运算,本质是将大整数按位拆分,用字符串(或数组)保存每一位的数字,完全模拟小学阶段学习的竖式加减乘除规则,逐位计算并处理进位、借位。

  • 存储方式:通常用字符串保存,高位在前、低位在后(和日常书写习惯一致),在计算时会反转字符串,让下标 0 对应最低位,方便对齐处理。
  • 核心操作:进位处理(加法、乘法)、借位处理(减法)、前导零去除、数值大小比较。

三、实现方法

1)加法

运算原理

加法完全模拟小学竖式加法规则:

  1. 将两个数右对齐,从最低位开始逐位相加;
  2. 每一位的和 = 第一个数当前位 \(a[i]\) + 第二个数当前位 \(b[i]\) + 低位传来的进位 \(carry\)
  3. 当前位结果 = 和 % 10,向高位的进位 = 和 / 10;
  4. 所有位计算完成后,如果还有剩余进位,要在最高位补上。
// 代码中先将字符串反转,是为了让数组下标 0 对应最低位,天然实现右对齐,无需手动补位
string add(string a, string b)
{
    if (a.length() < b.length()) swap(a, b); // 大的数字放“上面”,方便操作
    reverse(a.begin(), a.end()), reverse(b.begin(), b.end());
    string res;
    int carry = 0;
    for (size_t i = 0; i < a.length(); ++i)
    {
        int x = a[i] - '0';
        int y = i < b.length() ? b[i] - '0' : 0;
        int sum = x + y + carry;
        res += (sum % 10) + '0';
        carry = sum / 10;
    }
    if (carry) res += carry + '0';
    reverse(res.begin(), res.end());
    return res;
}

2)减法

运算原理

减法同样模拟小学竖式减法规则:

  1. 先比较两个数的大小,用大数减小数,保证结果非负,最后再补上负号;
  2. 将两个数右对齐,从最低位开始逐位相减;
  3. 如果当前位被减数 < 减数 + 借位,就向高位借 1 当 10,同时标记借位;
  4. 计算完成后,结果会产生多余的前导零,需要去除。
为什么要去除前导零?

减法计算完成后,高位可能会出现连续的 0。例如计算 \(1000 - 999 = 1\),按位计算后反转得到的中间结果会是 \(0001\),前面的三个 0 就是前导零。前导零不符合数字的书写规范,也会影响后续运算的长度判断,因此必须去除多余的前导零,仅保留最后一个有效数字(注意结果为 0 时保留一个 0)。

string sub(string a, string b)
{
    bool swaped = false;
    if (a.length() < b.length() || (a.length() == b.length() && a < b))
    {
        swap(a, b);
        swaped = true;
    }
    reverse(a.begin(), a.end()), reverse(b.begin(), b.end());
    string res;
    int carry = 0;
    for (size_t i = 0; i < a.length(); ++i)
    {
        int x = a[i] - '0';
        int y = i < b.length() ? b[i] - '0' : 0;
        int sum = x - y - carry;
        if (sum < 0)
        {
            sum += 10;
            carry++;
        } else carry = 0;
        res += sum + '0';
    }
    /*-----去除前导零-----*/
    while (res.length() > 1 && res.back() == '0') res.pop_back();
    if (swaped) res += '-';
    reverse(res.begin(), res.end());
    return res;
}

3)乘法

运算原理

乘法模拟小学竖式乘法的错位相加规则:

  1. 第一个数的第 i 位,乘以第二个数的第 j 位,结果会落在结果的第 i+j 位上(最低位从 0 开始计数);
  2. 先把所有位的乘积累加到位对应的位置,暂不处理进位;
  3. 全部相乘完成后,统一从低位到高位处理进位,每一位满十进一;
  4. 最后去除高位多余的前导零,再转换为字符串。

注意:两个长度分别为 La 和 Lb 的数相乘,结果长度最多为 La+Lb 位。

string mul(string a, string b)
{
    int La = a.size(), Lb = b.size();
    vector<int> res(La + Lb, 0);
    reverse(a.begin(),a.end()), reverse(b.begin(),b.end());

    for (size_t i = 0; i < a.length(); i++)
    {
        int x = a[i] - '0';
        for (size_t j = 0; j < b.length(); j++)
        {
            int y = b[j] - '0';
            res[i+j] += x * y; // 把每次每位的相乘结果放在数组中,先不进位
        }
    }

    /*-----进位-----*/
    for (int i = 0; i < (int)res.size(); i++)
        if (res[i] >= 10)
        {
            res[i+1] += res[i] / 10;
            res[i] %= 10;
        }

    /*-----去除前导零-----*/
    while (res.size() > 1 && res.back() == 0) res.pop_back();

    string ret;
    for (int i = res.size()-1; i >= 0; i--) ret += res[i] + '0';

    return ret;
}

4)除法

运算原理

这里实现的是高精度整数除以普通 int 整数(高精度除以高精度复杂度更高,日常场景中高精度除低精度最常用),同样模拟竖式除法规则:

  1. 从被除数的最高位开始,逐位向下计算;
  2. 当前位数值 = 上一位的余数 × 10 + 当前位数字;
  3. 当前位的商 = 当前位数值 ÷ 除数,新的余数 = 当前位数值 % 除数;
  4. 逐位计算完成后,去除商的前导零,同时可以返回余数。

代码中提供了两种写法:一种返回商和余数的 pair,一种只返回商,核心逻辑完全一致。

// 高精度除法:a / b,a为大整数(字符串),b为int,返回商和余数
pair<string, int> div(string a, int b)
{
    string res;
    int remainder = 0; // 余数
    for (size_t i = 0; i < a.length(); ++i)
    {
        remainder = remainder * 10 + (a[i] - '0');
        res += (remainder / b) + '0';
        remainder %= b;
    }
    // 去除前导零
    size_t pos = res.find_first_not_of('0');
    if (pos == string::npos) res = "0";
    else res = res.substr(pos);

    return {res, remainder};
}

string divide(string a, int b)
{
    if (b == 0) return "Error: Division by zero";
    
    vector<int> num;
    for (char c : a) num.push_back(c - '0');
    
    vector<int> result;
    int remainder = 0;
    
    for (size_t i = 0; i < num.size(); i++)
    {
        int current = remainder * 10 + num[i];
        result.push_back(current / b);
        remainder = current % b;
    }
    
    // 去除前导零
    size_t start = 0;
    while (start < result.size() && result[start] == 0) start++;
    
    if (start == result.size()) return "0"; // 等于0时保留一个0
    
    string ret = "";
    for (size_t i = start; i < result.size(); i++) ret += result[i] + '0';
    
    return ret;
}

四、总结

高精度四则运算的核心思想高度统一:用字符串/数组拆分存储大整数的每一位,完全模拟人工竖式计算的规则

  • 加法、乘法:从低位到高位计算,重点处理进位;
  • 减法:从低位到高位计算,重点处理借位,注意正负号判断;
  • 除法(低精度除数):从高位到低位计算,逐位落数,保留余数。

所有运算最终都需要处理前导零,保证输出格式符合常规数字的书写习惯。这套实现可以处理任意长度的正整数运算,是算法竞赛和大数计算场景的基础工具。

五、最后碎碎念

由于学期接近期末,最近都在准备期末考试,所以更新慢了许多,还请见谅!此外如果你发现我的博客有错误的地方、语义矛盾或者有更好的方法,欢迎你的建议与指导!