有没有一个 c++ gmp 库函数与 python gmpy2 库 divm(..） 函数相同?

Is there a c++ gmp library function that does the same as the python gmpy2 library divm(...) function?

本文关键字：divm 函数 gmpy2 c++ 有一个 gmp 库函数 python      更新时间：2023-10-16

正如标题所说，我试图在C++ gmp库中找到一个与gmpy2 python库的divm(...(函数相同的函数。

我无法想象 gmpy2 有这种方法，而 gmp C++ 库没有任何东西可以做同样的计算。

如果这不存在，我将不胜感激任何关于从头开始制作divm 函数的建议(仍然必须使用 gmp，因为我正在通过 mpz_class 处理大于标准整数的值(。

谢谢！

该函数的 C 源代码可以在这里找到，在GMPy_MPZ_Function_Divm中。如您所见，它并不是与底层 C 代码的一对一映射，但是，当您剥离所有 Python 引用计数内容时，它看起来非常基本：

// numz, denz, and modz are copies of the arguments.
// resz, gcdz ar temporaries.
if (mpz_invert(resz, denz, modz)) {    // inverse exists?
ok = 1;
} else {                               // num, den AND mod have a gcd > 1?
mpz_init(gcdz);
mpz_gcd(gcdz, numz, denz);
mpz_gcd(gcdz, gcdz, modz);
mpz_divexact(numz, numz, gcdz);
mpz_divexact(denz, denz, gcdz);
mpz_divexact(modz, modz, gcdz);
mpz_clear(gcdz);
ok = mpz_invert(resz, denz, modz);
}
if (ok) {
mpz_mul(resz, resz, numz);
mpz_mod(resz, resz, modz);
mpz_clear(numz);
mpz_clear(denz);
mpz_clear(modz);
return resz;
} else {
mpz_clear(numz);
mpz_clear(denz);
mpz_clear(modz);
return NULL;
}

在数学上，它只是计算divm(a, b, m)的表达式b-1* a % m，当然假设该表达式是有效的(例如，b不是零(。

相关文章：