なぜMSVC

Question

に、このポアソン分布の実装の戻り0は、次の実装は常にVisual Studioのコンパイラを使用して0を返し、なぜ私は好奇心旺盛ですが、gccと打ち鳴らす

int poissrand(double lambda){ 
    int k=0; 
    double L=exp(-lambda), p=1; 
    do { 
    ++k; 
    p *= rand()/(double)INT_MAX; 
    } while (p > L); 
    return --k; 
} 

例here（VCでコンパイルした場合、それは正常に動作しません）およびhere（gcc）

Answer 1

rand()の範囲は異なります。すべてのバージョンでがstdlib.hに指定されています。

あなたがでpを掛けた値が範囲0..1にすることになっている場合は、ループの前にメインで

p *= (double)rand()/RAND_MAX;