高精度乘法

标签：

（转）

1. /*
3. 高精度乘法
4. 输入：两行，每行表示一个非负整数（不超过10000位）
5. 输出：两数的乘积。
6. */
8. #include <stdio.h>
9. #include <string.h>
10. #include <stdlib.h>
11. #include <memory.h>
13. #define MAX 10001
15. int high_precision(int *sum,int *a,int *b,int a_len,int b_len)
16. {
17.     int i,j;
18.     memset(sum,0,sizeof(sum));
19.     int sum_len = 0 ;
20.     for(i=1;i<=a_len;i++) /*用数组模拟运算*/
21.         for(j=1,sum_len=i-1;j<=b_len;j++)
22.             sum[++sum_len] += b[j]*a[i];
23.     for(i=1;i<=sum_len;i++)/*进位处理*/
24.         if(sum[i] >= 10)
25.         {
26.             if(sum[sum_len] >= 10)
27.                 sum_len++;
28.             sum[i+1] += sum[i]/10;
29.             sum[i] %= 10;
30.         }
32.     return sum_len;
33. }
35. void multiply(char* a, char* b, char* c)
36. {
37.     int i, j, ca, cb, * s;
38.     ca = strlen(a);
39.     cb = strlen(b);
40.     s = (int*)malloc(sizeof(int) * (ca + cb));
41.     for (i = 0; i < ca + cb; i++)
42.         s[i] = 0;
43.     for (i = 0; i < ca; i++)
44.         for (j = 0; j < cb; j++)
45.             s[i+j+1] += (a[i] - ‘0‘) * (b[j] - ‘0‘);
46.     for (i = ca + cb - 1; i >= 0; i--)
47.         if (s[i] >= 10)
48.         {
49.             s[i-1] += s[i] / 10;
50.             s[i] %= 10;
51.         }
52.     i = 0;
53.     while (s[i] == 0)
54.         i++;
55.     for (j = 0; i < ca + cb; i++, j++)
56.         c[j] = s[i] + ‘0‘;
57.     c[j] = ‘\0‘;
58.     free(s);
59. }
61. int main(int argc,char *argv[])
62. {
63.     int a[MAX]= {0},b[MAX]= {0},sum[MAX*2]= {0};
64.     int a_len=0,b_len=0,sum_len=0;
65.     int i,j;
66.     char c_a[MAX],c_b[MAX],c_sum[MAX*2];
68.     scanf("%s%s",c_a,c_b);
70.     a_len = strlen(c_a);
71.     b_len = strlen(c_b);
73.     for(i=1,j=a_len-1;i<= a_len;i++,j--)
74.         a[i] = c_a[j] - ‘0‘;
75.     for(i=1,j=b_len-1;i<= b_len;i++,j--)
76.         b[i] = c_b[j] - ‘0‘;
78.     sum_len = high_precision(sum,a,b,a_len,b_len) ;
80.     for(i=sum_len;i>=1;i--)
81.         printf("%d",sum[i]);
82.     putchar(‘\n‘);
83.     multiply(c_a,c_b,c_sum);
84.     printf("%s\n",c_sum);
86.     return 0 ;
87. }

高精度乘法

标签：

原文地址：http://www.cnblogs.com/yskyskyer123/p/4530110.html