斯特林公式可以求出n!的近似值,但是如果需要求精确值的话,就要采取另外的办法了。‘
当n<=1000的时候,可以直接模拟求阶乘,用一个数组保存阶乘的每一位,大概3000的数组就行。程序如下:
#include<stdio.h>
#include<string.h>
#define N 3000
int f[N];//保存阶乘的位数
int main()
{
int n,i,j,c;
while(~scanf("%d",&n))
{
memset(f,0,sizeof(f));
f[0]=1; //f[0] 必须初始为1
for(i=2;i<=n;i++)
{
c=0; //进位
for(j=0;j<N;j++)
{
f[j]=f[j]*i+c;
c=f[j]/10;
f[j]=f[j]%10;
}
}
for(i=N-1;i>=0;i--) if(f[i]) break; //忽略前导0
for(j=i;j>=0;j--) printf("%d",f[j]);
printf("\n");
}
return 0;
}
但是当n太大的时候,比如10000,上述算法会超时,所以 还可以采用一种万进位的方法,就是数组中的每一个数都保存4位数,超过四位数才进位,形式跟上面类似,但是减少了循环次数,节省时间。
#include<stdio.h>
#include<string.h>
#define N 10000
int f[N]; //每一个数保存阶乘中的四位数
int main()
{
int n,i,j,p,c; //p表示已经用了数组的几个数 c是进位
while(~scanf("%d",&n))
{
memset(f,0,sizeof(f));
f[0]=1;
p=1;
for(i=2;i<=n;i++)
{
c=0;
for(j=0;j<p;j++)
{
f[j]=f[j]*i+c;
c=f[j]/N;
f[j]=f[j]%N;
}
if(c)
{
p++;
f[p-1]+=c;
}
}
printf("%d",f[p-1]);
for(i=p-2;i>=0;i--)
printf("%04d",f[i]); //不足4位前面补0
printf("\n");
}
return 0;
}