随机数发生器碰撞测试中碰撞太多

Question

我期待lkolizji变量在128左右，但是对于大量的生成数字和“盒子”来说，它要高得多。较小数字的结果很好。我不知道为什么会发生这种情况。这里是我的代码，带有错误答案的示例参数。好结果（约128）的示例是int lPrzedzialow = 1000000; int iLiczb = 16000;

#include <iostream> 
#include <gsl/gsl_rng.h> 
#include <stdlib.h> 
#include<cmath> 
#include <algorithm> 
using namespace std; 
int main (void) 
{ 
//Random number 
    unsigned int seed=2596524; 
    gsl_rng * r=gsl_rng_alloc (gsl_rng_mt19937); 
    gsl_rng_set(r,seed); 
    gsl_rng_env_setup(); 
//Parameters 
    int lPrzedzialow=10000000000;//number of boxes 
    int iLiczb = 1600000;//number of random numbers 
    int z,lKolizji=0;//lKolizji holds collision number 
    vector<int> lwKomorkach(iLiczb);//number of boxes of random numbers 
    long double dlPrzedzialu=1./(lPrzedzialow); 
//number of box of a random number 
    for (int i = 0; i < iLiczb; i++) 
    { 
     lwKomorkach[i] = floor((gsl_rng_uniform (r)/dlPrzedzialu)); 
    } 
//sorting 
    sort(lwKomorkach.begin(), lwKomorkach.end()); 
//how many collisions 
    for(z=0;z<=iLiczb-1;z++) 
    { 
     if(lwKomorkach[z+1]==lwKomorkach[z]){lKolizji++;} 
    } 
    double pdf[lKolizji]; 
    pdf[0]=exp(-128); 
    double spdf=exp(-128); 
    for(int h=1;h<lKolizji;h++){ 
    pdf[h]=pdf[h-1]*128./(h); 
    spdf+=pdf[h]; 
    } 
    double pwyzsze=1.-spdf; 
    cout<<endl<<lKolizji<<" "<<spdf<<" "<<pwyzsze<<endl; 
    gsl_rng_free (r); 
    return 0; 
    } 

Answer 1

这个数字：10000000000，对于32位整数来说太大了。事实上，它相当于1,410,065,408，大约相信它的大小的1/7。