【算法模板】高精度问题

说在前面，高精度问题仅针对c；javap与python自带这个函数。

将大整数视作一个数组，先存储低位数。

一、大整数加法

#include <iostream>
#include <vector>

using namespace std;
const int N = 1e6+10;

//C=A+B
vector<int>add(vector<int> &A,vector<int> &B)            //这里的&是为了提高速度，减少一次数组拷贝
{
    vector<int> C;
    
    int t=0;//t是进位
    for (int i = 0; i < A.size()||i < B.size(); i ++ )   //  
    {
        if(i<A.size())      t+=A[i];                      //t=A[i]+B[i]+t
        if(i<B.size())      t+=B[i];
        C.push_back(t%10);                                //t取余是C[i]
        t/=10;                                            //看有没有进位
        
    }
    if(t)   C.push_back(1);                               //为了c的最高位，若循环结束t不等于0，还要进一位
    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);                                  //auto是自动判断，就是vector<int> C
    
    for(int i =C.size()-1;i>=0;i--)         cout<<C[i];
    
    
    return 0;
}

问题一：为什么a[i]-'0' 点击这里

问题二: vector 点击这里

二、大整数减法

#include <iostream>
#include <vector>


using namespace std;

bool cmp(vector<int>& A, vector<int>& B)                        //判断是否A >= 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 1;
}

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;                    //t=A[i]-B[i]-t,因为B位数少，所以判断B[i]是否还有位数再减
        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) //去掉前导0，如123-120=003；C.size()>1为123-123=000
        C.pop_back();
        
    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');

    if (cmp(A, B))                                     //判断是否A>B,若不是则写个负号再算B-A
    {
        auto C = sub(A, B);
        
        for (int i = C.size() - 1; i >= 0; i--)          printf("%d", C[i]);

    }
    else
    {
        auto C = sub(B,A);
        
        printf("-");
        
        for (int i = C.size() - 1; i >= 0; i--)          printf("%d", C[i]);
    }

    return 0;

}

C.push_back((t + 10) % 10)是由于 t>=0 t=t

t<0 t= t+10

21-15；个位1-5=-4；t=-4那么t<0就要加10，得到6，21-15的个位才是6。从而总结为这个公式 (t + 10) % 10

三、大整数乘法

#include<iostream>
#include<vector>
const int N = 1e6+10;
using namespace std;

//C=A*b
vector<int> mul(vector <int> &A,int b)
{
    vector<int> C;
    int  t=0;  //进位
    for(int i = 0; i < A.size() || t ;i++)    //当i没有循环完或进位没有进完，就一直做
    {
        t = A[i] * b +t ;                //实现乘法
        C.push_back(t%10);              //取每一位
        t /= 10;                        //得到进位
    }
    
    return C;
}
int main()
{
    string a;
    int b;
    vector<int> A;
    cin>>a>>b;
    for(int i = a.size() - 1 ; i >= 0; i --)            A.push_back(a[i] -'0');
    
    if(b==0)    printf("0");
    else    
    {
        auto C = mul(A,b);
    
        for(int i = C.size() - 1; i >= 0; i -- )            printf("%d",C[i]);
    }
            return 0;
}

四、大整数除法

#include<iostream>
#include<vector>
#include<algorithm>

using namespace std;

//C=A/b   余数是r
vector<int> divi (vector<int> &A,int b,int &r)//r是引用
{
    vector<int> C;
    
    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());            //除法是从高位开始push，所以反过来
    while(C.size() > 1 && C.back()==0)      C.pop_back(); //去掉前导0
 
    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  = divi(A,b,r);
    
    for(int i  = C.size()-1;i >= 0;i--)         printf("%d",C[i]);
    cout<<endl<<r<<endl;
    
    return 0;
}