堆结构的实现和应用

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,即没有父亲的时候交换结束。

算法思想(小堆):

  1. 将目标节点与父节点比较。
  2. 如果目标节点的值比父节点的值小,则交换目标节点与父节点的位置,并把父节点当作新的目标节点继续向上交换,当目标节点的值比父亲节点大,或者目标节点变成根节点,则停止交换,此时,数组就可以看作一个小堆了。

在这里插入图片描述

//假如是小堆
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;
		}
	}
}

堆的删除

删除堆是删除堆顶的数据,因为删除堆尾的数据意义不大

  1. 挪动覆盖删除堆顶元素,重新建堆。
    挪动覆盖,不能保证还是堆,父子关系全乱了,只能重新建堆,代价太大。
  2. 将堆顶的数据根最后一个数据一换,然后删除数组最后一个数据,再进行向下调整算法。
    在这里插入图片描述
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);
}

注意:堆为空的时候不能删除

向下调整算法

算法思想(小堆):

  1. 从根节点开始,选出左右孩子中值较小的孩子
  2. 如果较小的孩子的值比父节点的值小,则交换该孩子节点与父节点的位置,并把该孩子点当作新的父亲节点继续向下交换,当较小孩子节点的值比父亲节点大,或者孩子节点出了数据有效范围,则停止交换,此时,数组就可以看作一个小堆了。

向下调整算法有一个前提:左右子树必须是一个堆,才能调整。
在这里插入图片描述

代码如下:

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 方法二

  1. 利用数组进行向上调整建堆或向上调整建堆的方法进行建堆。
  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问题,能想到的最简单直接的方式就是排序,但是:如果数据量非常大,排序就不太可取了(可能数据都不能一下子全部加载到内存中)。最佳的方式就是用堆来解决,基本思路如下:

  1. 用数据集合中前K个元素来建堆
    • 前k个最大的元素,则建小堆
    • 前k个最小的元素,则建大堆
  2. 用剩余的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;
}

我们创造好文件之后可以,可以手动修改文件内容。

在这里插入图片描述

运行结果如图:
在这里插入图片描述

  • 17
    点赞
  • 28
    收藏
    觉得还不错? 一键收藏
  • 2
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值