一、希尔排序算法
基本思想
希尔排序通过选择一个间隔序列来将原始数组分成若干个子数组,对这些子数组进行插入排序,然后再逐步减少间隔直至1,最终完成整个数组的排序。
步骤
- 选择间隔 :首先选择一个初始间隔
gap
,通常是数组长度的一半。 - 分组并排序 :将数组按照
gap
分隔成若干个子数组,对这些子数组进行插入排序。 - 缩小间隔 :每次迭代都将
gap
缩小一半,重复步骤2,直到gap
变为1。
示例
假设我们有一个数组 [5, 2, 4, 6, 1, 3]
:
- 选择间隔 :
gap = 3
(数组长度为6,取一半) - 分组并排序 :
[5, 4]
->[4, 5]
[2, 1]
->[1, 2]
[6, 3]
->[3, 6]
- 经过第一轮排序后数组变为
[4, 1, 5, 3, 2, 6]
- 缩小间隔 :
gap = 1
(数组长度为6,取一半) - 分组并排序 :
- 进行插入排序,得到最终排序结果
[1, 2, 3, 4, 5, 6]
- 进行插入排序,得到最终排序结果
性能分析
- 时间复杂度:平均情况下接近 O(n^(3/2)),具体取决于所选的间隔序列;最坏情况为 O(n^2)。
- 空间复杂度:O(1)(原地排序)。
二、代码
cpp
#include <stdlib.h>
#include <stdio.h>
#include <time.h>
// 函数声明
int* create_and_generate_random_array(int size);
void print_array(int *array, int size);
void shell_sort(int *array, int size);
int generate_random_size();
int main() {
int size = generate_random_size(); // 随机生成数组大小
int *array = create_and_generate_random_array(size);
if (array == NULL) {
// 如果内存分配失败
printf("Memory allocation failed\n");
free(array);
return 1;
}
// 打印原始数组(如果需要,可以取消注释)
// printf("Original array:\n");
// print_array(array, size);
// 获取开始时间
clock_t start_time = clock();
// 对数组进行希尔排序
shell_sort(array, size);
// 获取结束时间
clock_t end_time = clock();
// 计算时间差并转换为毫秒
double execution_time = ((double)(end_time - start_time) / CLOCKS_PER_SEC) * 1000;
// 打印排序后的数组(如果需要,可以取消注释)
// printf("Sorted array:\n");
// print_array(array, size);
printf("array_size = %d\n", size);
// 打印执行时间
printf("Execution time: %.2f ms\n", execution_time);
// 释放分配的内存
free(array);
return 0;
}
// 生成随机数组大小
int generate_random_size() {
srand(time(NULL));
return rand() % 9000 + 1000; // 生成1000到9999之间的随机数
}
// 创建并生成随机数组
int* create_and_generate_random_array(int size) {
int *array = (int *)malloc(sizeof(int) * size);
if (array == NULL) {
// 如果内存分配失败
return NULL;
}
// 使用当前时间作为随机数种子
srand(time(NULL));
for (int i = 0; i < size; i++) {
array[i] = rand() % 1000; // 生成0到999之间的随机数
}
return array;
}
// 打印数组
void print_array(int *array, int size) {
for (int i = 0; i < size; i++) {
printf("%d ", array[i]);
}
printf("\n");
}
// 希尔排序
void shell_sort(int *array, int size) {
int gap, i, j, temp;
for (gap = size / 2; gap > 0; gap /= 2) {
for (i = gap; i < size; i++) {
temp = array[i];
for (j = i; j >= gap && array[j - gap] > temp; j -= gap) {
array[j] = array[j - gap];
}
array[j] = temp;
}
}
}