算法导论
5.1.2 题目这样描述: 给出RANDOM(a,b)函数的一种实现,它只调用RANDOM(0,1)。 并给出程序的期望运行时间。
想法:构造出一个样本集合U,其包含的元素个数为b-a+1,这样可以实现与a~b的一一映射,那么问题转化为,从集合U中等概率的取出一个元素。
然而手头只有RANDOM(0,1)函数这一个资源,它能做的事情就是可以等概率的获得0或1,很自然的想到用它来构造一个只包含0和1的二进制序列集合。
构造出的二进制集合很容易的和十进制的数字产生一一映射。如调用两次RANDOM(0,1),可以获得集合{00,01,10,11}对应{0,1,2,3}。
这样一来,调用k次RANDOM(0,1)可以获得的二进制集合对应到集合{0,1,2,3.....2^k-1}.
a. 基于以上分析,我们可以令2^k-1 >= b-a+1, 从而求得一个最小的整数k, 这样获得的集合将能够包含所有要输出的随机数值。
invoke RANDOM(0,1) k times ---> Binary string ----> convert to int x (consider as unsigned Binary)
if x < b-a+1 Return a+x;
else go back and redo;
这样可以实现等概率的获得a~b中的随机数, 弊端是,当k过大时会造成很大一部分的CPU时钟浪费,效率低下。令A={0,1,2...b-a+1},当U-A的元素数量很多的
时候,程序很大一部分时间在重复获得2进制串。并且获得的随机数虽是等概率的,但概率却不等于1/(b-a+1).
b. 为了改进上面的算法,可以将获得的集合{0,1,2....2^k}均分为b-a+1个区间, 生成的二进制串对应的十进制数值落在那个区间即输出区间对应的随机数。
这样就可以得到等概率且概率为1/(b+a-1)的随机数了。