已知有个rand7()的函数,返回1到7随机天然数,利用这个rand7()构造rand10() 随机1~10。bash
要保证rand10()在整数1-10的均匀分布,能够构造一个1-10n的均匀分布的随机整数区间(n为任何正整数)。 假设x是这个1-10n区间上的一个随机整数,那么x%10+1就是均匀分布在1-10区间上的整数。 (rand7()-1)*7+rand7()能够构造出均匀分布在1-49的随机数:dom
rand7()-1 {0,1,2,3,4,5,6}
(rand7()-1)*7 A={0,7,14,21,28,35,42}
rand7() B={1,2,3,4,5,6,7}
(rand7()-1)*7+rand7() 生成的整数均匀分布在1-49之间
复制代码
因为(rand7()-1)*7+rand7()能够构造出均匀分布在1-49的随机数(缘由见下面的说明),能够将41~49这样的随机数剔除掉,获得的数1-40仍然是均匀分布在1-40的,这是由于每一个数均可以当作一个独立事件。 代码以下:函数
int rand_10()
{
int x = 0;
do
{
x = 7 * (rand7() - 1) + rand7();
}while(x > 40);
return x % 10 + 1;
复制代码
已知random3()这个随机数产生器生成[1, 3]范围的随机数,请用random3()构造random5()函数,生成[1, 5]的随机数?ui
int random5()
{
int val = 0;
do
{
val = 3 * (random3() - 1) + random3();
}while(val > 5);
return val;
复制代码
如何从[1, 9]范围的数生成[1, 5]的数呢? 就是 rejection sampling 方法,即生成[1, 9]的随机数,若是数的范围不在[1, 5]内,则从新取样spa
已知random_m()随机数生成器的范围是[1, m] 求random_n()生成[1, n]范围的函数,m < n && n <= m *m 通常解法:code
int random_n()
{
int val = 0;
int t; //t为n的最大倍数,且知足t<m*m
do
{
val = m * (random_m() - 1) + random_m();
}while(val > t);
return val;
复制代码