排序算法之分配排序
罗朝辉(http://www.cppblog.com/kesalin)
转载请注明出处
排序是数据处理中经常使用的一种重要运算,在计算机及其应用系统中,花费在排序上的时间在系统运行时间中占有很大比重,其重要性勿需多言。下文将介绍常用的如下排序方法,对它们进行简单的分析和比较,并提供 C 语言实现。
所谓排序,就是要将一堆记录,使之按关键字递增(或递减)次序排列起来。根据排序所采用的策略,可以分为如下五种:
1、插入排序(直接插入排序、希尔排序)
2、交换排序(冒泡排序、快速排序)
3、选择排序(直接选择排序、堆排序)
4、归并排序
5、桶排序(桶排序,基数排序);
---------------------------------------------------------------------------------
前面讲了插入排序,交换排序,选择排序,归并排序,下面接着来讲桶排序,基数排序。
桶排序和基数排序均属于分配排序。分配排序的基本思想:排序过程无须比较关键字,而是通过用额外的空间来"分配"和"收集"来实现排序,它们的时间复杂度可达到线性阶:O(n)。简言之就是:用空间换时间,所以性能与基于比较的排序才有数量级的提高!
桶排序(Bucket Sort),也称箱排序
基本思想:设置若干个箱子,依次扫描待排序的记录 array[0],array[1],…,array[n - 1],把关键字等于 k 的记录全都装入到第 k 个箱子里(分配),然后按序号依次将各非空的箱子里的记录收集起来,从而完成排序。
桶排序所需要的额外空间取决于关键字的个数,若 array[0..n - 1] 中关键字的取值范围是 0 到 m - 1 的整数,则必须设置 m 个箱子。因此箱排序要求关键字的类型是有限类型,否则可能要无限个箱子。一般情况下每个箱子中存放多少个关键字相同的记录是无法预料的,故箱子的类型应设计成链表为宜。
在实现上,桶排序一般采用另外的变种。即:把 [0, m) 划分为 m 个大小相同的子区间,每一子区间是一个桶。然后将 n 个记录分配到各个桶中。因为关键字序列是均匀分布在[0,m)上的,所以一般不会有很多个记录落入同一个桶中。由于同一桶中的记录其关键字不尽相同,所以必须采用关键字比较的排序方法(通常用插入排序)对各个桶进行排序,然后依次将各非空桶中的记录收集起来即可。
C 代码实现:
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif)
struct bucket_node
{
int key;
bucket_node* next;
};
![](http://www.cppblog.com/Images/OutliningIndicators/None.gif)
// 取得数组中最大数的位数
int get_max_digital_count(int* array, int length)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ContractedBlock.gif)
{
assert(array && length > 0);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
int i = 0;
int max = array[0];
int maxDigitalCount = 1;
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 1; i < length; i++)
{
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (max < array[i])
{
max = array[i];
}
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
while ((max / 10) > 0)
{
max %= 10;
++maxDigitalCount;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
return maxDigitalCount;
}
![](http://www.cppblog.com/Images/OutliningIndicators/None.gif)
// 取得数 num 中从低到高第 n 位上的数字
int get_ditital_at(int num, int n)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ContractedBlock.gif)
{
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
while (--n > 0)
{
num /= 10;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
return (num % 10);
}
![](http://www.cppblog.com/Images/OutliningIndicators/None.gif)
// 箱/桶排序
//
void bucket_sort(int* array, int length)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ContractedBlock.gif)
{
assert(array && length >= 0);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (length <= 1)
{
return;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
int i, index;
bucket_node* temp = NULL;
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
bucket_node bucket[10] =
{0, }; // 根据数字个数 0 ~ 9 建立 10 个桶
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
int count = get_max_digital_count(array, length);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 建立数据节点
bucket_node* data = (bucket_node*)malloc(length * sizeof(bucket_node));
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (!data)
{
printf("Error: out of memory!\n");
return;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 0; i < length; i++)
{
data[i].key = array[i];
data[i].next = NULL;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 分配
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 0; i < length; i++)
{
index = get_ditital_at(data[i].key, count);
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (bucket[index].next == NULL)
{
bucket[index].next = &data[i];
}
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
else
{
temp = &bucket[index];
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
while (temp->next != NULL && temp->next->key < data[i].key)
{
temp = temp->next;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
data[i].next = temp->next;
temp->next = &data[i];
}
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 收集
index = 0;
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 0; i < 10; i++)
{
temp = bucket[i].next;
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
while (temp != NULL)
{
array[index++] = temp->key;
temp = temp->next;
}
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
free(data);
}
时间复杂度分析:
桶排序的平均时间复杂度是线性的,即 O(n)。但最坏情况仍有可能是 O(n ^ 2)。
空间复杂度分析:
桶排序只适用于关键字取值范围较小的情况,否则所需箱子的数目 m 太多导致浪费存储空间和计算时间。
基数排序(Radix Sort)
基本思想:基数排序是对桶排序的改进和推广。如果说桶排序是一维的基于桶的排序,那么基数排序就是多维的基于桶的排序。我这么说,可能还不是太清楚。比方说:用桶排序对 [0, 30] 之间的数进行排序,那么需要 31 个桶,分配一次,收集一次,完成排序;那么基数排序则只需要 0 - 9 总共 10 个桶(即关键字为数字 0 - 9),依次进行个位和十位的分配和收集从而完成排序。
C 代码实现:
// 基数排序
//
void radix_sort(int* array, int length)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedBlockStart.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ContractedBlock.gif)
{
assert(array && length >= 0);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (length <= 1)
{
return;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
const int buffer_size = length * sizeof(int);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
int i, k, count, index;
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
int bucket[10] =
{0, }; // 根据数字个数 0 ~ 9 建立 10 个桶
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
int* temp = (int*)malloc(buffer_size);
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
if (!temp)
{
printf("Error: out of memory!\n");
return;
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
count = get_max_digital_count(array, length);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (k = 1; k <= count; ++k)
{
memset(bucket, 0, 10 * sizeof(int));
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 统计各桶中元素的个数
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 0; i < length; ++i)
{
index = get_ditital_at(array[i], k);
++bucket[index];
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 为每个记录创建索引下标
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 1; i < 10; ++i)
{
bucket[i] += bucket[i - 1];
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 按索引下标顺序排列
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = length - 1; i >= 0; --i)
{
index = get_ditital_at(array[i], k);
assert(bucket[index] - 1 >= 0);
temp[--bucket[index]] = array[i];
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
// 一趟桶排序完毕,拷贝结果
memcpy(array, temp, buffer_size);
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
#ifdef DEBUG_SORT
debug_print(" 第 %d 趟排序:", k);
![](http://www.cppblog.com/Images/OutliningIndicators/ExpandedSubBlockStart.gif)
for (i = 0; i < length; ++i)
{
debug_print("%d ", array[i]);
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
debug_print("\n");
#endif
}
![](http://www.cppblog.com/Images/OutliningIndicators/InBlock.gif)
free(temp);
}
![](http://www.cppblog.com/Images/OutliningIndicators/None.gif)
时间复杂度分析:
基数排序的时间负责度为 O(n)。
空间复杂度:
基数排序所需的辅助存储空间为 O(n + r * d),其中 r 为记录中关键字分量的最大个数,d 为关键字的个数。比如说:待排序为 0 - 999,那么分量的最大个数为 3,关键字的个数为 10(0 - 9)。
补充:
基数排序是稳定的。
若排序文件不是以数组 array 形式给出,而是以单链表形式给出(此时称为链式的基数排序),则可通过修改出队和入队函数使表示箱子的链队列无须分配结点空间,而使用原链表的结点空间。人队出队操作亦无需移动记录而仅需修改指针。虽然这样一来节省了一定的时间和空间,但算法要复杂得多,且时空复杂度就其数量级而言并未得到改观。
=======================================================================================
测试:
在前文《排序算法之插入排序》测试代码的基础上添加两行代码即可:
{"桶/箱排序", bucket_sort},
{"基数排序", radix_sort},
测试结果:
=== 桶/箱排序 ===
original: 65 32 49 10 8 72 27 42 18 58 91
sorted: 8 10 18 27 32 42 49 58 65 72 91
original: 10 9 8 7 6 5 4 3 2 1 0
sorted: 0 1 2 3 4 5 6 7 8 9 10
=== 基数排序 ===
original: 65 32 49 10 8 72 27 42 18 58 91
第 1 趟排序:10 91 32 72 42 65 27 8 18 58 49
第 2 趟排序:8 10 18 27 32 42 49 58 65 72 91
sorted: 8 10 18 27 32 42 49 58 65 72 91
original: 10 9 8 7 6 5 4 3 2 1 0
第 1 趟排序:10 0 1 2 3 4 5 6 7 8 9
第 2 趟排序:0 1 2 3 4 5 6 7 8 9 10
sorted: 0 1 2 3 4 5 6 7 8 9 10