高精度
高精度适用情况
大数的长度最长为 10^6位
- 高精度加法通过对数位分别储存运算,并维护一个用于处理进位的变量,从而在不造成溢出的前提下实现大数相加。
- Java和Python中虽然不需要考虑大整数相加的溢出问题,但这种算法的思想很值得借鉴。
高精度储存
- 从高位对应着 数据中表示也是高位。
- 当然输出的时候需要倒叙输出
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-HGw3PIPX-1641802022731)(…/…/…/图库/1590506552684.png)]
int main()
{
// 使用字符串读入
string a, b;
vector<int> A, B;
cin >> a >> b; //a = "123456"
// 使用vector逆序读入,变成整数需要减去偏移量0
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');
// 相当于vector<int>
vector<int> C = add(A, B);
// 倒序输出
for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
return 0;
}
高精度加法
很简单。就是模拟竖式相加。注意:
- t为进位
- 从低到高位计算 c[i] = (A[i] + B[i] + t) % 10, 进位 t = (A[i] + B[i] + t) / 10;
- 计算完之后如果 t大于0,记得进位加到 c数组最后。
// C = A + B, A >= 0, B >= 0
vector<int> add(vector<int> &A, vector<int> &B)
{
if (A.size() < B.size()) return add(B, A); //始终保持长的数 + 短的数
vector<int> C; // vector C 储存和
int t = 0;
for (int i = 0; i < A.size(); i ++ )
{
t += A[i];//进位累加
if (i < B.size()) t += B[i]; // 数B没了
C.push_back(t % 10);
t /= 10;
}
if (t) C.push_back(t);//最高位如果进位
return C;
}
例题
https://www.acwing.com/problem/content/submission/code_detail/1607027/
#include <bits/stdc++.h>
using namespace std;
// C = A + B, A >= 0, B >= 0
vector<int> add(vector<int> &A, vector<int> &B)
{
if (A.size() < B.size()) return add(B, A); //始终保持长的数 + 短的数
vector<int> C; // vector C 储存和
int t = 0;
for (int i = 0; i < A.size(); i ++ )
{
t += A[i];//进位累加
if (i < B.size()) t += B[i]; // 数B没了
C.push_back(t % 10);
t /= 10;
}
if (t) C.push_back(t);//最高位如果进位
return C;
}
int main()
{
// 使用字符串读入
string a, b;
vector<int> A, B;
cin >> a >> b; //a = "123456"
// 使用vector逆序读入,变成整数需要减去偏移量0
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');
// 相当于vector<int>
vector<int> C = add(A, B);
// 倒序输出
for (int i = C.size() - 1; i >= 0; i--) printf("%d", C[i]);
return 0;
}
高精度减法,A - B(注意需要保证 A >= B)
A-B中两种情况:
A和B对位相减的话,两种情况:
- 从低位到高位计算 c[i] = (A[i] - B[i] - t) >= 0 ? (A[i] - B[i] - t) % 10 : ((A[i] - B[i] - t) + 10) % 10, 借位 t = (A[i] - B[i] - t) >= 0 ? 0 : 1
- 注意 : 1、答案有可b能是负数
2、最后存储在 vector 中的字符串可能会存在前导 0
#include <iostream>
#include <vector>
using namespace std;
const int N = 1000010;
//判断A是否比B大
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;
}
//去前导零
void trimZero(vector<int> &A)
{
while (A.back() == 0 && A.size() > 1) A.pop_back();
}
vector<int> sub(vector<int> &A, vector<int> &B)
{
vector<int> C;//C为A-B的差
int t = 0;
for (int i = 0; i < A.size(); i++)
{
t = A[i] - t;//减去借位
if (i < B.size()) t -= B[i]; //前两行实际是 t = A[i] - t - B[i]
C.push_back((t + 10) % 10);//加上借位,如果t为负,那需要借位的
if (t < 0) t = 1; //向前借位
else t = 0;
}
trimZero(C);//去除前导0
return C;
}
int main()
{
string a, b;
cin >> a >> b;
vector<int> A, B, C;
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 = sub(A, B); //若B比A大,将A-B变为B-A,加上负号
else {
C = sub(B, A);
printf("-");
}
for (int i = C.size() - 1; i >= 0; i--) cout << C[i];//倒序输出
return 0;
}
高精度乘法
以123 x 12 为例子
- 123的每一位都乘以 12
- c为每一位要输出的数字
- t为进位
- 可能有前导0,需要去除
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-dgIdxush-1641802022737)(…/…/…/图库/1590450263226.png)]
高精度乘以高精度
#include <iostream>
#include <vector>
using namespace std;
vector<int> mul(vector<int> &A, vector<int> &B) {
// 初始化为 0,且999*99最多 5 位
vector<int> C(A.size() + B.size(), 0);
for (int i = 0; i < A.size(); i++)
for (int j = 0; j < B.size(); j++)
C[i + j] += A[i] * B[j];
int t = 0;
for (int i = 0; i < C.size(); i++) {
// i = C.size() - 1时 t 一定小于 10
t += C[i];
C[i] = t % 10;
t /= 10;
}
// 必须要去前导 0,因为最高位很可能是 0
while (C.size() > 1 && C.back() == 0) C.pop_back();
return C;
}
int main() {
string a, b;
cin >> a >> b; // a = "1222323", b = "2323423423"
vector<int> A, B;
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');
vector<int> C = mul(A, B);
for (int i = C.size() - 1; i >= 0; i--)
cout << C[i]; //倒序输出
return 0;
}
高精度乘以低精度
#include <iostream>
#include <vector>
using namespace std;
vector <int> mul(vector <int> & A, int b) {
vector <int> C;
int t = 0;
for (int i = 0; i < A.size(); i ++) {
t += A[i] * b; // 例如:t + A[i] * b = 7218
C.push_back(t % 10); // 只取个位 8
t /= 10; // 721 看作 进位
}
while (t) { // 处理最后剩余的 t
C.push_back(t % 10);
t /= 10;
}
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');
vector <int> C = mul(A, b);
for (int i = C.size() - 1; i >= 0; i --) {
cout << C[i];
}
return 0;
}
高精度除法
需要带一个余数 r
下一次除的被除数就是 r*10 + A[i]
r = (r*10 + A[i])% b
C[i] =(r*10 + A[i])/ b
从高位开始除
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-iXXvXAe5-1641802022738)(…/…/…/图库/1590451533408.png)]
// A / b = C ... r, A >= 0, b > 0
vector<int> div(vector<int>& A, int& b, int& r){ // r引用
vector<int> C;
r = 0;//r为余数
for(int i = A.size() - 1; i >= 0; i --){ // 从高位除起
r = r * 10 + A[i];
C.push_back(r / b);//商
r %= b;
}
reverse(C.begin(), C.end()); // 翻转
while(C.size() > 1 && C.back() == 0) // 去除前导0
C.pop_back();
return C;
}
// 高精度除法
#include <iostream>
#include<vector>
#include<algorithm>
using namespace std;
// A被除数 b为除数 r为余数
vector<int> div(vector<int> A, int b, int &r) {
vector<int> C;
r = 0; //余数
for (int i = A.size() - 1; i >= 0; i--) { // 除法是从最高位开始算的
r = r * 10 + A[i];// 余数*10+当前这一位
C.push_back(r / b);
r %= b;
}
reverse(C.begin(), C.end()); //倒过来是因为 除法是从高位到地位
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 r; //余数
vector<int> C = div(A, b, r);
for (int i = C.size() - 1; i >= 0; i--) {
cout << C[i];
}
cout << endl << r << endl;
return 0;
}