使用Blitz++生成随机数的例子

#include <ctime>
#include <iostream>
#include <random/uniform.h>    //均匀分布
#include <random/normal.h>    //正态分布
#include <random/exponential.h> //指数分布
#include <random/discrete-uniform.h>//离散均匀分布
#include <random/beta.h> //Beta 分布
#include <random/gamma.h>//Gamma 分布
#include <random/chisquare.h>//χ^2 分布
#include <random/F.h> //F分布
using namespace ranlib;
using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{

// 连续均匀随机分布[0,1)
Uniform<double> run;
run.seed((unsigned int)time(0));
cout<<run.random()<<endl;

//    连续均匀随机分布[0,1]
UniformClosed<double> run1;
run1.seed((unsigned int)time(0));
cout<<run1.random()<<endl;

//    连续均匀随机分布(0,1)
UniformOpen<double> run2;
run2.seed((unsigned int)time(0));
cout<<run2.random()<<endl;

//    连续均匀随机分布(0,1]
UniformOpenClosed<double> run3;
run3.seed((unsigned int)time(0));
cout<<run3.random()<<endl;

//  1到99的离散均匀分布的随机数

DiscreteUniform<long> rdu(100);
rdu.seed((unsigned int)time(0));
cout<<rdu.random()<<endl;

// 正态分布
Normal<double> rnor(10,2);
rnor.seed((unsigned int)time(0));
cout<<rnor.random()<<endl;
// Beta分布有问题
//     Beta<double> rba;
//     cout<<rba.random()<<endl;
//     cout<<numeric_limits<double>::min();

// 指数分布

Exponential<double> rex(0.5);
cout<<rex.random();
//..
return 0;
}