算法导论第5章 概率分析和随机算法（5.1）

看看概率和你的代码和数学有什么关系，换个（概率的）角度来看看代码，也许有所收获。擅长百度、谷歌代码完成功能的码农看过来。

第5章  概率分析和随机算法

5.1  雇佣问题

概率分析：在问题的分析中应用概率技术。

随机算法：如果一个算法的行为，不只由输入决定，也由随机数生成器决定，那么称这个算法为随机算法。

1. 习题解答

5.1-1

如果可以知道任意两个人的排名的大小，自然等价于得到了总排名。

5.1-2

解：这题等价于二进制表示任意数，其中一次random相当于一位二进制，其取值为0或1.

例如a=3,b=5;共3个数：3,4,5。等价于0,1,2.只需要2个random即可。00,01,10分别表示3,4,5。

分析其期望运行时间，理论上，N位数表示的是偶数个数，所以有可能多一个数。故一次生成的概率为X=(b-a)/(b-a+1);根据引理5.2，EX=(b-a)/(b-a+1)（次）；

再分析运行一次的时间复杂度。位数为w=（[lg(high-low+1)]的向上取整），故需要循环w次依次取随机数。

综上所述：期望运行时间为EX*w=

附C语言下实现如下：

//统计上概率地返回在low和high之间的数

#define random(x) (rand()%x)

int t5_1_2(int low,int high)

{

    int i=0;

    int count=0;

    double j=log(high-low+1)/log(2);

    if(j-(int)j>0) {

        j++;

    }

    while(1) {

        count=0;

        for(i=0; i<j; i++) {

            count+=random(2)*pow(2,i);

        }

        if(low+count<=high) {

            break;

        }

    }

    return low+count;

}

测试代码：

#include <stdio.h>

#include <stdlib.h>

#include <math.h>

#define random(x) (rand()%x)

int main()

{

    int j=0;

    for(; j<11; j++) {

        int a1=t5_1_2(1,4);

        int a2=t5_1_2(11,14);

        printf("between %d and %d : %d\n",a1,a2,t5_1_2(a1,a2));

    }

    return 0;

}

 

PS:是不是很神奇~~~~。下一题更有趣。

5.1-3

解：根据上一题的启发。一位表示一次random。

根据离散型随机变量的相关知识，我们可以猜出来：

只要找到两个等概率的取值就好了。例如(0,1)和(1,0)。他们的概率分别是p(1-p)和(1-p)p。我们如上题一样舍去(0,0)和(1,1)即可。

取值

00

01

10

11

P

(1-p)(1-p)

(1-p)p

p(1-p)

pp

为了复习一下概率，我把它当做概率题来做好了。

设X表示两位random的值。

设

Y=1 当X==01或10；

Y=0 当X==11 或00

故

取值

1

0

P

2p(1-p)

Pp+(1-p)(1-p)

计算其期望：

EX=2p(1-p).

再分析其单次运行时间为O(1)；

所以其运行时间期望为EX=2p(1-p).

附C语言实现代码：

//有偏输出0或1

int biasedRandom()

{

    int i=random(3);

    if(i>=1) return 1;

    else return 0;

}

int t5_1_3( )

{

    int j1=0,j2=0;

    while (1) {

        j1 =biasedRandom();

        j2 =biasedRandom();

        if(j1+j2!=1) {

            continue;

        }

        if(j1==1) return 1;

        else return 0;

    }

}

测试代码：

    int j=0;

    int i0=0;

    int i1=0;

    for(; j<500; j++) {

        if( t5_1_3( )==0) i0++;

        else i1++;

    }

    printf("i0=%d,i1=%d",i0,i1);

    return 0;

欢迎加入QQ群181539446