模拟堆——模板

关于堆的操作：

1. 插入一个数：heap[++ size] = x; up(size);
2. 求集合当中的最小值：heap[1]
3. 删除最小值：heap[1] = heap[size]; size--; down(1);
4. 删除任意一个元素（stl 实现不了）：heap[k] = heap[size]; size--; down(k); up(k);
5. 修改任意一个元素 （stl 实现不了）：heap[k] = x; down(k); up(k);

// h[N]存储堆中的值, h[1]是堆顶，x的左儿子是2x, 右儿子是2x + 1
// ph[k]存储第k个插入的点在堆中的位置
// hp[k]存储堆中下标是k的点是第几个插入的
int h[N], ph[N], hp[N], size;

// 交换两个点，及其映射关系
void heap_swap(int a, int b)
{
    swap(ph[hp[a]],ph[hp[b]]);
    swap(hp[a], hp[b]);
    swap(h[a], h[b]);
}

void down(int u)
{
    int t = u;
    if (u * 2 <= size && h[u * 2] < h[t]) t = u * 2;
    if (u * 2 + 1 <= size && h[u * 2 + 1] < h[t]) t = u * 2 + 1;
    if (u != t)
    {
        heap_swap(u, t);
        down(t);
    }
}

void up(int u)
{
    while (u / 2 && h[u] < h[u / 2])
    {
        heap_swap(u, u / 2);
        u >>= 1; // u = u >> 1
    }
}

// O(n)建堆
for (int i = n / 2; i; i -- ) down(i);

分析

i为什么从n/2开始down？

首先要明确要进行down操作时必须满足左儿子和右儿子已经是个堆。

开始创建堆的时候，元素是随机插入的，所以不能从根节点开始down，而是要找到满足下面三个性质的结点：

1. 左右儿子满足堆的性质。
2. 下标最大（因为要往上遍历）
3. 不是叶结点（叶节点一定满足堆的性质）

那这个点为什么时n/2？看图。

---

堆排序：

输入一个长度为n的整数数列，从小到大输出前m小的数。

#include <iostream>
using namespace std;

const int N = 100010;

int n, m, h[N], cnt;

void down(int u)
{
    int t = u; // 记一下，原来的下标
    // 首先判断左右儿子分别满足<= cnt，然后取它们之间比q[u]最小的元素下标
    if (2 * u <= cnt && h[2 * u] < h[t]) t = u * 2;
    if (2 * u + 1 <= cnt && h[2 * u + 1] < h[t]) t = u * 2 + 1;
    if (u != t)
    {
    	// 作交换，继续往下调整
        swap(h[u], h[t]);
        down(t);
    }
}

int main()
{
    scanf("%d%d", &n, &m);
    cnt = n;
    for (int i = 1; i <= n; i++) scanf("%d", &h[i]);
    // O(n) 进行堆处理
    for (int i = n / 2; i; i--) down(i);

    while (m--)
    {
        printf("%d ", h[1]);
        // 用堆尾元素覆盖堆顶元素，个数减一，向下调整堆
        h[1] = h[cnt], cnt--, down(1);
    }
    return 0;
}

---

相关例题：

堆排序
模拟堆