大数相乘问题

        我们在现实的计算中难免会遇到益处问题，加法溢出，乘法溢出等，因为我们的类型变量表示的范围是有限的。这里我们讨论下有溢出的乘法问题的解决。这种方法我们可以通过两个数组来存储乘数和被乘数来解决问题。

代码1：

#include <iostream> using namespace std; void Multi(int *num1,int size1,int *num2,int size2) { int size = size1 + size2 - 1; int *ret = new int[size]; memset(ret,0,sizeof(int) * size); for (int i = size2-1; i >= 0; --i) { int k = size2-i-1; for (int j = size1-1; j >= 0; --j) { int temp = num2[i] * num1[j] + ret[k]; if (temp > 10) { ret[k] = temp % 10; ret[k+1] += temp / 10; } else { ret[k] = temp; } ++k; } } cout << "两数相乘后的结果为: " << endl; for (i = size-1; i >= 0; --i) { cout << ret[i]; } cout << endl; } int main() { int a[5] = {1,2,3,4,5}; int b[3] = {4,5,6}; Multi(a,5,b,3); return 0; }

结果运行正确，但是代码中有点小错误，不知道大家能发现吗？

好的，现在我们分析一下上述代码:

int size = size1 + size2 - 1;

定义新数组的大小等于两个数组大小之和减1；至于这个1是怎么来的呢？我们可以从两数乘法中归纳出来。我们可以用我们初中或是小学学的竖式相乘来观看：

         1  2  3  4  5

        *       4  5  6

 ------------------------

         7 4   0  7  0

      6 1 7   2  5

  4  9 3 8   0

-------------------------

  5  6 2 9   3  2  0

 

我们可以看出从第二行开始依次往左移一位。移动的位数恰为size2 - 1; 而本来的行数等于size1;

所以新数组的长度就是 size1 + size2 - 1;

哈哈，错误就在这里了，在这里我有点想当然了，这个例子用到的中间数组的长度确实是size1 + size2 - 1；

然而若是我们的书改变一下呢。即若是最高位的数字相乘大于10呢，此时又会有进位产生。所以真正的数组长度应该为:

size = size1 + size2;

接着看代码:

int *ret = new int[size]; memset(ret,0,sizeof(int) * size);

动态分配一个数组，并将它各个空间初始化为0；

for (int i = size2-1; i >= 0; --i) { int k = size2-i-1; for (int j = size1-1; j >= 0; --j) { int temp = num2[i] * num1[j] + ret[k]; if (temp > 10) { ret[k] = temp % 10; ret[k+1] += temp / 10; } else { ret[k] = temp; } ++k; } }

从乘数的最低位开始依次乘被乘数的每一位，然后与当前对应的数组空间的值相加。

然后判断这个值是否大于10，若大于10，我们将其余数放入当前空间，商+=于高位值并把结果放入该高位。这样一次进行下去，最后得到的结果便是我们的最后结果了。数组该数组即可。

代码2：

#include <iostream> using namespace std; void Multi(int *num1,int size1,int *num2,int size2) { int size = size1 + size2 - 1; int *ret = new int[size]; memset(ret,0,sizeof(int) * size); /* 遍历乘数的每一个数，让它成语被乘数的每一位然后存到新的数组中. 然后判断新数组的每一位时候大于10，然后除10进位，数组留在原位。 */ for (int i = size2-1; i >= 0; --i) { int k = size2 - i - 1; for (int j = size1 - 1; j >= 0; --j) { ret[k++] += num2[i] * num1[j]; } } for (i = 0; i < size; ++i) { if (ret[i] >= 10) { ret[i+1] += ret[i] / 10; ret[i] %= 10; } } cout << "两数相乘后的结果为: " << endl; for (i = size-1; i >= 0; --i) { cout << ret[i]; } cout << endl; } int main() { int a[5] = {1,2,3,4,5}; int b[3] = {4,5,6}; Multi(a,5,b,3); return 0; }

这个代码的大概思想和上面那个差不多。就不在解释了。。呵呵。。