堆结构的实现和应用
- 1. 堆结构的代码实现
- 2. 向上调整算法和向下调整算法时间复杂度
- 3. 堆排序
- 3.1 方法一
- 3.2 方法二
- 3.3 建堆时间复杂度
- 4. TOP-K问题
1. 堆结构的代码实现
首先,我们先定义一个堆的结构,这里和顺序表的结构一样
typedef int HPDataType;typedef struct Heap
{HPDataType* a;int size;int capacity;
}HP;
初始化堆
void HeapInit(HP* hp)
{hp->a = NULL;hp->capacity = 0;hp->size = 0;
}
销毁堆
void HeapDestroy(HP* hp)
{free(hp->a);hp->a = NULL;hp->size = 0;hp->capacity = 0;
}
因为堆的空间是动态开辟的,所以在堆使用结束后,需要把动态开辟的空间释放,避免造成内存泄漏。
堆的插入
在堆尾插入数据后并不能保证插入之后的数据是堆,所以,在插入之后,要进行向上调整算法,直到满足堆。
void HeapPush(HP* hp, HPDataType x)
{assert(hp);HeapCheckCapacity(hp);hp->a[hp->size] = x;hp->size++;HeapUp(hp->a, hp->size - 1);
}
这里也可以不把扩容的代码单独写成一个函数,因为这里只有这一个地方会用到扩容操作,但是还是建议封装成一个函数,一方面可以提高代码的可读性,另一方面平时我们写代码的时候这个操作一般都是封装成一个函数,所以这样写可以在一定程度上避免代码出错。
检查是否需要扩容
void HeapCheckCapacity(HP* hp)
{if (hp->size == hp->capacity){int newcapacity = hp->capacity == 0 ? 4 : hp->capacity * 2;HPDataType* ptr = (HPDataType*)realloc(hp->a, sizeof(HPDataType) * newcapacity);if (ptr == NULL){perror("realloc fail");return NULL;}hp->a = ptr;hp->capacity = newcapacity;}
}
向上调整算法
要求插入之前这个数组是大堆或者小堆,如果时小堆,当插入的孩子比父亲小就交换,直到孩子等于0,即没有父亲的时候交换结束。
算法思想(小堆):
- 将目标节点与父节点比较。
- 如果目标节点的值比父节点的值小,则交换目标节点与父节点的位置,并把父节点当作新的目标节点继续向上交换,当目标节点的值比父亲节点大,或者目标节点变成根节点,则停止交换,此时,数组就可以看作一个小堆了。
//假如是小堆
void AdjustUp(HPDataType* a, int child)
{while (child > 0)//调整到根节点处停止{int parent = (child - 1) / 2;if (a[child] < a[parent])//当孩子比父亲小就交换{HPDataType tmp = a[child];a[child] = a[parent];a[parent] = tmp;}else//如果孩子比父亲大就变成堆了,停止交换。{break;}}
}
堆的删除
删除堆是删除堆顶的数据,因为删除堆尾的数据意义不大
- 挪动覆盖删除堆顶元素,重新建堆。
挪动覆盖,不能保证还是堆,父子关系全乱了,只能重新建堆,代价太大。 - 将堆顶的数据根最后一个数据一换,然后删除数组最后一个数据,再进行向下调整算法。
void HeapPop(HP* hp)
{assert(hp);assert(!HeapEmpty(hp));Swap(&hp->a[0], &hp->a[hp->size - 1]);hp->size--;AdjustDown(hp->a,hp->size,0);
}
注意:堆为空的时候不能删除
向下调整算法
算法思想(小堆):
- 从根节点开始,选出左右孩子中值较小的孩子
- 如果较小的孩子的值比父节点的值小,则交换该孩子节点与父节点的位置,并把该孩子点当作新的父亲节点继续向下交换,当较小孩子节点的值比父亲节点大,或者孩子节点出了数据有效范围,则停止交换,此时,数组就可以看作一个小堆了。
向下调整算法有一个前提:左右子树必须是一个堆,才能调整。
代码如下:
void AdjustDown(HPDataType* a, int n, int parent)
{int child = parent * 2 + 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 = parent * 2 + 1;}else{break;//当较小的孩子比父亲大,则停止交换}}
}
交换元素
因为多处需要用到交换函数,所以我们单独封装成一个函数
void Swap(HPDataType* p1, HPDataType* p2)
{HPDataType tmp = *p1;*p1 = *p2;*p2 = tmp;
}
取出堆顶元素
HPDataType HeapTop(HP* hp)
{assert(hp);assert(!HeapEmpty(hp));return hp->a[0];
}
检测堆是否为空,如果为空返回非零结果,如果不为空返回0
bool HeapEmpty(HP* hp)
{return hp->size == 0;
}
获取堆中有效元素个数
int HeapSize(HP* hp)
{assert(hp);return hp->size;
}
全部代码实现
Heap.h
#include <stdio.h>
#include <stdlib.h>
#include <stdbool.h>
#include <assert.h>typedef int HPDataType;typedef struct Heap
{HPDataType* a;int size;int capacity;
}HP;void HeapInit(HP* hp);void HeapDestroy(HP* hp);void HeapPush(HP* hp,HPDataType x);void HeapPop(HP* hp);HPDataType HeapTop(HP * hp);bool HeapEmpty(HP* hp);int HeapSize(HP* hp);
Heap.c
#include "Heap.h"void HeapInit(HP* hp)
{assert(hp);hp->a = NULL;hp->capacity = 0;hp->size = 0;
}void HeapDestroy(HP* hp)
{assert(hp);free(hp->a);hp->a = NULL;hp->capacity = 0;hp->size = 0;
}void HeapCheckCapacity(HP* hp)
{assert(hp);if (hp->size == hp->capacity){int newcapacity = hp->capacity == 0 ? 4 : hp->capacity * 2;HPDataType* ptr = (HPDataType*)realloc(hp->a, sizeof(HPDataType) * newcapacity);if (ptr == NULL){perror("realloc fail");exit(-1);}hp->a = ptr;hp->capacity = newcapacity;}
}void Swap(HPDataType* p1, HPDataType* p2)
{HPDataType tmp = *p1;*p1 = *p2;*p2 = tmp;
}void AdjustUp(HPDataType* a, int child)
{while (child > 0){int parent = (child - 1) / 2;if (a[child] < a[parent]){Swap(&a[child], &a[parent]);}else{break;}child = parent;}
}void HeapPush(HP* hp, HPDataType x)
{assert(hp);HeapCheckCapacity(hp);hp->a[hp->size] = x;hp->size++;AdjustUp(hp->a, hp->size - 1);
}void AdjustDown(HPDataType* a, int n, int parent)
{int child = parent * 2 + 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 = parent * 2 + 1;}else{break;}}
}bool HeapEmpty(HP* hp)
{return hp->size == 0;
}void HeapPop(HP* hp)
{assert(hp);assert(!HeapEmpty(hp));Swap(&hp->a[0], &hp->a[hp->size - 1]);hp->size--;AdjustDown(hp->a,hp->size,0);
}HPDataType HeapTop(HP* hp)
{assert(hp);assert(!HeapEmpty(hp));return hp->a[0];
}int HeapSize(HP* hp)
{assert(hp);return hp->size;
}
Test.c
#include "Heap.h"int main()
{HP hp;HeapInit(&hp);int a[] = { 65,100,70,32,50,60 };int sz = sizeof(a) / sizeof(a[0]);int i = 0;for (i = 0; i < sz; i++){HeapPush(&hp, a[i]);}for (i = 0; i < sz; i++){printf("%d\n", HeapTop(&hp));HeapPop(&hp);}HeapDestroy(&hp);return 0;
}
运行结果如图:
这里只是把数据从小到大打印出来,并没有改变数组的内容,所以不是堆排序。
2. 向上调整算法和向下调整算法时间复杂度
假如树的高度为h,那么向上调整算法和向下调整算法,调整的次数最多为h-1次,也可以认为h次,因为,我们之前已经推出高度h和节点个数N的关系,那么我们很容易推导处下图的公式。
所以,向上调整算法和向下调整算法时间复杂度用大O的渐进法表示就是O(logN)。
3. 堆排序
3.1 方法一
#include "Heap.h"void HeapSort(int* a, int n)
{HP hp;HeapInit(&hp);int i = 0;for (i = 0; i < n; i++){HeapPush(&hp, a[i]);}i = 0;while (!HeapEmpty(&hp)){int top = HeapTop(&hp);a[i++] = top;HeapPop(&hp);}HeapDestroy(&hp);
}int main()
{int a[] = { 3,5,1,4 };int sz = sizeof(a) / sizeof(a[0]);HeapSort(a, sz);int i = 0;for (i = 0; i < sz; i++){printf("%d ", a[i]);}return 0;
}
运行结果如图:
这样虽然也可以实现堆排序,但是有弊端:必须先有一个堆,太麻烦,需要拷贝数据空间复制度比较高。
3.2 方法二
- 利用数组进行向上调整建堆或向上调整建堆的方法进行建堆。
- 小堆,选出最小的,首尾交换,最小的放到最后的位置,把最后一个数据,不看作堆里面的,向下调整后就可以选出次小的,最终,排出来的数据就是降序。
向上调整建堆:
void HeapSort(int* a, int n)
{//O(N * log(N))//建堆 -- 向上调整建堆int i = 0;for (i = 1; i < n; i++){// 升序 -- 建大堆// 降序 -- 建小堆AdjustUp(a, i);}O(N * logN)int end = n - 1;while (end > 0){Swap(&a[0], &a[end]);//在调整,选出次小的数AdjustDown(a, end, 0);end--;}
}int main()
{int a[] = { 3,5,1,4 };int sz = sizeof(a) / sizeof(a[0]);HeapSort(a, sz);int i = 0;for (i = 0; i < sz; i++){printf("%d ", a[i]);}return 0;
}
运行结果如图:
向下调整建堆
void HeapSort(int* a, int n)
{//O(N)int i = 0;for (i = (n - 1 - 1) / 2; i >= 0; i--){AdjustDown(a, n, i);}O(N * logN)int end = n - 1;while (end > 0){Swap(&a[0], &a[end]);AdjustDown(a, end, 0);end--;}
}int main()
{int a[] = { 3,5,1,4 };int sz = sizeof(a) / sizeof(a[0]);HeapSort(a, sz);int i = 0;for (i = 0; i < sz; i++){printf("%d ", a[i]);}return 0;
}
运行结果如图:
3.3 建堆时间复杂度
因为堆是完全二叉树,而满二叉树也是完全二叉树,此处为了简化使用满二叉树来证明(时间复杂度本来看的就是近似值,多几个节点不影响最终结果):
4. TOP-K问题
TOP-K问题:即求数据结合中前K个最大的元素或者最小的元素,一般情况下数据量都比较大。
比如:专业前10名、世界500强、富豪榜、游戏中前100的活跃玩家等。
对于Top-K问题,能想到的最简单直接的方式就是排序,但是:如果数据量非常大,排序就不太可取了(可能数据都不能一下子全部加载到内存中)。最佳的方式就是用堆来解决,基本思路如下:
- 用数据集合中前K个元素来建堆
- 前k个最大的元素,则建小堆
- 前k个最小的元素,则建大堆
- 用剩余的N-K个元素依次与堆顶元素来比较,不满足则替换堆顶元素
将剩余N-K个元素依次与堆顶元素比完之后,堆中剩余的K个元素就是所求的前K个最小或者最大的元素。
#include <stdio.h>
#include <stdlib.h>
#include <time.h>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;while (child < n){if (child + 1 < n && a[child] > a[child + 1]){child++;}if (a[child] < a[parent]){Swap(&a[child], &a[parent]);parent = child;child = parent * 2 + 1;}else{break;}}
}void CreateNDate()
{// 造数据int n = 100;srand(time(NULL));const char* file = "data.txt";FILE* fin = fopen(file, "w");if (fin == NULL){perror("fopen error");return;}for (size_t i = 0; i < n; ++i){int x = rand() % 100;fprintf(fin, "%d\n", x);}fclose(fin);
}void PrintTopK(int k)
{const char* file = "data.txt";FILE* fout = fopen(file, "r");if (fout == NULL){perror("fopen error");return;}int* kminheap = (int*)malloc(sizeof(int) * k);if (kminheap == NULL){perror("malloc error");return;}for (int i = 0; i < k; i++){fscanf(fout, "%d", &kminheap[i]);}// 建小堆for (int i = (k - 1 - 1) / 2; i >= 0; i--){AdjustDown(kminheap, k, i);}int val = 0;while (!feof(fout)){fscanf(fout, "%d", &val);if (val > kminheap[0]){kminheap[0] = val;AdjustDown(kminheap, k, 0);}}for (int i = 0; i < k; i++){printf("%d ", kminheap[i]);}printf("\n");
}int main()
{//CreateNDate();PrintTopK(5);return 0;
}
我们创造好文件之后可以,可以手动修改文件内容。
运行结果如图: