高精度算法模板
高精度一般是对于大于 10 的 ±6 次方的计算,其数据范围大大超过了标准数据类型能表示的范围的运算,使用高精度来帮助运算。
高精度运算类似于 竖式运算 ,将字符串转化为 数组 ,其数组的每一位元素对应着十进制数,满足十进制的 “满十进一” 的运算法则。
需要注意的是,高精度数组存储采用 倒序存储 ,从低位向高位逐位来计算,输出时也是采用倒序输出。
一、高精度整数加法
给定两个正整数(不含前导0),计算它们的和。
(C = A + B,A >= 0, B >= 0)
代码如下:
#include <iostream>
#include <vector>
using namespace std;
const int N = 1e6 + 10;
// C = A + B
vector<int> add(vector<int> &A, vector<int> &B) //加上&直接搜索数组A和B,不用全部遍历,节省时间
{
vector<int> C;
int t = 0; //进位,低位满10向高位进位,低位变为0
for (int i = 0; i < A.size() || i < B.size(); i++){
//两个if把两个数组相同位上的数相加
if (i < A.size())
t += A[i];
if (i < B.size())
t += B[i];
//把相加后的结果除以10求余,压入C数组中
C.push_back(t % 10);
//t再除以10,放入高位
t /= 10;
}
//如果最高位有数,则压入C数组中
if (t)
C.push_back(t);
//最后返回C数组
return C;
}
int main(){
string a, b;
vector<int> A, B;
cin >> a >> b; //例如 a = "123456"
for (int i = a.size() - 1; i >= 0; i--)
A.push_back(a[i] - '0'); //倒序存放, 数组A = [6, 5, 4, 3, 2, 1]
for (int i = b.size() - 1; i >= 0; i--)
B.push_back(b[i] - '0'); //同理, 数组B也是从低位(个位)存储, 从小到大
//使用auto编译器会自动判断数据是什么类型
auto C = add(A, B);
for (int i = C.size() - 1; i >= 0; i--)
printf("%d", C[i]);
return 0;
}
二、高精度整数减法
给定两个正整数(不含前导0),计算它们的差,计算结果可能为负数。
1、前提:
C = A - B,满足 A >= B,A >= 0,B >= 0 (如果 B >= A,则将其转化为 -(B - A) )
2、代码段:(只考虑正数的情况)
#include <iostream>
#include <vector>
using namespace std;
//判断是否有 A >= B
bool cmp(vector<int> &A, vector<int> &B)
{
//A长度不等于B,若A大于B,返回true,否则返回false
if (A.size() != B.size())
return A.size() > B.size();
//因为是倒序数组,高位在后面,所以从后遍历
//如果A > B, 返回true,否则返回flase
for (int i = A.size() - 1; i >= 0; i--){
if (A[i] != B[i])
return A[i] > B[i];
}
return true;
}
// C = A - B
vector<int> sub(vector<int> &A, vector<int> &B) //加上&直接搜索数组A和B,不用全部遍历,节省时间
{
vector<int> C;
//定义一个借位的情况
int t = 0;
for (int i = 0; i < A.size(); i++){
//t = A的这一位数
t = A[i] - t;
//判断B在这一位上有没有数,如果有,则减去B[i]
if (i < B.size())
t -= B[i];
//这种写法包含了两种情况:如果t >= 0, 则直接减;若t < 0, 则向高位借位(即t + 10)后再减
//将这一位相减后的结果压入C中
C.push_back((t + 10) % 10);
if (t < 0) //t < 0,需要借位,标记为1
t = 1;
else //t >= 0,不需要借位,标记为0
t = 0;
}
//删除前导0
while (C.size() > 1 && C.back() == 0)
C.pop_back();
//最后返回C数组
return C;
}
int main(){
string a, b;
vector<int> A, B;
cin >> a >> b; //例如 a = "123456"
for (int i = a.size() - 1; i >= 0; i--)
A.push_back(a[i] - '0'); //倒序存放, 数组A = [6, 5, 4, 3, 2, 1]
for (int i = b.size() - 1; i >= 0; i--)
B.push_back(b[i] - '0'); //同理, 数组B也是从低位(个位)存储, 从小到大
//如果A > B,返回A - B
if (cmp(A, B)){
auto C = sub(A, B);
for (int i = C.size(); i >= 0; i--)
printf("%d", C[i]);
}
//否则返回 -(B - A)
else {
auto C = sub(B, A);
printf("-");
for (int i = C.size(); i >= 0; i--)
printf("%d", c[i]);
}
return 0;
}
三、高精度整数乘法
给定两个非负整数(不含前导0)A 和 B,要求计算 A × B 的值
一般是 高精 × 低精,用 A × b 表示
(C = A * b,A >= 0,b >= 0)
代码如下:
#include <iostream>
#include <vector>
using namespace std;
//C = A * b
vector<int> mul(vector<int> &A, int b)
{
vector<int> C;
int t = 0;
for (int i = 0; i < A.size() || t; i++)
{
if (i < A.size())
t += A[i] * b;
//同加法一样处理进位,逐位压入数组中
C.push_back(t % 10);
t /= 10;
}
//删除前导0
while (C.size() > 1 && C.back() == 0)
C.pop_back();
return C;
}
int main()
{
string a;
int b;
cin >> a >> b;
vector<int> A;
for (int i = a.size() - 1; i >= 0; i--) //倒序插入
A.push_back(a[i] - '0');
auto C = mul(A, b);
for (int i = C.size() - 1; i >= 0; i--) //倒序输出
printf("%d", C[i]);
return 0;
}
四、高精度整数除法
给定两个非负整数(不含前导0)A 和 B,要求计算 A / B 的商和余数 (第一行输出所求的商,第二行输出所求余数)
一般是 高精 ÷ 低精,用 A ÷ b 表示
(A / b = C ··· r,A >= 0,b > 0)
代码如下:
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
vector<int> div(vector<int> &A, int b, int &t) //传入t的地址,便于直接对余数进行处理
{
vector<int> C;
t = 0;
for (int i = A.size() - 1; i >= 0; i--)
{
//将上次的余数×10再加上当前位的数字,得到该位的被除数
t = t * 10 + A[i];
//所得即为商在这一位的数字
C.push_back(t / b);
t %= b;
}
//由于在除法运算中,从高位到低位运算,因此前导0在数组前面,所以需要将其翻转,将前导0置于尾部,从而便于删除前导0
reverse(C.begin(), C.end());
//删除前导0
while (C.size() > 1 && C.back() == 0)
C.pop_back();
return C;
}
int main()
{
string a;
int B;
cin >> a >> B;
vector<int> A;
for (int i = a.size() - 1; i >= 0; i--)
A.push_back(a[i] - '0');
int t; //t为余数
auto C = div(A, B, t);
for (int i = C.size() - 1; i >= 0; i--)
cout << C[i];
cout << endl << t << endl;
return 0;
}