快速中值求取算法

中值，顾名思义，就是指一个从小到大的序列的中间的那一个数。一般的讲，中值比平均值还要更加稳定。如一个序列中的某一个值被误乘以了100，平均值则会有很大的波动，但是中位数则不会发生太大的变化。

但是如果对数据先排序，然后再进行取中值，则比较耗时。以下是一种快速提取中位数的算法，算法描述如下：
1：取队首，队尾和中间的三个数，通过交换数据，使得队尾最大，中间的数据最小，队首的数据为哨兵。
2：交换中间和第2个数据，通过变换数据，使存在一个位置A，在该位置前的数据都小于哨兵，在该位置后的数据都大于或等于哨兵。
3：如果A的位置在队中之后，则更新队列为1~A，否则为A~end。并继续调用这个过程。

举例如下：
数组：3，6，2，1，9，8，0，4，5，7
经过步骤1并交换中间和第2个数的值之后，则为：
7，3，2，1，6，8，0，4，5，9；此序列中，7为哨兵
经过第一次迭代后，序列为：
4，3，2，1，6，5，0，7，8，9，此序列中，在数字7的位置是正确的。
第二次迭代则取1~7个数进行求解中位数，其中，哨兵为1。迭代之后如下：
0，1，2，3，6，5，4，7，8，9，同样的，此序列中，数字1的位置也是正确的。
第三次迭代则以第3~7个数进行求解，哨兵为5，此时，5即为中位数。

以下为非递归式调用的C代码。

int CQuickMedian(int * pnData, const int knLength){    int nLow = 0;    int nHigh = 0;    int nMiddle = 0;    int nMedian = 0;    int nLTmp = 0;    int nHTmp = 0;    nMedian = (knLength - 1) >> 1;    nHigh = knLength - 1;    while(1)    {        if(nHigh == nLow)        {            return pnData[nHigh];        }        if(nHigh == nLow + 1)        {            return pnData[nHigh] > pnData[nLow] ? pnData[nLow] : pnData[nHigh];        }        nMiddle = (nHigh + nLow)>>1;        if(pnData[nLow] > pnData[nHigh])        {            CSwap(&pnData[nHigh], &pnData[nLow]);        }        if(pnData[nMiddle] > pnData[nHigh])        {            CSwap(&pnData[nMiddle], &pnData[nHigh]);        }        if(pnData[nMiddle] > pnData[nLow])        {            CSwap(&pnData[nMiddle], &pnData[nLow]);        }        CSwap(&pnData[nMiddle], &pnData[nLow+1]);        nLTmp = nLow + 2;        nHTmp = nHigh - 1;        while(1)        {            while(pnData[nLTmp] <= pnData[nLow])            {                nLTmp ++;            }            while(pnData[nHTmp] >= pnData[nLow])            {                nHTmp --;            }            if(nLTmp > nHTmp)            {                CSwap(&pnData[nHTmp], &pnData[nLow]);                if(nHTmp > nMedian)                {                    nHigh = nHTmp - 1;                }                else                {                    nLow = nLTmp - 1;                }                break;            }            else            {                CSwap(&pnData[nLTmp], &pnData[nHTmp]);                nLTmp++;                nHTmp--;            }                   }           }}