快速排序:
写快排的注意事项
1.单趟排序hoare
2.不写优化只说优化就行
理想的情况下:每次排序都是二分,直到二分到最后,那就相当于递归高度次(logN),每一层单趟排都是O(N),时间复杂度O(NlogN) 空间复杂度就是递归的深度 O(logN)
最坏的情况:数组有序(或者接近有序)是最坏情况。假如数组是升序,key是left,要排升序,有N层,则每层都排(N,N-1,N-2…,3,2,1),时间复杂度是O(NN)。因此,优化最差情况1.三数取中 2.小区间排序优化
相遇有两种情况:right找left或者left找right
如果是right先走,最坏在遇见left的时候停止,而这个位置是left的位置,值肯定是比key小的,但left先走的话,最坏情况遇见right,此时值是比key大的,交换会出问题。因此,left是key需要让right先走。
void QuickSort(int* a,int n)
{
int left = 0, right = n-1;
int keyi = left;
while(left < right)
{
// 找小
// 如果没有left < right可能会导致right一直--,越界访问eg:[1,2,3,4,5]
while(left < right && a[right] >= a[keyi])
--right;
// 找大
//如果没有a[left] <= a[keyi]的=符号时,可能会导致死循环eg:[5,1,2,5,6];
while(left < right && a[left] <= a[keyi])
++left;
Swap(&a[left], &a[right]);
}
Swap(&a[keyi], &a[left]);
}
第二阶段:分治(递归),对每次的key的左右区间看成子问题,传入递归一直分治,最终完成。递归结束标志:当这个子区间不存在或者只有单个值时返回。
void QuickSort(int* a, int begin, int end)
{
if(begin >= end) // 迭代的终点:左区间>=右区间时,迭代结束。
return;
int left = begin, right = end-1;
int keyi = left;
while(left < right)
{
// 找小
while(left < right && a[right] >= a[keyi])
--right;
// 找大
while(left < right && a[left] <= a[keyi])
++left;
Swap(&a[left], &a[right]);
}
int meeti = left;
Swap(&a[keyi], &a[meeti]);
//左:[begin,meeti-1] meeti 右:[meeti+1,end]
QuickSort(a, begin, meeti-1); // 传参就是[0, n-1]
QuickSort(a, meeti+1, end);
}
挖坑法(代码在后面):将key的位置看作是一个坑(hole),然后利用part1的思路,找到大的将值放进坑中,然后此时right就是hole,再找小的值,将小的值放进hole中,此时left是坑…,直到相遇,把最早的key值放进相遇的坑中。
挖坑法的初次排序后的值可能和hoare的值不同。
前后指针法(代码在后面):cur找比keyi位置小的值,找到后prev++,再交换prev和cur的值,直到cur走到数组尾,然后将keyi和prev的位置的值交换。
快排最差情况的优化1:(防止最坏的情况)
三数取中:对快排影响最大的是key,如果key取值越接近中位数则越接近二分,效率越高
三数取中是在最左边,左右边和最中间的值,选择他们三个不是最大也不是最小的数。
三数取中,将需要的那个中间的数换到最左边的位置,这样就不用操心key位置变换的问题了。
// 快速排序的三数取中 优化快排的最差情况
int GetMidIndex(int* a, int left, int right)
{
int mid = (left + right) >> 1; //就是(left + right)/2
// left, mid, right
if(a[left] < a[mid])
{
if(a[mid] < a[right])
return mid;
else if(a[left] > a[right])
return left;
else
return right;
}
else // a[left] > a[mid]
{
if(a[mid] > a[right])
return mid;
else if(a[left] < a[right])
return left;
else
return right;
}
}
// 快速排序hoare版本 -- 最原始的 左右指针法
int PartSort1(int* a, int left, int right)
{
int midIndex = GetMidIndex(a, left, right); //找到中间数
Swap(&a[left], &a[midIndex]); //再交换到left位置
int keyi = left;
while(left < right)
{
// 找小
while(left < right && a[right] >= a[keyi])
--right;
// 找大
while(left < right && a[left] <= a[keyi])
++left;
Swap(&a[left], &a[right]);
}
Swap(&a[keyi], &a[left]);
return left;
}
// 迭代过程
void QuickSort(int* a, int begin, int end)
{
if(begin >= end) // 迭代的终点:左区间>=右区间时,迭代结束。
return;
int keyi = PartSort1(a, begin, end);
//左:[begin,keyi -1] keyi 右:[keyi +1,end]
QuickSort(a, begin, keyi -1); // 传参就是[0, n-1]
QuickSort(a, keyi +1, end);
}
快排最差情况的优化2:(防止最后递归太多)
小区间优化:对于最后面几个数的排序不再进行递归,直接选择插入排序对最后几个数进行排序。(作用不大)
// 迭代过程
void QuickSort(int* a, int begin, int end)
{
if(begin >= end) // 迭代的终点:左区间>=右区间时,迭代结束。
return;
// 如果子区间数据较多继续选key单趟,分割子区间分治递归
if(end - begin > 20)
{
int keyi = PartSort1(a, begin, end);
//左:[begin,keyi -1] keyi 右:[keyi +1,end]
QuickSort(a, begin, keyi -1); // 传参就是[0, n-1]
QuickSort(a, keyi +1, end);
}
else //如果子区间数很少了还要进行递归就很不划算。
{
//数据少了直接使用插入排序算法。 插入排序传的是数组和元素个数。
InsortSort(a+begin, end-begin+1);
}
}
将快排中的单趟排序独立出来
// 快速排序hoare版本 -- 最原始的 左右指针法
int PartSort1(int* a, int left, int right)
{
int keyi = left;
while(left < right)
{
// 找小
while(left < right && a[right] >= a[keyi])
--right;
// 找大
while(left < right && a[left] <= a[keyi])
++left;
Swap(&a[left], &a[right]);
}
Swap(&a[keyi], &a[left]);
return left;
}
//挖坑法:
int PartSort2(int* a, int left, int right)
{
int midIndex = GetMidIndex(a, left, right); //找到中间数
Swap(&a[left], &a[midIndex]); //再交换到left位置
int key = a[left];
while(left < right)
{
//找小
while(left < right && a[right] >= key )
right--;
// 放到左边的坑位,右边就是新的坑
a[left] = a[right];
//找大
while(left < right && a[left] <= key )
left++;
//放到右边的坑位,左边就是新的坑
a[right] = a[left];
}
a[left] = key;
return left;
}
// 快速排序前后指针法
int PartSort3(int* a, int left, int right)
{
int midIndex = GetMidIndex(a, left, right); //找到中间数
Swap(&a[left], &a[midIndex]); //再交换到left位置
int keyi = left;
int prev = left, cur = prev+1;
while(cur <= right)
{
if(a[cur] < a[key])
{
prev++;
if(prev != cur)
{
Swap(&a[cur], &a[prev]);
}
}
cur++
}
Swap(&a[cur], &a[keyi]);
return prev;
}
// 迭代过程
void QuickSort(int* a, int begin, int end)
{
if(begin >= end) // 迭代的终点:左区间>=右区间时,迭代结束。
return;
int keyi = PartSort1(a, begin, end);
//左:[begin,keyi -1] keyi 右:[keyi +1,end]
QuickSort(a, begin, keyi -1); // 传参就是[0, n-1]
QuickSort(a, keyi +1, end);
}
非递归法:递归最大问题是如果递归的太深,栈空间不够导致溢出,只能改成非递归
1.改循环—斐波那契
2.树遍历的非递归和快排非递归,只能用stack(栈)存储数据模逆递归
非递归法将数组区间一直入栈,分割后先入左边再入右边,这样栈可以每次先出右边再出左边,。直到分割的区间只有一个数或者没有数了这半边就不用入栈了。
栈里面存的是数组的下标。
// 实现支持动态增长的栈
typedef int STDataType;
typedef struct Stack
{
STDataType* a;
int top; // 栈顶
int capacity; // 容量
}Stack;
// 初始化栈
void StackInit(Stack* pst)
{
assert(pst);
pst->a = (STDataType*)malloc(sizeof(STDataType)*4);
pst->top = 0;
pst->capacity = 4;
}
// 销毁栈
void StackDestroy(Stack* pst)
{
assert(pst);
free(pst->a);
pst->a = NULL;
pst->top = 0;
pst->capacity = 0;
}
// 入栈
void StackPush(Stack* pst, STDataType data)
{
assert(pst);
if(pst->top == pst->capacity)
{
STDataType* tmp = (STDataType*)realloc(pst->a, sizeof(STDataType)*pst->capacity*2);
if(tmp == NULL)
{
printf("realloc fail\n");
exit(-1); //结束掉整个程序
}
pst->a = tmp;
pst->capacity *= 2;
}
pst->a[pst->top] = data;
pst->top++;
}
// 出栈
void StackPop(Stack* pst)
{
assert(pst);
assert(!StackEmpty(pst)); //如果栈已经空了就不要top--了
pst->top--;
}
// 获取栈顶元素
STDataType StackTop(Stack* pst)
{
assert(pst);
assert(!StackEmpty(pst)); //如果栈已经空了就不要top--了
return pst->a[pst->top - 1];
}
// 检测栈是否为空,如果为空返回非零结果,如果不为空返回0
int StackEmpty(Stack* pst)
{
assert(pst);
return pst->top == 0;
}
int StackSize(Stack* pst)
{
assert(pst);
return pst->top;
}
// 快速排序 非递归实现
void QuickSortNonR(int* a, int left, int right)
{
Stack st;
StackInit(&st);
//先入左再入右
StackPush(&st, begin);
StackPush(&st, end);
while(!StackEmpty(&st))
{
int left, right;
right = StackTop(&st);
StackPop(&st);
left = StackTop(&st);
StackPop(&st);
int keyi = PartSort1(a, left, right);
if(left < keyi - 1) //如果左值小于keyi-1说明左边还有区间
{
StackPush(&st, left);
StackPush(&st, keyi-1);
}
if(keyi+1 < right)
{
StackPush(&st, keyi+1);
StackPush(&st, right);
}
}
StackDestory(&st);
}