大数阶乘

// 阶乘.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; ? }