已知有個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;
複製代碼