问题描述一
- 假设你希望以 1 / 2 1/2 1/2的概率输出0与1。你可以自由使用一个输出0或1的过程BIASED-RANDOM。它以某概率 p p p 输出1,概率 1 − p 1-p 1−p 输出0,其中 0 < p < 1 0<p<1 0<p<1,但是 p p p 的值未知。请给出一个利用BIASED-RANDOM作为子程序的算法,返回一个无偏的结果,能以概率 1 / 2 1/2 1/2返回0,以概率 1 / 2 1/2 1/2返回1。
问题求解一
- 我们通过调用两次BIASED-RANDOM,对其结果进行分析。
- 我们记第一次BIASED-RANDOM的结果为x,第二次BIASED-RANDOM的结果为y。
x , y ∈ { 0 , 1 } x,y\in \left\{0,1\right\} x,y∈{0,1}
那么,这两个随机事件的结果的概率分布如下表
x=0 | x=1 | |
---|---|---|
y=0 | ( 1 − p ) 2 (1-p)^2 (1−p)2 | p ( 1 − p ) p(1-p) p(1−p) |
y=1 | p ( 1 − p ) p(1-p) p(1−p) | p 2 p^2 p2 |
我们注意到,{ x=1 and y=0 } 与 { x=0 and y=1 }这两个事件发生的概率相等,故可以以这两个事件为基础,构造等概率生成0,1的算法,即可满足要求。
算法以伪码的形式给出:
UNBIASED-RANDOM
while true
x=BIASED-RANDOM
y=BIASED-RANDOM
if x!=y
return x
问题描述二
- 调用RANDOM ( a , b ) (a,b) (a,b)将返回一个介于a和b之间的整数,并且每个整数以等概率出现。例如,RANDOM ( 0 , 1 ) (0,1) (0,1)产生0的概率是1/2,产生1的概率也是1/2。调用RANDOM ( 3 , 7 ) (3,7) (3,7)将返回3、4、5、6或7,每个出现的概率都是1/5。每次RANDOM返回的整数独立于前面调用的返回值。
- 现在问题为:请描述RANDOM ( a , b ) (a,b) (a,b)过程的一种实现,它只调用RANDOM ( 0 , 1 ) (0,1) (0,1)。
问题求解二
Step1.计算数(b-a)的二进制位数n
Step2.调用RANDOM(0,1)n次,并将n次的结果以二进制数的形式进行保存为m
Step3.将二进制数m转换为十进制数m‘
Step4.计算a+m’,若a+m’ ∈ [ a , b ] \in [a,b] ∈[a,b],则返回a+m’;若a+m’ ∉ [ a , b ] \notin [a,b] ∈/[a,b]则回到Step2,重新生成。
问题来源
《算法导论》(原书第3版) —— Thomas H.Cormen / Charles E.Leiserson / Ronald L.Rivest / Clifford Stein