如何使boost :: exponential_distribution工作？

Question

我一直在玩Boost.Random一天，而boost::uniform_int_distribution<>效果很好，我遇到了問題boost::exponential_distribution<>。

一個簡單的程序是勝過千言萬語：

#include <iostream> 
#include <boost/random/mersenne_twister.hpp> 
#include <boost/random/exponential_distribution.hpp> 

int main() { 
    boost::mt19937 gen; 
    boost::exponential_distribution<> dis; 
    std::cout << dis(gen) << "\n"; 
    return 0; 
} 

編譯鏗鏘3.0，使用Boost 1.39.1（不，我不能升級加速）。

輸出總是相同的：nan。

我找不到任何報告的問題，所以我想這是我沒有正確使用庫...任何線索將不勝感激。

Answer 1

給定隨機數r，均勻分佈在[0,1）上，-log(r)分佈在（0，無窮遠），分佈爲exp(-x)。 前者是boost::uniform_01()。

如果你需要一個分配p exp(-px)，那麼它的-(1/p)log(r)。

在這兩種情況下，這裏的log(x)都是自然對數（基數爲e）。

UPD：使用boost::variate_generator似乎爲我工作（提高1.43）：

#include <iostream> 
#include<boost/random.hpp> 

int main() { 
    boost::mt19937 rng(11u); 
    boost::variate_generator< boost::mt19937&, boost::exponential_distribution<> > rndm(rng, boost::exponential_distribution<>()) ; 
    std::cout<<rndm()<<"\n"; 
} 

我沒有，雖然檢查的分佈。