文章目录
前言
复习acwing算法基础课的内容,本篇为讲解基础算法:高精度运算,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上
一、高精度运算概述
我们都知道,我们定义一个数,不论是整数还是浮点数,它都是有范围的都是有限的,就拿unsigned long long
为例,它的最大值也就1e19,意思是最多只有19位,而我们想要进行40位,50位甚至更多位数的运算,显然就不能简单的定义int
,long long
等类型进行计算了,这时候我们就可以运用到高精度算法,把每一位上的数都存到数组之中,本博客的高精度算法中,用到了STL中的vector
,关于vector
,详见:vector
二、高精度算法
1. 高精度加法
AcWing791. 高精度加法
本题链接:AcWing791. 高精度加法
本博客给出题目截图:
高精度加法模板
本模板来自ACWing算法基础课
vector<int> add(vector<int> &A, vector<int> &B)
{
if (A.size() < B.size()) return add(B, A); //保证A的长度始终 ≥ B的长度,因为下面的for循环是通过遍历A的每个元素
vector<int> C;
//C中存储的就是A + B后的结果
int t = 0;
for (int i = 0; i < A.size(); i ++ )
{
t += A[i];
if (i < B.size()) t += B[i];
C.push_back(t % 10);
t /= 10;
}
if (t) C.push_back(t); //有可能最高位加法存在进位,如果t存在,就加入到C
return C;
}
AC代码
#include <iostream>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
vector<int> add(vector<int> &A, vector<int> &B)
{
if (A.size() < B.size()) return add(B, A);
vector<int> C;
int t = 0;
for (int i = 0; i < A.size(); i ++ )
{
t += A[i];
if (i < B.size()) t += B[i];
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;
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');
auto C = add(A, B);
for (int i = C.size() - 1; i >= 0; i -- ) cout << C[i];
cout << endl;
return 0;
}
2.高精度减法
AcWing 792. 高精度减法
本题链接:AcWing 792. 高精度减法
本博客给出题目截图:
高精度减法模板
本模板来自ACWing算法基础课
// C = A - B, 满足A >= B, A >= 0, B >= 0
vector<int> sub(vector<int> &A, vector<int> &B)
{
vector<int> C;
for (int i = 0, t = 0; i < A.size(); i ++ )
{
t = A[i] - t;
if (i < B.size()) t -= B[i];
C.push_back((t + 10) % 10); //因为t可能是负的
//举个例子,我们计算10 - 8,对于个位而言就是-8,我们通过(t + 10) % 10的方法就可以得到2
if (t < 0) t = 1; //证明存在借位,即像高位借1
else t = 0;
}
while (C.size() > 1 && C.back() == 0) C.pop_back(); //去除前导0
//比如111111118 - 111111117 = 000000001,而我们要输出1,通过这一步骤可以实现,要求C.size() > 1是因为有可能答案就是0
return C;
}
AC代码
保证始终是绝对值大的数去减绝对值小的数:
#include <iostream>
#include <cstring>
#include <vector>
using namespace std;
vector <int> A, B, C;
bool cmp(vector <int> &A, vector <int> &B) //比较A和B谁大
{
if (A.size() != B.size()) return A.size() > B.size();
else
{
for (int i = A.size() - 1; i >= 0; i -- )
if (A[i] != B[i]) return A[i] > B[i];
}
return true;
}
void sub(vector <int> &A, vector <int> &B)
{
int t = 0;
for (int i = 0; i < A.size(); i ++ )
{
t = A[i] - t;
if (i < B.size()) t -= B[i];
C.push_back((t + 10) % 10);
if (t < 0) t = 1;
else t = 0;
}
while (C.size() > 1 && C.back() == 0) C.pop_back();
}
int main()
{
string a, b;
cin >> 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');
if (cmp(A, B)) sub (A, B);
else
{
sub (B, A);
printf("-");
}
for (int i = C.size() - 1; i >= 0; i -- ) cout << C[i];
return 0;
}
3.高精度乘法
AcWing 793. 高精度乘法
本题链接:AcWing 793. 高精度乘法
本博客给出题目截图:
高精度乘法模板
本模板来自ACWing算法基础课
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;
}
while (C.size() > 1 && C.back() == 0) C.pop_back(); //去除前导0
//比如 123 * 0,按照这个运算的话是 0000,答案是0,故需要去除前导0
return C;
}
AC代码
#include <iostream>
#include <vector>
#include <cstring>
using namespace std;
vector <int> A, B;
void mul(vector <int> &A, int b)
{
int t = 0;
for (int i = 0; i < A.size() || t; i ++ )
{
t += A[i] * b;
B.push_back(t % 10);
t /= 10;
}
while (B.size() > 1 && B.back() == 0) B.pop_back();
}
int main()
{
string a;
int b;
cin >> a >> b;
for (int i = a.size() - 1; i >= 0; i -- ) A.push_back(a[i] - '0');
mul(A, b);
for (int i = B.size() - 1; i >= 0; i -- ) cout << B[i];
return 0;
}
4.高精度除法
AcWing 794. 高精度除法
本题链接:AcWing 794. 高精度除法
本博客给出本题截图:
高精度除法模板
本模板来自ACWing算法基础课
vector<int> div(vector<int> &A, int b, int &r) //传入的是r的地址,所以可以直接对r的值进行修改
{
vector<int> C;
r = 0;
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());
//除法和加法,减法,乘法运算都不一样的是除法运算是高位向低位进行,而其他三种运算是从低位到高位运算
//所以除法运算的前导0其实是在开头而不是尾部,而我们只能从尾部删除一个元素,所以我们要把C翻转
while (C.size() > 1 && C.back() == 0) C.pop_back(); //去除前导0
return C;
}
AC代码
#include <iostream>
#include <vector>
#include <cstring>
#include <algorithm>
using namespace std;
vector <int> A, C;
int r;
void div(vector <int> A, int b)
{
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) C.pop_back();
}
int main()
{
string a;
int b;
cin >> a >> b;
for (int i = a.size() - 1; i >= 0; i -- ) A.push_back(a[i] - '0');
div(A, b);
for (int i = C.size() - 1; i >= 0; i -- ) cout << C[i];
cout << endl << r;
return 0;
}
三、时间复杂度
关于高精度算法的时间复杂度以及证明,后续会给出详细的说明以及证明过程,目前先鸽了。