c++实现高精度

我们知道c++不同于python,python对于数值的运算是可以直接使用高精度的(因为python的代码不需要交代变量的类型),而c++对于每个变量的定义都需要限定类型,c++中存储空间最大的整数类型为long long int型,能表示最大的数字为2^{63}-1,但是对于某些题目来说,要表示的数字远远大于long long int能够表示的数字,所以这就需要实现高精度。

思想:小学的时候学习了采用列竖式的方法计算加法和乘法,所以我们利用c++来模拟这个过程,我们用一个较大的数组记录结果的每一位,数字的低位位于数组的低位,在计算的过程中,记录下乘数和被乘数的每一位相乘或相加后得到的结果,并将结果分为进位和当前位,将进位和当前位都加入之前的运算结果中,并用一个变量记录数字的当前位数。

以洛谷中阶乘之和为例子,在该例子中,既要实现高精度的加法又要实现高精度的加法,代码如下:

#include<bits/stdc++.h>
using namespace std;

//result数组用来存储最终的结果,factorial用来存储阶乘的结果,tmp存储中间结果,下标零的位置不记录
int result[1005],factorial[1005],tmp[1005];

int main(void){
    //carry用来记录进位数字相加或相乘中的进位,length记录阶乘计算结果的位数,ansLength记录最终结果的位数
    //x记录当前要乘的数,mid记录中间结果,t记录从哪个位置开始相加
    int n,carry=0,ansLength=1,length=1,tmplength=1,x,mid,t;
    scanf("%d",&n);
    result[ansLength]=1;//初始化为存有1!
    factorial[length]=1;
    for(int i=2;i<=n;i++){
        //在(i-1)!基础上计算i!
        memset(tmp,0,sizeof(tmp));
        x=i;
        t=0;
        tmplength=length;
        while (x)
        {
            carry=0;
            mid = x%10;
            x /=10;
            //将mid和阶乘的各位相乘
            for(int j=1;j<=length;j++){
                int k = mid * factorial[j];
                tmp[j+t] +=(k+carry)%10;//更新结果
                carry=(k+carry)/10;//更新进位
            }
            tmplength=length+t;
            if(carry){
                ++tmplength;
                tmp[tmplength] +=carry;
                carry=0;
            }
            t++;
                
        }

        length=tmplength;
        for(int j=1;j<=length;j++){
            factorial[j]=tmp[j];
        }
        // for(int j=length;j>0;j--)
        //     printf("%d",factorial[j]);
        // printf("\n");
        int addtmp,addcarry=0;
        //将阶乘结果相加
        for(int j=1;j<=length;j++){
            addtmp=result[j]+factorial[j];
            result[j]=(addtmp+addcarry)%10;
            addcarry=(addtmp+addcarry)/10;
        }
        ansLength=length;
        if(addcarry)
            result[++ansLength]=addcarry;
        
    }
    for(int i=ansLength;i>0;i--){
        printf("%d",result[i]);
    }
    printf("\n");
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值