查找第k小（中位数）

中位数就是一列数按从小到大排序后中间的那个数，如有n个数，如果n是奇数那么中位数就是n div 2+1的那个数，偶数的话有下中位数和上中位数分别为第n div 2和n div 2+1。查找中位数就是查找一列数中第n div 2或n div 2+1小的数。

    查找第k小的数最直观的方法就是先排序再输出数组中对应下标的值，这样做的时间复杂度取决于排序的时间复杂度，如果是快速排序的话是O(nlogn)，其 实有更快捷的算法，形式上像是快速排序，实质上是把快速排序的过程简化了一下，去掉与此无关的操作就得到了该算法，因此，该算法时间复杂度同为 O(nlogn)但省去不少步骤，因而比快速排序的算法要快。为方便介绍先简单复习一下快速排序，快速排序是一种二分的算法首先选定一个数作为基准，再把 数组按该基准分成左右两块，分别是小于和大于该基准的所有数，再递归对这两块分别进行上述划分过程直到区间被分成单个的数。在这里不是每次都要递归划分， 如果小于基准的那块包含g-1个数那么该基准就是第g小，个数可以由前后座标的差加上1得到，如果g=k这时已经找到答案可以退出程序了。如果g> k，那么在小于基准那一块中的第k小的数就是全部的第k小，如果g<k，那么就在大于基准的那块中查找第k-g小，这样的查找已经缩小了数据规模， 也就快捷的多了。基准的选取可以有很多方式，一般选取最后一个或第一个，我的程序中选最后一个。但如果是已经有序的数列，那么这种基准的划分结果就很不平 衡，因此采用随机挑选基准的方法就会好很多。下面程序中rs函数就是随机选定基准进行划分。divided就是划分的过程。

   在divided中有两个指针i和j，所有下标小于i的数组元素均小于基准，在i与j之间的元素均大于基准，在程序结束时i的位置就是基准在排序后所在的位置。

下面是pascal的程序，是我从c语言翻译过来的，偶然发现在lazarus中使用c语言的行注释方法并没有报错，于是没有改过来

program k_min; type d=array[0..10000]of integer; var n, i, k:integer; a:d; function divide(var a:d;p,r:integer): integer; {把当前区间的数组以最后一个元素分为两块，左边的部分小于该元素小于右边的元素 } var x, i, j, t:integer; begin x:= a[r]; {取最后一个元素} i:= p; {划分点，程序运行到末尾时，当下标小于i时元素小于x，大于i时元素大于x } if p < r then begin {扫描，当a[j]<x时放到i的位置，并将i向后移动一位， 这样将使j前所有小于x的全部放到i的左侧 } for j:= p to r-1 do begin if a[j] < x then begin if i <> j then begin t:= a[i]; a[i]:= a[j]; a[j]:= t; end; inc(i); end; end; if i <> r then begin //将x放到i的位置。 t := a[i]; a[i] := a[r]; a[r] := t; end; end; divide:=i; //返回i，此时以a[i]为标准，分成大于a[i]和小于a[i]; end; function rs(var a:d; p,r:integer):integer; //random selection var x,t:integer; begin x:= random(r - p)+p; //查找一个随机位置，并与最后一个元素掉换 t:= a[r]; a[r]:= a[x]; a[x]:= t; rs:= divide(a, p, r); //返回划分结果是第几个元素 end; function kmin(var a:d; p,r,k:integer): integer; //第k小的元素 var q, g:integer; begin if p < r then begin q:= rs(a, p, r); //随机选取一个数划分 g:= q - p + 1; //q是该区间中第g个元素 if g = k then kmin:= a[q] //由于a[q]前面的元素全部小于a[q]，若q是第k个元素就是第k小 else if g > k then kmin:= kmin(a, p, q - 1, k) //若k在g的前面则只需在a[q]前面找第k小 else kmin:= kmin(a, q + 1, r, k - g); //否则在a[q]后找k－g小 end else kmin:= a[p]; end; begin randomize; write( '输入有几个数:'); readln(n); writeln('输入每个数'); for i:= 0 to n-1 do read( a[i]); writeln( '您要查找第几个数'); read(k); writeln( kmin(a, 0, n - 1, k)); read(k); end.