Question

我试图实现Karatsuba算法进行乘法运算。我有点关注这个wiki page中的伪代码。但我总是得到这个错误：

由信号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;
}

Answer 1

我发现问题的原因。这是因为这些方面：

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相同。因为它会导致递归无限。

Answer 2

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

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

确保将int pow(int, int);放在顶部。

在递归函数中由信号SIGSEGV（地址边界错误）终止

2 个答案: