一、算法思想
快速排序是Hoare于1962年提出的一种二叉树结构的交换排序方法,其基本思想是:任意选取待排序序列中的一个元素作为基准值,按照该排序码将待排序集合分割成两个子序列,左子序列中所有元素均小于基准值,右子序列中所有元素均大于基准值,然后再对左右子序列重复该过程,直到所有元素都排列在相应位置上为止。
注:后续实现均以取序列最右侧元素为基准值为例。
其中,将区间按照基准值划分为左右两部分的常见方式有三种:
1.hoare版
通过两个标志位,分别从左往右找比基准值大的元素和从右向左找比基准值小的元素,然后交换两个元素,重复上述过程直到两个标志位相遇,最后再将基准值交换到相遇位置,然后再对左右子序列重复上述过程,直到整个序列完成排序。
图示:
2.挖坑法
挖坑法大体思路与hoare法思路相同,其基本思想是:将基准值保存到标记位中,这样最右侧位置就形成了一个坑位,然后左侧标注位往右遍历找比基准值大的元素,找到后将该元素填入右侧坑位中,该位置就又形成了一个新的坑位,然后从右侧往左遍历找比基准值小的元素,将该元素填入到坑位中,重复上述过程直到左右标志位相遇,最后将基准值放入最后的坑位中,最对左右子序列重复上述过程直到整个序列排序完成。
图示基准值取的为最左侧元素:
3.前后指针版
初始设置cur、prev两个标志指针,cur标志序列第一个元素,prev标志cur前一个位置,cur位置的元素若大于基准值,cur向前前进,若小于基准值,则先对prev进行加一,然后交换cur和prev标记位置的元素,这样就能保证cur与prev之间的元素都大于基准值,prev之前的元素都小于基准值,重复上述过程,直到cur超过序列最右侧位置,最后进行一次判断,若prev标记位置不是序列最后一个位置,则将基准值交换到prev交换位置,即完成左右子序列划分,再对左右子序列重复上述过程,直到整个序列完成排序。
图示:
二、算法缺陷与优化
1.算法缺陷
1.1基准值取值
若所取基准值为序列中的最大或最小值,那么每趟划分后,基准值的一侧将会出现没有元素的情况,就相当于每趟只完成了一个元素的排序,那么快速排序的时间复杂度此时将达到,效率较低。
1.2递归超限
算法递归深度会随着元素的增多而加深,若序列元素元素量非常巨大,可能会造成递归超限。
2.优化方法
2.1三位取中法
我们可以通过选取序列中开始、中间、末尾三个位置中元素大小居中的元素作为基准值,此方法可以大大降低基准值取到序列内最大或最小元素的概率。
注:为了不改变代码的逻辑实现,若取得的基准值的位置不是基准值默认位置,可先提前将基准值交换到所取基准值的默认位置处。
2.2设置阈值
我们可以发现,随着划分次数的增加,子序列内的元素个数会不断减少,而且元素数量较少时,各类排序算法的效率差距其实是可以忽略不记的,而比较适合元素个数比较少的排序就是直接插入排序,所以我们可以在递归的过程中设计一个阈值,当子序列中元素个数不超过阈值时,递归不再进行,而是直接对当前子序列进行直接插入排序,然后返回。
当然,这种方法只能进行一定程度的缓解,并不能完全解决递归超限的问题。
所以,我们可以提前估算一下递归的深度N,再在递归深度上设置一个阈值count,若N小于count直接进行快速排序;若N大于count,先通过快速排序递归count次,然后再对每个分组中进行其他排序(比如:堆排序)。
这样即不会对算法效率造成太大影响,也避免了递归超限的问题。
2.3循环实现
使用循环的方式实现快速排序,我们会发现直接转循环并不容易实现,结合递归调用的特性,后调用的先结束退出,先调用的后结束,符合栈后进先出的特点,所有可以利用栈来循环实现快速排序。
每次对序列进行分割后,利用栈依次压入左右子序列的右边界和左边界,这样下一次循环就会依次拿到一个子序列的左右边界,然后再对此序列进行分割,循环进行上述操作,直到栈为空即完成序列排序。(可加上优化方法,为子序列元素个数加上阈值)
三、接口实现
1.快速排序
void QuickSort(int* array, int left, int right) {//快速排序 if (right - left <= 16) {//阈值设置为16 InsertSort(array + left, right - left);//达到阈值,进行直接插入排序 } else { if (right - left <= 1) {//区间内少于等于1个元素,不需要再排 return; } int div = PartSort2(array, left, right);//找一个基准值对区间元素进行划分,分割完成后返回基准值位置 QuickSort(array, left, div);//对基准值左侧进行快排 QuickSort(array, div + 1, right);//对基准值右侧进行快排 } }
2.hoare版
int PartSort1(int* array, int left, int right) {//快速排序hoare版 int begin = left; int end = right - 1; int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != end) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[end]); } int key = array[end];//取最后一个元素为基准值 while (begin < end) {//begin与end没有相遇 while (begin < end && array[begin] <= key) {//1.从前往后找比基准值key大的元素 begin++; } while (begin < end && array[end] >= key) {//2.从后往前找比基准值key小的元素 end--; } if (begin != end) {//3.交换两个元素 Swap(&array[begin], &array[end]); } } Swap(&array[begin], &array[right - 1]);//4.将基准值交换到区间中间位置 return begin;//5.返回基准值下标 }
3.挖坑法
int PartSort2(int* array, int left, int right) {//快速排序挖坑法 int begin = left; int end = right - 1; int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != end) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[end]); } int key = array[end];//1.取最后一个元素为基准值,并挖坑end while (begin < end) {//begin与end没有相遇 while (begin < end && array[begin] <= key) {//2.从前往后找比基准值key大的元素 begin++; } if (begin != end) { array[end--] = array[begin];//3.填end坑,挖begin坑,end前移一位 } while (begin < end && array[end] >= key) {//4.从后往前找比基准值key小的元素 end--; } if (begin != end) { array[begin++] = array[end];//5.填begin坑,挖end坑,begin后移一位 } } array[end] = key;//6.基准值放入最后一个坑内 return end;//7.返回基准值下标 }
4.前后指针版
int PartSort3(int* array, int left, int right) {//快速排序前后指针版 int cur = left;//标记第一个元素 int prev = cur - 1;//位于cur之后一个位置 int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != right - 1) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[right - 1]); } int key = array[right - 1];//取最后一个元素为基准值 while (cur < right) { if (array[cur] < key && ++prev != cur) {//始终保存cur与prev之间的元素都大于基准值key Swap(&array[cur], &array[prev]); } cur++; } if (++prev != right - 1) {//将基准值交换到prev的下一个位置 Swap(&array[right - 1], &array[prev]); } return prev; }
5.非递归版(循环版)
void QuickSortNonR(int* array, int size) {//快速排序非递归版 Stack s; StackInit(&s); int left = 0; int right = size; StackPush(&s, right);//依次压入右边界和左边界 StackPush(&s, left); while (!StackEmpty(&s)) {//栈不为空,循环进行 left = StackTop(&s);//获取栈顶,左边界 StackPop(&s);//左边界出栈 right = StackTop(&s);//获取栈顶,右边界 StackPop(&s);//右边界出栈 if (right - left <= 16) {//子序列达到阈值 InsertSort(array + left, right - left);//进行直接插入排序 } else { int div = PartSort2(array, left, right); //基准值左侧(left,div);右侧(div+1,right) //先压右侧子序列边界(先压右边界,再压左边界) StackPush(&s, right); StackPush(&s, div + 1); //再压左侧子序列边界 StackPush(&s, div); StackPush(&s, left); } } }
四、接口测试
1.测试用例
void TestPartSort() {//快速排序测试函数 int array[] = { 5,4,8,1,9,7,3,2,6,0 }; int length = sizeof(array) / sizeof(array[0]); printf("排序前:"); PrintArray(array, length); QuickSort(array, 0, length);//递归方式 //QuickSortNonR(array, length);//非递归方式 printf("\n排序后:"); PrintArray(array, length); }
2.测试结果
递归方式:
非递归方式:
五、性能分析
1.时间复杂度
最坏情况
如果序列初始已经是有序状态,那么每次都只能排好一个元素,且剩余元素都需要逐个比较一次,此时快速排序的时间性能退化比较严重,所有最坏情况下时间复杂度为。
最好情况
根据算法思想可知,如果每次划分的基准值都恰好处于中间大小位置,那么每次都能将其划分为两个等长的子序列,此时快速排序的时间性能达到最佳,所有最好情况下时间复杂度为
平均情况
综合两种情况,快速排序的时间复杂度为。
2.空间复杂度
由于递归过程中,需要保存相关信息,且大小与递归次数有关,所有其空间复杂度为。
3.稳定性
根据算法思想可知,在划分过程中会改变相同元素的相对位置,所有快速排序是不稳定的。
六、完整代码
1.QuickSort.c
#include"Stack.h" void QuickSort(int* array, int left, int right);//快速排序 int PartSort1(int* array, int left, int right);//快速排序hoare版 int PartSort2(int* array, int left, int right);//快速排序挖坑法 int PartSort3(int* array, int left, int right);//快速排序前后指针版 int GetMiddleIndex(int* array, int left, int right);//三位取中法(选取元素大小相对靠中的为基准值) void QuickSortNonR(int* array, int size);//快速排序非递归版 void InsertSort(int* array, int size);//直接插入排序 void PrintArray(int* array, int size);//数组打印 void Swap(int* num1, int* num2);//整数交换 void TestPartSort();//快速排序测试函数 int main() { TestPartSort(); return 0; } void QuickSort(int* array, int left, int right) {//快速排序 if (right - left <= 16) {//阈值设置为16 InsertSort(array + left, right - left);//达到阈值,进行直接插入排序 } else { if (right - left <= 1) {//区间内少于等于1个元素,不需要再排 return; } int div = PartSort2(array, left, right);//找一个基准值对区间元素进行划分,分割完成后返回基准值位置 QuickSort(array, left, div);//对基准值左侧进行快排 QuickSort(array, div + 1, right);//对基准值右侧进行快排 } } int PartSort1(int* array, int left, int right) {//快速排序hoare版 int begin = left; int end = right - 1; int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != end) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[end]); } int key = array[end];//取最后一个元素为基准值 while (begin < end) {//begin与end没有相遇 while (begin < end && array[begin] <= key) {//1.从前往后找比基准值key大的元素 begin++; } while (begin < end && array[end] >= key) {//2.从后往前找比基准值key小的元素 end--; } if (begin != end) {//3.交换两个元素 Swap(&array[begin], &array[end]); } } Swap(&array[begin], &array[right - 1]);//4.将基准值交换到区间中间位置 return begin;//5.返回基准值下标 } int PartSort2(int* array, int left, int right) {//快速排序挖坑法 int begin = left; int end = right - 1; int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != end) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[end]); } int key = array[end];//1.取最后一个元素为基准值,并挖坑end while (begin < end) {//begin与end没有相遇 while (begin < end && array[begin] <= key) {//2.从前往后找比基准值key大的元素 begin++; } if (begin != end) { array[end--] = array[begin];//3.填end坑,挖begin坑,end前移一位 } while (begin < end && array[end] >= key) {//4.从后往前找比基准值key小的元素 end--; } if (begin != end) { array[begin++] = array[end];//5.填begin坑,挖end坑,begin后移一位 } } array[end] = key;//6.基准值放入最后一个坑内 return end;//7.返回基准值下标 } int PartSort3(int* array, int left, int right) {//快速排序前后指针版 int cur = left;//标记第一个元素 int prev = cur - 1;//位于cur之后一个位置 int midPos = GetMiddleIndex(array, left, right);//获取元素大小相对靠中的元素下标 if (midPos != right - 1) {//将该元素交换到末尾,从而不影响后续代码的逻辑实现 Swap(&array[midPos], &array[right - 1]); } int key = array[right - 1];//取最后一个元素为基准值 while (cur < right) { if (array[cur] < key && ++prev != cur) {//始终保存cur与prev之间的元素都大于基准值key Swap(&array[cur], &array[prev]); } cur++; } if (++prev != right - 1) {//将基准值交换到prev的下一个位置 Swap(&array[right - 1], &array[prev]); } return prev; } void QuickSortNonR(int* array, int size) {//快速排序非递归版 Stack s; StackInit(&s); int left = 0; int right = size; StackPush(&s, right);//依次压入右边界和左边界 StackPush(&s, left); while (!StackEmpty(&s)) {//栈不为空,循环进行 left = StackTop(&s);//获取栈顶,左边界 StackPop(&s);//左边界出栈 right = StackTop(&s);//获取栈顶,右边界 StackPop(&s);//右边界出栈 if (right - left <= 16) {//子序列达到阈值 InsertSort(array + left, right - left);//进行直接插入排序 } else { int div = PartSort2(array, left, right); //基准值左侧(left,div);右侧(div+1,right) //先压右侧子序列边界(先压右边界,再压左边界) StackPush(&s, right); StackPush(&s, div + 1); //再压左侧子序列边界 StackPush(&s, div); StackPush(&s, left); } } } int GetMiddleIndex(int* array, int left, int right) {//三位取中法(选取元素大小相对靠中的为基准值) int mid = (left + right) / 2; //返回基准值的下标 if (array[left] < array[right - 1]) { if (array[mid] < array[left]) { return left; } else if (array[mid] > array[right - 1]) { return right - 1; } else { return mid; } } else {//array[left]>=array[right-1] if (array[mid] > array[left]) { return left; } else if (array[mid] < array[right - 1]) { return right - 1; } else { return mid; } } } void InsertSort(int* array, int size) {//直接插入排序 for (int i = 1; i < size; i++) {//从1开始循环,默认数组内第一个元素为有序序列 int end = i - 1;//标记已排序序列最后位置下标 int key = array[i];//依次拿取数组内元素 while (end >= 0 && key < array[end]) {//key从前往后比较:小于当前元素,继续往前走 array[end + 1] = array[end];//将当前元素往后移一个位置 end--; } array[end + 1] = key;//key大于等于当前元素,插入到当前位置之后 } } void PrintArray(int* array, int size) {//数组打印 for (int i = 0; i < size; i++) { printf("%d ", array[i]); } } void Swap(int* num1, int* num2) {//整数交换 int temp = *num1; *num1 = *num2; *num2 = temp; } void TestPartSort() {//快速排序测试函数 int array[] = { 5,4,8,1,9,7,3,2,6,0 }; int length = sizeof(array) / sizeof(array[0]); printf("排序前:"); PrintArray(array, length); //QuickSort(array, 0, length);//递归方式 QuickSortNonR(array, length);//非递归方式 printf("\n排序后:"); PrintArray(array, length); }
2.Stack.h
#pragma once #include<stdio.h> #include<assert.h> #include<string.h> #include<malloc.h> // 支持动态增长的栈 typedef int STDataType; typedef struct Stack { STDataType* array; int top; // 栈顶 int capacity; // 容量 }Stack; /*void StackInit(Stack* ps, int capacity);// 初始化栈 void StackPush(Stack* ps, STDataType data);// 入栈 void StackPop(Stack* ps);// 出栈 STDataType StackTop(Stack* ps);// 获取栈顶元素 int StackSize(Stack* ps);// 获取栈中有效元素个数 int StackEmpty(Stack* ps);// 检测栈是否为空,为空返回1,非空返回0 void StackDestroy(Stack* ps);// 销毁栈 void Stack_Test();//功能测试函数 */ // 初始化栈 void StackInit(Stack* ps) { assert(ps); ps->array = (STDataType*)malloc(sizeof(STDataType) * 3); if (ps->array == NULL) {//申请失败 assert(0); return; } ps->capacity = 3; ps->top = 0; } void StackCheckCapacity(Stack* ps) { assert(ps); if (ps->top == ps->capacity) {//栈满 int newCapacity = ps->capacity * 2;//每次扩容扩大两倍 STDataType* temp = (STDataType*)malloc(sizeof(STDataType) * newCapacity);//开辟新空间 if (temp == NULL) { assert(0); return; } memcpy(temp, ps->array, sizeof(STDataType) * ps->capacity);//将元数据拷贝到新空间中 free(ps->array);//释放旧空间,使用新空间 ps->array = temp; ps->capacity = newCapacity; } } void StackPush(Stack* ps, STDataType data) {//入栈 StackCheckCapacity(ps);//判断是否栈满,栈满则扩容 ps->array[ps->top] = data; ps->top++; } int StackEmpty(Stack* ps) {//栈判空 assert(ps); return 0 == ps->top;//空返回1,非空返回0 } void StackPop(Stack* ps) {//出栈 if (StackEmpty(ps)) {//判空 return; } ps->top--; } STDataType StackTop(Stack* ps) {//获取栈顶元素 assert(ps); return ps->array[ps->top - 1]; } int StackSize(Stack* ps) {//获取栈中有效元素个数 assert(ps); return ps->top; } void StackDestroy(Stack* ps) {//栈销毁 assert(ps); if (ps->array) { free(ps->array); ps->capacity = 0; ps->top = 1; } }