0

我正在尝试更改函数 extended_gcd 中 mpz_t 类型的 d 的值，然后将其传递回我的主函数，当我执行程序时出现分段错误。

    int extended_gcd(mpz_t c, mpz_t a, mpz_t b)
    {

      mpz_t x, lastx, y, lasty, temp, quotient, temp2;

      mpz_init(temp); 
      mpz_init(temp2); 
      mpz_init(a);
      mpz_init(b);
      mpz_init(c);
      mpz_init(quotient); 
      mpz_init_set_ui(x,0); 
      mpz_init_set_ui(lastx,1);
      mpz_init_set_ui(lasty,0);
      mpz_init_set_ui(y,1);

while(!(mpz_cmp(b, 0)==0))
{
    mpz_set(temp, b);
    mpz_cdiv_q(quotient, a, b);
    mpz_mod(b, a, b);
mpz_set(a, temp);
    mpz_set(temp, x);
mpz_mul(temp2, quotient, x);
mpz_sub(x, lastx, temp2);
    mpz_set(lastx, temp);
    mpz_set(temp, y);
    mpz_mul(temp2, quotient, y);
mpz_sub(y, lasty, temp2);
    mpz_set(lasty, temp);                 
}

mpz_set(c, lastx);        
mpz_clear(x);
mpz_clear(y);
mpz_clear(lastx);
mpz_clear(lasty);
mpz_clear(temp);
mpz_clear(temp2);
mpz_clear(quotient);
 return 0;

}

这是我主要的代码：

    extended_gcd(d,e,phi);  
    printf("d=");
    mpz_out_str(stdout, 10, d);
    printf("\n");

知道可能出了什么问题，我该如何解决？感谢您的时间。

4

1 回答 1

1

我看到一个问题：您不应该在extended_gcd() 的参数上调用mpz_init()。您应该只对在extended_gcd() 中声明的mpz_t 变量调用mpz_init()/mpz_clear()。

于 2012-06-14T06:26:47.380 回答