高精度原理介绍及代码实现

news2024/10/6 16:23:16

目录

高精度

引入

使用场景

实现原理

高精度加法

数据存储

加法实现

总代码

高精度减法

与加法的不同点:

总代码

高精度乘法

总代码

高精度除法

总结

总注意点

减法注意点

高精度

引入

所谓高精度并不是很高级难懂的东西,只是对传统的加减法模拟实现

使用场景

高精度算法(High Accuracy Algorithm)的出现是为了处理超大数据的数学计算问题。在一般的科学计算中,我们可能会遇到需要计算小数点后几百位甚至更多的数字,或者处理几千亿、几百亿这样的大数字。这些数字超出了标准数据类型(如整型、实型)能够表示的范围,因此无法直接在计算机中正常存储和计算

实现原理

在高精度算法中,我们并不会把数据用int,long long,double这种数据结构来存储,而是用数组来存储,接下来我就讲讲其原理

高精度加法

数据存储

先从简单的加减运算说起,我们拿 23 + 25这个例子来举例,根据小学所学知识,我们知道当我们计算时应该列出这样子的式子:

如果我们要进行加法,那么我们将从右加到左,一一进位

而我们的算法将这里的27处理为:A[2] = {7,2};25处理为:B[2] = {5,2};

他们的和仍然是一个数组:C[2] = {2, 5};

你可能好奇,为什么是倒着存?🧐🤨

俺在前面也说过:从右加到左,因此,当我倒着存时,最低位(27的最低位为7)就在A[0],而当A[0] + B[0]要进位时,进位的1就被加到了C[1],综上所述,倒着存是为了方便进位(i进的位直接存在了i + 1的位置)(也就是7 + 5进的位存在了C[1])

倒着存应该怎么写呢?🤨🤨🤨

很简单,我们要注意读入数据时不能把存数据的变量设置为int,long long了,而是要设置为string

//a + b = c
​
int A[505];//第一个加数的数组
int B[505];//第二个加数的数组
int C[506];//和的数组
int len_c;
int main()
{
    string a, b;//注意用字符串读入
    cin >> a >> b;
​
    int len_a = a.size(), len_b = b.size();//len_a表示a有几个数字(若a = 1,则len_a = 1,若a = 10,则len_a = 2)
    int len_c = max(len_a, len_b);//a 和 b的和c最初设置为a 和 b的最大值,因为要确保每位都参与加法,就得让位数多的每位都加上
    
    //逆序存
    for (int i = 0; i < len_a; i++)
        A[len_a - 1 - i] = a[i] - '0';//注意 - ‘0’,因为是字符读入,且逆序存
    for (int i = 0; i < len_b; i++)
        B[len_b - 1 - i] = b[i] - '0';//注意 - ‘0’
    Add();
    
    return 0;
}

接下来就是重要的加法实现了

加法实现

先看代码

void Add()
{
    for (int i = 0; i < len_c; i++)
    {
        C[i] += A[i] + B[i];
        C[i + 1] += C[i] / 10;
        C[i] %= 10;
    }
​
    if (C[len_c])
        len_c++;
}
  • for循环:for循环是从0开始的,这是因为我们存储的时候0位放的就是加数的最低位(27的7)

  • A[i] + B[i]:模仿传统加法(我们例子中就是:7 + 5)

  • C[i + 1] += C[i] / 10;:模仿进位(我们例子中就是:(7 + 5) / 10 = 1,因此在C[1]上进位了1)

  • C[i] %= 10; : 本位进位完后就要变为只有个位(我们例子中就是:(7 + 5) /%10 = 2,因此在C[0]上只有2)

  • if (C[len_c]) len_c++;:如果结果中最高位下一位不为0,说明len_c的长度要++

    如同:

    9 7 + 3,在len_c = max(len_a, len_b);后,len_c的长度为2,但是当计算完97 + 3 = 100后:之前的最高位(第二位)的下一位(第三位)不为0(100的1),长度变为3,因此要更新len_c

没想到吧,到这里就把高精度讲完啦,你看到这里很棒了哦,给你一个赞 d=====( ̄▽ ̄*)b

总代码

//a + b = c
​
int A[505];//第一个加数的数组
int B[505];//第二个加数的数组
int C[506];//和的数组
int len_c;
​
void Add()
{
    for (int i = 0; i < len_c; i++)
    {
        C[i] += A[i] + B[i];
        C[i + 1] += C[i] / 10;
        C[i] %= 10;
    }
​
    if (C[len_c])
        len_c++;
}
​
int main()
{
    string a, b;//注意用字符串读入
    cin >> a >> b;
​
    int len_a = a.size(), len_b = b.size();
    len_c = max(len_a, len_b);
    for (int i = 0; i < len_a; i++)
        A[len_a - 1 - i] = a[i] - '0';//注意 - ‘0’,且逆序存
    for (int i = 0; i < len_b; i++)
        B[len_b - 1 - i] = b[i] - '0';//注意 - ‘0’
​
    Add();
​
    for (int i = 0; i < len_c; i++)
        cout << C[len_c - 1 - i];//注意逆序输出,因为存的时候是逆序的
    return 0;
}

高精度减法

这也是类似的,高精度通用的就是:

  • 数据用字符串读入,用数组存数

  • 要倒着存储数据

  • 注意进位

  • 注意逆序输出

与加法的不同点:

敲重点了哦~ψ(`∇´)ψ

  • 要看两个数据谁更大,要把大的作为a,小的作为b

  • 若起初a < b,用标记标记该答案为负数,以便最后得出答案时输出符号

总代码

//与加法相似
//加了调转字符串的操作
int main()
{
    string str1, str2;
    cin >> str1;
    cin >> str2;
​
    //两个数字最大为10^10086,因此数组设为10088就可(咱图个吉利)
    int a[10088] = { 0 }, b[10088] = { 0 }, c[10088] = { 0 };
    
    int flag = 0;
    //标记是否进行了调转字符串
​
    //与加法不同的是减法要调转字符串,把位数多的放前面,方便计算
    //至于负号,可以先标记,最后再进行处理
​
    if (str2.size() > str1.size()/*如果起初的b > a,---->要调转a, b, 把原先的a 变为 b*/||(str2.size() == str1.size()&&str1<str2))
    {
        /*!!!!!!!!!!!!!!!!!!!!注意或者后面的情况,若只有或者前面部分的,若是3-4这种情况,则无法算出正确答案*/
        //str1<str2:前者的数字小于后面的数,虽然这是字符串,但仍然可以这样比较
        flag = 1;
        //标记进行了调换:方便最后输出 - 号
        swap(str1, str2);
        //调转的函数,具体的大家可以去自行了解
    }
​
    for (int i = 0; i < str1.size(); i++)
    {
        a[i] = str1[str1.size() - 1 - i] - '0';
        //减'0'别忘了
    }
​
    for (int i = 0; i < str2.size(); i++)
    {
        b[i] = str2[str2.size() - 1 - i] - '0';
    }
​
    //因为一开始就进行了调换最长的数放在str1,所以这里不需要取str1和str2中最长的赋值给len
    for (int i = 0; i < str1.size(); i++)
    {
        c[i] = a[i] - b[i];
        if (c[i] < 0)
        {
            a[i + 1]--;
            c[i] = a[i] + 10 - b[i];
            //记得给a[i]加10,就是模拟实际的计算
        }
    }
​
    //因为一开始就进行了调换最长的数放在str1,所以这里不需要取str1和str2中最长的赋值给len,叶不用进行len++
    //直接判断是否要进行len--
    int len = str1.size();
​
    /*是while,不是if*/
​
    while(c[len - 1] == 0 && len > 1)
        len--;
​
​
    if (flag == 1)//输出负号
        printf("-");
    for (int i = 0; i < len; i++)//仍然是倒着输出
    {
        printf("%d", c[len - 1 - i]);
    }
    printf("\n");
​
    return 0;
}

高精度乘法

总代码

int c[5000];
int len_a, len_b, len_c;
void mul(int a[], int b[])
{
    for (int i = 0; i < len_a; i++)
    {
        for (int j = 0; j < len_b; j++)
        {
            c[i + j] += a[i] * b[j];//注意是 +=
            c[i + j + 1] += c[i + j] / 10;
            c[i + j] = c[i + j] % 10;
        }
    }
​
    while (len_c > 1 && c[len_c - 1] == 0)
    //注意要len_c > 1,如果乘积为0,就会一直减下去,最后什么都不输出,这是错误的
    //注意是c[len_c - 1] == 0
        len_c--;
}
​
int main()
{
    string a, b;
    cin >> a >> b;
    len_a = a.size(), len_b = b.size();
    len_c = len_a + len_b;
    int aarr[2005], barr[2005];
​
    //倒着存
    for (int i = 0; i < len_a; i++)
        aarr[i] = a[len_a - 1 - i] - '0';//记得 - '0'
    for (int i = 0; i < len_b; i++)
        barr[i] = b[len_b - 1 - i] - '0';
    
    mul(aarr, barr);
​
    //cout << len_c << endl;
    //倒着输出
    for (int i = 0; i < len_c; i++)
        cout << c[len_c - 1 - i];
    return 0;
}

高精度除法

脑子不好的小菜鸟正在学习中::>_<::

总结

总注意点

综上所述,所谓的高精度就是模仿传统的运算法则,注意总点如下:

  • 不用数据类型存储数据,因为计算机存不下那么大的数字

  • 数组存储,且逆序存储

  • 注意进位,还有前置0

    (如:100 - 90 = 010,这个最高位的0是不输出的,要用如下句子处理掉:

    while (len_c > 1 && c[len_c - 1] == 0)
        //注意要len_c > 1,如果答案为0,就会一直减下去,最后什么都不输出,这是错误的
        //注意是c[len_c - 1] == 0
            len_c--;

减法注意点

  • 要看两个数据谁更大,要把大的作为a,小的作为b

  • 若起初a < b,要用标记标记该答案为负数,以便最后得出答案时输出符号

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

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

相关文章

改进YOLOv5,YOLOv5+CBAM注意力机制

目录 1. 目标检测模型 2. YOLOv5s 3. YOLOv5s融合注意力机制 4. 修改yolov5.yaml文件 5. ChannelAttentionModule.py 6. 修改yolo.py 1. 目标检测模型 目标检测算法现在已经在实际中广泛应用&#xff0c;其目的是找出图像中感兴趣的对象&#xff0c;并确定对象的类别和位…

一文读懂设计模式-单例模式

单例模式&#xff08;Singleton Pattern&#xff09;提供了一种创建对象的最佳方式 单例模式涉及到一个单一的类&#xff0c;该类负责创建自己的对象&#xff0c;同时确保只有单个对象被创建&#xff0c;这个类提供了一种访问其唯一的对象的方式&#xff0c;可以直接访问&…

Redis—图文详解高可用原因

本文不会讲解Redis的用途&#xff0c;关于用途会发另一片文章讲解&#xff0c;本文主要讲的是高可用的原理。 Redis高可用主要有以下三个原因&#xff1a;主从模式(上一篇讲Kafka的文章里有涉及到)&#xff0c;哨兵模式&#xff0c;Redis-Cluster(Redis集群)。 什么是主从模式…

mysql集群NDBcluster引擎在写入数据时报错 (1114, “The table ‘ads‘ is full“)

问题描述&#xff1a;mysql集群在写入数据时&#xff0c;出现上述报错 问题原因&#xff1a;表数据已满&#xff0c;一般是在集群的管理节点设置里面datamemory的值太小&#xff0c;当数据量超过该值时就会出现该问题 解决方案&#xff1a; 修改集群管理节点的config.ini里面…

【Linux 网络】网络编程套接字 -- 详解

⚪ 预备知识 1、理解源 IP 地址和目的 IP 地址 举例理解&#xff1a;&#xff08;唐僧西天取经&#xff09; 在 IP 数据包头部中 有两个 IP 地址&#xff0c; 分别叫做源 IP 地址 和目的 IP 地址。 如果我们的台式机或者笔记本没有 IP 地址就无法上网&#xff0c;而因为…

操作系统原理与实验——实验十四多道批处理作业调度算法应用二

样例 样例一&#xff1a; 1 100 0 5 JOB1 10:06 42分钟 55K JOB2 10:20 30分钟 40K JOB3 10:30 24分钟 35K JOB4 10:36 15分钟 25K JOB5 10:42 12分钟 20K 2 100 0 5 JOB1 10:06 42分钟 55K JOB2 10:20 30分钟 40K JOB3 10:30 24分钟 35K JOB4 10:36 15分钟 25K JOB5 10…

开发一款相亲交友小程序

uni-app框架&#xff1a;使用Vue.js开发跨平台应用的前端框架&#xff0c;编写一套代码&#xff0c;可编译到Android、小程序等平台。 框架支持:springboot/Ssm/thinkphp/django/flask/express均支持 前端开发:vue.js 可选语言&#xff1a;pythonjavanode.jsphp均支持 运行软件…

探索共享内存:解锁并发编程的潜力

文章目录 序言shm 原理对shm的理解通过代码认识shm调用shmget方法实现 序言 system V版本 指定的一种规则(俗话说一流公司定规则,二流公司重服务,三流公司重技术).这个规则虽然有很多种(消息队列,共享内存等只是比较出名的几个).但是在内核的相关技术解决上是类似的,因为都是基…

NIUKE SQL:进阶挑战 (下)

05 窗口函数 专用窗口函数 SQL136 每类试卷得分前3名 找到每类试卷得分的前3名&#xff0c;如果两人最大分数相同&#xff0c;选择最小分数大者&#xff0c;如果还相同&#xff0c;选择uid大者 SELECT tag,uid,ranking FROM(SELECTtag,ER.uid,ROW_NUMBER() OVER (PARTITION …

jdk8的新特征

1&#xff1a; jdk8中新增的方法 在jdk8中对接口进行了增强&#xff0c;在jdk8之前 interface 接口名{ 静态常量&#xff1a; 抽象方法&#xff1a; } 在jdk8之后 interface 接口名{ 静态常量&#xff1a; 抽象方法&#xff1a; 默认方法&#xff1a; 静态方法&#xff1a; } 2…

【IMX6ULL项目】IMX6ULL下Linux实现产测工具框架

电子产品量产测试与烧写工具。这是一套软件&#xff0c;用在我们的实际生产中&#xff0c; 有如下特点&#xff1a; 1.简单易用&#xff1a; 把这套软件烧写在 SD 卡上&#xff0c;插到 IMX6ULL 板子里并启动&#xff0c;它就会自动测试各个模块、烧写 EMMC 系统。 工人只要按…

【算法】dfs

快乐的流畅&#xff1a;个人主页 个人专栏&#xff1a;《算法神殿》《数据结构世界》《进击的C》 远方有一堆篝火&#xff0c;在为久候之人燃烧&#xff01; 文章目录 引言一、全排列1. 决策树2. 设计代码1. 全局变量2. dfs函数3. 细节问题 二、子集解法一1. 决策树2. 设计代码…

Spring Boot:让微服务开发像搭积木一样简单!

带你一探 Spring Boot 的自动配置和 Starter POMs 的神奇之处&#xff0c;展示如何通过几个简单的步骤就能让你的微服务应用在云端翱翔&#xff01; 文章目录 1. 引言1.1 简述Spring框架的起源与重要性1.2 阐述文章目的&#xff1a;深入解析Spring核心功能与应用实践2. 背景介绍…

C语言 | Leetcode C语言题解之第85题最大矩形

题目&#xff1a; 题解&#xff1a; int maximalRectangle(char** matrix, int matrixSize, int* matrixColSize) {int m matrixSize;if (m 0) {return 0;}int n matrixColSize[0];int left[m][n];memset(left, 0, sizeof(left));for (int i 0; i < m; i) {for (int j …

251 基于matlab的动态粒子群算法

基于matlab的动态粒子群算法。普通粒子群算法无法感知外界环境的变化&#xff0c;在外界环境发生改变时无法实时进行响应&#xff0c;因而缺乏动态环境寻优能力。在普通粒子群算法基本上通过增加敏感粒子得到一种动态粒子群算法&#xff0c;该算法通过实时计算敏感粒子的适应度…

Redis集群安装

将Redis安装包分别上传到3个文件夹&#xff0c;并解压缩 #编译并指定安装目录 cd /root/usr/local/redis-cluster/redis-7001/redis-6.2.6/ make make PREFIX/root/usr/local/redis-cluster/redis-7001 install # cd /root/usr/local/redis-cluster/redis-7002/redis-6.2.6/ m…

iZotope RX 11 for Mac 激活版:让您的音频焕发生机!

在追求音频完美的道路上&#xff0c;iZotope RX 11 for Mac是您的得力助手。它凭借先进的音频修复技术和丰富的音频增强工具&#xff0c;让您的音频作品焕发出前所未有的生机与活力。无论您是专业的音频工程师&#xff0c;还是业余的音乐爱好者&#xff0c;都能在这款工具中找到…

栈的讲解

栈的概念及结构 栈&#xff1a;一种特殊的线性表&#xff0c;其只允许在固定的一端进行插入和删除元素操作。 进行数据插入和删除操作的一端称为栈顶&#xff0c;另一端称为栈底&#xff08;因为先进后出&#xff09;。栈中的数据元素遵守后进先出LIFO&#xff08;Last In Firs…

【Threejs进阶教程-算法篇】1.常用坐标系介绍与2d/3d随机点位算法

2d/3d随机算法 学习ThreeJS的捷径坐标系简介平面直角坐标系和极坐标系空间直角坐标系圆柱坐标系球坐标系球坐标系与直角坐标系的转换 基于坐标系系统的随机点位算法平面直角坐标系随机平面直角坐标系随机的变形 空间直角坐标系随机二维极坐标系随机圆柱坐标系随机基于Cylinderc…

Python 整数类型(int)详解:无限范围与多种进制

引言 在编程中&#xff0c;整数是最基本的数据类型之一。不同编程语言对整数的处理方式各不相同&#xff0c;这往往影响到程序的性能和开发者的选择。本文将深入探讨 Python 中的整数类型&#xff08;int&#xff09;&#xff0c;其独特的处理方式&#xff0c;以及它在日常编程…