高精度运算

一、前言

本篇笔记以图解形式为主，望理解（其实是今天想摸一下鱼）。除法基本上只有被除数会用到高精度，除数为整数，在了解了前面几种高精度的实现方式后就很容易写出来。

二、高精度算法原理

1、存储方式

堆栈和数组都可以，实现原理大同小异，我还是习惯用堆栈写高精度加法。

2、算法原理（高精度加法）

数组简单的进位加法（如下图1）

注意：用数组存储的时候低位放前面，高位放后面。下图是堆栈实现的原理。

代码实现：

#include<iostream>
#include<vector>
using namespace std;

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;
    cin>>a>>b;
    
    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=add(A,B);
    
    for(int i=C.size()-1;i>=0;i--)cout<<C[i];
    return 0;
}

3、算法原理（高精度减法）

代码实现

#include<iostream>
#include<vector>
#include<algorithm>
using namespace  std;

vector<int>mul(vector<int>A,vector<int>B)
{
    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)%10);
        if(t<0)t=-1;
        else t=0;
    }
    while(!C.back()&&C.size()>1)C.pop_back();
    
    return C;
}

bool cmp(vector<int >A,vector<int>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;
}

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=mul(A,B);
    }
    else 
    {
        cout<<"-";
        C=mul(B,A);
    }
    
    for(int i=C.size()-1;i>=0;i--)cout<<C[i];
    
    return 0;
}

3、算法原理（高精度乘法）

注意数组存储的下标从1开始。

代码实现

/*数组实现高精度乘法*/
#include<iostream>
#include<vector>
using namespace std;
const int N=100010;

int a1[N],b1[N];
int c[N];

int main()
{
    string a,b;
    cin>>a>>b;
    
    int la=a.size(),lb=b.size();
    for(int i=la-1;i>=0;i--)a1[la-i]=a[i]-'0';
    for(int i=lb-1;i>=0;i--)b1[lb-i]=b[i]-'0';
    
    for(int i=1;i<=la;i++)
    {
        for(int j=1;j<=lb;j++)
        {
           c[i+j-1]+=a1[i]*b1[j];
         }
    }
    
    int lc=la+lb;
    for(int i=1;i<=lc;i++)
    {
        if(c[i]>9)
        {
            c[i+1]+=c[i]/10;
            c[i]%=10;
        }
    }
    
    int len=lc;
    while(!c[len]&&len>1)len--;
    
    for(int i=len;i>=1;i--)cout<<c[i];
    return 0;
}

三、小结

好感慨啊，两个半月前还只能看着别的dailao写的高精度题解现在也能发表一些本蒟蒻对于高精度的见解了(*^▽^*)，还是有点开心的。希望自己在将来也能能够取得更大的进步，加油。