几道有趣的概率题

这几天在刷题，遇到几道有趣的概率题，需要设计算法，编程实现，整理一下：

给定N张扑克牌和一个随机函数，设计一个洗牌算法

void fun(vector<int>& v)                                                                                                                                                                
{      
     srand(time(0));
     for (int n = v.size(); n > 0; n--)
     {  
         int index = rand() % n;
         swap(v[n-1], v[index]);
     }  
 }  

---

如何等概率地从n个数中随机抽出m个数？

//假设这n个数的序号依次为0,1,2,...,n-1,数组名为num
void knuth1(int* pNum, int m, int n)
{
    srand((unsigned int)time(0));
    for (int i=0; i<n; i++)
    {
        if (rand()%(n-i) < m) //rand()%(n-i)的取值范围是[0, n-i）
        {
            cout << pNum[i] << endl;
            m--;
        }
    }
}

或者维护一个set，每次取一个元素删除下次继续取，也是可以的，但是有一点不好，当set中元素很少时候，rand（）函数可能要很久才产生这个数，效率不高。

如果n未知，可参考 从未知大小的n个数中取m个数，使各数被取出的概率相等

---

给定一个能够生成0,1两个数的等概率随机数生成器”，如何生成⼀个产生0,1,2,3的等概率随机数生成器?
和上题类似，如何用rand7生成rand9?

int rand0_3()
{
    return rand0_1()+rand0_1();
}

int rand9()
{
    int n=0;
    do{
        n=(rand()7-1)*7+rand()7;
    }while(n>45)
    return n%9+1;
}
//rand7 --》 rand5
int Rand5(){
    int x = ~(1<<31); // max int
    while(x > 5)
        x = Rand7();
    return x;
}

---

有一枚硬币,以p的概率产生正面，以1-p的概率产生背面，如何利用它产生个0.5概率的生成器？

制作 1 2 发生概率都是 1 / 2 的发生器，连续发生2次，则发生00,11的概率为p*p,(1-p)(1-p),发生10,01的概率都为p(1-p),在发生10时返回1，发生01时返回2，则发生1，2的概率相等

制作 1 2 3 发生概率都是 1 / 3的发生器，连续发生3次，则发生001,010,100的概率都为p*p*(1-P),或者是110,101,011概率都为p*(1-p)*(1-p),则用001,010,100分别对应1,2,3返回，即可使得发生1,2,3的概率都为1/3