mr 的大數階乘

1. #include <stdio.h>
2. #define N 100 //n的值
3. #define BITS 200 //不小於結果的位元
4. int main(void)
5. {
6.     int s[(BITS+11)/6]={1,1},k=2,i=1,a=0;
7.     
9. for(;i<=*s||(k++<N&&(a=0,i=1));a+=s[i]*k,s[i]=a%1000000,a/=1000000,*s==i++&&a&&(*s)++);
10.     for(printf("%d",s[*s]);--*s;printf("%06d",s[*s]));
11.     return 0;
12. }
16. char *mul(char *s, char *a, char *b)
17. {
18.     size_t alen=strlen(a), blen=strlen(b), slen=alen+blen, i=slen, j, k;
19.     if ((a[0]=='0' && a[1]=='/0') || (b[0]=='0' && b[1]=='/0'))
20.     {
21.         s[0]='0';
22.         s[1]='/0';
23.         return s;
24.     }
25.     while (i) s[--i]=0;
26.     for (; i<alen; i++)
27.     {
28.         for (j=0; j<blen; j++)
29.         {
30.             k=-1;
31.             s[i+j+1] += (a[i]&15)*(b[j]&15);
32.             while (s[i+j-k]>9)
33.             {
34.                 s[i+j-k-1] += s[i+j-k]/10;
35.                 s[i+j-k] %= 10;
36.                 k++;
37.             }
38.         }
39.     }
40.     for (i=0; i<slen; i++) s[i]|='0';
41.     s[i]='/0';
42.     if (s[0]=='0')
43.         for (i=0; s[i]=s[i+1]; i++);
44.     return s;
45. }