问题:rand7能生成1-7的整数随机数。要求利用rand7生成1-10的整数随机数。最初发这个文的时候有错误,现在改正过来。另外提出并解决一个更实际的问题:如何按照正态分布模拟生成一个班的学生成绩。感慨一下知识的浩瀚和自己的渺小,并感谢一下楼主minglinag,hpsmouse,以及众多网友。可以参看原帖.
1. 舍去法
rand7生成7个整数,没有办法均匀的映射成10个整数,但是运行两次rand7可以生成49个数字,如果这49个数字是均匀分布的,舍去多余的9个,剩下的40个正好可以用模10运算映射到10个整数上。
代码1
int i;
do
{
i = 7 * (rand7() - 1) + rand7(); // it is now uniformly random between 1 and 49
} while(i > 40); // it is now uniformly random between 1 and 40
return i % 10 + 1; // result is now uniformly random between 1 and 10
这个算法做到了从40个数字均匀映射到1到10,这个说法有些抽象,进一步形象的来说明,考虑如下的种子矩阵:
int seed7[7][7] = {
{1 , 2 , 3 , 4 , 5 , 6 , 7},
{8 , 9 , 10, 1 , 2 , 3 , 4},
{5 , 6 , 7 , 8 , 9 , 10, 1},
{2 , 3 , 4 , 5 , 6 , 7 , 8},
{9 , 10, 1 , 2 , 3 , 4 , 5},
{6 , 7 , 8 , 9 , 10, 0 , 0},
{0 , 0 , 0 , 0 , 0 , 0 , 0}
};
如果用x=0...6,y=1...7,则
变换i=7x+y与矩阵中每个元素位置与(x,y)唯一对应,也就是x选择行,y选择列,如果x和y都是均匀分布,那么这49个
位置有相同的被选中的
概率=1/49。下面这行代码实现了这个变换:
i = 7 * (rand7() - 1) + rand7();
其中7*(rand7()-1)相当于选择种子矩阵中的行,第二个rand7相当于选择列,而最后的模10+1运算,就恰好生成了矩阵中每个元素的值,但是(i>40)这个循环条件把最后的9个值变为了0。因此
代码1等价于在种子矩阵中做选择。
再来看,在这个矩阵中1-10数字一共出现了40次,每个数字出现4次,0出现了9次。我们在这个表里筛选,如果选中了0则重新选一次,如果非0则返回结果,那么1-10恰好可以得到 平均机会。(更理论也更准确地说,这是一个 条件概率。我们看一下在不选0的前提下选择到1的概率,令B=不选择零,A=选择到1,那么 P(A)=4/49,P(B)=40/49. 因为只要选择到1,B就必然成立,所以P(B|A)&