C++11随机数发生器

前言

一直知道所谓的"随机数"都是伪随机，事实上也是满足某种规则生成的。有些程序测试时通常需要一个随机数源，但在新标准出现之前，C++都是依赖简单的C库函数rand来生成随机数的。最近突然看到了C++11中的随机数发生器，简直如获珍宝，下面会谈一谈这两者的区别。

如何产生随机数

利用C库函数

头文件<stdlib.h>，但是注意在linux下stdlib.h包含srandom 和random ，但在VC下stdlib.h包含的是srand和rand

x = rand()%11; /*产生1~10之间的随机整数*/

y = rand()%51 - 25; /*产生-25 ~ 25之间的随机整数*/

z = ((double)rand()/RAND_MAX)*(b-a) + a;/*产生区间[a,b]上的随机数*/

整理一下常见的产生随机数的通用表达公式为:

取得(0,x)的随机整数：rand()%x；

取得(a,b)的随机整数：rand()%(b-a)；

取得[a,b)的随机整数：rand()%(b-a)+a；

取得[a,b]的随机整数：rand()%(b-a+1)+a；

取得(a,b]的随机整数：rand()%(b-a)+a+1；

取得0-1之间的浮点数：rand()/double(RAND_MAX)

比如说随便写个函数：

#define random(x) (rand()%x)
int main()
{
	for (int i = 0; i < 10; ++i)
		cout << random(11) << " ";
	cout << endl;
	system("pause");
}

你会发现两次运行这一段程序输出结果是一样的，说明rand这个函数所生成的随机数是一次性的，为了改进这一点，我们可以利用srand()函数，srand()是用来设置rand()产生随机数时的随机数种子。在调用rand()函数产生随机数前，必须先利用srand()设好随机数种子（seed）, 如果未设随机数种子, rand()在调用时会自动设随机数种子为1。上面的两个例子就是因为没有设置随机数种子，每次随机数种子都自动设成相同值1 ，进而导致rand()所产生的随机数值都一样。

而我们一般会用时间来做为种子:

#include<iostream>
#include<cstdlib>
#include<ctime>
using namespace std;

#define random(a,b) (((double)rand()/RAND_MAX)*(b-a)+a)

void main() {
	srand((int)time(0));

	for (int i = 0; i < 100; i++) {
		cout << random(0, 10) << " ";
	}
	cout << endl;

	system("pause");
}

C++随机数发生器

随机数发生器=随机数引擎+随机数分布类

由于C库函数生成的为均匀分布的伪随机整数，但是如果我们需要非均匀分布的数的话，就会手工转换rand生成的随机数的范围、类型和分布，这样又引入了非随机性，不妥。所以我们可以用定义在头文件random中的随机数库通过一组写作的类来解决这些问题:随机数引擎和随机数分布类。引擎是用来生成随机unsigned整数序列的，分布是使用引擎返回服从特定概率分布的随机数。

首先利用引擎来产生随机数:

#include<iostream>
#include<random>
#include<ctime>
using namespace std;

void main() {
	default_random_engine e;   //生成随机无符号数
	for (size_t i = 0; i < 10; ++i)
		cout << e() << " ";

	cout << endl;
	cout << "seed Random:" << endl;

	e.seed(int(time(0)));
	for (size_t i = 0; i < 10; ++i)
		cout << e() << " ";
	system("pause");
}

//OUTPUT:
3499211612 581869302 3890346734 3586334585 545404204 4161255391 3922919429 949333985 2715962298 1323567403
seed Random:
1209134335 963543660 450936573 838635935 555864927 2629301182 2752494759 772597830 1658101593 1483430385 

但是一般来说，随机数引擎的输出是不能直接使用的，因为生成的随机数的值范围通常与我们需要的不符，所以还要加上分布类型，我们再给个完整的例子体会一下:

#include<iostream>
#include<random>
#include<string>
using namespace std;

vector<unsigned> randVec()
{
	static default_random_engine e;
	static uniform_int_distribution<unsigned> u(0, 9);
	vector<unsigned> ret;
	for (size_t i = 0; i < 20; ++i)
		ret.push_back(u(e));
	return ret;
}

vector<float> normalVec()
{
	cout << "正态分布:" << endl;
	static default_random_engine e;
	static normal_distribution<float> n(4, 1.5);//均值4，方差1.5
	vector<float> res;
	vector<unsigned> vals(9);
	for (size_t i = 0; i != 20; ++i)
	{
		unsigned tmp = n(e);
		unsigned v = lround(tmp);
		if (v < vals.size())
			++vals[v];
		res.push_back(tmp);

	}
	for (size_t j = 0; j != vals.size(); ++j)
		cout << j << ": " << string(vals[j], '*') << endl;
	return res;
}

template<typename T>
void display(vector<T> a)
{
	for (int i = 0; i < a.size(); ++i)
	{
		cout << a[i] << " ";
	}
	cout << endl;
}

void main() {
	vector<unsigned> res;
	res=randVec();
	cout << "均匀分布:" << endl;
	display(res);
	vector<float> res2;
	res2 = normalVec();
	display(res2);
	system("pause");
}

结果如下：

更多的分布参数设置见下图: