CLRS 5.1雇佣问题

5.1-1
由于我们可以比较任意两个应聘者，这说明我们可以有一个应聘者的优劣顺序，因此我们知道应聘者的全部次序。

5.1-2
设 $n=b-a$，我们需要找到最少的 $c$ 位二进制数，使得 $2^c\ge n$，即 $c=\lfloor \lg n\rfloor+1$，则：
1) 调用 $c$ 次RANDOM(0,1)，产生 $c$ 位二进制数 $r$；
2) 若 $r$ 大于 $n$，重新 1)；
3) 返回 $a+r$。

生成某个随机数的概率为$(\frac12)^c$，若此数大于 $n$ 则需要重复第 2 步，设第一次生成的数小于等于 $n$ 的概率为 $p=\frac{n}{2^c}$ ，否则需要第二次生成 $r$ 。这是一个几何分布，我们要 $1/p$ 的概率生成不大于 $n$ 的数，即 $\frac{2^c}{n}$。我们调用 $c$ 次RANDOM(0,1)，因此有：

$$O\bigg(\frac{2^c}{n} c\bigg) = O\bigg(\frac{\lg(b-a)2^{\lg(b-a)}}{b-a}\bigg)= O(\lg(b-a))$$

#include <iostream>
#include <random>
#include <ctime>
using namespace std;

int random_01()
{
    static default_random_engine e(time(0));
    static bernoulli_distribution b;//伯努利分布
    if (b(e))
        return 1;
    else return 0;
}

int random_ab(int a, int b)
{
    int n = b - a;
    int c = log(static_cast<double>(n)) / log(static_cast<double>(2)) + 1;
    while (true)
    {
        int res = 0;
        for (int i = 0; i != c; ++i)
            res += random_01() << i;
        if (res + a <= b)
            return res + a;
    }
}

int main()
{
    for (int i = 0; i < 10; ++i)
        cout << random_ab(0, 100) << ' ';
    cout << endl;
    return 0;
}

5.1-3
使用两次BIASED-RANDOM()产生两个随机数 $x$，$y$（$x,y\in (0,1)$），若 $x=y$ 则返回 $x$，否则重新产生 $x,y$。
出现两个1的概率为 $p^2$，出现两个0的概率为 $(1-p)^2$，出现一个1和一个0的概率是 $2p(1-p)$。

$$\eqalign { E(T) &= \sum_{i=1}^{\infty} {2p(1-p)(1-2p(1-p))^{i-1})} \\ & = O(\frac {1}{2p(1-p)})}$$