「C++」论高精度
大家好,我是Charzie。在编程领域,高精度计算是一个常见的问题。当标准的整型或浮点型无法满足我们的计算需求时,高精度计算就显得尤为重要。在C++中,虽然标准库没有直接提供高精度数据类型,但我们可以通过一些技巧和工具类来实现高精度计算。
为什么需要高精度?
在编程中,我们经常会遇到一些大数计算的问题,比如大数的加减乘除、大素数判断等。这些问题的标准数据类型(如int、long long等)往往无法处理,因为它们有数据范围的限制。此时,我们就需要使用高精度计算来解决问题。
C++中实现高精度的方法
在C++中,实现高精度计算的方法主要有两种:使用数组模拟和使用高精度库。
使用数组模拟
我们可以使用数组来模拟高精度数。每个数组元素存储高精度数的一位,然后通过模拟手工计算的过程来实现高精度计算。这种方法虽然比较繁琐,但能够很好地理解高精度计算的原理。
例如,对于高精度加法,我们可以从低位到高位依次相加,注意处理进位问题。对于高精度乘法,我们可以使用竖式乘法的思想,将两个高精度数相乘的结果逐位计算并存储。
使用高精度库
为了简化高精度计算的实现,我们可以使用一些现成的高精度库,如GMP(GNU Multiple Precision Arithmetic Library)或Boost.Multiprecision库。这些库提供了丰富的高精度数据类型和运算函数,可以大大简化高精度计算的实现。
例如,使用GMP库,我们可以直接定义高精度整数类型,并进行加减乘除等运算,而无需关心底层的实现细节。
高精度计算的应用场景
高精度计算在很多领域都有广泛的应用,比如密码学、大数据分析、科学计算等。在密码学中,高精度计算可以用于实现公钥密码算法中的大数运算;在大数据分析中,高精度计算可以帮助我们处理海量的数据;在科学计算中,高精度计算可以提高计算的精度和可靠性。
高精度结构体:Bigint
这算是送给大家的一个礼物,是使用C++进行编写的结构体Bigint,集加、减、乘、除、输入、输出为一体。
#intclude<bits/stdc++.h>
struct Bigint {
std::vector<int> digits; // 存储每一位数字,低位在前
bool isNegative; // 标记是否为负数
Bigint() : isNegative(false) {}
// 构造函数,从字符串构造Bigint
Bigint(const std::string& str) : isNegative(false) {
if (str[0] == '-') {
isNegative = true;
for (size_t i = 1; i < str.length(); ++i) {
digits.push_back(str[i] - '0');
}
} else {
for (size_t i = 0; i < str.length(); ++i) {
digits.push_back(str[i] - '0');
}
}
std::reverse(digits.begin(), digits.end()); // 将数字反转,使得高位在后,方便计算
}
// 加法
Bigint operator+(const Bigint& b) const {
// 假设两个数都是正数,负数情况可以通过取反来处理
Bigint result;
int carry = 0;
size_t maxLength = std::max(digits.size(), b.digits.size());
for (size_t i = 0; i < maxLength || carry; ++i) {
int aDigit = i < digits.size() ? digits[i] : 0;
int bDigit = i < b.digits.size() ? b.digits[i] : 0;
int sum = aDigit + bDigit + carry;
carry = sum / 10;
result.digits.push_back(sum % 10);
}
while (result.digits.size() > 1 && result.digits.back() == 0) {
result.digits.pop_back(); // 去除前导零
}
return result;
}
// 减法(假设this >= b)
Bigint operator-(const Bigint& b) const {
Bigint result;
int borrow = 0;
for (size_t i = 0; i < digits.size(); ++i) {
int aDigit = digits[i];
int bDigit = i < b.digits.size() ? b.digits[i] : 0;
int diff = aDigit - bDigit - borrow;
if (diff < 0) {
diff += 10;
borrow = 1;
} else {
borrow = 0;
}
result.digits.push_back(diff);
}
while (result.digits.size() > 1 && result.digits.back() == 0) {
result.digits.pop_back(); // 去除前导零
}
return result;
}
// 乘法(简化版,未优化)
Bigint operator*(const Bigint& b) const {
Bigint result;
result.digits.resize(digits.size() + b.digits.size(), 0);
for (size_t i = 0; i < digits.size(); ++i) {
for (size_t j = 0; j < b.digits.size(); ++j) {
int product = digits[i] * b.digits[j];
result.digits[i + j] += product / 10; // 进位部分加到更高位上
result.digits[i + j + 1] += product % 10; // 个位数加到下一位上
}
}
for (size_t i = 0; i < result.digits.size() - 1; ++i) { // 处理进位
if (result.digits[i] >= 10) {
result.digits[i + 1] += result.digits[i] / 10;
result.digits[i] %= 10;
}
}
while (result.digits.size() > 1 && result.digits.back() == 0) {
result.digits.pop_back(); // 去除前导零
}
return result;
}
// 输入输出
friend std::istream& operator>>(std::istream& in, Bigint& b) {
std::string s;
in >> s;
b = Bigint(s);
return in;
}
friend std::ostream& operator<<(std::ostream& out, const Bigint& b) {
if (b.isNegative) out << '-';
for (auto it = b.digits.rbegin(); it != b.digits.rend(); ++it) {
out << *it;
}
return out;
}
};
// 示例使用
int main() {
Bigint a, b, sum, diff, prod;
std::cin >> a >> b;
sum = a + b;
diff = a - b; // 假设a >= b
prod = a * b;
std::cout << "Sum: " << sum << std::endl;
std::cout << "Difference: " << diff << std::endl;
std::cout << "Product: " << prod << std::endl;
return 0;
}
结语
高精度计算是编程领域的一个重要主题,它可以帮助我们解决标准数据类型无法处理的大数计算问题。在C++中,我们可以通过数组模拟或使用高精度库来实现高精度计算。随着技术的发展,高精度计算将在更多的领域得到应用,为我们的生活和工作带来更多的便利和可能性。
最后,我想说的是,高精度计算虽然有一定的难度,但只要我们掌握了基本的原理和方法,并多加练习和实践,就一定能够熟练掌握这一技能。希望这篇博客能够对大家有所帮助,让我们一起探索高精度计算的奥秘吧!
「C++」论高精度的更多相关文章
- 把虚拟教练带回家,「EuMotus」想用AI实现高精度运动反馈系统
https://36kr.com/p/5089139.html 无需穿戴设备,只需一个红外摄像头和+已安装好EuMotus专利软件的手提电脑 由政府主导的高达2200亿美金的健身与运动支出,15%的健 ...
- 「译」JUnit 5 系列:条件测试
原文地址:http://blog.codefx.org/libraries/junit-5-conditions/ 原文日期:08, May, 2016 译文首发:Linesh 的博客:「译」JUni ...
- 「译」JUnit 5 系列:扩展模型(Extension Model)
原文地址:http://blog.codefx.org/design/architecture/junit-5-extension-model/ 原文日期:11, Apr, 2016 译文首发:Lin ...
- JavaScript OOP 之「创建对象」
工厂模式 工厂模式是软件工程领域一种广为人知的设计模式,这种模式抽象了创建具体对象的过程.工厂模式虽然解决了创建多个相似对象的问题,但却没有解决对象识别的问题. function createPers ...
- 「C++」理解智能指针
维基百科上面对于「智能指针」是这样描述的: 智能指针(英语:Smart pointer)是一种抽象的数据类型.在程序设计中,它通常是经由类型模板(class template)来实做,借由模板(tem ...
- 「JavaScript」四种跨域方式详解
超详细并且带 Demo 的 JavaScript 跨域指南来了! 本文基于你了解 JavaScript 的同源策略,并且了解使用跨域跨域的理由. 1. JSONP 首先要介绍的跨域方法必然是 JSON ...
- 「2014-5-31」Z-Stack - Modification of Zigbee Device Object for better network access management
写一份赏心悦目的工程文档,是很困难的事情.若想写得完善,不仅得用对工具(use the right tools),注重文笔,还得投入大把时间,真心是一件难度颇高的事情.但,若是真写好了,也是善莫大焉: ...
- 「2014-3-18」multi-pattern string match using aho-corasick
我是擅(倾)长(向)把一篇文章写成杂文的.毕竟,写博客记录生活点滴,比不得发 paper,要求字斟句酌八股结构到位:风格偏杂文一点,也是没人拒稿的.这么说来,arxiv 就好比是 paper 世界的博 ...
- 「2014-3-17」C pointer again …
记录一个比较基础的东东-- C 语言的指针,一直让人又爱又恨,爱它的人觉得它既灵活又强大,恨它的人觉得它太过于灵活太过于强大以至于容易将人绕晕.最早接触 C 语言,还是在刚进入大学的时候,算起来有好些 ...
- 「2014-3-13」Javascript Engine, Java VM, Python interpreter, PyPy – a glance
提要: url anchor (ajax) => javascript engine (1~4 articles) => java VM vs. python interpreter =& ...
随机推荐
- 基础 IO (Linux学习笔记)
基础IO 1.重谈文件 空文件在磁盘也要占据空间 文件 = 内容 + 属性 文件操作 = 对文件内容+对属性 or 对文件内容加属性 标定一个文件,必须使用文件路径加文件名[唯一性] 如果没有指明对应 ...
- JavaScript中字符串小知识
1. 字符串是不可变的 字符串一旦创建就是不可变的,后续的修改都是新建一个新的字符串而不是在原有的字符串上修改 // 在内存中开辟 可以存放五个字母的空间 str指向该空间 let str = 'fi ...
- 力扣661(java)-图片平滑器(简单)
题目: 图像平滑器 是大小为 3 x 3 的过滤器,用于对图像的每个单元格平滑处理,平滑处理后单元格的值为该单元格的平均灰度. 每个单元格的 平均灰度 定义为:该单元格自身及其周围的 8 个单元格的 ...
- 一文详解用eBPF观测HTTP
简介: 随着eBPF推出,由于具有高性能.高扩展.安全性等优势,目前已经在网络.安全.可观察等领域广泛应用,同时也诞生了许多优秀的开源项目,如Cilium.Pixie等,而iLogtail 作为阿里内 ...
- 为了让你在“口袋奇兵”聊遍全球,Serverless 做了什么?
简介: 江娱互动是一家新兴的游戏企业,自 2018 年成立伊始,江娱互动就面向广阔的全球游戏市场,通过创造有趣的游戏体验,在竞争激烈的游戏市场占得一席之地.仅仅 2 年的时间,江娱互动就凭借 Topw ...
- 一文详解 | 开放搜索兼容Elasticsearch做召回引擎
简介:开放搜索发布开源兼容版,支持阿里云Elasticsearch做搜索召回引擎,本文详细介绍阿里云ES用户如何通过接入开放搜索兼容版丰富行业分词库,提升查询语义理解能力,无需开发.算法投入,即可获 ...
- [Trading] 买卖如何移动 ( 影响 ) 市场价格
大多数人都知道市场价格的变化是因为买卖行为,但却没有多少人了解买卖行为是如何影响市场价格的. 乍一看,这可能令人困惑,因为每一笔市场交易都要求总有一个买家和一个卖家. 首先,重要的是要明白市场上总是有 ...
- Raft 共识算法2-领导者选举
Raft 共识算法2-领导者选举 Raft算法中译版地址:https://object.redisant.com/doc/raft中译版-2023年4月23日.pdf 英原论文地址:https://r ...
- Java根据URL截图的4种方式
方案选择 XHTMLRenderer(不要用) PhantomJs(三方库,已停更) Puppeteer(Chrome团队开发和维护) Selenium(支持多浏览器.多语言,服务器需要安谷歌浏览器) ...
- Java开启异步的两种方式
二.Java开启异步的两种方式 1.注解开启:@Async 1.1.配置异步的线程池 必须配置异步线程池,否则异步不会生效. @EnableAsync 注解:指定异步线程池.不指定默认使用:Simpl ...