算法学习系列(十五):最小堆、堆排序

目录

引言

这个堆排序的话,考的还挺多的,主要是构建最小堆,并且在很多情况下某些东西还用得着它来优化,比如说迪杰斯特拉算法可以用最小堆优化,然后面试和考研用的也是挺多的,总之开始吧。

一、最小堆概念

本文只讲述最小堆,其一这个用的最多,而且跟最大堆来说其实都是差不多的,就一个小于一个大于

最小堆:首先是一个完全二叉树,然后每个结点都小于或等于其两个儿子,性质:根结点是整个堆的最小值。因为父亲是最小的,然后儿子又作为其儿子的父亲,也是其最小的,所以推出堆根是最小的。

存储方式 :是用数组来存的,i 号下标的儿子为 2 * i,2 * i + 1 ,i 号下标的父亲为 i / 2

STL:优先级队列就是最小堆

二、堆排序模板(最小堆)

整体思路:先构建一个最小堆,然后输出堆根,再把堆根删了,再次构建,重复往返

删除堆根:用 h[size] 覆盖 h[1] ,size-- ,down(1)

这个模板我们用例题来说明

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

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

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

数据范围
1≤m≤n≤105,1≤数列中元素≤109

输入样例:
5 3
4 5 1 3 2

输出样例:
1 2 3
c 复制代码
#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int N = 1e5+10;

int h[N];  //h[i]代表第堆中的i号下标对应的值
int n, m, cnt;  //cnt代表堆中的数量

//元素变小了,就up
void up(int u)
{
    while(u / 2 && h[u / 2] > h[u])  //若u>=1 并且比父亲大就交换,然后u变成父亲继续判断
    {
        swap(h[u], h[u/2]);
        u >>= 1;
    }
}

//元素变大了,就down
void down(int u)  //将下标为u的元素下移
{
    int t = u;
    if(u * 2 <= cnt && h[2 * u] < h[t]) t = 2 * u;  
    if(u * 2 + 1 <= cnt && h[2 * u + 1] < h[t]) t = 2 * u + 1;  //判断最小的值t
    
    if(t != u)  //若有一个儿子比自己小
    {
        swap(h[u], h[t]);  //交换值
        down(t);  //再次判断这个儿子
    }
}

int main()
{
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= n; ++i) scanf("%d", &h[i]);
    cnt = n;
    
    for(int i = cnt / 2; i; --i) down(i);  //从倒数第二层的元素开始down到根,可以构建一个最小堆
    
    while(m--)
    {
        printf("%d ", h[1]);  //每次都是堆根的元素最小
        swap(h[1], h[cnt]);
        cnt--;
        down(1);
    }
    
    return 0;
}

三、模拟堆

我们还是用例题来说明

然后这个有一个要求就是要对第k个插入的元素进行操作,但是我们又不知道第k个元素是谁,只知道下标,所以得维护两个数组ph[i] ,hp[i],代表ph[k] == i,hp[i] == k,h[i] == a

维护一个集合,初始时集合为空,支持如下几种操作:
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
c 复制代码
#include <cstdio>
#include <cstring>
#include <iostream>

using namespace std;

const int N = 1e5+10;

int h[N], hp[N], ph[N];
int n, cnt, idx;

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

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

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_heap(t,u);
        down(t);
    }
}

int main()
{
    scanf("%d", &n);
    
    while(n--)
    {
        char op[5];
        int k, x;
        
        scanf("%s", op);
        if(!strcmp(op,"I"))
        {
            scanf("%d", &x);
            cnt++, idx++;
            h[cnt] = x;
            hp[cnt] = idx, ph[idx] = cnt;
            up(cnt);
        }
        else if(!strcmp(op,"PM"))
        {
            printf("%d\n", h[1]);
        }
        else if(!strcmp(op,"DM"))
        {
            swap_heap(1,cnt);
            cnt--;
            down(1);
        }
        else if(!strcmp(op,"D"))
        {
            scanf("%d", &k);
            k = ph[k];
            swap_heap(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;
}
相关推荐
shymoy21 分钟前
Radix Sorts
数据结构·算法·排序算法
风影小子30 分钟前
注册登录学生管理系统小项目
算法
黑龙江亿林等保32 分钟前
深入探索哈尔滨二级等保下的负载均衡SLB及其核心算法
运维·算法·负载均衡
lucy1530275107935 分钟前
【青牛科技】GC5931:工业风扇驱动芯片的卓越替代者
人工智能·科技·单片机·嵌入式硬件·算法·机器学习
杜杜的man1 小时前
【go从零单排】迭代器(Iterators)
开发语言·算法·golang
怀旧6661 小时前
spring boot 项目配置https服务
java·spring boot·后端·学习·个人开发·1024程序员节
小沈熬夜秃头中୧⍤⃝1 小时前
【贪心算法】No.1---贪心算法(1)
算法·贪心算法
木向2 小时前
leetcode92:反转链表||
数据结构·c++·算法·leetcode·链表
阿阿越2 小时前
算法每日练 -- 双指针篇(持续更新中)
数据结构·c++·算法
skaiuijing2 小时前
Sparrow系列拓展篇:对调度层进行抽象并引入IPC机制信号量
c语言·算法·操作系统·调度算法·操作系统内核