算法实现求n的阶乘（防止溢出）



    

求大整数n阶乘，在找工作笔试和面试的过程中，不止一次遇到这个问题，用一个for循环迭代出的结果肯定是不行的，即直接用int，默认是32位，它能表示的最大值为2，147，483，647,但是12的阶乘为479，001，600,13的阶乘为6，227，020，800，所以当n为13的时候已经溢出了。所以当n为更大的值时，需要采用巧妙的方法来防止溢出。

• 编辑
• 删除

我们可以用数组的方式来存储每一位数据。

具体代码如下：

1. #include <iostream>    
2. #include <cstdio>    
3. #define MAX 40000    
4. using namespace std;    
5.     
6. int main(void)    
7. {    
8.     int n;    
9.     while(scanf("%d",&n)!=EOF&&n>=0)    
10.     {    
11.         int i,j;    
12.         int a[MAX];      //存数运算结果    
13.         int p,h;           //p存储当前结果的位数,h为进位    
14.         a[1]=1;    
15.         p=1;    
16.         for(i=2;i<=n;i++)   //循环与2,3,4.....n相乘    
17.         {    
18.             for(j=1,h=0;j<=p;j++)    //让a[]的每位与i相乘    
19.             {    
20.                 a[j] = a[j] * i + h ;    
21.                 h = a[j] / 10 ;    
22.                 a[j] = a[j] % 10 ;    
23.             }    
24.             while(h>0)         //如果h不为0    
25.             {    
26.                 a[j] = h % 10 ;    
27.                 h = h / 10 ;    
28.                 j++;    
29.             }    
30.             p = j - 1 ;            //将当前的位数赋给p    
31.         }    
32.         for(i=p;i>=2;i--)    
33.         {    
34.             printf("%d",a[i]);    
35.         }    
36.         printf("%d\n",a[i]);    
37.     }    
38.     return 0;    
39. }    

#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.

这样就可以计算任意数的阶乘，而不用担心溢出，自己的一点总结，积累知识，温故而知新,为自己加油！