普通的二叉树是不适合用数组来存储的,因为可能会存在大量的空间浪费。而完全二叉树更适合使用顺序结构存储。现实中我们通常把堆(一种二叉树)使用顺序结构的数组来存储,需要注意的是这里的堆和操作系统虚拟进程地址空间中的堆是两回事,一个是数据结构,一个是操作系统中管理内存的一块区域分段。
小堆:将根结点最小的堆叫做小堆,也叫最小堆或小根堆。
大堆:将根结点最大的堆叫做大堆,也叫最大堆或大根堆。
堆的性质:
现在我们给出一个数组,逻辑上看作一棵完全二叉树。我们通过从根节点开始的向下调整算法可以把它调整成一个小堆。
但是,使用向下调整算法需要满足一个前提:
向下调整算法的基本思想(以建小堆为例):
代码如下:
//交换函数
void Swap(int* x, int* y){
int tmp = *x;
*x = *y;
*y = tmp;
}
//堆的向下调整(小堆)
void AdjustDown(int *a, int n, int parent) {
//child记录左右孩子中值较小的孩子的下标
int child = 2 * parent + 1;//先默认其左孩子的值较小
while (child < n) {
if (child + 1 < n && a[child + 1] < a[child]){//右孩子存在并且右孩子比左孩子还小
child++;//较小的孩子改为右孩子
}
if (a[child] < a[parent]){//左右孩子中较小孩子的值比父结点还小
//将父结点与较小的子结点交换
Swap(&a[child], &a[parent]);
//继续向下进行调整
parent = child;
child = 2 * parent + 1;
} else{//已成堆
break;
}
}
}
使用堆的向下调整算法,最坏的情况下(即一直需要交换结点),需要循环的次数为:h - 1次(h为树的高度)。而h = log2(N+1)(N为树的总结点数)。所以堆的向下调整算法的时间复杂度为:O(logN) 。
上面说到,使用堆的向下调整算法需要满足其根结点的左右子树均为大堆或是小堆才行,那么如何才能将一个任意树调整为堆?
答案很简单,我们只需要从倒数第一个非叶子结点开始,从后往前,按下标,依次作为根去向下调整即可。
代码:
for (int i = (n - 1 - 1) / 2; i >= 0; i--) {
AdjustDown(php->a, n, i);
}
那么建堆的时间复杂度又是多少呢?
当结点数无穷大时,完全二叉树与其层数相同的满二叉树相比较来说,它们相差的结点数可以忽略不计,所以计算时间复杂度的时候我们可以将完全二叉树看作与其层数相同的满二叉树来进行计算。
总结一下:
当我们在一个堆的末尾插入一个数据后,需要对堆进行调整,使其仍然是一个堆,这时需要用到堆的向上调整算法。
向上调整算法的基本思想(以建小堆为例):
代码如下:
//交换函数
void Swap(HPDataType *x, HPDataType *y) {
HPDataType tmp = *x;
*x = *y;
*y = tmp;
}
//堆的向上调整(小堆)
void AdjustUp(HPDataType *a, int child) {
int parent = (child - 1) / 2;
while (child > 0) { //调整到根结点的位置截止
if (a[child] < a[parent]) {//孩子结点的值小于父结点的值
//将父结点与孩子结点交换
Swap(&a[child], &a[parent]);
//继续向上进行调整
child = parent;
parent = (child - 1) / 2;
} else {//已成堆
break;
}
}
}
首先,必须创建一个堆类型,该类型中需包含堆的基本信息:存储数据的数组、堆中元素的个数以及当前堆的最大容量。
typedef int HPDataType;
// 堆的结构 - 顺序表
typedef struct Heap {
HPDataType *a;
int size;
int capacity;
} Heap;
然后我们需要一个建堆函数,对刚创建的堆进行初始化,注意在初始化期间要将传入数据建堆。
// 堆的创建 - 建堆算法
void HeapCreate(Heap *php, HPDataType *a, int n) {
assert(php);
php->a = (HPDataType *) realloc(php->a, sizeof(HPDataType) * n);
if (php->a == NULL) {
perror("realloc fail");
exit(-1);
}
// 内存复制 数组复制到php->a中
memcpy(php->a, a, sizeof(HPDataType) * n);
php->size = php->capacity = n;
// 建堆算法 - 从最后一个叶子节点的父亲节点开始向下调整,然后从父亲节点的前所有节点都向下调整一次
// 最后节点的父亲的坐标是(n-1-1)/2 n-1因为n是节点个数,坐标从0开始
for (int i = (n - 1 - 1) / 2; i >= 0; i--) {
AdjustDown(php->a, n, i);
}
}
为了避免内存泄漏,使用完动态开辟的内存空间后都要及时释放该空间,所以,一个用于释放内存空间的函数是必不可少的。
// 堆的销毁
void HeapDestroy(Heap *php) {
assert(php);
free(php->a);
php->a = NULL;
php->size = php->capacity = 0;
}
// 堆的打印
void HeapPrint(Heap *php) {
assert(php);
for (int i = 0; i < php->size; i++) {
printf("%d ", php->a[i]);
}
printf("\n");
}
数据插入时是插入到数组的末尾,即树形结构的最后一层的最后一个结点,所以插入数据后我们需要运用堆的向上调整算法对堆进行调整,使其在插入数据后仍然保持堆的结构。
// 堆的插入
void HeapPush(Heap *php, HPDataType x) {
assert(php);
// 内存满了,扩容
if (php->size == php->capacity) {
int newcapacity = php->capacity == 0 ? 4 : php->capacity * 2;
HPDataType *tmp = (HPDataType *) realloc(php->a, sizeof(HPDataType) * newcapacity);
if (tmp == NULL) {
perror("realloc fail");
exit(-1);
}
php->a = tmp;
php->capacity = newcapacity;
}
// 插入
php->a[php->size] = x;
php->size++;
// 插入完,开始向上调整建堆
// 将数组和child的位置传过去
AdjustUp(php->a, php->size - 1);
}
堆的删除,删除的是堆顶的元素,但是这个删除过程可并不是直接删除堆顶的数据,而是先将堆顶的数据与最后一个结点的位置交换,然后再删除最后一个结点,再对堆进行一次向下调整。
原因:我们若是直接删除堆顶的数据,那么原堆后面数据的父子关系就全部打乱了,需要全体重新建堆,时间复杂度为O(N) 。若是用上述方法,那么只需要对堆进行一次向下调整即可,因为此时根结点的左右子树都是小堆,我们只需要在根结点处进行一次向下调整即可,时间复杂度为O(log N)
// 堆的删除
void HeapPop(Heap *php) {
assert(php);
assert(php->size > 0);
// 堆的删除,因为不能破坏堆的结构,所以将堆顶,和堆底最后一个数据交换,然后删除堆底,堆顶再向下调整,保持堆的结构
// 1.交换堆顶和堆底
Swap(&php->a[0], &php->a[php->size - 1]);
// 2.删除堆底
php->size--;
// 3.向下调整
AdjustDown(php->a, php->size, 0);
}
获取堆顶的数据,即返回数组下标为0的数据。
// 取堆顶
HPDataType HeapTop(Heap *php) {
assert(php);
assert(php->size > 0);
return php->a[0];
}
获取堆的长度,即返回堆结构体中的size变量。
// 求堆的长度
size_t HeapSize(Heap *php) {
assert(php);
return php->size;
}
堆的判空,即判断堆结构体中的size变量是否为0。
// 堆的判空
bool HeapEmpty(Heap *php) {
assert(php);
return php->size == 0;
}
#pragma once
#include
#include
#include
#include
#include
typedef int HPDataType;
// 堆的结构 - 顺序表
typedef struct Heap {
HPDataType *a;
int size;
int capacity;
} Heap;
// 堆的初始化
void HeapInit(Heap *php) {
assert(php);
php->a = NULL;
php->size = php->capacity = 0;
}
// 堆的打印
void HeapPrint(Heap *php) {
assert(php);
for (int i = 0; i < php->size; i++) {
printf("%d ", php->a[i]);
}
printf("\n");
}
// 堆的销毁
void HeapDestroy(Heap *php) {
assert(php);
free(php->a);
php->a = NULL;
php->size = php->capacity = 0;
}
// 向上调整 - 大堆
void AdjustUp(HPDataType *a, int child) {
// 1.找父亲
int parent = (child - 1) / 2;
// 2.跟父亲比大小,如果是大堆,知道父亲大于孩子循环结束,如果一直小于孩子,一直交换,然后循环结束条件是child==0
while (child > 0) {
// 孩子大于父亲则交换
if (a[child] > a[parent]) {
Swap(&a[child], &a[parent]);
child = parent;
parent = (child - 1) / 2;
} else {
break;
}
}
}
// 向下调整 - 大堆
void AdjustDown(HPDataType *a, int n, int parent) {
// 如果是大堆,先找父亲的孩子中的大的,然后跟他交换
// 先假设左孩子是大的,如果不是,重新设置为右孩子是大的
int child = parent * 2 + 1;
// child的值不会越界,所以循环条件是child < n
while (child < n) {
// 重新设置最大的孩子,如果右孩子大,就++child。特殊情况:最后的节点,只有左孩子,没有右孩子,所以还要加条判断,左孩子+1
if (child + 1 < n && a[child] < a[child + 1]) {
child++;
}
// 1.父亲小于孩子,交换,继续向下调整
// 2.父亲大于孩子,跳出
if (a[parent] < a[child]) {
Swap(&a[parent], &a[child]);
// 交换后,重新设置parent,找下一个位置开始向下调整
parent = child;
child = parent * 2 + 1;
} else {
break;
}
}
}
// 堆的插入
void HeapPush(Heap *php, HPDataType x) {
assert(php);
// 内存满了,扩容
if (php->size == php->capacity) {
int newcapacity = php->capacity == 0 ? 4 : php->capacity * 2;
HPDataType *tmp = (HPDataType *) realloc(php->a, sizeof(HPDataType) * newcapacity);
if (tmp == NULL) {
perror("realloc fail");
exit(-1);
}
php->a = tmp;
php->capacity = newcapacity;
}
// 插入
php->a[php->size] = x;
php->size++;
// 插入完,开始向上调整建堆
// 将数组和child的位置传过去
AdjustUp(php->a, php->size - 1);
}
// 堆的删除
void HeapPop(Heap *php) {
assert(php);
assert(php->size > 0);
// 堆的删除,因为不能破坏堆的结构,所以将堆顶,和堆底最后一个数据交换,然后删除堆底,堆顶再向下调整,保持堆的结构
// 1.交换堆顶和堆底
Swap(&php->a[0], &php->a[php->size - 1]);
// 2.删除堆底
php->size--;
// 3.向下调整
AdjustDown(php->a, php->size, 0);
}
// 取堆顶
HPDataType HeapTop(Heap *php) {
assert(php);
assert(php->size > 0);
return php->a[0];
}
// 求堆的长度
size_t HeapSize(Heap *php) {
assert(php);
return php->size;
}
// 堆的判空
bool HeapEmpty(Heap *php) {
assert(php);
return php->size == 0;
}
// 交换
void Swap(HPDataType *p1, HPDataType *p2) {
HPDataType tmp = *p1;
*p1 = *p2;
*p2 = tmp;
}
// 堆的创建 - 建堆算法
void HeapCreate(Heap *php, HPDataType *a, int n) {
assert(php);
php->a = (HPDataType *) realloc(php->a, sizeof(HPDataType) * n);
if (php->a == NULL) {
perror("realloc fail");
exit(-1);
}
// 内存复制 数组复制到php->a中
memcpy(php->a, a, sizeof(HPDataType) * n);
php->size = php->capacity = n;
// 建堆算法 - 从最后一个叶子节点的父亲节点开始向下调整,然后从父亲节点的前所有节点都向下调整一次
// 最后节点的父亲的坐标是(n-1-1)/2 n-1因为n是节点个数,坐标从0开始
for (int i = (n - 1 - 1) / 2; i >= 0; i--) {
AdjustDown(php->a, n, i);
}
}
堆排序即利用堆的思想来进行排序,总共分为两个步骤:
代码如下:
#include
// 交换
void Swap(int *p1, int *p2) {
int tmp = *p1;
*p1 = *p2;
*p2 = tmp;
}
// 向下调整 - 大堆
void AdjustDown(int *a, int n, int parent) {
// 如果是大堆,先找父亲的孩子中的大的,然后跟他交换
// 先假设左孩子是大的,如果不是,重新设置为右孩子是大的
int child = parent * 2 + 1;
// child的值不会越界,所以循环条件是child < n
while (child < n) {
// 重新设置最大的孩子,如果右孩子大,就++child。特殊情况:最后的节点,只有左孩子,没有右孩子,所以还要加条判断,左孩子+1
if (child + 1 < n && a[child] < a[child + 1]) {
child++;
}
// 1.父亲小于孩子,交换,继续向下调整
// 2.父亲大于孩子,跳出
if (a[parent] < a[child]) {
Swap(&a[parent], &a[child]);
// 交换后,重新设置parent,找下一个位置开始向下调整
parent = child;
child = parent * 2 + 1;
} else {
break;
}
}
}
// 对数组进行堆排序
void HeapSort(int *a, int n) {
// 向上调整建堆 -- N*logN
/*for (int i = 1; i < n; ++i)
{
AdjustUp(a, i);
}*/
// 向下调整建堆 - 大堆 O(N) _
for (int i = (n - 1 - 1) / 2; i >= 0; i--) {
AdjustDown(a, n, i);
}
int end = n - 1;
while (end > 0) {
// 第一个和最后一个交换,除了最后一个,剩下的进行建堆调整,把最大的调整到堆顶
Swap(&a[0], &a[end]);
// end为坐标,坐标比个数多一个,所以下面的end是剩余的个数
AdjustDown(a, end, 0);
end--;
}
}
int main() {
int arr[10] = {32, 43, 56, 76, 84, 12, 45, 67, 43, 37};
HeapSort(arr, sizeof(arr) / sizeof(int));
for (int i = 0; i < sizeof(arr) / sizeof(int); i++) {
printf("%d ", arr[i]);
}
}
TOP-K问题:即求数据结合中前K个最大的元素或者最小的元素,一般情况下数据量都比较大。
比如:专业前10名、世界500强、富豪榜、游戏中前100的活跃玩家等。
对于Top-K问题,能想到的最简单直接的方式就是排序,但是:如果数据量非常大,排序就不太可取了(可能数据都不能一下子全部加载到内存中)。最佳的方式就是用堆来解决,基本思路如下:
将剩余N-K个元素依次与堆顶元素比完之后,堆中剩余的K个元素就是所求的前K个最小或者最大的元素。
//交换函数
void Swap(int *x, int *y) {
int tmp = *x;
*x = *y;
*y = tmp;
}
//堆的向下调整(小堆)
void AdjustDown(int *a, int n, int parent) {
//child记录左右孩子中值较小的孩子的下标
int child = 2 * parent + 1;//先默认其左孩子的值较小
while (child < n) {
if (child + 1 < n && a[child + 1] < a[child]) {//右孩子存在并且右孩子比左孩子还小
child++; //较小的孩子改为右孩子
}
if (a[child] < a[parent]) {//左右孩子中较小孩子的值比父结点还小
//将父结点与较小的子结点交换
Swap(&a[child], &a[parent]);
//继续向下进行调整
parent = child;
child = 2 * parent + 1;
} else {//已成堆
break;
}
}
}
int *getLeastNumbers(int *arr, int arrSize, int k, int *returnSize) {
*returnSize = k;
if (k == 0)
return NULL;
//用数组的前K个数建小堆
int i = 0;
int *retArr = (int *) malloc(sizeof(int) * k);
for (i = 0; i < k; i++) {
retArr[i] = arr[i];
}
for (i = (k - 1 - 1) / 2; i >= 0; i--) {
AdjustDown(retArr, k, i);
}
//剩下的N-k个数依次与堆顶数据比较
for (i = k; i < arrSize; i++) {
if (arr[i] > retArr[0]) {
retArr[0] = arr[i];//堆顶数据替换
}
AdjustDown(retArr, k, 0);//进行一次向下调整
}
return retArr;//返回最大的k个数
}
时间复杂度:O(k + N * log k)空间复杂度:O(k)