博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
快速排序实现及其pivot的选取
阅读量:5010 次
发布时间:2019-06-12

本文共 3903 字,大约阅读时间需要 13 分钟。

coursera上斯坦福的算法专项在讲到快速排序时,称其为最优雅的算法之一。快速排序确实是一种比较有效的排序算法,很多类库中也都采用了这种排序算法,其最坏时间复杂度为$O(n^2)$,平均时间复杂度为$O(nlogn)$,且其不需要额外的存储空间。

基本步骤

快速排序主要使用了分治的思想,通过选取一个pivot,将一个数组划分为两个子数组。其步骤为:

1.从数组中选择一个元素作为pivot
2.重新排列数组,小于pivot的在pivot的左边,大于pivot的在其右边。
3.递归地对划分后的左右两部分重复上述步骤。

简单的伪代码如下:

qsort%E9%80%89%E5%8C%BA_001.png

其中最主要的就是partition划分过程了。

划分过程

partition过程需要首先选择一个pivot,然后将小于pivot的元素放到左半部分,大于pivot的放到右半部分,并且最终pivot的位置及为其在排序好的数组中的最终位置。

这里使用第一个元素作为pivot,若选择其他元素作为pivot,则将其交换到第一个元素,这样可以保证代码的一致性及容易实现。示意图如下:

qsort%E9%80%89%E5%8C%BA_002.png

这里使用i和j,i和j最初为p+1的位置,在遍历的过程中i始终指向>p的第一个元素,j始终指向当前待遍历的元素,若a[j] < p,则将其与a[i]进行交换。相关过程如下:

qsort%E9%80%89%E5%8C%BA_003.png

基本实现如下:

/**    * a[l+1],...,a[i-1] < p    * a[i],...,a[j-1] > p    */    private static int partition(int[] a, int l, int r) {        int p = a[l];        int i = l + 1;        for (int j=l+1; j<=r; j++) {            if (a[j] < p) {                swap(a, j, i);                i++;            }        }        swap(a, l, i-1);        return i-1;    }

基本实现

public class QuickSort {    public static void qSort(int[] a) {        if (a == null || a.length <= 1) {            return;        }            qSort(a, 0, a.length-1);    }    private static void qSort(int[] a, int l, int r)    {        if (l >= r) {            return;        }        int pos = partition(a, l, r);        qSort(a, l, pos - 1);        qSort(a, pos + 1, r);    }    /**    * a[l+1],...,a[i-1] < p    * a[i],...,a[j-1] > p    */    private static int partition(int[] a, int l, int r) {        int p = a[l];        int i = l + 1;        for (int j=l+1; j<=r; j++) {            if (a[j] < p) {                swap(a, j, i);                i++;            }        }        swap(a, l, i-1);        return i-1;    }    //返回pivot下标 选择第一个元素    private static int choosePivotFirst(int[] a, int l, int r) {        return l;    }       private static void swap(int[] a, int x, int y) {        int temp = a[x];        a[x] = a[y];        a[y] = temp;    }

pivot的选取

根据斯坦福算法专项课,然我们实现三种不同的pivot选取方式,并计算相应比较次数,分别为choose first, choose last, median of three, 还可以进行随机选取,这也是快速排序为什么是一种随机化算法。

pivot的选取决定了快速排序的运行时间,下面对几种特殊情况进行分析:

1.最坏情况

假设我们始终选取第一个元素作为pivot, 并且输入数组是有序的,那么每次划分后面所有元素都大于pivot, 每次只能将问题规模减少1,所以运行时间为$n+n-1+n-2+...+1$ = $O(n^2)$.

2.最好情况

最好情况为每次选取的pivot都能将数组平均地划分为两部分,由于划分的过程为$O(n)$,所以总的运行时间为$$T(n) = 2T(n/2) + O(n)$$根据主方法,时间复杂度为O(nlogn)。

3.随机选取

每次运行过程中,随机选取pivot, 通常能得到比较好的结果。

选取方式及实现

斯坦福算法专项课上让我们实现三种不同的选取方式,选取第一个,最后一个,以及三数取中。

1.choose first

该种方式最为简单,只需返回子数组的第一个元素下标即可,下面为其实现:

//返回pivot下标 选择第一个元素private static int choosePivotFirst(int[] a, int l, int r) {    return l;}

2.choose last

选择最后一个元素,实现如下:

//选择最后一个元素作为pivotprivate static int choosePivotLast(int[] a, int l, int r) {        return r;}

3.median-of-three

选取第一个、最后一个以及中间的元素的中位数,如4 5 6 7, 第一个4, 最后一个7, 中间的为5, 这三个数的中位数为5, 所以选择5作为pivot,8 2 5 4 7, 三个元素分别为8 5 7, 中位数为7, 所以选择最后一个元素7作为pivot,其实现如下:

//median-of-three pivot ruleprivate static int choosePivotMedianOfThree(int[] a, int l, int r) {        int mid = 0;    if ((r-l+1) % 2 == 0) {        mid = l + (r-l+1)/2 - 1;    } else {        mid = l + (r-l+1)/2;    }    //只需要找出中位数即可,不需要交换    //有的版本也可以进行交换    if (((a[l]-a[mid]) * (a[l]-a[r])) <= 0) {        return l;    } else if (((a[mid]-a[l]) * (a[mid]-a[r])) <= 0)    {        return mid;    } else {        return r;    }}

最后的划分过程如下:

private static int partition(int[] a, int l, int r) {    //pivot选择方式    //int pi = choosePivotFirst(a, l, r);    //int pi = choosePivotLast(a, l, r);    int pi = choosePivotMedianOfThree(a, l, r);    //始终将第一个元素作为pivot, 若不是, 则与之交换    if (pi != l) {        swap(a, pi, l);    }    int p = a[l];    int i = l + 1;    for (int j=l+1; j<=r; j++) {        if (a[j] < p) {            swap(a, j, i);            i++;        }    }    swap(a, l, i-1);    return i-1;}

注意最后的划分过程相比于之前增加的pivot的选取方式,而不是单纯地将第一个元素作为pivot, 可以看到,若第一个元素不是pivot, 需要将pivot与第一个元素进行交换,这样保证代码的统一性。

总结与感想

1.学会体会这些算法背后的思想,为什么要这样设计

2.对于比较复杂的算法,学会使用特殊情况进行分析

参考资料:

(1) coursera斯坦福算法专项课part1

(2)

转载于:https://www.cnblogs.com/litexy/p/9744355.html

你可能感兴趣的文章
关于C语言中return的一些总结
查看>>
Codeforces Round #278 (Div. 2)
查看>>
51. N-Queens
查看>>
Linux 命令 - 文件搜索命令 locate
查看>>
CTFMON。exe
查看>>
spark
查看>>
[Angular] Angular Global Keyboard Handling With EventManager
查看>>
[Python] Execute a Python Script
查看>>
[Angular 2] Using Promise to Http
查看>>
[Grunt] grunt.template
查看>>
一、基础篇--1.1Java基础-Object类中常见的方法,为什么wait notify会放在Object里边...
查看>>
UVa 10079 - Pizza Cutting
查看>>
Ubuntu最小化桌面快捷键Super+D不生效解决
查看>>
Cookie&Session会话跟踪技术
查看>>
UNIX环境高级编程 第17章 高级进程间通信
查看>>
ES的Zen发现机制
查看>>
【hibernate】1、Hibernate的一个注解 @Transient
查看>>
flex入门----基础知识
查看>>
HihoCoder 1877 - Approximate Matching
查看>>
【转】C++多继承的细节
查看>>