c - 如何修复适用于某些值但不适用于其他值的代码

Question

当我输入某些值时，我得到答案，但其他值则没有，系统保持阻塞状态

int pgcdfon(mpz_t a, mpz_t b, mpz_t *x, mpz_t *y) {
    if (mpz_cmp_ui(a, 0) == 0) {
        mpz_set_ui(*x, 0);
        mpz_set_ui(*y, 1);
        return b;
    }
    mpz_t x1, y1, res, h;
    mpz_init(x1);
    mpz_init(y1);
    mpz_init(res);
    mpz_t o;
    mpz_init(o);
    mpz_init(h);
    mpz_mod(res, b, a);
    int pgcd = pgcdfon(res, a, &x1, &y1);
    mpz_div(o, b, a);
    mpz_mul(h, o, x1);
    mpz_sub(*x, y1, h);
    mpz_set(*y, x1);
    return pgcd;
}

void main() {
    mpz_t x, y, q, w, g;
    mpz_init(x);
    mpz_init(y);
    mpz_init(q);
    mpz_init(w);
    mpz_init(g);
    gmp_printf("donner q \n ");
    gmp_scanf("%Zd", &q);
    gmp_printf("donner w \n");
    gmp_scanf("%Zd", &w);
    mpz_set(g, (pgcdfon(q, w, &x, &y)));
    gmp_printf("\n\n pgcd  de (%Zd, %Zd) = %Zd  \n\n", q, w, g);
}

Answer 1

类型mpz_t是一个(伪装的)数组类型，因此不能从函数返回。您必须更改函数的原型(prototype)以返回 out 参数中的值，而不是返回它。这意味着，像这样的函数

mpz_t myfun (...) {
     mpz_t ret;
     mpz_init(ret);

     // fill ret with someting

     return ret;
}

无法工作，但必须重写为:

void myfun (..., mpz_t *out) {
     mpz_t ret;
     mpz_init(ret);

     // fill ret with stuff
     mpz_init(*out); // if not initialized by caller
     mpz_set(*out, ret);
}

将此技术应用于您的函数会产生

void pgcdfon(mpz_t a, mpz_t b, mpz_t *x, mpz_t *y, mpz_t *ret) {
    if (mpz_cmp_ui(a, 0) == 0) {
        mpz_set_ui(*x, 0);
        mpz_set_ui(*y, 1);
        mpz_init(*ret);
        mpz_set(*ret, b);
        return;
    }
    mpz_t x1, y1, res, h;
    mpz_init(x1);
    mpz_init(y1);
    mpz_init(res);
    mpz_t o;
    mpz_init(o);
    mpz_init(h);
    mpz_mod(res, b, a);
    mpz_t pgcd;
    pgcdfon(res, a, &x1, &y1, &pgcd);
    mpz_div(o, b, a);
    mpz_mul(h, o, x1);
    mpz_sub(*x, y1, h);
    mpz_set(*y, x1);
    mpz_init(*ret);
    mpz_set(*ret, pgcd);
}

int main() {
    mpz_t x, y, q, w, g;
    mpz_init(x);
    mpz_init(y);
    mpz_init(q);
    mpz_init(w);
    mpz_init(g);
    gmp_printf("donner q \n ");
    gmp_scanf("%Zd", &q);
    gmp_printf("donner w \n");
    gmp_scanf("%Zd", &w);
    pgcdfon(q, w, &x, &y, &g);
    gmp_printf("\n\n pgcd  de (%Zd, %Zd) = %Zd  \n\n", q, w, g);
}

它应该按您的预期工作。