大数阶乘算法
“10000的阶乘怎么算?”“10000”这个数字太大了,无论用什么数据类型保存结果都会溢出。使用数组来模拟数字,这样无论结果数字有多大,只要数组的长度够长就能表示出来。
首先定义一个足够长的数组。
拿10000的阶乘为例,最后的结果长度是35660位,所以我们定义一个40000个成员的数组就可以了。
int result[40000];
其核心思想就是把计算结果每一位上的数字保存到一个数组成员中,例如:
把124保存至数组中,保存结果应该是
result[0] 4
result[1] 2
result[2] 1
这样肯定是没有问题的,一个int型数据存放一个小于10的数是绝对不会溢出。
把整个数组看成一个数字,这个数字和一个数相乘的时候,需要每一位都和这个乘数进行相乘运算还需要把前一为的进位加上。运算方法和小学数学是一样的,乘积的个位是当前位上应该表示的数字,10位以上的需要进位。
int 结果 = result[x] * 乘数 + 进位;
每一位的计算结果有了,把这个结果的个位数拿出来放到这个数组元素上:
result[x] = 结果%10;
接下来的工作就是计算出进位:
进位 = 结果 / 10;
这样一位一位的把整个数组计算一遍,最后可能还有进位,用同样的方法,把进位的数值拆成单个数字,放到相应的数组元素中。
最后输出一下结果,从最高位吧数字打印一遍就OK了。
代码:
#include <stdio.h>
int main()
{
int carry,n;
int a[40001];
int digit;// 记录有多少位了;
int temp; // 记录每次乘积;
while(scanf("%d",&n)!=EOF)
{
a[0]=1;
digit=1;
for(int i=2; i<=n; i++)
{
carry = 0;
for(int j = 0; j < digit; ++j)
{
temp=a[j]*i+carry; //每一位都乘以 i,加上上一位多出来的数;
a[j]=temp%10;//每一位存小于10的数;
carry=temp/10; // 还剩下多少;
}
while(carry) //如果没有存完,剩下的继续往下存;
{
a[digit++]=carry%10;
carry/=10;
}
}
for(int k = digit - 1; k >= 0; --k)
printf("%d",a[k]);
printf("\n");
printf("length=%d\n",digit);
}
return 0;
}