高精度算阶乘

1. /*
2. 习题 92：高精度算阶乘★★
3. 题目描述：
4. 算出n!的完整结果，其中n!=1*2*3*...*n
6. 输入：
7. 多组测试数据，一行一组，每行仅一个数n
8. 其中0<=n<=10000
10. 输出：
11. 输出相应的n!的结果，必须精确到个位
13. 样例输入：
14. 10
15. 20
16. 100
18. 样例输出：
19. 3628800
20. 2432902008176640000
21. 933262154439441526816992388562667004907159682643816214685929
22. 638952175999932299156089414639761565182862536979208272237582
23. 51185210916864000000000000000000000000
25. 其它信息：
26. 最后一个100!的结果由于过长，故拆分成三行，每行60字符，请见谅
28. 难度：Easy
29. */
30. #define MAX 10000
31. #include<stdio.h>
32. #include<string.h>
33. int main()
34. {
35.     int n;
36.     int i,j,k;
37.     char result[MAX];
38.     char tmp[MAX];
39.     while(scanf("%d",&n)!=EOF){
40.         if(n==0){
41.             printf("1/n");
42.             continue;
43.         }
44.         memset(result,'0',MAX);
45.         memset(tmp,'0',MAX);
46.         tmp[0]='1';
47.         for(i=1;i<=n;i++){
48.             memset(result,'0',MAX);
49.             for(j=0;j<i;j++){
50.                 for(k=0;k<MAX;k++){
51.                     result[k]=result[k]+(tmp[k]-48);
52.                     if(result[k]>57){
53.                         result[k+1]=result[k+1]+(result[k]-48)/10;
54.                         result[k]=result[k]-(result[k]-48)/10*10;
55.                     }
56.                 }
57.             }
58.             memcpy(tmp,result,MAX);
59.         }
60.         for(i=MAX-1;i>=0;i--)
61.             if(result[i]!='0')break;
62.         for(;i>=0;i--)printf("%c",result[i]);
63.         printf("/n");
64.     }
65. }

<script src="http://pagead2.googlesyndication.com/pagead/show_ads.js" type=text/javascript> </script>