C++高精度算法的简单实现

news2026/5/18 22:49:58
一、基本原理1、存储方式采用数字记录高精度数字数组的第一个元素存储数据长度比如记录数字为1024示例如下2、计算方式采用模拟立竖式计算比如加法的计算流程如下图所示10249000这里只给出加法的计算说明其他的以此类推减法与加法基本一致。乘法和除法略有不同通过示例图表示也复杂还不如通过代码去理解本质的方法就是模拟笔算的立竖式计算。二、辅助方法1、字符串转高精度长度记录在数组第一个元素中123456789101112/// summary/// 通过字符串初始化/// /summary/// param namea[in]高精度数组/param/// param namevalue[in]字符串首地址/paramstaticvoidloadStr(int* a,constchar* value){//记录长度a[0] strlen(value);for(inti 1; i a[0]; i)a[i] value[a[0] - i] -0;}2、整型转高精度12345678910111213141516171819/// summary/// 通过无符号整型初始化/// /summary/// param namea[in]高精度数组/param/// param namevalue[in]整型值/paramstaticvoidloadInt(int* a, uint64_t value){for(size_ti 1; i 8096; i){a[i] value % 10;value / 10;if(!value){//记录长度a[0] i;return;}}}3、比较123456789101112131415/// summary/// 比较两个高精度数的大小/// /summary/// param namea[in]第一个数/param/// param nameb[in]第二个数/param/// returns1是ab,0是ab,-1是ab/returnsstaticintcompare(int* a,int* b){if(a[0] b[0])return1;if(a[0] b[0])return-1;for(inti a[0]; i 0; i--)if(a[i] b[i])return1;elseif(a[i] b[i])return-1;return0;}4、打印123456789/// summary/// 打印输出结果/// /summarystaticvoidprint(int* a) {if(!a[0])printf(0);for(inti a[0]; i 0; i--)printf(%d, a[i]);}三、算法实现原理就不做具体介绍了四种计算的核心都是模拟立竖式计算。1、加法为了保证代码相对简单当b长度较小时可能会做一些多余的计算不影响结果。123456789101112131415161718/// summary/// 加法(累加)///结果会保存在a中/// /summary/// param namea[in]被加数/param/// param nameb[in]加数/paramstaticvoidacc(int* a,int* b){intlen a[0] b[0] ? a[0] : b[0];memset(a a[0] 1, 0, (len - a[0] 1) *sizeof(int));memset(b b[0] 1, 0, (len - b[0] 1) *sizeof(int));for(inti 1; i len; i) {inttemp a[i] b[i];a[i] temp % 10;a[i 1] temp / 10;}if(a[len 1])a[0];}2、减法12345678910111213141516171819202122232425262728/// summary/// 减法(累减)///结果会保存在a中/// /summary/// param namea[in]被减数被减数必须大于等于减数/param/// param nameb[in]减数/paramstaticvoidsubc(int* a,int* b) {memset(b b[0] 1, 0, (a[0] - b[0]) *sizeof(int));for(inti 1; i a[0]; i){inttemp a[i] - b[i];a[i] temp;if(temp 0){//借位a[i 1] - 1;a[i] 10;}}//记录长度for(inti a[0]; i 0; i--)if(a[i]){a[0] i;return;}a[0] 0;}3、乘法123456789101112131415161718192021222324252627282930313233/// summary/// 乘法/// /summary/// param namea[in]被乘数/param/// param nameb[in]乘数/param/// param namec[out]结果数组长度必须大于等于aLenbLen1/paramstaticvoidmul(int* a,int* b,intc[]) {c[a[0] b[0]] 0;memset(c, 0,sizeof(int) * (a[0] b[0] 1));for(inti 1; i a[0]; i){intj;intd 0;//被乘数的一位去乘以乘数的每一位for(j 1; j b[0]; j){inttemp a[i] * b[j] c[j i - 1] d;c[j i - 1] temp % 10;d temp / 10;}if(d){c[j i - 1] d;}}//记录长度for(inti a[0] b[0]; i 0; i--)if(c[i]){c[0] i;return;}}4、除法采用了升阶减法实现123456789101112131415161718192021222324252627282930313233343536/// summary/// 除法/// 依赖减法subc/// /summary/// param namea[in]被除数,被除数必须大于除数/param/// param nameb[in]除数/param/// param namec[out]商数组长度大于等于aLen-bLen1/param/// param namemod[out]余数数组长度大于等于aLen/param/// param nametemp[in]临时缓冲区由外部提供以提高性能数组长度大于等于aLen-bLen1/paramstaticvoiddivi(int* a,int* b,int* c,int* mod,int* temp) {//相差的阶数intdigit a[0] - b[0] 1;memcpy(mod, a, (a[0] 1) *sizeof(int));memset(c, 0,sizeof(int) * (digit 1));memset(temp, 0,sizeof(int) * digit);while(digit){//升阶memcpy(temp digit, b 1,sizeof(int) * b[0]);temp[0] b[0] digit - 1;//减法while(compare(mod, temp) ! -1){subc(mod, temp);c[digit];}digit--;}//记录长度for(inti a[0] - b[0] 1; i 0; i--)if(c[i]){c[0] i;return;}}四、使用示例1、加法计算累加1234567891011121314intmain() {int64_t n;intnum[1024];intnum2[1024];std::cin n;loadInt(num, 0);for(int64_t i 1; i n; i){loadInt(num2, i);acc(num, num2);}print(num);return0;}结果2、减法两个任意n位数的减法数字1大于数字2。1234567891011intmain(){inta1[8096], a2[8096];std::string s1, s2;std::cin s1 s2;loadStr(a1, s1.c_str());loadStr(a2, s2.c_str());subc(a1, a2);print(a1);return0;}结果#数字1752425289999999999999652142141414141414146666676667677682324000001302461646520#数字2587891851201874512000000000154515100202121555555555555555555555545477910232111#计算结果1645334387981254879996521419868990412120251111211121221267684444558245514144093、乘法计算阶乘1234567891011121314151617181920intmain() {int64_t n;intnum[8192];intnum2[8192];intnum3[8192];int* p1 num;int* p2 num3;std::cin n;loadInt(num, 1);for(int64_t i 1; i n; i){loadInt(num2, i);mul(p1, num2, p2);int* temp p1;p1 p2;p2 temp;}print(p1);return0;}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2623082.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

SpringBoot-17-MyBatis动态SQL标签之常用标签

文章目录 1 代码1.1 实体User.java1.2 接口UserMapper.java1.3 映射UserMapper.xml1.3.1 标签if1.3.2 标签if和where1.3.3 标签choose和when和otherwise1.4 UserController.java2 常用动态SQL标签2.1 标签set2.1.1 UserMapper.java2.1.2 UserMapper.xml2.1.3 UserController.ja…

wordpress后台更新后 前端没变化的解决方法

使用siteground主机的wordpress网站,会出现更新了网站内容和修改了php模板文件、js文件、css文件、图片文件后,网站没有变化的情况。 不熟悉siteground主机的新手,遇到这个问题,就很抓狂,明明是哪都没操作错误&#x…

网络编程(Modbus进阶)

思维导图 Modbus RTU(先学一点理论) 概念 Modbus RTU 是工业自动化领域 最广泛应用的串行通信协议,由 Modicon 公司(现施耐德电气)于 1979 年推出。它以 高效率、强健性、易实现的特点成为工业控制系统的通信标准。 包…

UE5 学习系列(二)用户操作界面及介绍

这篇博客是 UE5 学习系列博客的第二篇,在第一篇的基础上展开这篇内容。博客参考的 B 站视频资料和第一篇的链接如下: 【Note】:如果你已经完成安装等操作,可以只执行第一篇博客中 2. 新建一个空白游戏项目 章节操作,重…

IDEA运行Tomcat出现乱码问题解决汇总

最近正值期末周,有很多同学在写期末Java web作业时,运行tomcat出现乱码问题,经过多次解决与研究,我做了如下整理: 原因: IDEA本身编码与tomcat的编码与Windows编码不同导致,Windows 系统控制台…

利用最小二乘法找圆心和半径

#include <iostream> #include <vector> #include <cmath> #include <Eigen/Dense> // 需安装Eigen库用于矩阵运算 // 定义点结构 struct Point { double x, y; Point(double x_, double y_) : x(x_), y(y_) {} }; // 最小二乘法求圆心和半径 …

使用docker在3台服务器上搭建基于redis 6.x的一主两从三台均是哨兵模式

一、环境及版本说明 如果服务器已经安装了docker,则忽略此步骤,如果没有安装,则可以按照一下方式安装: 1. 在线安装(有互联网环境): 请看我这篇文章 传送阵>> 点我查看 2. 离线安装(内网环境):请看我这篇文章 传送阵>> 点我查看 说明&#xff1a;假设每台服务器已…

XML Group端口详解

在XML数据映射过程中&#xff0c;经常需要对数据进行分组聚合操作。例如&#xff0c;当处理包含多个物料明细的XML文件时&#xff0c;可能需要将相同物料号的明细归为一组&#xff0c;或对相同物料号的数量进行求和计算。传统实现方式通常需要编写脚本代码&#xff0c;增加了开…

LBE-LEX系列工业语音播放器|预警播报器|喇叭蜂鸣器的上位机配置操作说明

LBE-LEX系列工业语音播放器|预警播报器|喇叭蜂鸣器专为工业环境精心打造&#xff0c;完美适配AGV和无人叉车。同时&#xff0c;集成以太网与语音合成技术&#xff0c;为各类高级系统&#xff08;如MES、调度系统、库位管理、立库等&#xff09;提供高效便捷的语音交互体验。 L…

(LeetCode 每日一题) 3442. 奇偶频次间的最大差值 I (哈希、字符串)

题目&#xff1a;3442. 奇偶频次间的最大差值 I 思路 &#xff1a;哈希&#xff0c;时间复杂度0(n)。 用哈希表来记录每个字符串中字符的分布情况&#xff0c;哈希表这里用数组即可实现。 C版本&#xff1a; class Solution { public:int maxDifference(string s) {int a[26]…

【大模型RAG】拍照搜题技术架构速览:三层管道、两级检索、兜底大模型

摘要 拍照搜题系统采用“三层管道&#xff08;多模态 OCR → 语义检索 → 答案渲染&#xff09;、两级检索&#xff08;倒排 BM25 向量 HNSW&#xff09;并以大语言模型兜底”的整体框架&#xff1a; 多模态 OCR 层 将题目图片经过超分、去噪、倾斜校正后&#xff0c;分别用…

【Axure高保真原型】引导弹窗

今天和大家中分享引导弹窗的原型模板&#xff0c;载入页面后&#xff0c;会显示引导弹窗&#xff0c;适用于引导用户使用页面&#xff0c;点击完成后&#xff0c;会显示下一个引导弹窗&#xff0c;直至最后一个引导弹窗完成后进入首页。具体效果可以点击下方视频观看或打开下方…

接口测试中缓存处理策略

在接口测试中&#xff0c;缓存处理策略是一个关键环节&#xff0c;直接影响测试结果的准确性和可靠性。合理的缓存处理策略能够确保测试环境的一致性&#xff0c;避免因缓存数据导致的测试偏差。以下是接口测试中常见的缓存处理策略及其详细说明&#xff1a; 一、缓存处理的核…

龙虎榜——20250610

上证指数放量收阴线&#xff0c;个股多数下跌&#xff0c;盘中受消息影响大幅波动。 深证指数放量收阴线形成顶分型&#xff0c;指数短线有调整的需求&#xff0c;大概需要一两天。 2025年6月10日龙虎榜行业方向分析 1. 金融科技 代表标的&#xff1a;御银股份、雄帝科技 驱动…

观成科技:隐蔽隧道工具Ligolo-ng加密流量分析

1.工具介绍 Ligolo-ng是一款由go编写的高效隧道工具&#xff0c;该工具基于TUN接口实现其功能&#xff0c;利用反向TCP/TLS连接建立一条隐蔽的通信信道&#xff0c;支持使用Let’s Encrypt自动生成证书。Ligolo-ng的通信隐蔽性体现在其支持多种连接方式&#xff0c;适应复杂网…

铭豹扩展坞 USB转网口 突然无法识别解决方法

当 USB 转网口扩展坞在一台笔记本上无法识别,但在其他电脑上正常工作时,问题通常出在笔记本自身或其与扩展坞的兼容性上。以下是系统化的定位思路和排查步骤,帮助你快速找到故障原因: 背景: 一个M-pard(铭豹)扩展坞的网卡突然无法识别了,扩展出来的三个USB接口正常。…

未来机器人的大脑:如何用神经网络模拟器实现更智能的决策?

编辑&#xff1a;陈萍萍的公主一点人工一点智能 未来机器人的大脑&#xff1a;如何用神经网络模拟器实现更智能的决策&#xff1f;RWM通过双自回归机制有效解决了复合误差、部分可观测性和随机动力学等关键挑战&#xff0c;在不依赖领域特定归纳偏见的条件下实现了卓越的预测准…

Linux应用开发之网络套接字编程(实例篇)

服务端与客户端单连接 服务端代码 #include <sys/socket.h> #include <sys/types.h> #include <netinet/in.h> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <arpa/inet.h> #include <pthread.h> …

华为云AI开发平台ModelArts

华为云ModelArts&#xff1a;重塑AI开发流程的“智能引擎”与“创新加速器”&#xff01; 在人工智能浪潮席卷全球的2025年&#xff0c;企业拥抱AI的意愿空前高涨&#xff0c;但技术门槛高、流程复杂、资源投入巨大的现实&#xff0c;却让许多创新构想止步于实验室。数据科学家…

深度学习在微纳光子学中的应用

深度学习在微纳光子学中的主要应用方向 深度学习与微纳光子学的结合主要集中在以下几个方向&#xff1a; 逆向设计 通过神经网络快速预测微纳结构的光学响应&#xff0c;替代传统耗时的数值模拟方法。例如设计超表面、光子晶体等结构。 特征提取与优化 从复杂的光学数据中自…