数据结构：堆（堆排序和模拟堆）

1.堆排序

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

输入格式
第一行包含整数n和m。
第二行包含n个整数，表示整数数列。

输出格式
共一行，包含m个整数，表示整数数列中前m小的数。

数据范围
1≤m≤n≤10^5，
1≤数列中元素≤10^9
输入样例：
5 3
4 5 1 3 2
输出样例：
1 2 3

算法的基本思想：
按完全二叉树存放数据，则数据中x的左二子为2x,右儿子为2x+1的关系，利用down（向下移动）、up（向上移动）函数实现题目要求。
求最小值： heap[1];
插入：heap[++size] = x; up[size];
删除最小值： heap[1] = heap[size]; size–; down(1);
删除任意值：heap[k] = heap[size]; size–; down(k); up(k);
修改任意值：heap[k] = x; down(k); up(k);

程序代码：

#include <iostream>
using namespace std;

const int N = 100010;
int n, m, s, h[N];

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

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

int main()
{
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; i ++)  scanf("%d", &h[i]);
    s = n;
    
    for(int i = n / 2; i ; i --)  down(i);
    
    while(m --)
    {
        printf("%d ", h[1]);
        h[1] = h[s];
        s --;
        down(1);
    }
    return 0;
}

2.模拟堆

维护一个集合，初始时集合为空，支持如下几种操作：
“I x”，插入一个数x；
“PM”，输出当前集合中的最小值；
“DM”，删除当前集合中的最小值（数据保证此时的最小值唯一）；
“D k”，删除第k个插入的数；
“C k x”，修改第k个插入的数，将其变为x；
现在要进行N次操作，对于所有第2个操作，输出当前集合的最小值。

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

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

数据范围
1≤N≤10⁵
−10⁹≤x≤10⁹
数据保证合法。

输入样例：
8
I -10
PM
I -10
D 1
C 2 8
I 6
PM
DM
输出样例：
-10
6

程序代码：

#include <iostream>
#include <algorithm>
#include <string.h>
using namespace std;

const int N = 100010;
int h[N], ph[N], hp[N], s;//ph存放第k个插入点的下标, 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]);
}

void down(int u)
{
    int t = u;
    if(u * 2 <= s && h[u * 2] < h[t])  t = u * 2;
    if(u * 2 + 1 <= s && 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 / 2] > h[u])
    {
        heap_swap(u, u / 2);
        u /= 2;
    }
}

int main()
{
    int n, m = 0;
    scanf("%d", &n);
    while(n --)
    {
        char op[10];
        int k, x;
        
        scanf("%s", op);
        if(!strcmp(op, "I"))
        {
            scanf("%d", &x);
            s ++, m ++;
            ph[m] = s, hp[s] = m;
            h[s] = x;
            up(s);
        }
        else if(!strcmp(op, "PM"))  printf("%d\n", h[1]);
        else if(!strcmp(op, "DM"))
        {
            heap_swap(1, s);
            s --;
            down(1);
        }
        else if(!strcmp(op, "D"))
        {
            scanf("%d", &k);
            k = ph[k];
            heap_swap(k, s);
            s --;
            down(k), up(k);
        }
        else
        {
            scanf("%d%d", &k, &x);
            k = ph[k];
            h[k] = x;
            down(k), up(k);
        }
    }
    return 0;
}