随机化算法实验(Sherwood型线性时间选择)精要.doc

算法分析与设计实验报告 第 八 次实验 姓名 学号 班级 时间 12.26上午 地点 工训楼309 实验名称 随机化算法实验(Sherwood型线性时间选择) 实验目的 通过上机实验，要求掌握 Sherwood 型线性时间选择算法的问题描述、算法设计思想、程序设计。 使用舍伍德型选择算法， 根据不同的输入用例， 能准确的输出用例中的中值， 并计算出程序运行所需要的时间。 实验原理 问题描述： 给定任意几组数据，利用舍伍德型选择算法，找出数组中的中值并输出（若 数组为奇数个则输出中值，若数组为偶数个则输出第 n/2 小元素）。 基本思想： ?设A是一个确定性算法，当它的输入实例为x时所需的计算时间记为tA(x)。设Xn是算法A的输入规模为n的实例的全体，则当问题的输入规模为n时，算法A所需的平均时间为。这显然不能排除存在xXn使得的可能性。希望获得一个随机化算法B，使得对问题的输入规模为n的每一个实例均有。这就是舍伍德算法设计的基本思想。当s(n)与tA(n)相比可忽略时，舍伍德算法可获得很好的平均性能。 实验步骤 （1）先判断是否需要进行随机划分即（k∈(1,n)?n1?）； （2）产生随机数j,选择划分基准，将a[j]与a[l]交换； （3）以划分基准为轴做元素交换，使得一侧数组小于基准值，另一侧数组值大于基准值； （4）判断基准值是否就是所需选择的数，若是，则输出；若不是对子数组重复步骤(2)(3)(4)。 关键代码 //计算a[l:r]中第k小元素 templateclass Type Type select(Type a[],int l,int r,int k) { while(true) { if(l=r) //此种情况表示只有一个元素 { return a[l]; } int i = l; int j = r+1; Type pivot = a[l]; //轴值每次都选择最左边的元素 //以划基准为轴做元素交换 while(true) { //将pivot的元素换到轴值左边 //将pivot的元素换到轴值右边 while(a[++i]pivot); while(a[--j]pivot); if(i=j) { break; } Swap(a[i],a[j]); } if(j-l+1 == k)//如果轴值是要选的元素即第k小 { return pivot; } //如果轴值不是要选的元素 //a[j]必然小于等于pivot,做最后一次交换，满足左侧比pivot小，右侧比pivot大 a[l] = a[j]; a[j] = pivot; //对子组重复划分过程 if(j-l+1k) //第k小在pivot的右边部分 { k = k-j+l-1;//右侧：k-(j-l+1)=k-j+l-1 l = j + 1; } else //第k小在pvot的左边部分 { r = j - 1; } } } 输入较小的序列数： 输入较大的序列数： 不输出数组数只输出结果比较： 实验分析 在Sherwood型线性时间选择算法中，有两种实现方法，一种是随机划分基准；一种是随机洗牌预处理。 随机划分基准处理如下： 随机洗牌预处理如下： if(l=r) templateclass Type { return a[l]; } void Shuffle(Type a[],int n) int i = l, { //随机数类 static