大数运算：计算20000的阶乘

fucks

1. #include <stdio.h>
   
2. #include <stdlib.h>
   
3. 
   
4. char* real_mult(const char* a, const char* b) {
   
5.         int i, j, k, m, n, len_a, len_b;
   
6.         len_a = strlen(a);
   
7.         len_b = strlen(b);
   
8.         int* res_p;
   
9.         char* c;
   
10.         res_p = (int*)malloc((len_a + len_b) * sizeof(int));
    
11.         if (!res_p) {
    
12.         printf((const char*)"malloc int[%d] fail!\n", len_a + len_b);
    
13.         for ( ; ; );
    
14.         }
    
15.         for (i = len_a + len_b - 1; i >= 0; i--) {
    
16.                 res_p[i] = 0;
    
17.         }
    
18.         for (i = 0; i < len_a; i++) {
    
19.                 for (j = 0; j < len_b; j++) {
    
20.                         k = i + j + 1;
    
21.                         res_p[k] += (a[i] - '0') * (b[j] - '0');
    
22.                 }
    
23.         }
    
24.         res_p[0] = 0;
    
25.         for ( ; k > 0; k--) {
    
26.                 res_p[k - 1] += res_p[k] / 10;
    
27.                 res_p[k] %= 10;
    
28.         }
    
29.         for (j = 0; j < len_a + len_b; j++) {
    
30.         if (res_p[j] != 0) {
    
31.             break;
    
32.         }
    
33.         }
    
34.         c = (char*)malloc((len_a + len_b + 1 - j) * sizeof(char));
    
35.         if (!c) {
    
36.         printf((const char*)"malloc char[%d] fail!\n", len_a + len_b + 1 - j);
    
37.         for ( ; ; );
    
38.         }
    
39.         for (i = 0; i < len_a + len_b - j; i++) {
    
40.                 c[i] = (char)(res_p[i + j] + '0');
    
41.         }
    
42.         free(res_p);
    
43.         c[i] = '\0';
    
44.         return c;
    
45. }
    
46. 
    
47. char* mult1(char* a, const char* b) {
    
48.     char* c = real_mult((const char*)a, b);
    
49.     free(a);
    
50.     return c;
    
51. }
    
52. 
    
53. char* order(int n) {
    
54.     char temp[16];
    
55.     char* res = (char*)malloc(32);
    
56.     for ( ; n > 1; n--) {
    
57.         sprintf((char*)temp, "%d", n);
    
58.         res = mult1(res, (const char*)temp);
    
59.     }
    
60.     return res;
    
61. }
    
62. 
    
63. int main(int argc, const char* argv[]) {
    
64.     char* c;
    
65.     c = order(20000);
    
66.     printf("%s\n", c);
    
67.     free(c);
    
68.     return 0;
    
69. }

复制代码