Visual C++2008 C++ TR1 随机数编程笔记

喜欢ヅ旅行 2022-09-30 14:49 196阅读 0赞

At the core of any pseudorandom number generation software is a routine for generating uniformly distributed random integers.  
In C++ TR1 you have your choice of several core generators that it calls “engines.” The following four engine classes are supported in the Visual Studio 2008 feature pack

(微软已经发布了Visual Studio 2008的Services Pack 1，它包含了此前发布的feature pack，以及完整的TR1支持，后来又发布了一个修正：[VC 2008 SP1: Problems with STL/TR1 after installing VS2008 SP1][VC 2008 SP1_ Problems with STL_TR1 after installing VS2008 SP1]，关于：[VC9 SP1 Hotfix For The vector<function<FT>> Crash][VC9 SP1 Hotfix For The vector_function_FT_ Crash]，关于文档：[微软TR1文档][TR1]).

*  **linear\_congruential** uses a recurrence of the form x(i) = (A \* x(i-1) + C) mod M
 *  **mersenne\_twister** implements the famous Mersenne Twister algorithm
 *  **subtract\_with\_carry** uses a recurrence of the form x(i) = (x(i - R) - x(i - S) - cy(i - 1)) mod M in integer arithmetic
 *  **subract\_with\_carry\_01** uses a recurrence of the form x(i) = (x(i - R) - x(i - S) - cy(i - 1)) mod 1 in floating point arithmetic

Each engine has a seed() method that accepts an unsigned long argument to specify the random number generation seed. It is also possible to set the seed in more detail using template parameters unique to each engine.

微软的C++ TR1可以生成以下分布的随机数：

<table style="width: 420px;">    
 <tbody> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a Bernoulli distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a binomial distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates an exponential distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a gamma distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a geometric distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a normal distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a Poisson distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a uniform integer distribution.</span></td> 
  </tr> 
  <tr> 
   <td><span style="font-family: Verdana; font-size: x-small;">Generates a uniform floating-point distribution.</span></td> 
  </tr> 
 </tbody> 
</table>

**试验一：在VS2008中先建一空的VC++项目文件，然后添加新建项cpp文件如下：**

\#include <random>  
\#include <iostream>

void main()\{  
std::tr1::mt19937 eng; // a core engine class：Mersenne Twister generator  
std::tr1::normal\_distribution<double> dist;  
std::tr1::uniform\_int<int> unif(1, 52);  
  
for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数  
std::cout << dist(eng)<<std::endl;  
  
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数  
std::cout << unif(eng) << std::endl;  
\}

在循环中，每循环一次，就调用mt19937 eng一次，产生一个随机数输出。

**试验二：关于种子seed**

\#include <random>  
\#include <iostream>  
\#include <time.h>

void main()\{  
std::tr1::mt19937 eng; // a core engine class：Mersenne Twister generator  
std::tr1::normal\_distribution<double> dist;  
std::tr1::uniform\_int<int> unif(1, 52);  
eng.seed((unsigned int)time(NULL)); // reseed base engine 设置种子用\#include <time.h>, 不能用\#include <time>  
  
for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数  
std::cout << dist(eng)<<std::endl;  
//eng.seed(); // reseed base engine  
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数  
std::cout << unif(eng) << std::endl;  
\}

**试验三：随机数写入文件**

\#include <random>  
\#include <iostream>  
\#include <fstream>  
\#include <time.h>

using namespace std;  
using namespace std::tr1;

void main()  
\{  
mt19937 eng; // a core engine class：Mersenne Twister generator  
normal\_distribution<double> dist;  
uniform\_int<int> unif(1, 52); //uniform\_int 类以相同的概率在一个范围内抽取整数，它的构造函数有两个参数，分别表示抽取范围的最大值和最小值，需要注意的是抽取范围是个闭区间，也就是这两个值也可能被抽取到。

eng.seed((unsigned int)time(NULL)); // 设置种子用\#include <time.h>, 不能用\#include <time>

for (int i = 0; i < 10; ++i) //产生正态分布的10个随机数  
cout << dist(eng)<<endl;  
  
 ofstream fileout("fileout.dat");  
for(int i = 0; i < 5; ++i) //产生均匀分布的在1到52之间的五个整数随机数(含1和52)，即\[1,52\]闭区间  
fileout << unif(eng)<< endl;  
  
 fileout.close();  
\}

**试验四：第三方"Mersenne Twister"随机数生成程序使用试验**(程序来源：Agner Fog[http://www.agner.org/random/][http_www.agner.org_random])

// 使用说明:从网站下载压缩包，[http://www.agner.org/random/randomc.zip][http_www.agner.org_random_randomc.zip]  
// 展开后,将其中的randomc.h头文件及mersenne.cpp文件Copy到项目文件夹,  
// 并将它们加入到项目中，其中包括"Mersenne Twister"的实现

\#include <iostream>  
\#include <time.h>  
\#include "randomc.h" // define classes for random number generators

using namespace std;

void main()  
\{

int seed = (int)time(0); // random seed

// choose one of the random number generators:  
CRandomMersenne RanGen(seed); // make instance of random number generator  
cout<<"/n/nRandom integers in interval from 0 to 99:/n";  
for (int i = 0; i < 40; i++) \{  
int ir = RanGen.IRandom(0,99);  
 cout<<ir<<" ";  
\}  
  
cout <<endl;

cout<<"/n/n/n/nRandom floating point numbers in interval from 0 to 1:/n";  
for (int i = 0; i < 40; i++) \{  
float fr = RanGen.Random();  
 cout<<fr<<" ";  
\}  
  
cout <<endl;  
\}

**试验五：第三方"Mother-Of-All"随机数生成程序使用试验**(程序来源：Agner Fog[http://www.agner.org/random/][http_www.agner.org_random])

// 使用说明:从网站下载压缩包，[http://www.agner.org/random/randomc.zip][http_www.agner.org_random_randomc.zip]  
// 展开后,将其中的randomc.h头文件及mother.cpp文件Copy到项目文件夹,  
// 并将它们加入到项目中，其中包括"Mother-Of-All" generator invented by George Marsaglia 的实现

\#include <iostream>  
\#include <time.h>  
\#include "randomc.h" // define classes for random number generators

using namespace std;

void main()  
\{

int seed = (int)time(0); // random seed

// choose one of the random number generators:  
CRandomMother RanGen(seed); // make instance of random number generator  
cout<<"/n/nRandom integers in interval from 0 to 99:/n";  
for (int i = 0; i < 40; i++) \{  
int ir = RanGen.IRandom(0,99);  
cout<<ir<<" ";  
\}  
  
cout <<endl;

cout<<"/n/n/n/nRandom floating point numbers in interval from 0 to 1:/n";  
for (int i = 0; i < 40; i++) \{  
float fr = RanGen.Random();  
cout<<fr<<" ";  
\}  
  
cout <<endl;  
\}

**试验六：第三方"SFMT"随机数生成程序使用试验**(程序来源：Agner Fog[http://www.agner.org/random/][http_www.agner.org_random])

重要提示：在编译前，可以修改头文件sfmt.h中的\#define MEXP以及下面相应的宏代码：Choose one of the possible Mersenne exponents. Higher values give longer cycle length and use more memory。SFMT利用了SSE2指令，速度最快，但只适合intel系列的部分芯片。

// 使用说明:从网站下载压缩包，[http://www.agner.org/random/randomc.zip][http_www.agner.org_random_randomc.zip]  
// 展开后,将其中的头文件及sfmt.cpp文件Copy到项目文件夹,  
// 并将它们加入到项目中，其中包括"SFMT" 的实现

\#include <iostream>  
\#include <time.h>  
\#include "sfmt.h" // define classes for random number generators

using namespace std;

void main()  
\{

int seed = (int)time(0); // random seed

// choose one of the random number generators:  
CRandomSFMT1 RanGen(seed); //注意可以是CRandomSFMT，是CRandomSFMT0，或CRandomSFMT1  
cout<<"/n/nRandom integers in interval from 0 to 99:/n";  
for (int i = 0; i < 40; i++) \{  
int ir = RanGen.IRandomX(0,99);  
cout<<ir<<" ";  
\}  
  
cout <<endl;

cout<<"/n/n/n/nRandom floating point numbers in interval from 0 to 1:/n";  
for (int i = 0; i < 40; i++) \{  
float fr = RanGen.Random();  
cout<<fr<<" ";  
\}  
  
cout <<endl;  
\}

[VC 2008 SP1_ Problems with STL_TR1 after installing VS2008 SP1]: https://connect.microsoft.com/VisualStudio/Downloads/DownloadDetails.aspx?DownloadID=15359
[VC9 SP1 Hotfix For The vector_function_FT_ Crash]: http://blogs.msdn.com/vcblog/archive/2008/12/17/vc9-sp1-hotfix-for-the-vector-function-ft-crash.aspx
[TR1]: http://msdn.microsoft.com/en-us/library/bb982198.aspx
[http_www.agner.org_random]: http://www.agner.org/random/
[http_www.agner.org_random_randomc.zip]: http://www.agner.org/random/randomc.zip