第10章概率算法全解.pptVIP

考虑快速排序中的划分过程，选定一个轴值将序列ri～rj进行划分，使得比轴值小的元素都位于轴值的左侧，比轴值大的元素都位于轴值的右侧，假定轴值的最终位置是s，则： （1）若k=s，则rs就是第k小元素； （2）若ks，则第k小元素一定在序列ri～rs-1中； （3）若ks，则第k小元素一定在序列rs+1～rj中； 整数因子分解问题可以归结为整数因子划分和素数测试：假设对整数n进行因子分解，首先进行素数测试，如果n是素数，则分解完成；否则，再进行因子划分，找到n的一个非平凡因子m，并递归地对m和n/m进行因子分解。 下面讨论整数因子划分问题。 bool Prime2(int n) { composite=false; a=Random(2, n-1); //产生[2, n-1]的随机数 result=Power(a, n-1, n); if (composite | | (result!=1)) return false; else return true; } 算法Factor是对范围在1～ 的所有整数进行了试除而得到了n的最小素数因子，其时间复杂性是O( )。对于一个正整数n，其位数为 ，则算法Factor的时间复杂性是O(10m/2)。假定每次循环只需要1纳秒，它也需要花费1 000年的时间来分解一个40位左右的坚固的（Hard）合数。 “坚固的”合数是指这个数是两个规模相当的素数的乘积。 到目前为止，还没有找到求解整数因子划分问题的多项式时间算法。 求解整数因子划分问题的拉斯维加斯型概率算法在开始时选取0～n-1范围内的随机数x1，然后递归地由下式产生无穷序列x1, x2, …, xk, …。 对于i=2k（k=0, 1, 2, …），以及2k＜j≤2k+1，计算xj-xi与n的最大公因子d，如果d大于1，则d即是n的非平凡因子，算法实现对n的一次划分。 计算xj-xi与n的最大公因子可以采用欧几里德算法，具体算法如下： 算法10.6——最大公因子 int Gcd(int a, int b) { r=a % b while (r!=0) { a=b; b=r; r=a % b; } return b; } 算法10.7——整数因子划分 int Pollard(int n) { i=1; k=2; x=Random(0, n-1); // x为[0, n-1]区间的随机整数 y=x; while (true) { i++; x=(x*x-1) % n; d=Gcd(y-x, n); if (d1) return d; //若y-x与n存在最大公约数d，则d即为n的非平凡因子 if (i= =k) { y=x; k*=2; } } } 算法Pollard中的while循环执行约 次后，会得到n的一个素数因子p。由于n的最小素数因子 ，故算法Pollard可在O(n1/4)的时间内找到n的一个素数因子。 10.4 蒙特卡罗(Monte Carlo)型概率算法 10.4.1 主元素问题 10.4.2 素数测试问题 10.4.1 主元素问题 设T[n]是一个含有n个元素的数组，x是数组T的一个元素，如果数组中有一半以上的元素与x相同，则称元素x是数组T的主元素（Major Element）。 例如，在数组T[7]={3, 2, 3, 2, 3, 3, 5}中，元素3就是主元素。 蒙特卡罗型概率算法求解主元素问题可以随机地选择数组中的一个元素T[i]进行统计，如果该元素出现的次数大于n/2，则该元素就是数组的主元素，算法返回true；否则随机选择的这个元素T[i]不是主元素，算法返回false。此时，数组中可能有主元素也可能没有主元素。如果数组中存在主元素，则非主元素