高精度算法的意义
先来看一组数据:
int类型的范围:-2147483648~2147483647
long long类型的范围:-9223372036854775808~9223372036854775807
一般地,我们对于整数的存储都在这范围内,但是在某些情况(比如斐波那契数列数列递推到后面数值过于庞大)下还是无法满足。而C++又不会像Python一样自带高精度,这个时候,便需要我们自行编写高精度算法。高精度算法通常有高精度与高精度(或单精度)的加法、减法、乘法。除法在我的博客里不做讲解,一般用不到.
高精度算法的通用思想
1.字符串 倒序 存储至数组
2.按照一般的加、减、乘法进行计算
3.处理进位或退位
4.处理前导零及因进退位引起的长度变化
5.再次倒序输出数组元素(即数字的正序)
今天讲述高精度加单精度
数据准备
const int MAXN = 205; //数字的位数(可更改)
string num1; //高精度数字
int num2, len; //单精度数字;高精度数字的长度
int a[MAXN]; //用于高精度数字存储并计算结果
倒序存储
a数组的正序下标 | 高精度数字的倒序下标 |
---|---|
0 | len - 1 |
1 | len - 2 |
2 | len - 3 |
3 | len - 4 |
到这里我们就可以直接推导出高精度数字的倒序下标等于len - i - 1
for (int i = 0; i < len; i++) {
a[i] = num1[len - i- 1] - '0';
}
一般加法
没错,就直接在a[0](也就是竖式的个位)加上单精度数
a[0] += num2;
处理进位
如图,是一个竖式过程,9 + 9 = 18
,很明显,需要进位,我们将18中的8写在了当前位,然后为更高一位进了1,,这是我们的竖式,如何用代码实现呢?显然当前位的8是a[i] %= 10
,而往前进的一位则是a[i] / 10
for (int i = 0; i < len; i++) {
a[i + 1] += a[i] / 10;
a[i] %= 10;
}
处理进位引起的长度变化
由于可能出现最高位进位的情况,但我们的进位处理只针对一般情况而处理到len - 1
位。我们需要通过while
循环来进行最高位及往上的进位。从第len
位开始。
while (a[len]) {
a[len + 1] += a[len] / 10; //和上面的进位原理一样
a[len] %= 10;
len++; //最高位进位引起了长度变化,要加1
}
输出
见代码,无难处。
for (int i = len - 1;i >= 0; i--) {
cout << a[i];
}
完整代码
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
//高精度加单精度
const int MAXN = 205;
string num1;
int num2, len;
int a[MAXN];
int main() {
cin >> num1 >> num2;
len = num1.size();
for (int i = 0;i < len; i++) {
a[i] = num1[len - i- 1] - '0';
}
a[0] += num2;
for (int i = 0; i < len; i++) {
a[i + 1] += a[i] / 10;
a[i] %= 10;
}
while (a[len]) {
a[len + 1] += a[len] / 10;
a[len] %= 10;
len++;
}
for (int i = len - 1;i >= 0; i--) {
cout << a[i];
}
return 0;
}