高精度N的阶乘-N!

题目:输入一个正整数N(0<N<=30),求N!
代码:

#include<stdlib.h>
#include<cstdio>
#include<string>
#include<iostream>
using namespace std;
int a[2005]={0};
int main()
{
    int n,i,j,t=0,sum=0,b=0;//t represent carry bit
    cin>>n;
    a[1]=1;
    
    for(i=2;i<=n;i++)
    {
        for(j=1;j<=2000;j++)//enough bit to store the number
        {
            sum=a[j]*i+t;//each bit x i 
            a[j]=sum%10;//finished bit number 
            t=sum/10;//turn to the next bit
        }
    }
    
    if(n==1)cout<<"1";
    else
    {   for(j=2000;j>=1;j--)
        {
            if(a[j]!=0)b=1;
        
            if(b==1)cout<<a[j];
        }
    }

    return 0;
}
/*
   24
 x  5
 ----
  120
  abc
*/

说句实在的,这题我解决了很久,思路并没有像我想的那样快速而又灵活。
定义一个大数组a存储位数,这题很明显如果使用普通的递归函数调用或者用一个for循环暴力解决会使最终的数据超过int的32位整数范围。
那么思路这时就比较清晰:用数组存储每一位的位数,最后输出这个数组。
使用了两个for循环来模拟N!

第一个for循环:for(i=2;i<=n;i++)模拟与数i相乘。
第二个for循环: for(j=1;j<=2000;j++) 遍历数组中存储的每一位,在我第一次想到他的时候,问题在于取定j的范围,也就是最后结果的N!的位数。
这样去做太过麻烦,我做了两三次尝试无果,最后发现只需定义一个比较大足以容得下最大的N!的位数即可。

第二个for循环中的语句:

                for(j=1;j<=2000;j++)//enough bit to store the number
        {
            sum=a[j]*i+t;//each bit x i 
            a[j]=sum%10;//finished bit number 
            t=sum/10;//turn to the next bit
        }

这里来说明几个变量的意义:1.sum 是位数和数i的乘积与进位数t之和 2.a[j]存储位数 3.t代表进位数
sum=a[j]*i+t; 在阶乘的过程中,从小的位数到大的位数依次乘以此时的数i,sum存储此时第j位位数与i相乘之后再加上进位数t的和。
a[j]=sum%10;位数保留sum的个位数。
t=sum/10;进位数t=sum/10。

举个例子:24*5 先对24的个位4进行处理,4乘以5等于20,此时sum=20,留下的位数为0,进位数t=2;再对十位数2进行处理,2乘以5等于10,此时的sum=2x5+t=10+2=12,保留个位数字2为位数,进制数t=12/10=1;那么下一次执行的时候,sum=0+t=1,a[j]=1,t=0。实际上对数字位数的遍历操作已经结束了。

此时保证了在遍历数组元素(位数)的过程中,不会发生越界的现象,我在刚刚开始的时候对所有数组的元素赋值0,保证不会出现奇怪的错误和多位的情况。
最后输出的时候,用一个状态变量b,判断是否遇到了我们存储的“最大的”位数,是的话开始输出。

在C语言中计算高精度阶乘(n!)是一项挑战,因为阶乘的结果随着n的增加会非常快地变得非常大,可能会超出int或long long等标准数据类型的范围。为此,我们通常需要借助数组或者字符串来存储中间结果,并使用迭代方法来计算。 以下是一个简单的高精度阶乘算法示例,使用了循环和动态内存管理: ```c #include <stdio.h> #include <stdlib.h> // 定义结构体存放整数的数组,作为高精度数值 typedef struct { int digits[100]; // 可以根据需要调整数组大小 int size; // 数组实际长度减一 } BigInt; // 初始化一个BigInt为零 BigInt zero = { {0}, 0 }; // 阶乘函数 BigInt factorial(int n) { BigInt result = zero; for (int i = 1; i <= n; ++i) { result.digits[result.size++] = i; // 将每个数字加到结果中 if (result.size > 1 && result.digits[result.size - 1] == 0) { // 压缩结果 --result.size; } carry(result); // 进位操作 } return result; } // 进位操作,将当前最大位上的值累加到下一位 void carry(BigInt &num) { int carry = num.digits[num.size - 1]; num.digits[num.size - 1] = 0; while (carry > 0 && num.size > 0) { num.digits[num.size - 1]++; carry -= num.digits[num.size - 1]; ++num.size; } if (carry > 0) { num.digits = realloc(num.digits, sizeof(int) * (num.size + 1)); num.digits[num.size] = carry; ++num.size; } } // 输出高精度阶乘结果 void print(BigInt num) { for (int i = num.size - 1; i >= 0; --i) { printf("%d", num.digits[i]); } printf("\n"); } int main() { int n; printf("请输入一个整数: "); scanf("%d", &n); print(factorial(n)); return 0; }
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值