大数阶乘

// 阶乘.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h" #include <stdio.h> #include <math.h> #include <stdlib.h> /* 求N!的位数公式 log10(1)+log10(2)+···+log10(n) 取整加1 */ int wei(int n) {  int i;  double sum=0;  for(i=1;i<=n;i++)   sum+=log10((double)i); /* 以万为进制，一位可以表示4个数，减少存储空间 */  return (int)((sum+1)/4+1); } /* 高精度阶乘核心 */ int _tmain(int argc,_TCHAR* argv[]) {  int i,j,n,k,m=0,jinwei,weishu=1,temp;  unsigned int *x;  scanf("%d",&n); /* 依据阶乘位数申请动态数组 */  x=(unsigned int*) malloc(wei(n)*sizeof(int));  x[0]=1;  for(i=2;i<=n;i++)  {   jinwei=0;      for(j=1;j<=weishu;j++)   {        temp=x[j-1]*i+jinwei;        if (temp>=1)    { /* 以万为进制，提高运算速度 */     x[j-1]=temp%10000;     jinwei=temp/10000;    }   }     while(jinwei)   {    weishu++;    x[weishu-1]=jinwei%10000;    jinwei/=10000;        }  } /* 先输出第一个数，防止首位出现0 */  printf("%d",x[weishu-1]);     for(j=weishu-2;j>=0;j--) printf("%04d",x[j]); /* 输出其余的数，因为万进制，需要补0 */  system("pause");  return 0;   }