搜索
1、舍伍德(Sherwood)算法设A是一个确定性算法，当它的输入实例为x时所需的计算时间记为tA(x)。设Xn是算法A的输入规模为n的实例的全体，则当问题的输入规模为n时，算法A所需的平均时间为。这显然不能排除存在x∈Xn使得的可能性。希望获得一个随机化算法B，使得对问题的输入规模为n的每一个实例均有。这就是舍伍德算法设计的基本思想。当s(n)与tA(n)相比可忽略时，舍伍德算法可获得很好的平均性能。2、线性时间选择算法1)随机划分选择基准对于选择问题而言,用拟中位数作为划分基准可以保证在最坏的情况下用线性时间完成选择。如果只简单地用待划分数组的第一个元素作为划分基准，则算法的平均性能较好，而在最坏的情况下需要O(n^2)计算时间。舍伍德选择算法则随机地选择一个数组元素作为划分基准，这样既保证算法的线性时间平均性能，又避免了计算拟中位数的麻烦。非递归的舍伍德型选择算法如下：[cpp]viewplaincopy1.//随机化算法线性时间选择随机划分选择基准2.#includestdafx.h3.#includeRandomNumber.h4.#includeiostream5.usingnamespacestd;6.7.templateclassType8.Typeselect(Typea[],intl,intr,intk);9.10.templateclassType11.Typeselect(Typea[],intn,intk);12.13.templateclassType14.inlinevoidSwap(Type&a,Type&b);15.16.intmain()17.{18.inta[]={5,7,3,4,8,6,9,1,2};19.cout原数组为：endl;20.for(inti=0;i9;i++)21.{22.couta[i];23.}24.coutendl;25.cout所给数组第7小元素为：select(a,9,7)endl;26.return0;27.}28.29.//计算a[0:n-1]中第k小元素30.//假设a[n]是一个键值无穷大的元素31.templateclassType32.Typeselect(Typea[],intn,intk)33.{34.if(k1||kn)35.{36.cout请输入正确的k!endl;37.return0;38.}39.returnselect(a,0,n-1,k);40.}41.42.//计算a[l:r]中第k小元素43.templateclassType44.Typeselect(Typea[],intl,intr,intk)45.{46.staticRandomNumberrnd;47.while(true)48.{49.if(l=r)50.{51.returna[l];52.}53.54.inti=l,55.j=l+rnd.Random(r-l+1);//随机选择划分基准56.57.Swap(a[i],a[j]);58.59.j=r+1;60.Typepivot=a[l];61.62.//以划分基准为轴做元素交换63.while(true)64.{65.while(a[++i]pivot);66.while(a[--j]pivot);67.if(i=j)68.{69.break;70.}71.Swap(a[i],a[j]);72.}73.74.if(j-l+1==k)//第k小75.{76.returnpivot;77.}78.79.//a[j]必然小于pivot,做最后一次交换，满足左侧比pivot小，右侧比pivot大80.a[l]=a[j];81.a[j]=pivot;82.83.//对子数组重复划分过程84.if(j-l+1k)85.{86.k=k-j+l-1;//右侧：k-(j-l+1)=k-j+l-187.l=j+1;88.}89.else90.{91.r=j-1;92.}93.}94.}95.96.templateclassType97.inlinevoidSwap(Type&a,Type&b)98.{99.Typetemp=a;100.a=b;101.b=temp;102.}程序运行结果如图：2)随机洗牌预处理有时也会遇到这样的情况，即所给的确定性算法无法直接改造成舍伍德型算法。此时可借助于随机预处理技术，不改变原有的确定性算法，仅对其输入进行随机洗牌，同样可收到舍伍德算法的效果。例如，对于确定性选择算法，可以用下面的洗牌算法shuffle将数组a中元素随机排列，然后用确定性选择算法求解。这样做所收到的效果与舍伍德型算法的效果是一样的。[cpp]viewplaincopy1.//随机化算法线性时间选择输入预处理，洗牌2.#includestdafx.h3.#includeRandomNumber.h4.#includeiostream5.usingnamespacestd;6.7.templateclassType8.Typeselect(Typea[],intl,intr,intk);9.10.templateclassType11.Typeselect(Typea[],intn,intk);12.13.templateclassType14.voidShuffle(Typea[],intn);15.16.templateclassType17.inlinevoidSwap(Type&a,Type&b);18.19.intmain()20.{21.inta[]={5,7,3,4,8,6,9,1,2};22.cout原数组为：endl;23.for(inti=0;i9;i++)24.{25.couta[i];26.}27.coutendl;28.Shuffle(a,9);//洗牌29.cout洗牌后数组为：endl;30.for(inti=0;i9;i++)31.{32.couta[i];33.}34.coutendl;35.cout所给数组第7小元素为：select(a,9,7)endl;36.return0;37.}38.39.//计算a[0:n-1]中第k小元素40.//假设a[n]是一个键值无穷大的元素41.templateclassType42.Typeselect(Typea[],intn,intk)43.{44.if(k1||kn)45.{46.cout请输入正确的k!endl;47.return0;48.}49.returnselect(a,0,n-1,k);50.}51.52.//计算a[l:r]中第k小元素53.templateclassType54.Typeselect(Typea[],intl,intr,intk)55.{56.while(true)57.{58.if(l=r)59.{60.returna[l];61.}62.inti=l;63.intj=r+1;64.Typepivot=a[l];65.66.//以划分基准为轴做元素交换67.while(true)68.{69.while(a[++i]pivot);70.while(a[--j]pivot);71.if(i=j)72.{73.break;74.}75.Swap(a[i],a[j]);76.}77.78.if(j-l+1==k)//第k小79.{80.returnpivot;81.}82.83.//a[j]必然小于pivot,做最后一次交换，满足左侧比pivot小，右侧比pivot大84.a[l]=a[j];85.a[j]=pivot;86.87.//对子数组重复划分过程88.if(j-l+1k)89.{90.k=k-j+l-1;//右侧：k-(j-l+1)=k-j+l-191.l=j+1;92.}93.else94.{95.r=j-1;96.}97.}98.}99.100.templateclassType101.inlinevoidSwap(Type&a,Type&b)102.{103.Typetemp=a;104.a=b;105.b=temp;106.}107.108.//随机洗牌算法109.templateclassType110.voidShuffle(Typea[],intn)111.{112.staticRandomNumberrnd;113.for(inti=0;in;i++)114.{115.intj=rnd.Random(n-i)+i;116.Swap(a[i],a[j]);117.}118.}程序运行结果如图：