AcWing 839 模拟堆

题目描述：

维护一个集合，初始时集合为空，支持如下几种操作：

1. “I x”，插入一个数x；
2. “PM”，输出当前集合中的最小值；
3. “DM”，删除当前集合中的最小值（当最小值不唯一时，删除最早插入的最小值）；
4. “D k”，删除第k个插入的数；
5. “C k x”，修改第k个插入的数，将其变为x；

现在要进行N次操作，对于所有第2个操作，输出当前集合的最小值。

输入格式

第一行包含整数N。

接下来N行，每行包含一个操作指令，操作指令为”I x”，”PM”，”DM”，”D k”或”C k x”中的一种。

输出格式

对于每个输出指令“PM”，输出一个结果，表示当前集合中的最小值。

每个结果占一行。

数据范围

1≤N≤10^5
−10^9≤x≤10^9
数据保证合法。

输入样例：

8
I -10
PM
I -10
D 1
C 2 8
I 6
PM
DM

输出样例：

-10
6

分析：

本题要求模拟的堆能够做到删除和修改第k个插入的元素，这是普通的堆做不到的。要想修改第k个插入的元素，我们需要知道它的下标是多少，所以用数组ph[k]存储第k个插入的元素在堆h中的下标。另外，在上滤和下滤时，需要交换两个节点，比如要交换某节点和其父节点，ph数组的值也得交换，但是如何知道要交换的两个节点是第几个插入的数呢？只好再定义一个hp数组存储某节点是第几个插入的数。这样在交换节点时，我们需要交换三个数组。

void heap_swap(int a, int b)
{
    swap(ph[hp[a]],ph[hp[b]]);
    swap(hp[a], hp[b]);
    swap(h[a], h[b]);
}

要交换在h中下标为a，b的两个节点，需要先交换ph数组的值，a是在第hp[a]次被插入的，b是在第hp[b]次被插入的，所以需要先交换ph[hp[a]]和ph[hp[b]]，使用过hp数组原本的值后才能开始交换hp数组的值，最后交换元素的值。

上滤下滤操作较为简单，唯一要注意的地方是，删除某节点的步骤，调用上述的交换节点的函数交换待删除的节点和堆的末节点，然后是该上滤还是下滤呢？若值被改变的节点值变大了，就需要下滤，否则就需要上滤，不妨都调用一下，上滤和下滤的函数顶多会执行一次。对于修改第k个插入节点值之后的操作也是要么上滤要么下滤，如出一辙。

总的代码： 

#include <iostream>
#include <algorithm>
#include <cstring>
using namespace std;
const int N = 100010;
int h[N], ph[N], hp[N], size;//ph[k]-第k个插入点的下标,hp[k]-堆里第k个点是第几个插入的点
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;
    }
}

int main()
{
    int n, m = 0;
    scanf("%d", &n);
    while (n -- )
    {
        char op[5];
        int k, x;
        scanf("%s", op);
        if (!strcmp(op, "I"))//插入元素
        {
            scanf("%d", &x);
            size ++ ;
            m ++ ;
            ph[m] = size, hp[size] = m;//初始化两个指针数组
            h[size] = x;//插入末尾然后上滤
            up(size);
        }
        else if (!strcmp(op, "PM")) printf("%d\n", h[1]);//输出堆顶的值
        else if (!strcmp(op, "DM"))//删除最小值
        {
            heap_swap(1, size);
            size -- ;
            down(1);
        }
        else if (!strcmp(op, "D"))//删除第k个插入的元素
        {
            scanf("%d", &k);
            k = ph[k];
            heap_swap(k, size);
            size -- ;
            up(k);//尝试上滤
            down(k);//尝试下滤
        }
        else//修改第k个插入的元素
        {
            scanf("%d%d", &k, &x);
            k = ph[k];
            h[k] = x;
            up(k);
            down(k);
        }
    }

    return 0;
}