1

我正在尝试实现 Karatsuba 乘法算法。我有点遵循这个 wiki页面中的伪代码。但我总是收到这个错误：

由信号 SIGSEGV 终止（地址边界错误）

当我用其他东西替换导致递归发生的行时：

z0 = multiply(a, c);
z1 = multiply(b, d);
z2 = multiply(a+b, c+d);

错误消失了。

这是我的代码：

#include <iostream>
#include <math.h>

long int multiply(int x, int y);
int get_length(int val);

int main()
{
  int x = 0, y = 0;
  long int result = 0;

  std::cout << "Enter x: ";
  std::cin >> x;
  std::cout << "Enter y: ";
  std::cin >> y;

  result = multiply(x, y);
  std::cout << "Result: " << result << std::endl;
  return 0;
}

long int multiply(int x, int y)
{
  if(x < 10 || y < 10) {
    return x * y;
  }

  int x_len = get_length(x);
  int y_len = get_length(y);

  long int z0 = 0 , z1 = 0, z2 = 0;
  int a = 0, b = 0, c = 0, d = 0;

  a = x / pow(10, x_len);
  b = x - (a * pow(10, x_len));
  c = y / pow(10, y_len);
  d = y - (c * pow(10, y_len));

  z0 = multiply(a, c);
  z1 = multiply(b, d);
  z2 = multiply(a+b, c+d);

  return (pow(10, x_len) * z0) + (pow(10, x_len/2) * (z2 - z1 - z0)) + z1;
}

int get_length(int val)
{
  int count = 0;
  while(val > 0) {
    count++;
    val /= 10;
  }
  return count;
}
4

2 回答 2

2

我找到了问题的原因。这是因为这些行：

a = x / pow(10, x_len);
b = x - (a * pow(10, x_len));
c = y / pow(10, y_len);
d = y - (c * pow(10, y_len));

它应该是x_len / 2代替x_len和相同的y_len。因为它导致递归是无限的。

于 2016-08-22T21:04:43.027 回答
0

您正在使用该pow函数进行整数幂。它不是整数函数。pow编写适合您应用程序的您自己的函数。例如：

int pow(int v, int q)
{
    int ret = 1;
    while (q > 1)
    {
        ret*=v;
        q--;
    }
    return ret;
}

确保int pow(int, int);在顶部放置一个。

于 2016-08-22T20:56:27.317 回答