C++ 高精度算法

目录

1.A+B

2.A-B

3.A*b

4.A/b

---

---

1.A+B

（1）若采用数组存储数，下标从起始到结束依次存放数的低位到高位，方便在数组最后的位置进位，反之，则需要从起始位置插入进位，相对麻烦。

（2）预设数组内存时，假设数据在0~10000，一般扩容为1e4+C，C为常数，目的是防止边界问题。

 此处采用vercor容器，相对数组能更好的操作数的相加运算和进位，模板如下：

//传入引用：在函数调用时在内存中不会生成副本，速度更快
vector<int>add(vector<int>&A,vector<int>&B)
{
    vector<int>C;
    int t = 0;
    for(int i=0;i<A.size()||i<B.size();i++)
    {
       if(i<A.size())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;
}

完整代码（输入A，B，输出A+B）：

#include`<iostream>`
#include`<vector>`

>using namespace std;

>const int M =1e6+10;

>vector<int>add(vector<int>&A,vector<int>&B)
{
   vector<int>C;
   int t = 0;
    for(int i=0;i<A.size()||i<B.size();i++)
    {
       if(i<A.size())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];
   }
   return 0;
}

---

2.A-B

和加法类似，但是要判断结果是否为负数，即A是否大于b，若小于输出要加上‘-’号，在最后的结果中，由于存在借位，最高位可能为零，需要去除前导0的情况，例如我们需要输出12345，而不是012345。

 模板如下：

//判断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;
}

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);
        if(t<0)t=1;
       else t  = 0;
    }
    
   while(C.size()>1&&C.back()==0)C.pop_back();//去除前导0的情况
   return C;//返回计算结果，即C容器
}

完整代码（输入A，B，输出A-B）如下：

#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()||t;i++)
   {
      if(i<A.size())t+=A[i]*b;
       C.push_back(t%10);
       t/=10;
   }
   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);
   
   if(C[C.size()-1]==0)
       cout<<0;
   else
       for(int i=C.size()-1;i>=0;i--)cout<<C[i];
   return 0;
}

---

3.A*b

计算A*b不能单纯用人工的计算方式，而是类似一位数乘上一个小于10的整数运算，即把b当作一个整体去与A的每一位相乘，例如56*12

（1）6*12=72，72/10=7，72%10=2，进位7，留下2
（2）5*12+7=67，67/10=6，67%10=7，进位6，留下7
（3）高位0加上进位6，完成最后计算，结果为6 7 2 ，也就是672

 模板如下：

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;
   }
   return C;
}

完整代码如下（输入A，b输出A*b的结果）：

#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()||t;i++)
   {
      if(i<A.size())t+=A[i]*b;
       C.push_back(t%10);
       t/=10;
   }
   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);
   
   if(C[C.size()-1]==0)
       cout<<0;
   else
       for(int i=C.size()-1;i>=0;i--)cout<<C[i];
   return 0;
}

---

4.A/b

计算A/b和A*b类似，但是从最高位开始，由一位开始除运算，留下的余数乘上10在加上下一位开始除运算，进位则是余数对10取模，也会存在前导0情况，而且对于计算的结果是从高位到低位尾插到vector中的，需要逆序，才是最后的结果。

模板如下：

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];
      C.push_back(r/b);
      r%=b;
  }
  reverse(C.begin(),C.end());
  
  while(C.size()>1&&C.back()==0)C.pop_back();
  
  return C;
}

完整代码（输入A，b，输出A/b和余数）：

#include<vector>

#include<algorithm>

using namespace std;

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];
      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=0;
   auto C = div(A,b,r);
  
  for(int i=C.size()-1;i>=0;i--)cout<<C[i];
  cout<<endl<<r<<endl;
  return 0;
}