快速排序

当元素较少时，用插入排序速度更快（N<10时有这一结论）。而且快速排序首先需要分割，所以至少要有4个元素。

#include "stdafx.h"

#define CutOff (3) //分割参数，至少要有3个元素以供分割

void InsertionSort(int A[], int N)
{
	int i = 0;
	int j = 0;
	int tmp;
	for (i=1; i<N; i++)
	{
		tmp = A[i];
		for (j=i; j>0 && A[j-1]>tmp; j--)
		{
			A[j] = A[j-1];
		}
		A[j] = tmp;
	}
}

void Swap(int *x,int *y)
{
	int tmp;
	tmp = *x;
	*x = *y;
	*y = tmp;
}

//三数中值分割法
int Median3(int A[], int Left, int Right)
{
	int center = (Left + Right) / 2;
	int B[3] = {A[Left],A[center],A[Right]};
	InsertionSort(B,3);
	A[Left] = B[0];
	A[center] = B[1];
	A[Right] = B[2];

	//排序后Center作为枢纽元，有Left<Center，可以为j游标提供警戒标记（j从右向左探寻）
	//将中间值放入倒数第二位置，可以为i游标提供警戒标记（i从左向右探寻）
		Swap(&A[center],&A[Right-1]);

	return A[Right-1];
}

void Qsort(int A[], int Left, int Right)
{
	int i = 0;
	int j = 0;
	int Pivot = 0;

	if (Left + CutOff > Right)
		InsertionSort(A + Left,Right - Left + 1);
	else
	{
		Pivot = Median3(A,Left,Right);
		i = Left;
		j = Right - 1;
		for ( ; ; )
		{
			while (A[++i] < Pivot){}//如果i找到比枢纽元大的值，则停止
			while (A[--j] > Pivot){}//如果j找到比枢纽元小的值，则停止
			if (i < j) //如果还需要继续探寻，则交换两元素使其位于正确的分组中
				Swap(&A[i],&A[j]);
			else
				break;//一旦i和j交叠，则一趟分组结束
		}
		Swap(&A[i],&A[Right - 1]);//循环结束后，i停止在比枢纽元大的元素上

		Qsort(A,Left,i-1);//注意位置i上是枢纽元，已经在它正确的位置上了，不用动
		Qsort(A,i+1,Right);
	}
}

void QuickSort(int A[], int N)
{
	Qsort(A,0,N-1);
}

int _tmain(int argc, _TCHAR* argv[])
{
	int A[10] = {3,5,9,7,1,4,2,0,8,6};
	QuickSort(A,10);
	return 0;
}

如果不使用三数中值分割法，可以得到比较简洁版本的快速排序算法，程序如下：

#include "stdafx.h"

void Swap(int *x,int *y)
{
	int tmp;
	tmp = *x;
	*x = *y;
	*y = tmp;
}

void Qsort (int A[], int Left, int Right)
{
	int pivot = A[(Left+Right)/2];
	int i = Left;
	int j = Right;
	do 
	{
		while (A[i] < pivot)
			++i;
		while (A[j] > pivot)
			--j;
		if (i <= j)
			Swap(&A[i++],&A[j--]);
	} while (i <= j);

	if (j > Left)
		Qsort(A,Left,j);
	if (i < Right)
		Qsort(A,i,Right);
}

void QuickSort (int A[], int N)
{
	Qsort(A,0,N-1);
}

int _tmain(int argc, _TCHAR* argv[])
{
	int A[10] = {3,5,9,7,1,4,2,0,8,6};
	QuickSort(A,10);
	return 0;
}