算法:计算n!

来源:互联网 发布:晨曦软件视频教程 编辑:程序博客网 时间:2024/06/06 19:13

基本思想:设n!是m为正整数,则10^(m-1)<=n!<=10^m-1<10^m,则m-1<=lg(n!)<m,[lgn!](表示取整)=m-1,所以

[lg1+lg2+lg3+lg4+......lgn]=m-1,m=[lg1+lg2+lg3+lg4+......lgn]+1。即确定了n!的位数。


#include<stdio.h>

#include<math.h>
void main()
{
int j,k,m,n,a[40000];
long g,t;
double s;
printf("请输入正整数n(n<10000):"); 
scanf("%d",&n);  //输入n
s=0;
for(k=2;k<=n;k++)
s+=log10(k);    //对数累加,确定n!的位数
m=(int)s+1;
for(k=1;k<=m;k++)
a[k]=0;   //数组清零
a[1]=1;g=0;
for(k=2;k<=n;k++)
for(j=1;j<=m;j++)
{
t=a[j]*k+g;   //数组累乘并进位
a[j]=t%10;
g=t/10;
}
printf(" %d!=",n);
for(j=m;j>=1;j--)
printf("%d",a[j]);
printf("\n  共%d位.\n",m);
}