我正在编写Median of Medians算法的C代码,以便在最坏情况线性时间内找到第k个最小元素。 我已经检查了我的代码以进行快速排序,交换等。 一切看起来都不错,但每次都不能正常工作。
给出的输入 -
n=12 kth=7
A[]=53 22 65 18 89 45 42 63 99 11 36 55
输出 -
Smallest at k=7 is 89
但输出需要
Smallest at k=8 is 53
函数调用 -
med_of_medians(A,0,n-1,kth);
代码 -
int med_of_medians(int A[], int a, int b, int kth)
{
if(a==b)
return 0;
int n=(b-a+1),median,pos,rank;
int i,med[(n+4)/5];
for(i=0;i<n/5;i++)
med[i]=find_median(A,(i*5)+a,((i+1)*5)-1);
if(n%5>0)
med[i++]=find_median(A,(n/5)*5+a,b);
median=(i==1)?med[0]:find_median(med,0,i-1);
pos=partition(A,a,b,median);
rank=pos-a+1;
if(rank==kth)
return A[pos];
else if(rank>kth)
return med_of_medians(A,a,pos-1,kth);
else
return med_of_medians(A,pos+1,b,kth-pos-1);
}
答案 0 :(得分:0)
您可能希望在第一个A[a]
中返回if
:
if (a == b)
return A[a];
a
电话中的第二个索引中缺少find_median
:
med[i]=find_median(A,(i*5)+a, a +((i+1)*5)-1);
并且还应添加a
以进行新的排名计算(应该是kth - rank
):
return med_of_medians(A,pos+1,b,kth-pos-1 + a);