堆

大雪菜的课（笔记）

数据结构（二）

3.堆

(1).模板(堆 —— 模板题 AcWing 838. 堆排序, AcWing 839. 模拟堆)

// 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;
    }
}

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

AcWing838. 堆排序

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

输入格式
第一行包含整数n和m。

第二行包含n个整数，表示整数数列。

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

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

#include <iostream>
using namespace std;
const int N=100010;
int h[N],cnt;
void down(int u){
    int t=u;
    if(u*2<=cnt&&h[u*2]<h[t])   t=u*2;
    if(u*2+1<=cnt&&h[u*2+1]<h[t])   t=u*2+1;
    if(t!=u){
        swap(h[t],h[u]);
        down(t);
    }
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)   scanf("%d",&h[i]);
    cnt=n;
    for(int i=n/2;i>=1;i--) down(i);
    while(m--){
        printf("%d ",h[1]);
        h[1]=h[cnt--];
        down(1);
    }
}

AcWing839. 模拟堆

题目
提交记录
讨论
题解
视频讲解

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

“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≤105
−109≤x≤109
数据保证合法。

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

#include <iostream>
#include <string.h>
using namespace std;
const int N=100010;
int h[N],ph[N],hp[N],cnt;
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<=cnt&&h[u*2]<h[t])   t=u*2;
    if(u*2+1<=cnt&&h[u*2+1]<h[t])   t=u*2+1;
    if(t!=u){
        heap_swap(t,u);
        down(t);
    }
}
void up(int u){
    while(u/2&&h[u/2]>h[u]){
        heap_swap(u/2,u);
        u>>=1;
    }
}
int main()
{
    int n,m=0;
    scanf("%d",&n);
    while(n--){
        char s[5];
        int k,x;
        scanf("%s",s);
        if(!strcmp(s,"I")){
            scanf("%d",&x);
            cnt++;
            m++;
            h[cnt]=x;
            ph[m]=cnt,hp[cnt]=m;
            up(cnt);
        }else if(!strcmp(s,"PM"))   printf("%d\n",h[1]);
        else if(!strcmp(s,"DM")){
            heap_swap(1,cnt);
            cnt--;
            down(1);
        }else if(!strcmp(s,"D")){
            scanf("%d",&k);
            k=ph[k];
            heap_swap(k,cnt);
            cnt--;
            up(k);
            down(k);
        }else{
            scanf("%d%d",&k,&x);
            k=ph[k];
            h[k]=x;
            up(k);
            down(k);
        }
    }
    return 0;
}