【转载】十大排序

转载自十大排序 并进行了优化补充。

0、算法概述

0.1 算法分类

十种常见排序算法可以分为两大类：

• 比较类排序：通过比较来决定元素间的相对次序，由于其时间复杂度不能突破O(nlogn)，因此也称为非线性时间比较类排序。
• 非比较类排序：不通过比较来决定元素间的相对次序，它可以突破基于比较排序的时间下界，以线性时间运行，因此也称为线性时间非比较类排序。

(如果洛谷打不开请点击链接)

0.2 算法复杂度

排序方法	时间复杂度(平均)	时间复杂度(最坏)	时间复杂度(最好)	空间复杂度	稳定性
插入排序	\(O(n^2)\)	\(O(n^2)\)	\(O(n)\)	\(O(1)\)	稳定
希尔排序	\(O(n^{1.3})\)	\(O(n^2)\)	\(O(n)\)	\(O(1)\)	不稳定
选择排序	\(O(n^2)\)	\(O(n^2)\)	\(O(n^2)\)	\(O(1)\)	不稳定
堆排序	\(O(nlog_2n)\)	\(O(nlog_2n)\)	\(O(nlog_2n)\)	\(O(1)\)	不稳定
冒泡排序	\(O(n^2)\)	\(O(n^2)\)	\(O(n)\)	\(O(1)\)	稳定
快速排序	\(O(nlog_2n)\)	\(O(n^2)\)	\(O(nlog_2n)\)	\(O(nlog_2n)\)	不稳定
归并排序	\(O(nlog_2n)\)	\(O(nlog_2n)\)	\(O(nlog_2n)\)	\(O(n)\)	稳定
计数排序	\(O(n+k)\)	\(O(n+k)\)	\(O(n+k)\)	\(O(n+k)\)	稳定
桶排序	\(O(n+k)\)	\(O(n^2)\)	\(O(n)\)	\(O(n+k)\)	稳定
基数排序	\(O(n \times k)\)	\(O(n \times k)\)	\(O(n \times k)\)	\(O(n+k)\)	稳定

0.3 相关概念

• 稳定：如果\(a\)原本在\(b\)前面，而\(a=b\)，排序之后\(a\)仍然在\(b\)的前面。
• 不稳定：如果\(a\)原本在\(b\)的前面，而\(a=b\)，排序之后$ a $可能会出现在 \(b\) 的后面。
• 时间复杂度：对排序数据的总的操作次数。反映当\(n\)变化时，操作次数呈现什么规律。
• 空间复杂度：是指算法在计算机

内执行时所需存储空间的度量，它也是数据规模\(n\)的函数。

1、冒泡排序（Bubble Sort）

冒泡排序是一种简单的排序算法。它重复地走访过要排序的数列，一次比较两个元素，如果它们的顺序错误就把它们交换过来。走访数列的工作是重复地进行直到没有再需要交换，也就是说该数列已经排序完成。这个算法的名字由来是因为越小的元素会经由交换慢慢“浮”到数列的顶端。

1.1 算法描述

• 比较相邻的元素。如果第一个比第二个大，就交换它们两个；
• 对每一对相邻元素作同样的工作，从开始第一对到结尾的最后一对，这样在最后的元素应该会是最大的数；
• 针对所有的元素重复以上的步骤，除了最后一个；
• 重复步骤1~3，直到排序完成。

1.2 动图演示

1.3 代码实现

void bubble_sort(int a[], int n) {
	for(int i = 1; i < n; i++)
	    for(int j = 1; j < n - i; j++) 
	        if(a[j] > a[j + 1]) swap(a[j], a[j + 1]);
} 

2、选择排序（Selection Sort）

选择排序(Selection-sort)是一种简单直观的排序算法。它的工作原理：首先在未排序序列中找到最小（大）元素，存放到排序序列的起始位置，然后，再从剩余未排序元素中继续寻找最小（大）元素，然后放到已排序序列的末尾。以此类推，直到所有元素均排序完毕。

2.1 算法描述

\(n\)个记录的直接选择排序可经过\(n-1\)趟直接选择排序得到有序结果。具体算法描述如下：

• 初始状态：无序区为\(R[1..n]\)，有序区为空；
• 第\(i\)趟排序\((i=1,2,3…n-1)\)开始时，当前有序区和无序区分别为\(R[1..i-1]\)和\(R(i..n）\)。该趟排序从当前无序区中-选出关键字最小的记录 \(R[k]\)，将它与无序区的第\(1\)个记录\(R\)交换，使\(R[1..i]\)和\(R[i+1..n]\)分别变为记录个数增加\(1\)个的新有序区和记录个数减少\(1\)个的新无序区；
• \(n-1\)趟结束，数组有序化了。

2.2 动图演示

　　

2.3 代码实现

void selection_sort(int a[], int n) {
	for(int i = 1; i < n; i++) 
	    for(int j = i + 1, j <= n; j++)
	        if(a[i] > a[j]) swap(a[j], a[i]);
} 

2.4 算法分析

表现最稳定的排序算法之一，因为无论什么数据进去都是\(O(n^2)\)的时间复杂度，所以用到它的时候，数据规模越小越好。唯一的好处可能就是不占用额外的内存空间了吧。理论上讲，选择排序可能也是平时排序一般人想到的最多的排序方法了吧。

3、插入排序（Insertion Sort）

插入排序（Insertion-Sort）的算法描述是一种简单直观的排序算法。它的工作原理是通过构建有序序列，对于未排序数据，在已排序序列中从后向前扫描，找到相应位置并插入。

3.1 算法描述

一般来说，插入排序都采用in-place在数组上实现。具体算法描述如下：

• 从第一个元素开始，该元素可以认为已经被排序；
• 取出下一个元素，在已经排序的元素序列中从后向前扫描；
• 如果该元素（已排序）大于新元素，将该元素移到下一位置；
• 重复步骤\(3\)，直到找到已排序的元素小于或者等于新元素的位置；
• 将新元素插入到该位置后；
• 重复步骤\(2~5\)。

3.2 动图演示

3.2 代码实现

void insert_sort(int a[], int n) {
	for(int i = 2; i <= n; i++) 
	    for(j = i; j > 1; j--)
	        if(a[j] < a[j - 1]) swap(a[j], a[j - 1]); else break;
} 

3.4 算法分析

插入排序在实现上，通常采用in-place排序（即只需用到\(O(1)\)的额外空间的排序），因而在从后向前扫描过程中，需要反复把已排序元素逐步向后挪位，为最新元素提供插入空间。

4、希尔排序（Shell Sort）

1959年Shell发明，第一个突破\(O(n^2)\)的排序算法，是简单插入排序的改进版。它与插入排序的不同之处在于，它会优先比较距离较远的元素。希尔排序又叫缩小增量排序。

4.1 算法描述

先将整个待排序的记录序列分割成为若干子序列分别进行直接插入排序，具体算法描述：

• 选择一个增量序列\(t_1，t_2，…，t_k\)，其中\(t_i>t_j，t_k=1\)；
• 按增量序列个数k，对序列进行k 趟排序；
• 每趟排序，根据对应的增量\(t_i\)，将待排序列分割成若干长度为\(m\) 的子序列，分别对各子表进行直接插入排序。仅增量因子为\(1\) 时，整个序列作为一个表来处理，表长度即为整个序列的长度。

4.2 动图演示

4.3 代码实现

void shell_sort(int a[], int n) {
	for(int div = n / 2; div >= 1; div /= 2)
	    for(int k = 1; k <= div; k++)
	        for(int i = div + k; i <= n; i += div)
	            for(int j = i; j > k; j -= div)
	                if(a[j] < a[j - div]) swap(a[j], a[j - div]); else break;
}

4.4 算法分析

希尔排序的核心在于间隔序列的设定。既可以提前设定好间隔序列，也可以动态的定义间隔序列。动态定义间隔序列的算法是《算法（第4版）》的合著者Robert Sedgewick提出的。　

5、归并排序（Merge Sort）

归并排序是建立在归并操作上的一种有效的排序算法。该算法是采用分治法（Divide and Conquer）的一个非常典型的应用。将已有序的子序列合并，得到完全有序的序列；即先使每个子序列有序，再使子序列段间有序。若将两个有序表合并成一个有序表，称为\(2\)路归并。

5.1 算法描述

• 把长度为n的输入序列分成两个长度为\(n/2\)的子序列；
• 对这两个子序列分别采用归并排序；
• 将两个排序好的子序列合并成一个最终的排序序列。

5.2 动图演示

5.3 代码实现

版本1

void msort(int b, int e) {
	int (b == e) return;
	int mid = (b + e) / 2, i = b, j = mid + 1, k = b;
	msort(b, mid), msort(mid + 1, e);
	while(i <= mid && j <= e) {
		if(a[i] <= a[j]) c[k++] = a[i++];
		else c[k++] = a[j++];
	}
	while(i <= mid) c[k++] = a[i++];
	while(j <= e) c[k++] = a[j++];
	for(int l = b; l <= e; l++) a[l] = c[l];
} 

版本2

void msort(int b,int e) {
	if(b === e) return;
	int mid = (b + e) / 2, i = b, j = mid + 1, k = b;
	msort(b, mid), msort(mid + 1, e);
	while(k <= e) {
		if(j > e || i <= mid && a[i] <= a[j]) c[k++] = a[i++];
		else c[k++] = a[j++];
	}
	for(int l = b; l <= e; l++) {
		a[l] = c[l];
	}
}

5.4 算法分析

归并排序是一种稳定的排序方法。和选择排序一样，归并排序的性能不受输入数据的影响，但表现比选择排序好的多，因为始终都是\(O(nlog_2n）\)的时间复杂度。代价是需要额外的内存空间。

6、快速排序（Quick Sort）

快速排序的基本思想：通过一趟排序将待排记录分隔成独立的两部分，其中一部分记录的关键字均比另一部分的关键字小，则可分别对这两部分记录继续进行排序，以达到整个序列有序。

6.1 算法描述

快速排序使用分治法来把一个串（list）分为两个子串（sub-lists）。具体算法描述如下：

• 从数列中挑出一个元素，称为 “基准”（pivot）；
• 重新排序数列，所有元素比基准值小的摆放在基准前面，所有元素比基准值大的摆在基准的后面（相同的数可以到任一边）。在这个分区退出之后，该基准就处于数列的中间位置。这个称为分区（partition）操作；
• 递归地（recursive）把小于基准值元素的子数列和大于基准值元素的子数列排序。

6.2 动图演示

6.3 代码实现

void qsort(int l, int r) {//应用二分思想 
	int mid = a[(l + r) / 2];//中间数 
	int i = l, j = r;
	do{
		while(a[i] < mid) i++;//查找左半部分比中间数大的数 
		while(a[j] > mid) j--;//查找右半部分比中间数小的数 
		if(i <= j) {//如果有一组不满足排序条件（左小右大）的数 
			swap(a[i], a[j]);//交换 
			i++;
			j--;
		}
	}while(i <= j);//这里注意要有= 
	if(l < j) qsort(l, j);//递归搜索左半部分 
	if(i < r) qsort(i, r);//递归搜索右半部分 
}

7、堆排序（Heap Sort）

堆排序（Heapsort）是指利用堆这种数据结构所设计的一种排序算法。堆积是一个近似完全二叉树的结构，并同时满足堆积的性质：即子结点的键值或索引总是小于（或者大于）它的父节点。

7.1 算法描述

• 将初始待排序关键字序列\((R1,R2….Rn)\)构建成大顶堆，此堆为初始的无序区；
• 将堆顶元素\(R[1]\)与最后一个元素\(R[n]\)交换，此时得到新的无序区\((R1,R2,……Rn-1)\)和新的有序区\((Rn)\),且满足\(R[1,2…n-1]<=R[n]；\)
• 由于交换后新的堆顶\(R[1]\)可能违反堆的性质，因此需要对当前无序区\((R1,R2,……Rn-1)\)调整为新堆，然后再次将\(R[1]\)与无序区最后一个元素交换，得到新的无序区\((R1,R2….Rn-2)\)和新的有序区\((Rn-1,Rn)\)。不断重复此过程直到有序区的元素个数为\(n-1\)，则整个排序过程完成。

7.2 动图演示

7.3 代码实现

#include<iostream>
#include<vector>
using namespace std;
 
// 递归方式构建大根堆(len是arr的长度，index是第一个非叶子节点的下标)
void adjust(vector<int> &arr, int len, int index)
{
    int left = 2*index + 1; // index的左子节点
    int right = 2*index + 2;// index的右子节点
 
    int maxIdx = index;
    if(left<len && arr[left] > arr[maxIdx])     maxIdx = left;
    if(right<len && arr[right] > arr[maxIdx])     maxIdx = right;
 
    if(maxIdx != index)
    {
        swap(arr[maxIdx], arr[index]);
        adjust(arr, len, maxIdx);
    }
 
}
 
// 堆排序
void heapSort(vector<int> &arr, int size)
{
    // 构建大根堆（从最后一个非叶子节点向上）
    for(int i=size/2 - 1; i >= 0; i--)
    {
        adjust(arr, size, i);
    }
  // 调整大根堆
    for(int i = size - 1; i >= 1; i--)
    {
        swap(arr[0], arr[i]);           // 将当前最大的放置到数组末尾
        adjust(arr, i, 0);              // 将未完成排序的部分继续进行堆排序
    }
}
int main()
{
    vector<int> arr = {8, 1, 14, 3, 21, 5, 7, 10};
    heapSort(arr, arr.size());
    for(int i=0;i<arr.size();i++)
    {
        cout<<arr[i]<<endl;
    }
    return 0;
}

8、桶排序（Bucket Sort）

桶排序是计数排序的升级版。它利用了函数的映射关系，高效与否的关键就在于这个映射函数的确定。桶排序 (Bucket sort)的工作的原理：假设输入数据服从均匀分布，将数据分到有限数量的桶里，每个桶再分别排序（有可能再使用别的排序算法或是以递归方式继续使用桶排序进行排）。

8.1 算法描述

• 设置一个定量的数组当作空桶；
• 遍历输入数据，并且把数据一个一个放到对应的桶里去；
• 对每个不是空的桶进行排序；
• 从不是空的桶里把排好序的数据拼接起来。

8.2 图片演示

当然，如上图，桶也可以有范围。

8.3 代码实现

#include <bits/stdc++.h>
using namespace std;
int n, a[1005], c[1005]; //1≤a[i]≤10004; 
int main() {
	cin >> n;
	for(int i = 1; i <= n; i++) {
		cin >> a[i];
		c[a[i]]++;
	}
	for(int i = 1; i <= 1005; i++) {
		for(int j = 1; j <= c[i]; j++) cout << i << " ";
	}
	cout << endl;
	return 0;
}

这份代码的桶是以单个数为一个桶。

8.4 算法分析

桶排序最好情况下使用线性时间\(O(n)\)，桶排序的时间复杂度，取决与对各个桶之间数据进行排序的时间复杂度，因为其它部分的时间复杂度都为\(O(n)\)。很显然，桶划分的越小，各个桶之间的数据越少，排序所用的时间也会越少。但相应的空间消耗就会增大。

8.5 计数排序

因为计数排序与桶排序较为相似，所以就没放计数排序，这里有一张计数排序的动图演示

9、基数排序（Radix Sort）

基数排序是按照低位先排序，然后收集；再按照高位排序，然后再收集；依次类推，直到最高位。有时候有些属性是有优先级顺序的，先按低优先级排序，再按高优先级排序。最后的次序就是高优先级高的在前，高优先级相同的低优先级高的在前。

9.1 算法描述

• 取得数组中的最大数，并取得位数；
• arr为原始数组，从最低位开始取每个位组成radix数组；
• 对radix进行计数排序（利用计数排序适用于小范围数的特点）；

9.2 动图演示

9.3 代码实现

int maxbit(int data[], int n) {//辅助函数，求数据的最大位数 
	int d = 1;
	int p = 10;
	for(int i = 0; i < n; i++) {
		while(data[i] >= p) {
			p *= 10;
			d++;
		}
	}
	return d;
}
void radixsort(int data[], int n) {//基数排序 
	int d = maxbit(data, n);
	int *tmp = new int[n];
	int *count = new int[10];//计数器 
	int i, j, k;
	int radix = 1;
	for(i = 1; i <= d; i++) { //进行d次排序 
	    for(j = 0; j < 10; j++) count[j] = 0;
	    for(j = 0; j < n; j++) {
	    	k = (data[j] / radix) % 10;
	    	count[k]++;
		}
		for(j = 1; j < 10; j++)
		    count[j] = count[j - 1] + count[j];
		for(j = n - 1; j >= 0; j--) {
			k = (data[j] / radix) % 10;
			tmp[count[k] - 1] = data[j];
			count[k]--;
		}
		for(j = 0; j < n; j++) data[j] = tmp[j];
		radix = radix * 10;
	}
	delete[]tmp;
	delete[]count;
}

9.4 算法分析

基数排序基于分别排序，分别收集，所以是稳定的。但基数排序的性能比桶排序要略差，每一次关键字的桶分配都需要\(O(n)\)的时间复杂度，而且分配之后得到新的关键字序列又需要\(O(n)\)的时间复杂度。假如待排数据可以分为\(d\)个关键字，则基数排序的时间复杂度将是\(O(d \times 2n)\)，当然\(d\)要远远小于\(n\)，因此基本上还是线性级别的。

基数排序的空间复杂度为\(O(n+k)\)，其中k为桶的数量。一般来说\(n>>k\)，因此额外空间需要大概n个左右。

10、sort

没想到吧！sort！STL万岁！

void sortsort(int l, int r) {
    sort(a + l, a + r + 1);
}

sort是STL的函数，在algorithm头文件中，它可以排序任意数组，int数组排序时间复杂度为\(O(nlog_2n)\)

为何什么说任意数组呢？如果数据类型是结构体呢？

struct stu{
    int id, score;
}student[100];

当我们用到结构体时，用sort必须重载结构体比较大小的运算符operator，但还有一种方法，就是用另一个函数支持结构体的大小比较，代码如下。

bool cmp(stu x, stu y) {
    if(x.score == y.score) {
        return x.id < y.id;
    }
    return x.score > y.score;
}

接下来就是用sort排序啦！

sort(student + 1, student + n + 1, cmp);

11、参考文献

想了解更多资料？戳我！