AcWing 838 堆排序

题目描述：

输入一个长度为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

分析：

 本题考察堆排序。可以直接使用STL的优先级队列，使用大根堆的话，不断往堆里存入元素，当堆中元素超过m个就删除堆顶元素，到最后留下的m个元素即是最小的m个元素。使用小根堆更为简洁，将所有元素都放到堆里，不断取堆顶元素，取m次就是最小的m个元素了。

#include <iostream>
#include <queue>
#include <vector>
using namespace std;
int n,m,x;
priority_queue<int,vector<int>,greater<int> > pq;
int main(){
    cin>>n>>m;
    for(int i = 0;i < n;i++){
        cin>>x;
        pq.push(x);
    }
    while(m--){
        cout<<pq.top()<<" ";
        pq.pop();
    }
    return 0;
}

有时需要手写堆来提高效率和添加功能，数组模拟堆实现堆排序的代码如下：

#include<iostream>
using namespace std;
const int maxn = 100005;
int m,n,h[maxn],size;
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){
        swap(h[u],h[t]);
        down(t);
    }
}

int main(){
    cin>>n>>m;
    for(int i = 1;i <= n;i++)   cin>>h[i];
    size = n;
    for(int i = n / 2;i;i--)    down(i);//O（n）建堆
    while(m--){
        cout<<h[1]<<" ";
        h[1] = h[size--];
        down(1);
    }
    return 0;
}

首先说下堆的存储，堆是一棵完全二叉树，使用可以直接采用数组模拟，某节点i的左右孩子下标分别为2i和2i+1，i的父节点为i / 2.

首先介绍下下滤操作，将某节点与左右孩子比较，选择较小的元素作为新的根节点，具体代码见上面的代码。

再介绍下建堆的操作：for(int i = n / 2;i;i--)    down(i);//O（n）建堆，由于堆是完全二叉树，n / 2的节点必然位于倒数第二层。下滤操作最多执行一次，n/4的节点位于倒数第三层，下滤操作最多执行两次...故一共会执行n / 4 * 1 + n / 8 * 2 + n / 16 * 3+...次。

比如上图的堆一共有n = 9个节点，n / 2 = 4，第四个节点是倒数第二层第一个节点。既然是从倒数第二层开始下滤，只需要下滤有孩子节点的节点，编号超过n / 2的节点度均为0，不需要下滤，故循环从n / 2开始。最坏情况下是满二叉树，n / 2是倒数第二层最后一个节点，倒数第二层有n / 4个节点需要下滤，倒数第三层有n / 8个节点需要下滤，总操作次数S = n / 4 * 1 + n / 8 * 2 + n / 16 * 3+...。后面乘的数表示需要下滤的层数，倒数第二层仅需下滤一层，倒数第三层需要下滤两层，以此类推。S / 2 = n / 8 * 1 + n / 16 * 2 + n /32 * 3,两式相减得S / 2 = n / 4 + n / 8 + n / 16 +...< n,故总的操作次数是O（n）级别的，这也就是建堆操作的的时间复杂度。

最后再介绍删除操作，堆顶元素始终为h[1],删除h[1]只需将最后一个元素赋给h[1],再不断下滤即可。