高精度算法一

news2024/9/24 9:25:26

目录

1. 基础知识

2. 大整数 + 大整数

3. 大整数 - 大整数


1. 基础知识

利用计算机进行数值计算,有时会遇到这样的问题:有些计算要求精度高,希望计算的数的位数可达几十位甚至几百位,虽然计算机的计算精度也算较高了,但因受到硬件的限制,往往达不到实际问题所要求的精度。我们可以利用程序设计的方法去实现这样的高精度计算。

这里我们要介绍:大整数 + 大整数,大整数 - 大整数,大整数 * 小整数,大整数 / 小整数。

大整数显然用简单的数据类型是无法表示的。所以我们选用数组来存储大整数的每一位。于是,就产生了这么一个问题:数组下标为0的元素是存最高位还是最低位呢?

显然,我们在运算过程中可能会发生进位,这时,如果我们用第一种(下标为0存最高位),在运算过程中最高位发生进位时,就需要在数组的前面添加一个新的位置来存放产生的进位。这就需要将数组的元素整体向后移动,这很麻烦。所以我们选用第二种方式来存储大整数(下标为0存最低位),这样在产生进位时我们只需要在数组的后面添加一个位置即可。这是非常容易滴。

2. 大整数 + 大整数

我们可以从人类计算加法的方法中找出做加法的规律:

下面以 345678 + 6789 距离分析:我们将存储进位的变量 carry 初始化为0,然后让 carry 去与两个大整数的个位相加,即 0 + 8 + 9 = 17,将得到的结果 17 对 10 取模,得到相加后的个位:7,push_back到存放结果的数组,然后继续将相加后的结果17除以10,作为下一位相加的进位。依次类推,直到将两个大整数的每一位相加。

#include<iostream>
#include<vector>
using namespace std;

//传引用可以提高效率,不用copy一遍数组了
vector<int> add(vector<int>& A, vector<int>& B)
{
	//存放结果的数组
	vector<int> C;

	//用来记录进位的变量,初始化为0
	int carry = 0;

	//对两个大整数的对应位置的数相加
	for (int i = 0; i < A.size() || i < B.size(); i++)
	{
		//因为两个大整数的长度可能有差异,当一个大整数的每一位加完了,就不需要继续了
		if (i < A.size())
			carry += A[i];
		if (i < B.size())
			carry += B[i];

		C.push_back(carry % 10);
		carry /= 10;
	}

	//如果最高位相加有进位(只可能是1,push_back 1 也行),再将进位push_back
	if (carry)
	{
		C.push_back(carry);
	}

	return C;

}

void test01()
{
	//两个大整数以字符串的形式输入
	string a, b;
	cout << "请输入两个整数,以回车键结束" << endl;
	cin >> a >> b;

	//两个数组用来存放大整数的每一位
	vector<int> A, B;

	//将输入的大整数的每一位拆分后添加到数组
	//数字字符 - '0' 即可得到数字字符对应的数字
	for (int i = a.size() - 1; i >= 0; i--)
	{
		A.push_back(a[i] - '0');
	}
	for (int i = b.size() - 1; i >= 0; i--)
	{
		B.push_back(b[i] - '0');
	}

	//用一下大整数加法的接口,并且用C接收最终结果
	vector<int>C = add(A, B);


	cout << "相加的结果:";
	//输出结果
	for (int i = C.size() - 1; i >= 0; i--)
	{
		printf("%d", C[i]);
	}
	cout << endl;
}

int main()
{

	//大整数的加法
	test01();

	system("pause");
	return 0;
}

3. 大整数 - 大整数

大整数减去大整数用户可能输入的是一个较小的大整数减去一个较大的大整数,我们就需要先进行一个比较,确保我们在做减法时是大的数减去小的数,如果需要在最终的结果上加上负号即可。

与大整数加大整数类似,对两个大整数的每一位做减法,显然这时有两种情况:当减数小于被减数,此时就需要向前借一位,而当减数大于等于被减数则不需要进位。

 这里还有一个问题就是:相减的结果在打印时可能会有先导0,例如:12345 - 12300 在没有去掉先导0,得到的结果就是:00045,这显然不是我们想要的。我们需要对存放结果的数组的最后一个元素进行判断,如果该元素为0,就需要将其pop掉。但是如果相减的最终结果是0,这样做存放结果的数组里面将没有元素了,所以我们需要保证存放结果的数组中至少有一个元素。

bool cmp(vector<int>& A, vector<int>& B)
{
	//先比较两整数的长度
	if (A.size() != B.size())
	{
		return A.size() - B.size();
	}
	//从最高位开始比较
	for (int i = A.size() - 1; i >= 0; i--)
	{
		if (A[i] != B[i])
			return A[i] > B[i];
	}

	//这里代表两个数相等
	return true;
}


vector<int> sub(vector<int>& A, vector<int>& B)
{
	vector<int> C;
	

	int borrow = 0;
	for (int i = 0; i < A.size(); i++)
	{
		//减数减去借位
		borrow = A[i] - borrow;
		//当i不小于被减数的大小就不用做减法了
		if (i < B.size())
			borrow -= B[i];

		//利用一般规律push_back结果
		C.push_back((borrow + 10) % 10);

		//有借位borrow置为1
		if (borrow < 0)
			borrow = 1;
		else
			borrow = 0;
	}
	//去掉先导0
	while (C.size() > 1 && C.back() == 0)
	{
		C.pop_back();
	}

	return C;
}


void test02()
{
	//两个大整数以字符串的形式输入
	string a, b;
	cout << "请输入两个整数,以回车键结束" << endl;
	cin >> a >> b;

	//两个数组用来存放大整数的每一位
	vector<int> A, B;

	//将输入的大整数的每一位拆分后添加到数组
	//数字字符 - '0' 即可得到数字字符对应的数字
	for (int i = a.size() - 1; i >= 0; i--)
	{
		A.push_back(a[i] - '0');
	}
	for (int i = b.size() - 1; i >= 0; i--)
	{
		B.push_back(b[i] - '0');
	}

	//对两个大整数做判断,确保是大数 - 小数
	if (cmp(A, B))
	{
		//用一下大整数加法的接口,并且用C接收最终结果
		vector<int>C = sub(A, B);

		cout << "相减的结果:";
		//输出结果
		for (int i = C.size() - 1; i >= 0; i--)
		{
			printf("%d", C[i]);
		}
		cout << endl;
	}
	else
	{
		//用一下大整数加法的接口,并且用C接收最终结果
		//如果B > A,交换位置传参即可,最后输出数据的时候加上负号
		vector<int>C = sub(B, A);


		cout << "相减的结果:";
		//输出结果
		cout << "-";
		for (int i = C.size() - 1; i >= 0; i--)
		{
			printf("%d", C[i]);
		}
		cout << endl;
	}

}

int main()
{

	//大整数的减法
	test02();

	system("pause");
	return 0;
}

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

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

相关文章

大数据之---Nifi-处理器类别_详细了解各种处理器的含义用途---大数据之Nifi工作笔记0007

处理器有很多类这里主要分了,数据转换的可以看到上面的这些是数据转换常用的 然后: 这里这个controlrate,用来控制流的速率,来保护作用,防止崩掉 distributeLoad用来负载均衡,减轻压力 这里的路由是,根据是否满足对应的条件,来执行是否把数据路由到某个地方. 然后来看数据库…

数据库技术-函数依赖、键与约束、范式

一、函数依赖 给定一个x&#xff0c;能唯一确定一个Y&#xff0c;就称x确定Y&#xff0c;或者说Y依赖于x&#xff0c;例如YX*X函数。 函数依赖又可扩展以下两种规则: 部分函数依赖:A可确定C&#xff0c;(A,B)也可确定C,(A,B)中的一部分&#xff08;即A&#xff09;可以确定C&a…

Spring Cloud Nacos实战(一)- 下载和安装

Spring Cloud Alibaba Nacos下载和安装 Nacos介绍 ​ Nacos&#xff08;Naming Configuration Service&#xff09; 是一个易于使用的动态服务发现、配置和服务管理平台&#xff0c;用于构建云原生应用程序 ​ 服务发现是微服务架构中的关键组件之一。Nacos 致力于帮助您发现…

使用nginx单独部署Vben应用

前言 本文主要介绍Vben使用nginx单独部署的方式&#xff0c;其实前端发展到现在已经不是当年的jsp&#xff0c;asp必须要和后端一起部署了。单独部署调试的工具也很多&#xff0c;比如vue-cli-service 和 Vben中用到的vite &#xff0c;当然这些我们一般用在开发的工程中。正式…

SpringCloud(三)Hystrix断路器服务降级、服务熔断、服务监控案例详解

七、Hystrix断路器 7.1 简介 分布式系统面临的问题 复杂分布式体系结构中的应用程序有数十个依赖关系&#xff0c;每个依赖关系在某些时候将不可避免地失败。 多个微服务之间调用的时候&#xff0c;假设微服务A调用微服务B和微服务C&#xff0c;微服务B和微服务C又调用其它的微…

C# chart绘图 鼠标响应

1、图形自动滚动设置 chart1.ChartAreas[0].AxisX.Maximum 横坐标显示区域最大值 chart1.ChartAreas[0].AxisX.Minimum 横坐标显示区域最小值 显示宽度 chart1.ChartAreas[0].AxisX.Maximum - chart1.ChartAreas[0].AxisX.Minimum chart1.ChartAreas[0].AxisX.Maximum x_d…

ChatGPT从业影响之个人浅见

文章目录前言ChatGPT简介AI简介ChatGPT对从业者的影响参考资料前言 ChatGRT最近很火&#xff0c;网上有一些观点&#xff0c;他会取代很多职业&#xff0c;让产业产生颠覆性的变化。今天聊聊自己的想法&#xff0c;他的出现&#xff0c;是否会让我&#xff08;程序员&#xff…

qmt api 接口文档

http://qmt.ptradeapi.com 用户自行安装 Python 三方库 对于有经验的 Python 开发者来说&#xff0c;平台提供了自行安装第三方库的方式。为了引入额外的第三方库&#xff0c;用户需要做如下一些操作&#xff1a; 安装前注意事项&#xff1a; 三方库的安装有可能会引起系统错误…

JVM学习篇对象创建与内存分配机制深度剖析

1. 对象的创建 对象创建的主要流程: 1.1 类加载检查 虚拟机遇到一条new指令时&#xff0c;首先将去检查这个指令的参数是否能在常量池中定位到一个类的符号引用&#xff0c;并且检查这个符号引用代表的类是否已被加载、解析和初始化过。如果没有&#xff0c;那必须先执行相…

【Java 面试合集】内存中为什么要区分栈和堆

内存中为什么要区分栈和堆 1. 概述 嗨&#xff0c;大家好Java 面试合集又来了&#xff0c;今天我们分享的主题很大众化以及普通。无论是Java 语言本身还是别的语言都会有所涉及&#xff0c;但是今天我们从Java角度来讲下 2. 分析 今天我们会从多个方向来分享这个话题 2.1 栈是…

ETCD本地多成员集群部署

目录安装部署etcdctl 操作etcd使用http请求操作etcd本地多成员集群搭建python获取成员信息参考安装部署 按照官网文档&#xff0c;安装release版本 https://etcd.io/docs/v3.4/install/ [rootVM-33-162-centos /usr/local/bin]# etcd --version etcd Version: 3.4.16 Git SH…

点云配准方法原理(NDT、ICP)

配准是点云处理中的一个基础问题&#xff0c;众多学者此问题进行了广泛而深入的研究&#xff0c;也出现了一系列优秀成熟的算法&#xff0c;在三维建模、自动驾驶等领域发挥着重要的作用。 本文主要介绍粗配准NDT (Normal Distribution Transform) 与 精配准ICP (Iterative Cl…

最新最全中小微企业研究数据:海量创业公司信息与获取投资信息(1985-2021年)

一、企业获取投资名单&资方信息 数据来源&#xff1a;搜企网、企查查、天眼查 时间跨度&#xff1a;1985年8月-2021年9月 区域范围&#xff1a;全国范围 数据字段&#xff1a;企业名称、时间、获得投资金额以及投资方信息 部分数据&#xff1a; DateCompany_nameUnit…

摄影师没了?!生成式人工智能即将降维打击摄影行业

本文是Mixlab无界社区成员的投稿&#xff1a;滚石deepfacelab和deepfacelive项目组成员摄影师失业了&#xff1f;&#xff1f;怎么说&#xff1f;##你还以为AI绘画影响的只是插画师行业吗&#xff1f;错了&#xff0c;摄影行业也即将面临技术洗牌。话不多说&#xff0c;先看一下…

java并发编程原理2 (AQS, ReentrantLock,线程池)

一、AQS&#xff1a; 1.1 AQS是什么&#xff1f; AQS就是一个抽象队列同步器&#xff0c;abstract queued sychronizer&#xff0c;本质就是一个抽象类。 AQS中有一个核心属性state&#xff0c;其次还有一个双向链表以及一个单项链表。 首先state是基于volatile修饰&#x…

分享113个HTML艺术时尚模板,总有一款适合您

分享113个HTML艺术时尚模板&#xff0c;总有一款适合您 113个HTML艺术时尚模板下载链接&#xff1a;https://pan.baidu.com/s/1ReoPNIRjkYov-SjsPo0vhg?pwdjk4a 提取码&#xff1a;jk4a Python采集代码下载链接&#xff1a;采集代码.zip - 蓝奏云 女性化妆用品网页模板 粉…

【Linux】用户分类+权限管理+umask+粘滞位说明

目录 1.用户分类 su指令 2.认识Linux权限 2.1 文件访问者的分类 2.2 文件类型和访问权限 a. 文件类型 file指令 b. 访问权限 2.3 文件权值的表示方法 a. 字母表示法 b. 八进制表示法 3.如何修改文件访问者的权限及相关指令 1. chmod指令 2. chown指令 3. chgrp指…

Python语言零基础入门教程(二十七)

Python OS 文件/目录方法 Python语言零基础入门教程&#xff08;二十六&#xff09; 61、Python os.utime() 方法 概述 os.utime() 方法用于设置指定路径文件最后的修改和访问时间。 在Unix&#xff0c;Windows中有效。 语法 utime()方法语法格式如下&#xff1a; os.uti…

SortableJS/Sortable拖拽组件,使用详细(Sortablejs安装使用)

简述 作为一名前端开发人员&#xff0c;在工作中难免会遇到拖拽功能&#xff0c;分享一个github上一个不错的拖拽js库&#xff0c;能满足我们在项目开发中的需要&#xff0c;支持Vue和React&#xff0c;下面是SortableJS的使用详细&#xff1b; 这个是sortableJS中文官方文档&…

kafka-6-python单线程操作kafka

使用Python操作Kafka&#xff1a;KafkaProducer、KafkaConsumer Python kafka-python API的帮助文档 1 kafka tools连接 (1)/usr/local/kafka_2.13-3.4.0/config/server.properties listeners PLAINTEXT://myubuntu:9092 advertised.listenersPLAINTEXT://192.168.1.8:2909…