定义
高精度算法(High Accuracy Algorithm)是处理大数字的数学计算方法。在一般的科学计算中,会经常算到小数点后几百位或者更多,当然也可能是几千亿几百亿的大数字。一般这类数字我们统称为高精度数,高精度算法是用计算机对于超大数据的一种模拟加,减,乘,除,乘方,阶乘,开方等运算。对于非常庞大的数字无法在计算机中正常存储,于是,将这个数字拆开,拆成一位一位的,或者是四位四位的存储到一个数组中, 用一个数组去表示一个数字,这样这个数字就被称为是高精度数。高精度算法就是能处理高精度数各种运算的算法,但又因其特殊性,故从普通数的算法中分离,自成一家。
代码
本模板的数据输入
int main(){
//读入
string A,B;
int b;
vector<int> AI,BI;
cin>>A>>B;
for(int i = A.length()-1;i>=0;i--){
AI.push_back(A[i]-48);
}
for(int i = B.length()-1;i>=0;i--){
BI.push_back(B[i]-48);
}
// gjd_add 加法
vector<int> C = gjd_add(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_sub 减法
C = gjd_sub(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_mul 乘法
C = gjd_mul(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_div 除法
vector<int> r; // 高精度除法返回的余数
C = gjd_div(AI,BI,r);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
for(int i = r.size()-1;i>=0;i--){
cout<<r[i];
}
}
高精度数组的比大小
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;
}
加法
vector<int> gjd_add(vector<int> &A, vector<int> &B){
if (A.size() < B.size()) return gjd_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;
}
减法
vector<int> gjd_sub(vector<int> &A, vector<int> &B)
{
vector<int> C;
if (!cmp(A,B))
{
C = gjd_sub(B,A);
C[C.size() - 1] = -C[C.size() - 1];
return 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);
if (t < 0) t = 1;
else t = 0;
}
while (C.size() > 1 && C.back() == 0) C.pop_back();
return C;
}
乘法
vector<int> gjd_mul(vector<int> &A, vector<int> &B) {
vector<int> C(A.size() + B.size() + 7, 0); // 初始化为 0,C的size可以大一点
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;
}
while (C.size() > 1 && C.back() == 0) C.pop_back(); // 必须要去前导 0,因为最高位很可能是 0
return C;
}
除法
vector<int> gjd_div(vector<int> &A, vector<int> &B, vector<int> &r){
vector<int> C;
//判断A是否小于B 若小于则答案为0 求余数
if(!cmp(A,B)){
C.push_back(0);
r.assign(A.begin(),A.end());
return C;
}
int j = B.size();
r.assign(A.end()-j,A.end());
while(j<=A.size()){
int k=0;
while(cmp(r,B)){
vector<int> s = gjd_sub(r,B);
r.clear();
r.assign(s.begin(),s.end());
k++;
}
C.push_back(k);
if(j<A.size()) r.insert(r.begin(),A[A.size()-j-1]);
if(r.size()>1&&r.back()==0) r.pop_back();
j++;
}
reverse(C.begin(),C.end());
while(C.size()>1&&C.back()==0) C.pop_back();
return C;
}
最终代码(AC Code For 洛谷P1932 A+B A-B A*B A/B A%B Problem):
#include<bits/stdc++.h>
using namespace std;
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;
}
vector<int> gjd_add(vector<int> &A, vector<int> &B)
{
if (A.size() < B.size()) return gjd_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;
}
vector<int> gjd_mul(vector<int> &A, vector<int> &B) {
vector<int> C(A.size() + B.size() + 7, 0); // 初始化为 0,C的size可以大一点
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;
}
while (C.size() > 1 && C.back() == 0) C.pop_back(); // 必须要去前导 0,因为最高位很可能是 0
return C;
}
vector<int> gjd_sub(vector<int> &A, vector<int> &B)
{
vector<int> C;
if (!cmp(A,B))
{
C = gjd_sub(B,A);
C[C.size() - 1] = -C[C.size() - 1];
return 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);
if (t < 0) t = 1;
else t = 0;
}
while (C.size() > 1 && C.back() == 0) C.pop_back();
return C;
}
vector<int> gjd_div(vector<int> &A, vector<int> &B, vector<int> &r){
vector<int> C;
//判断A是否小于B 若小于则答案为0 求余数
if(!cmp(A,B)){
C.push_back(0);
r.assign(A.begin(),A.end());
return C;
}
int j = B.size();
r.assign(A.end()-j,A.end());
while(j<=A.size()){
int k=0;
while(cmp(r,B)){
vector<int> s = gjd_sub(r,B);
r.clear();
r.assign(s.begin(),s.end());
k++;
}
C.push_back(k);
if(j<A.size()) r.insert(r.begin(),A[A.size()-j-1]);
if(r.size()>1&&r.back()==0) r.pop_back();
j++;
}
reverse(C.begin(),C.end());
while(C.size()>1&&C.back()==0) C.pop_back();
return C;
}
int main(){
//读入
string A,B;
int b;
vector<int> AI,BI;
cin>>A>>B;
for(int i = A.length()-1;i>=0;i--){
AI.push_back(A[i]-48);
}
for(int i = B.length()-1;i>=0;i--){
BI.push_back(B[i]-48);
}
// gjd_add 加法
vector<int> C = gjd_add(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_sub 减法
C = gjd_sub(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_mul 乘法
C = gjd_mul(AI,BI);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
// gjd_div 除法
vector<int> r; // 高精度除法返回的余数
C = gjd_div(AI,BI,r);
for(int i = C.size()-1;i>=0;i--){
cout<<C[i];
}
cout<<endl;
for(int i = r.size()-1;i>=0;i--){
cout<<r[i];
}
}