第九章概率算法 (2).pptVIP

计算机算法设计与分析计算机算法设计与分析计算机算法设计与分析第九章概率算法*计算机算法设计与分析*第1页，共38页，星期日，2025年，2月5日概率计算概率计算就是在算法中可采用随机选择计算的步骤、元素或参数等。它的基本特征是计算具有不确定性。它的解也不一定是最优解。它在很大程度上能降低算法的复杂度。在非标准算法中普遍了应用概率方法，主要有：(1)直接用概率进行数值计算；(2)用概率/随机进行选择；(3)利用概率加速搜索或避免陷于局部最优。第2页，共38页，星期日，2025年，2月5日直接用概率进行数值计算设f(x)是[0,1]上的连续函数，求I=∫f(x)dx。01y=f(x)G假设向单位正方形内随机投入n个点(xi,yi)，若有m个点落入G中，则I≈m/n。doubleDarts(intn){doublex,y;intk=0;staticRandomNumberdart;for(inti=1;i=n;i++){x=dart.fRandom();y=dart.fRandom();if(y=f(x))k++;}returnk/double(n);}第3页，共38页，星期日，2025年，2月5日划分基准的随机选择在快速排序算法中，若用拟中位数作为划分标准，可保证在线性时间内完成。但是确定拟中位数要付出额外开销。若选用第一个元素为划分基准，最坏时的时间复杂性为O(n2)。若在算法中采用随机选择一个元素作为划分标准，便可既保证算法的线性时间平均性能，又避免了计算拟中位数的麻烦。也可先对数组进行“洗牌”，然后再进行确定的排序算法。这样依然可取得同样的效果。第4页，共38页，星期日，2025年，2月5日“洗牌”后的快速排序voidShuffle(Typea[],intn){//随机洗牌算法staticRandomNumbermd;for(inti=1;in;i++){intj=md.Random(n–i+1)+i;Swap(a[i],a[j]);}}VoidQuiksortByShuffle(Typea[],intn){Shuffle(a,n);//将数组a洗牌Quiksort(a,n);}第5页，共38页，星期日，2025年，2月5日随机抽样在n个元素的集合中随机抽取m(0m≤n)个无重复的元素。为简单起见，假定所有元素的值都位于1至n之间。第6页，共38页，星期日，2025年，2月5日随机抽样我们采用下面的方法进行选择：1、首先将n个元素都标记为“未选择”；2、重复下列步骤直到抽取了m个不同的元素⑴产生一个1至n间的随机数r；⑵如果r标记为“未选择”，将它标记为“已选择”，并加入到抽样中。第7页，共38页，星期日，2025年，2月5日随机抽样intRandomSampling(S[n],A[m],m){mark[1..n]=False;count=0;while(countm){r=random(1,n);if(mark[r]==False){count++;A[count]=S[r];mark[r]=True;}}}第8页，共38页，星期日，2025年，2月5日判定素数的概率算法判定自然数是否是素数，不仅有重要理论意义，而且在密码学中具有重要实用价值。最简单的素数判定方法是依次测定从2到n?中是否存在n的因子，该算法的复杂度为O(n?)。筛法：将小于n的合数预先筛掉，而不用判断其是否为n的因子。它虽然没有降低算法的复杂度，但实际运行速度比前者要快得多。概率算法，保证一定概率的前提下简单判断。第9页，共38页，星期日，2025年，2月5日Fermat素数测试法Fermat定理：若p是素数，则对任意整数a，gcd(a,p)=1，则有ap–1≡1(modp)。显然，对素数p有pp–1≡1(modp)。对于一般的整数n，满足nn–1≡1(modn)的数目很少。满足的称为伪素数。就用是否满足nn–1≡1(modn)来判断n是否为素数。第10页，共38页，星期日，2025年，2月5日Fermat素数测试法BoolFermat_Prime(intn){a=2;power