求大整数n阶乘,在找工作笔试和面试的过程中,不止一次遇到这个问题,用一个for循环迭代出的结果肯定是不行的,即直接用int,默认是32位,它能表示的最大值为2,147,483,647,但是12的阶乘为479,001,600,13的阶乘为6,227,020,800,所以当n为13的时候已经溢出了。所以当n为更大的值时,需要采用巧妙的方法来防止溢出。
我们可以用数组的方式来存储每一位数据。
具体代码如下:
- #include <iostream>
- #include <cstdio>
- #define MAX 40000
- using namespace std;
- int main(void)
- {
- int n;
- while(scanf("%d",&n)!=EOF&&n>=0)
- {
- int i,j;
- int a[MAX]; //存数运算结果
- int p,h; //p存储当前结果的位数,h为进位
- a[1]=1;
- p=1;
- for(i=2;i<=n;i++) //循环与2,3,4.....n相乘
- {
- for(j=1,h=0;j<=p;j++) //让a[]的每位与i相乘
- {
- a[j] = a[j] * i + h ;
- h = a[j] / 10 ;
- a[j] = a[j] % 10 ;
- }
- while(h>0) //如果h不为0
- {
- a[j] = h % 10 ;
- h = h / 10 ;
- j++;
- }
- p = j - 1 ; //将当前的位数赋给p
- }
- for(i=p;i>=2;i--)
- {
- printf("%d",a[i]);
- }
- printf("%d\n",a[i]);
- }
- return 0;
- }
#include <iostream>
#include <cstdio>
#define MAX 40000
using namespace std;
int main(void)
{
int n;
while(scanf("%d",&n)!=EOF&&n>=0)
{
int i,j;
int a[MAX]; //存数运算结果
int p,h; //p存储当前结果的位数,h为进位
a[1]=1;
p=1;
for(i=2;i<=n;i++) //循环与2,3,4.....n相乘
{
for(j=1,h=0;j<=p;j++) //让a[]的每位与i相乘
{
a[j] = a[j] * i + h ;
h = a[j] / 10 ;
a[j] = a[j] % 10 ;
}
while(h>0) //如果h不为0
{
a[j] = h % 10 ;
h = h / 10 ;
j++;
}
p = j - 1 ; //将当前的位数赋给p
}
for(i=p;i>=2;i--)
{
printf("%d",a[i]);
}
printf("%d\n",a[i]);
}
return 0;
}
直接看代码比较难理解,可以举个例子,比如说我们现在得到了4!= 24,在数组中是以倒序的方式存储,即42,a[1] = 4,a[2] = 2(数组下标从1开始),现在要计算5!,内层for循环对目前的两位4和2进行处理,即每一位都乘以5,首先a[1] = a[1]*5 + 0(h)= 20,h从0变为2,a[1] = 0;然后a[2] = a[2] * 5 + 2 = 12,h从2变为1,a[2] = 2;这时内层for循环执行结束。
总结1:内层for循环的作用是更新已经存在的位数中的每一位与i的乘积的结果.
接着执行while循环,目前h = 1,a[3] = 1,h 从1变为0,然后退出while循环。处理结束就得到了5! = 120.p的位数增加到3.
总结2:while循环的作用是计算p更新之前的位数之后的进位,比如说p = 2的时候,数组中为42,然后内层for循环把42更新为02,while循环增加进位1,整个数组变为021.
这样就可以计算任意数的阶乘,而不用担心溢出,自己的一点总结,积累知识,温故而知新,为自己加油!