高精度加减乘除算法C++实现

#include<bits/stdc++.h>
using namespace std;
int a[100001],b[100001],c[100001],ans[100001],la,lb,lans,i,j,k,flag;
string s1,s2;
char ch;
int main()
{
    cin>>s1>>s2;
    la=s1.size();
    lb=s2.size();
    for(i=1;i<=la;i++)
        a[i]=s1[i-1]-'0';
    for(i=1;i<=lb;i++)
        b[i]=s2[i-1]-'0';
    memset(ans,0,sizeof(ans));
    if(la>=lb)
    {
        for(i=2;i<=la+1;i++)
            ans[i]+=a[i-1];
        for(i=la+1,j=lb;i>=la-lb+2,j>=1;i--,j--)
            ans[i]+=b[j];
    }
    else
    {
        for(i=2;i<=lb+1;i++)
            ans[i]+=b[i-1];
        for(i=lb+1,j=la;i>=lb-la+2,j>=1;i--,j--)
            ans[i]+=a[j];
    }
    for(i=max(la,lb)+1;i>=2;i--)
        if(ans[i]>=10)
            ans[i]-=10,ans[i-1]++;
    cout<<"sum=";
    for(i=1;i<=max(la,lb)+1;i++)
        if(i==1&&ans[i]==0);
        else
            cout<<ans[i];
    cout<<endl;
    memset(ans,0,sizeof(ans));
    cout<<"difference=";
    if(la<lb)
    {
        for(i=1;i<=lb;i++)
            ans[i]=b[i];
        for(i=lb,j=la;i>=lb-la+1,j>=1;i--,j--)
            ans[i]-=a[j];
        for(i=lb;i>=1;i--)
            if(ans[i]<0)
                ans[i]+=10,ans[i-1]--;
        cout<<'-';
    }
    else
        if(la==lb)
        {
            for(i=1;i<=la;i++)
                if(a[i]>b[i])
                    i=la;
                else
                    if(a[i]<b[i])
                        break;
            if(i>la)
            {
                for(i=1;i<=la;i++)
                    ans[i]=a[i];
                for(i=1;i<=la;i++)
                    ans[i]-=b[i];
            }
            else
            {
                for(i=1;i<=la;i++)
                    ans[i]=b[i];
                for(i=1;i<=la;i++)
                    ans[i]-=a[i];
                cout<<'-';
            }
            for(i=la;i>=1;i--)
                if(ans[i]<0)
                    ans[i]+=10,ans[i-1]--;
        }
        else
        {
            for(i=1;i<=la;i++)
                ans[i]=a[i];
            for(i=la,j=lb;i>=la-lb+1,j>=1;i--,j--)
                ans[i]-=b[j];
            for(i=la;i>=1;i--)
                if(ans[i]<0)
                    ans[i]+=10,ans[i-1]--;
        }
    for(i=1;i<=max(la,lb);i++)
        if(ans[i]||flag)
            cout<<ans[i],flag=1;
    cout<<endl;
    cout<<"product=";
    flag=0;
    memset(ans,0,sizeof(ans));
    for(i=la;i>=1;i--)
    {
        memset(c,0,sizeof(c));
        for(j=1;j<=lb;j++)
            c[j]=a[i]*b[j];
        for(j=i+lb+1,k=lb;j>=i+1,k>=1;j--,k--)
            ans[j]+=c[k];
    }
    for(i=la+lb+1;i>=1;i--)
        if(ans[i]>=10)
            ans[i-1]+=(ans[i]/10)*10,ans[i]=ans[i]%10;
    for(i=1;i<=la+lb+1;i++)
        if(ans[i]||flag)
            cout<<ans[i],flag=1;
    cout<<endl;
    memset(ans,0,sizeof(ans));
    cout<<"quotient=";
    if(la<lb)
    {
        cout<<0;
        return 0;
    }
    if(la==lb)
        for(i=1;i<=la;i++)
            if(a[i]<b[i])
            {
                cout<<0;
                return 0;
            }
            else
                if(a[i]>b[i])
                    break;
    lans=la-lb+1;
    for(i=1;i<=lans;i++)
    {
        while(1)
        {
            for(j=i,k=1;j<=lb+i-1,k<=lb;j++,k++)
                if(a[j]>b[k])
                {
                    j=lb+1;
                    break;
                }
                else
                    if(a[j]<b[k])
                        break;
            if(j>lb)
            {
                for(j=lb+i-1,k=lb;j>=i,k>=1;j--,k--)
                    if(a[j]-b[k]>=0)
                        a[j]-=b[k];
                    else
                        a[j]+=10,a[j-1]--,a[j]-=b[k];
                ans[i]++;
            }
            else
            {
                a[i+1]+=a[i]*10;
                a[i]=0;
                break;
            }
        }
    }
    for(i=1;i<=lans;i++)
        if(i==1&&ans[i]==0);
        else
            cout<<ans[i];
    cout<<endl<<"计算结束！";
    return 0;
}