高精度减法的问题！

xiaozhu2007

下面是一个实现高精度减法的程序，a[]，b[]存放两个操作数，r[]存放结果。
算法是用两个指针p,q操作数组的，其中p始终指向大的操作数组，q指向小的操作数组。
程序运行的结果有问题，请高手帮忙看看，程序哪些地方有问题？谢谢

1. 
   
2. #include <stdio.h>
   
3. #include <stdlib.h>
   
4. 
   
5. int main(int argc, char** argv)
   
6. {
   
7.         int                     a[100], b[100], r[100];//数组p保存结果
   
8.         int                     i = 0, j = 0;
   
9.         int                     len1, len2, val_p, val_q, k, c, m;
   
10.         int                     *p = NULL, *q = NULL;//p指向大数数组，q指向小数数组
    
11. 
    
12. 
    
13.         for(k = 0; k < 100; k++){
    
14.                 a[k] = 0;
    
15.                 b[k] = 0;
    
16.                 r[k] = 0;
    
17.         }
    
18. 
    
19.         while((c = getchar()) != '\n'){
    
20.                 a[i] = c - 48;
    
21.                 i++;
    
22.         }//输入一个减数,i是a的长度
    
23. 
    
24.         while((c = getchar()) != '\n'){
    
25.                 b[j] = c - 48;//字符转化为整数
    
26.                 j++;
    
27.         }//输入被减数,j是b的长度
    
28. 
    
29.         i--;
    
30.         j--;
    
31.         m = (i > j) ? i : j;//m是i,j中的最大值
    
32.         if(i > j){
    
33.                 p = a;
    
34.                 q = b;
    
35.                 len1 = i;
    
36.                 len2 = j;
    
37.         }
    
38.         else{
    
39.                 p = b;
    
40.                 q = a;
    
41.                 len1 = j;
    
42.                 len2 = i;
    
43.         }//p始终指向长数组,len1始终表示P所表示数的长度
    
44. 
    
45.         for(k = m; k > 0; k--){
    
46.                 if((len1 >= 0) && (len2 >= 0)){
    
47.                         val_p = *(p+len1);
    
48.                         val_q = *(q+len2);
    
49.                 }
    
50.                 else if((len1 >= 0) && (len2 < 0)){
    
51.                         val_p = *(p+i);
    
52.                         val_q = 0;
    
53.                 }
    
54. 
    
55.                 if(val_p >= val_q){
    
56.                         r[k] = val_p - val_q;
    
57.                 }
    
58.                 else{//借位
    
59.                         (*(p+len1-1))--;
    
60.                         r[k] = 10 + val_p - val_q;
    
61.                 }
    
62.                 len1--;
    
63.                 len2--;
    
64.         }
    
65. 
    
66.         if((i <= j) && (a[0] < b[0])){//此条件满足，表示减数小于被减数，结果是负数
    
67.                 printf("-");
    
68.         }
    
69.         if(r[0] != 0){
    
70.                 printf("%d", r[0]);
    
71.         }
    
72.         for(k = 1; k <= m; k++){
    
73.                 printf("%d", r[k]);
    
74.         }
    
75.         printf("\n");
    
76. 
    
77.         exit(0);
    
78. }
    

复制代码

win_hate

>>if((i <= j) && (a[0] < b[0])){//此条件满足，表示减数小于被减数，结果是负数

这个逻辑是错误的， a[0], b[0] 是最高位，为什么要用 && 呢？

而且当 i=j 时，a<b 并不意味着 a[0]<b[0]，比如  11<12

如果用十进制，一般的处理是这样的，a, b 都看成是 n 个十进制位，其中 n 是 a,b 中长度较大者。

用一个 c 去记录借位，0 表示没有借位， -1 表示借了位

如果最后 c 是 0，则得到结果，如果 c 是负数，则表示 a<b, 计算的实际上是 10^n+a-b，所以再用 10^n 减去当前结果（10^n+a-b）得到 b-a, 结果为负的  b-a。

如果发现 c 最终小于 0 后，不想用 10^n 去减中间结果，则直接用 b 去减 a 也可以。

agaric

<c语言接口与实现>  里面有介绍任意精度算法 ，你可以看看

win_hate

事实上在长度相等时，你可以逐位比较，弄清楚谁大再减，这样比较简洁。

xiaozhu2007

1. 
   
2. #include <stdio.h>
   
3. #include <stdlib.h>
   
4. 
   
5. int main(int argc, char** argv)
   
6. {
   
7.         int                     a[100], b[100], r[100];//数组p保存结果
   
8.         int                     i = 0, j = 0;
   
9.         int                     len1, len2, val_p, val_q, k, c, ch, m;
   
10.         int                     *p = a, *q = b;//p指向大数数组，q指向小数数组
    
11. 
    
12.         for(k = 0; k < 100; k++){
    
13.                 a[k] = 0;
    
14.                 b[k] = 0;
    
15.                 r[k] = 0;
    
16.         }
    
17. 
    
18.         while((ch = getchar()) != '\n'){
    
19.                 a[i] = ch - 48;
    
20.                 i++;
    
21.         }//输入一个减数,i是a的长度
    
22. 
    
23.         while((ch = getchar()) != '\n'){
    
24.                 b[j] = ch - 48;//字符转化为整数
    
25.                 j++;
    
26.         }//输入被减数,j是b的长度
    
27. 
    
28.         i--;
    
29.         j--;
    
30.         m = (i > j) ? i : j;//m是i,j中的最大值
    
31. 
    
32. //判断谁是大的操作数，并用p指向它，len1表示它的长度
    
33.         if(i > j){
    
34.                 p = a;
    
35.                 q = b;
    
36.                 len1 = i;
    
37.                 len2 = j;
    
38.         }
    
39.         else if(i < j){
    
40.                 p = b;
    
41.                 q = a;
    
42.                 len1 = j;
    
43.                 len2 = i;
    
44.         }
    
45.         else if(i == j){
    
46.                 for(k = 0; k < i; k++){
    
47.                         if(a[k] > b[k]){
    
48.                                 p = a;
    
49.                                 q = b;
    
50.                                 len1 = i;
    
51.                                 len2 = j;
    
52.                                 break;
    
53.                         }
    
54.                         if(a[k] < b[k]){
    
55.                                 p = b;
    
56.                                 q = a;
    
57.                                 len1 = j;
    
58.                                 len2 = i;
    
59.                                 break;
    
60.                         }
    
61.                 }//for
    
62.         }
    
63. 
    
64.         for(k = m; k >= 0; k--){
    
65.                 if((len1 >= 0) && (len2 >= 0)){
    
66.                         val_p = *(p+len1);
    
67.                         val_q = *(q+len2);
    
68. //                      printf("val_p = %d, val_q = %d\n", val_p, val_q);
    
69.                 }
    
70.                 else if((len1 >= 0) && (len2 < 0)){
    
71.                         val_p = *(p+len1);
    
72.                         val_q = 0;
    
73.                 }
    
74. 
    
75.                 if(val_p >= val_q){
    
76.                         r[k] = val_p - val_q;
    
77.                 }
    
78.                 else{//借位
    
79.                         (*(p+len1-1))--;
    
80.                         r[k] = 10 + val_p - val_q;
    
81.                 }
    
82.                 len1--;
    
83.                 len2--;
    
84.         }
    
85. 
    
86.         if(p == b){
    
87.                 printf("-");
    
88.         }
    
89. 
    
90.         if(r[0] != 0){
    
91.                 printf("%d", r[0]);
    
92.         }
    
93. 
    
94.         for(k = 1; k <= m; k++){
    
95.                 printf("%d", r[k]);
    
96.         }
    
97.         printf("\n");
    
98. 
    
99.         exit(0);
    
100. }
     

复制代码