Question

我正在尝试使用GMP库在C ++中生成一个巨大的随机数，但是在查找语法时遇到了问题。这与我发现的其他例子略有不同，因为我需要为随机数设置一个楼层和天花板。这就是我需要做的事情：

mpz_class high, low;

low  = pow(2,199);
high = pow(2,210);

// code to use the high and low numbers to generate the random number

我知道这不是很多，但是，我不知道在这一点上语法甚至是什么，我已经尝试了几件事，但我找不到任何东西让我告诉GMP使用高低范围进行数字生成。

思想？

Answer 1

使用@Less提供的逻辑，我写了以下内容来解决我的问题：

void 
makeprime ()
{
    // *********************** VARIABLE DECLARATION *********************** //
    // initilize the variables as gmp class instances
    mpz_t l, rand;
    unsigned long seed;
    // perform inits to create variable pointers with 0 value
    mpz_inits(l, rand);
    //mpz_init(rand);

    // calculate the random number floor
    mpz_ui_pow_ui(l, 2, 199);

    // initilze the state object for the random generator functions
    gmp_randstate_t rstate;
    // initialize state for a Mersenne Twister algorithm. This algorithm is fast and has good randomness properties.
    gmp_randinit_mt(rstate);

    // create the generator seed for the random engine to reference 
    gmp_randseed_ui(rstate, seed);

    /*
    Function:
    int mpz_probab_prime_p (const mpz_t n, int reps)

    Determine whether n is prime. Return 2 if n is definitely prime, return 1 if n is probably prime (without being certain), 
    or return 0 if n is definitely composite.
    */
    do {
        // return a uniformly distributed random number in the range 0 to n-1, inclusive.
        mpz_urandomb(rand, rstate, 310);

        // add the random number to the low number, which will make sure the random number is between the low and high ranges
        mpz_add(rand, rand, l);

        gmp_printf("randomly generated number: %Zd\n", rand);

    } while ( !(mpz_probab_prime_p(rand, 25)) );        

    // *********************** GARBAGE COLLECTION *********************** //
    // empty the memory location for the random generator state
    gmp_randclear(rstate);
    // clear the memory locations for the variables used to avoid leaks
    mpz_clear(l);
    mpz_clear(rand);
}

感谢@Less的逻辑和帮助！

Answer 2

来自Gmp Lib Documentation

功能：void mpz_urandomb（mpz_t rop，gmp_randstate_t state，   mp_bitcnt_t n）

生成0到0范围内的均匀分布的随机整数   2 ^ n-1，包括

所以，取210 - 199，并将其用作n，生成一个随机数并将结果添加到pow（2,199）。

如果你想要比2个上限更精细的东西，这对你来说不会有用。您可以使用上述相同的技术尝试无符号的int大小的随机函数：

- 功能：unsigned long gmp_urandomm_ui（gmp_randstate_t state，unsigned long n）

返回0到n-1范围内的均匀分布的随机数。

在这里，您可以找到您的粒度范围并将其用于n。然后将随机数添加到较低的值。限制是n必须小于MAXUINT，通常是2 ^ 32 -1

C ++ GMP生成随机数

2 个答案: