• 数据结构之堆的结构与实现


    目录

    一、堆的概念及结构

    1.1堆的概念

     1.2堆的性质

    1.3堆的结构

    二、堆的实现

    2.1堆向下调整算法(父亲与孩子做比较)

     2.2堆的向上调整算法(孩子与父亲做比较)

    2.3堆的创建(向下建堆)

     2.4向下建堆的时间复杂度

    2.5堆的插入

    2.6堆的删除

    2.7堆的完整代码实现

    三、堆的应用

    3.1堆排序

    3.2TOP-K问题


     

    一、堆的概念及结构

    1.1堆的概念

     1.2堆的性质

    堆中某个节点的值总是不大于或不小于其父节点的值;

    堆总是一棵完全二叉树。

    1.3堆的结构

    二、堆的实现

    2.1堆向下调整算法(父亲与孩子做比较)

    我们给出一个数组,逻辑上看做一颗完全二叉树。我们通过从根节点开始的向下调整算法可以把它调整成一个小堆。向下调整算法有一个前提:左右子树必须是一个堆,才能调整。
    以下面图片为例:建小堆过程中父亲不断与较小的孩子交换

    用代码来实现:

    1. void AdjustDown(HPDataType* a, int n, int parent)//n是参与向下算法的元素的个数
    2. {
    3. int child = parent * 2 + 1;
    4. while (child < n)
    5. {
    6. //建小堆,找到两个孩子中较小的那一个
    7. if (child + 1 < n && a[child + 1] < a[child])
    8. {
    9. child++;
    10. }
    11. //如果父亲不比孩子大,就证明已经是小堆了,直接跳出循环;
    12. //如果比孩子大就一直交换
    13. if (a[child] < a[parent])
    14. {
    15. Swap(&a[child], &a[parent]);
    16. parent = child;
    17. child = parent * 2 + 1;
    18. }
    19. else
    20. break;
    21. }
    22. }

     2.2堆的向上调整算法(孩子与父亲做比较)

    代码实现如下:

    1. void AdjustUp(HPDataType* a, int child)
    2. {
    3. int parent = (child - 1) / 2;
    4. while (child > 0)
    5. {
    6. if (a[child] < a[parent])
    7. {
    8. Swap(&a[child], &a[parent]);
    9. child = parent;
    10. parent = (child - 1) / 2;
    11. }
    12. else
    13. break;
    14. }
    15. }

    2.3堆的创建(向下建堆)

    我们给出一个数组,这个数组逻辑上可以看做一颗完全二叉树,但是还不是一个堆,现在我们通过算法,把它构建成一个堆。根节点左右子树不是堆,我们怎么调整呢?这里我们从倒数的第一个非叶子节点的 子树开始调整(向下调整),一直调整到根节点的树,就可以调整成堆。
    假定有数组 int a [] = { 1 , 5 , 3 , 8 , 7 , 6 };

     2.4向下建堆的时间复杂度

    因为堆是完全二叉树,而满二叉树也是完全二叉树,此处为了简化使用满二叉树来证明 ( 时间复杂度本来看的就是近似值,多几个节点不影响最终结果 )

     因此:向下建堆的时间复杂度为O(N)

    既然谈到了向下建堆的时间复杂度,不妨就算一下向上建堆的时间复杂度:

     冲两张图中可以看到:向下调整建堆的效率略高于向上调整建堆的效率,所以我上面所讨论的也都是向下调整建堆的实现方法。

    2.5堆的插入

    先插入一个 10 到数组的尾上,再进行向上调整算法,直到满足堆。

     代码实现:

    1. void HeapPush(Heap* hp, HPDataType x)
    2. {
    3. assert(hp);
    4. //判满以及扩容
    5. if (hp->_capacity == hp->_size)
    6. {
    7. int newCapacity = hp->_capacity == 0 ? 4 : 2 * hp->_capacity;
    8. HPDataType* tmp = (HPDataType*)realloc(hp->_a, sizeof(HPDataType) * newCapacity);
    9. if (tmp == NULL)
    10. {
    11. perror("realloc fail");
    12. exit(-1);
    13. }
    14. hp->_a = tmp;
    15. hp->_capacity = newCapacity;
    16. }
    17. hp->_a[hp->_size] = x;
    18. hp->_size++;
    19. AdjustUp(hp->_a, hp->_size - 1);
    20. }

    2.6堆的删除

    删除堆是删除堆顶的数据,将堆顶的数据根最后一个数据一换,然后删除数组最后一个数据,再进行向下调整算法。

     代码实现:

    1. void HeapPop(Heap* hp)
    2. {
    3. assert(hp);
    4. assert(hp->_size > 0);
    5. Swap(&hp->_a[0], &hp->_a[hp->_size - 1]);
    6. hp->_size--;
    7. AdjustDown(hp->_a, hp->_size, 0);
    8. }

    2.7堆的完整代码实现

    1. //Heap.h
    2. #pragma once
    3. #include
    4. #include
    5. #include
    6. #include
    7. typedef int HPDataType;
    8. typedef struct Heap
    9. {
    10. HPDataType* _a;
    11. int _size;
    12. int _capacity;
    13. }Heap;
    14. //堆的初始化
    15. void HeapInit(Heap* hp);
    16. // 堆的构建
    17. void HeapCreate(Heap* hp, HPDataType* a, int n);
    18. //交换
    19. void Swap(HPDataType* a, HPDataType* b);
    20. //向上调整
    21. void AdjustUp(HPDataType* a, int child);
    22. //向下调整
    23. void AdjustDown(HPDataType* a, int n, int parent);
    24. //打印
    25. void HeapPrint(Heap* hp);
    26. // 堆的销毁
    27. void HeapDestory(Heap* hp);
    28. // 堆的插入
    29. void HeapPush(Heap* hp, HPDataType x);
    30. // 堆的删除
    31. void HeapPop(Heap* hp);
    32. // 取堆顶的数据
    33. HPDataType HeapTop(Heap* hp);
    34. // 堆的数据个数
    35. int HeapSize(Heap* hp);
    36. // 堆的判空
    37. int HeapEmpty(Heap* hp);
    1. //Heap.c
    2. #include "Heap.h"
    3. void HeapInit(Heap* hp)
    4. {
    5. assert(hp);
    6. hp->_a = NULL;
    7. hp->_capacity = 0;
    8. hp->_size = 0;
    9. }
    10. void HeapCreate(Heap* hp, HPDataType* a, int n)
    11. {
    12. assert(hp);
    13. assert(a);
    14. hp->_a = (HPDataType*)malloc(sizeof(HPDataType)*n);
    15. if (hp->_a == NULL)
    16. {
    17. perror("malloc fail");
    18. exit(-1);
    19. }
    20. hp->_capacity = n;
    21. hp->_size = n;
    22. memcpy(hp->_a, a, sizeof(HPDataType) * n);
    23. for (int i = 1; i < n; i++)
    24. {
    25. AdjustUp(hp->_a, i);
    26. }
    27. }
    28. void Swap(HPDataType* a, HPDataType* b)
    29. {
    30. HPDataType tmp = *a;
    31. *a = *b;
    32. *b = tmp;
    33. }
    34. void AdjustUp(HPDataType* a, int child)
    35. {
    36. int parent = (child - 1) / 2;
    37. while (child > 0)
    38. {
    39. if (a[child] < a[parent])
    40. {
    41. Swap(&a[child], &a[parent]);
    42. child = parent;
    43. parent = (child - 1) / 2;
    44. }
    45. else
    46. break;
    47. }
    48. }
    49. void AdjustDown(HPDataType* a, int n, int parent)//n是参与向下算法的元素的个数
    50. {
    51. int child = parent * 2 + 1;
    52. while (child < n)
    53. {
    54. //建小堆,找到两个孩子中较小的那一个
    55. if (child + 1 < n && a[child + 1] < a[child])
    56. {
    57. child++;
    58. }
    59. //如果父亲不比孩子大,就证明已经是小堆了,直接跳出循环;
    60. //如果比孩子大就一直交换
    61. if (a[child] < a[parent])
    62. {
    63. Swap(&a[child], &a[parent]);
    64. parent = child;
    65. child = parent * 2 + 1;
    66. }
    67. else
    68. break;
    69. }
    70. }
    71. void HeapDestory(Heap* hp)
    72. {
    73. assert(hp);
    74. free(hp->_a);
    75. hp->_capacity = 0;
    76. hp->_size = 0;
    77. }
    78. void HeapPush(Heap* hp, HPDataType x)
    79. {
    80. assert(hp);
    81. //判满以及扩容
    82. if (hp->_capacity == hp->_size)
    83. {
    84. int newCapacity = hp->_capacity == 0 ? 4 : 2 * hp->_capacity;
    85. HPDataType* tmp = (HPDataType*)realloc(hp->_a, sizeof(HPDataType) * newCapacity);
    86. if (tmp == NULL)
    87. {
    88. perror("realloc fail");
    89. exit(-1);
    90. }
    91. hp->_a = tmp;
    92. hp->_capacity = newCapacity;
    93. }
    94. hp->_a[hp->_size] = x;
    95. hp->_size++;
    96. AdjustUp(hp->_a, hp->_size - 1);
    97. }
    98. void HeapPrint(Heap* hp)
    99. {
    100. assert(hp);
    101. for (int i = 0; i < hp->_size; i++)
    102. {
    103. printf("%d ", hp->_a[i]);
    104. }
    105. printf("\n");
    106. }
    107. void HeapPop(Heap* hp)
    108. {
    109. assert(hp);
    110. assert(hp->_size > 0);
    111. Swap(&hp->_a[0], &hp->_a[hp->_size - 1]);
    112. hp->_size--;
    113. AdjustDown(hp->_a, hp->_size, 0);
    114. }
    115. HPDataType HeapTop(Heap* hp)
    116. {
    117. assert(hp);
    118. assert(hp->_size > 0);
    119. return hp->_a[0];
    120. }
    121. int HeapSize(Heap* hp)
    122. {
    123. return hp->_size;
    124. }
    125. int HeapEmpty(Heap* hp)
    126. {
    127. assert(hp);
    128. if (hp->_size == 0)
    129. return 0;
    130. else
    131. return 1;
    132. }

    三、堆的应用

    3.1堆排序

    堆排序即利用堆的思想来进行排序,总共分为两个步骤:
    1. 建堆:
    升序:建大堆,降序:建小堆。
    2. 利用堆删除思想来进行排序
    建堆和堆删除中都用到了向下调整,因此掌握了向下调整,就可以完成堆排序。
    具体实现代码如下:
    1. void HeapSort1(int* a, int n)
    2. {
    3. //向上调整建堆
    4. /*for (int i = 1; i < n; i++)
    5. {
    6. AdjustUp(a, i);
    7. }*/
    8. //向下调整建堆
    9. for (int i = (n - 1 - 1) / 2; i >= 0; i--)//从第一个非叶子节点开始向下调整
    10. {
    11. AdjustDown(a, n, i);
    12. }
    13. //排序
    14. int end = n - 1;
    15. while (end)
    16. {
    17. Swap(&a[0], &a[end]);
    18. AdjustDown(a, end, 0);
    19. end--;
    20. }
    21. }

    3.2TOP-K问题

    TOP-K问题:即求数据结合中前K个最大的元素或者最小的元素,一般情况下数据量都比较大

    对于 Top-K 问题,能想到的最简单直接的方式就是排序,但是:如果数据量非常大,排序就不太可取了 ( 可能数据都不能一下子全部加载到内存中 ) 。最佳的方式就是用堆来解决,基本思路如下:
    1. 用数据集合中前 K 个元素来建堆 :
    k 个最大的元素,则建小堆,前 k 个最小的元素,则建大堆 。
    2. 用剩余的 N-K 个元素依次与堆顶元素来比较,不满足则替换堆顶元素
    将剩余 N-K 个元素依次与堆顶元素比完之后,堆中剩余的 K 个元素就是所求的前 K 个最小或者最大的元素。
    具体实现代码如下:
    1. void CreatNData()
    2. {
    3. // 造数据
    4. int n = 10000000;
    5. srand(time(0));
    6. const char* file = "data.txt";
    7. FILE* fin = fopen(file, "w");
    8. if (fin == NULL)
    9. {
    10. perror("fopen error");
    11. return;
    12. }
    13. //将数据写入data文件中
    14. for (int i = 0; i < n; ++i)
    15. {
    16. int x = (rand() + i) % 10000000;
    17. fprintf(fin, "%d\n", x);
    18. }
    19. fclose(fin);
    20. }
    21. void PrintTopK(const char* filename, int k)
    22. {
    23. FILE* fout = fopen(filename, "r");
    24. if (fout == NULL)
    25. {
    26. perror("fopen fail");
    27. exit(-1);
    28. }
    29. int* minHeap = (int*)malloc(sizeof(int) * k);
    30. if (minHeap == NULL)
    31. {
    32. perror("malloc fail");
    33. return;
    34. }
    35. for (int i = 0; i < k; i++)
    36. {
    37. fscanf(fout, "%d", &minHeap[i]);
    38. }
    39. for (int i = (k - 1 - 1) / 2; i >= 0; i--)
    40. {
    41. AdjustDown(minHeap, k, i);
    42. }
    43. //将剩余的n-k各元素与堆顶的元素进行交换
    44. int x = 0;
    45. while (fscanf(fout, "%d", &x) != EOF)
    46. {
    47. if (x > minHeap[0])
    48. {
    49. minHeap[0] = x;
    50. AdjustDown(minHeap, k, 0);
    51. }
    52. }
    53. //排序
    54. int end = k - 1;
    55. while (end)
    56. {
    57. Swap(&minHeap[0], &minHeap[end]);
    58. AdjustDown(minHeap, end, 0);
    59. end--;
    60. }
    61. for (int i = 0; i < k; i++)
    62. {
    63. printf("%d ", minHeap[i]);
    64. }
    65. free(minHeap);
    66. fclose(fout);
    67. }

  • 相关阅读:
    Filter快速入门、Filter执行流程、Filter使用细节、Listener概念、分类、ServletContextListener使用
    LQ0203 排它平方数【枚举+进制】
    android wifi启动流程中驱动的确认
    OpenCV学习笔记(四)——对视频的读取操作
    GO语言最常用的语法
    关于TypeScript中的泛型,希望这次能让你彻底理解
    y137.第八章 Servless和Knative从入门到精通 -- Serverless概念和基础(一)
    【C++】vector的模拟实现不会怎么办?看过来
    论格斗游戏中“一滴血反杀”现象的心理机制与启示
    计算摄像技术02 - 颜色空间
  • 原文地址:https://blog.csdn.net/m0_74265792/article/details/132968573