中位数不正常的中位数

时间:2015-10-08 16:29:45

标签: c median-of-medians

我正在编写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);
}

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);