c++ - 如何使用 MPIR 库计算 GCD

Question

我在使用算术库方面相对缺乏经验。我需要尝试编写一些代码来计算最大公分母，使用MPIR库中的一个函数。我找到了几种方法，我不明白他们定义变量的方式。例如：关于：

void mpz_gcd (mpz_t rop, mpz_t op1, mpz_t op2) 

我无法理解变量类型以及如何使用它们。谁能向我提供一些简单的代码来向我澄清这一点？

Answer 1

但无法理解变量类型。

如果它是您关心的变量类型mpz，那么它就是任意精度整数（与实数相反）。

就变量本身而言，rop是返回变量，将设置为op1和的 GCD op2。

例如，这里有一个完整的程序来说明该函数的使用：

#include <stdio.h>
#include <mpir.h>

int main (void) {
    mpz_t a, b, c;
    mpz_init (a); mpz_init (b); mpz_init (c);

    mpz_set_str (a, "1024", 10);
    mpz_set_str (b, "768", 10);
    mpz_gcd (c, a, b);

    printf ("GCD of "); mpz_out_str (stdout, 10, a);
    printf (" and ");   mpz_out_str (stdout, 10, b);
    printf (" is ");    mpz_out_str (stdout, 10, c);
    putchar ('\n');

    return 0;
}

输出是：

GCD of 1024 and 768 is 256

---

我想生成长位的随机值。

但是，如果（正如您在评论中提到的，虽然很难看出 GCD 计算的适用性）您需要一个 100 位的随机数，您可以使用类似的东西：

#include <stdio.h>
#include <time.h>
#include <mpir.h>

int main (void) {
    mpz_t a;
    gmp_randstate_t rstate;

    mpz_init (a);
    gmp_randinit_default (rstate);
    gmp_randseed_ui (rstate, time (0));

    mpz_urandomb (a, rstate, 100);
    printf ("100 random bits gives: "); mpz_out_str (stdout, 10, a);
    putchar ('\n');

    return 0;
}