如何使用GMP

Question

我試圖做到這一點，如下所示產生一定範圍內的隨機數：

int main(){ 
    mpz_t v; 
    mpf_class u,w,x[3][200];  

    for(i=0;i<n;i++){ 
    for(j=0;j<3;j++){ 
    mpz_random(v,1000000); 
    mpf_div_ui(w.get_mpf_t(),v,1000000); 
    mpf_div(u.get_mpf_t(),1,2); 
    mpf_sub(w.get_mpf_t(),w.get_mpf_t(),u); 
    mpf_mul(x[i][j].get_mpf_t(),w.get_mpf_t(),2); 
    } 
    } 
} 

但使用整數和漂浮起來不起作用。基本上嘗試分配隨機整數，除以最大整數可以給一個數字（不一定是整數）在0和1之間。減一半並乘以2給出一個-1到1之間的隨機數，這是我想要的，但是就像我說的那樣由於類型導致問題。

Answer 1

如果我正確理解你的程序的意圖，這應該是或多或少你所追求的。

int main(){ 
    gmp_randstate_t state; 
    mpf_class v, w, x[3][200]; 

    gmp_randinit_mt(state); 
    for(int i=0;i<3;i++){ 
    for(int j=0;j<200;j++){ 
    mpf_urandomb(v.get_mpf_t(), state, 4); 
    mpf_urandomb(w.get_mpf_t(), state, 256); 
    if (v < .5) w *= -1; 
    x[i][j] = w; 
    } 
    } 
} 

的state變量用於告知隨機數例行你想使用你的隨機數生成器哪種算法。在上面的代碼中，我已經初始化它以使用Mersenne Twister算法。

你沒有正確申報您的i和j指數。

到mpf_urandomb()第一呼叫得到0和1之間，用來決定是否下一個數應該是負數或不是4比特的隨機值。

第二次調用mpf_urandomb()會得到一個介於0和1之間的256位隨機值。如果前一個數字小於.5，則該值將被取反。

該值被分配給x[i][j]。