NC119 最小的K个数

最小的K个数

题目描述

给定一个数组,找出其中最小的K个数。例如数组元素是4,5,1,6,2,7,3,8这8个数字,则最小的4个数字是1,2,3,4。如果K>数组的长度,那么返回一个空的数组。

题目传送门

示例

输入

[4,5,1,6,2,7,3,8],4

输出

[1,2,3,4]

题解

方法一:直接排序

选择一种排序方法,对输的数据进行一次升序排序,通过循环取出前边最小的K个数。

时间复杂度:O( n l o g n nlog^n nlogn)

间复杂度:O(1)

class Solution {
public:
    vector<int> GetLeastNumbers_Solution(vector<int> input, int k) {
        vector<int> result;
        if(k==0 || input.empty() || k>input.size())
            return result;
        sort(input.begin(), input.end());
        for(int i=0; i<k; i++)
            result.push_back(input[i]);
        return result;
    }
};

方法二:优先队列

普通的队列是先进先出的数据结构,元素在队列尾追加,从队列头删除。在优先队列中,元素被赋予优先级,最高优先级的元素最先删除。优先队列包含队列的基本操作,添加了一个内部的排序,使用堆来实现内部的排序。建立一个容量为k的大根堆的优先队列。遍历一边元素,如果队列的大小<k就直接入队,否则让当前元素与队顶元素相比,如果队顶元素大,则出队并将当前元素入队。

时间复杂度:O( n l o g k nlog^k nlogk), 插入容量为k的大根堆时间复杂度为O( l o g k log^k logk), 一共遍历n个元素。
空间复杂度:O(k)

class Solution {
public:
    vector<int> GetLeastNumbers_Solution(vector<int> input, int k) {
        vector<int> result;
        if (k==0 || input.empty() || k>input.size())
            return result;
        priority_queue<int, vector<int>> pq;
        for(const int val : input) {
            if (pq.size() < k) {
                pq.push(val);
            } else {
                if(val < pq.top()) {
                    pq.pop();
                    pq.push(val);
                }
            }
        }
        while(!pq.empty()) {
            result.push_back(pq.top());
            pq.pop();
        }
        return result;
    }
};

方法三:利用set数据结构

set数据结构对插入的元素可以自动排序,建立大小为k的multiset,时间复杂度为O( n l o g k nlog^k nlogk)。不能直接修改set容器内的数据,所以只能删除某个元素再插入要修改的元素。

class Solution {
public:
    vector<int> GetLeastNumbers_Solution(vector<int> input, int k) {
        vector<int> result;
        if (k==0 || input.empty() || k>input.size())
            return result;
        // 仿函数greater<T>模板,从大到小排序,首元素为最大值
        multiset<int, greater<int> > leastNums;
        vector<int>::iterator vec_it=input.begin();
        for(; vec_it!=input.end(); vec_it++) {
            // 将前k个元素插入set集合
            if(leastNums.size() < k)
                leastNums.insert(*vec_it);
            else {
                // 第一个元素是最大值
                multiset<int, greater<int> >::iterator greatest_it=leastNums.begin();
                // 如果后续元素小于第一个元素,删除第一个,插入当前元素
                if(*vec_it < *(leastNums.begin())) {
                    leastNums.erase(greatest_it);
                    leastNums.insert(*vec_it);
                }
            }
        }
        for(multiset<int>::iterator set_it=leastNums.begin(); set_it!=leastNums.end(); set_it++)
            result.push_back(*set_it);
        return result;
        
    }
};

方法四:最大堆

用最大堆实现,堆的大小为k,调用algorithm里边的make_heap等现成的函数。

class Solution {
public:
    vector<int> GetLeastNumbers_Solution(vector<int> input, int k) {
        if (k==0 || input.empty() || k>input.size())
            return vector<int>{};
        vector<int> result(input.begin(), input.begin()+k);
        // 建立最大堆
        make_heap(result.begin(), result.end());
        for(int i=k; i<input.size(); i++) {
            if(result[0] > input[i]) {
                // 把最大值移动到最后一个元素
                pop_heap(result.begin(), result.end());
                // 移除最后一个元素
                result.pop_back();
                result.push_back(input[i]);
                // 重新选出最大值
                push_heap(result.begin(),result.end());
            }
        }
        sort_heap(result.begin(), result.end());
        return result;
        
    }
};

方法五:借鉴快排的Partition

快速排序有一个基准点,这个基准点左边都是较小的,右边都是较大值。只要这个基准点的索引值为K-1就完成了。

时间复杂度:O(n)

class Solution {
public:
    int partition(vector<int> &input, int left, int right) {
        int pivot = input[right-1];
        int i = left;
        for (int j=left; j<right-1; ++j) {
            if (input[j] < pivot) {
                swap(input[i++], input[j]);
            }
        }
        swap(input[i], input[right-1]);
        return i;
    }
    vector<int> GetLeastNumbers_Solution(vector<int> input, int k) {
        vector<int> result;
        if (k==0 || input.empty() || k>input.size())
            return result;
        int left = 0, right = input.size();
        while (left < right) {
            int pivot = partition(input, left, right);
            if (pivot+1 == k) {
                for(int i=0; i<k; i++) {
                    result.push_back(input[i]);
                }
            }
            if (pivot+1 < k) {
                left = pivot + 1;
            }  
            else {
                right = pivot;
            }
        }
        return result;
    }
};
下面是一个简单的MATLAB脚本,实现了多新息最小二乘辨识算法(MI-LS): ```matlab clear; clf; format short g % 设置参数 p = 1; % 创新长度 p=1,5,10 length1 = 1000; % 数据长度 % 多项式A(z),B(z),D(z)的阶数 na = 2; nb = 2; nc = 2; nd = 2; nl = na + nb; n2 = nc + nd; n = na + nb + nc + nd; % A(z),B(z),D(z)的系数向量 a = [1, -0.56, 0.42]; b = [0, 0.9, 0.6]; c = [1, -0.3, 0.2]; d = [1, 0.3, 0, -0.20]; c1 = conv(a, c); theta = [a(2:na+1); b(2:nb+1); c(2:nc+1); d(2:nd+1)]; % 参数向量 theta n = length(theta); fprintf('\nMI-LS算法\n'); fprintf('创新长度p=%d\n', p); % 准备输入数据和噪声数据 rand('state', 15); % 设置随机变量的状态 u = (rand(length1, 1) - 0.5) * sqrt(12); % 输入随机序列 randn('state', 15); v = randn(length1, 1); % 噪声随机序列 % 生成输出数据 y = zeros(length1, 1); w = zeros(n, 1); for t = n:length1 w(t) = theta' * [-w(t-1:-1:t-nc); v(t-1:-1:t-nd)] + v(t); y(t) = theta(1:na)' * [-y(t-1:-1:t-na); u(t-1:-1:t-nb)] + w(t); end % MI-LS算法 Phi = zeros(length1-p+1, n*p); Y = zeros(length1-p+1, 1); for t = p:length1 for k = t:-1:t-p+1 Phi(k, :) = [-y(k-1:-1:k-na); u(k-1:-1:k-nb); -w(k-1:-1:k-nc); v(k-1:-1:k-nd)]; end Y(t) = y(t); end theta_hat = Phi \ Y; % 参数估计 % 输出结果 fprintf('真实参数向量:\n'); disp(theta); fprintf('估计参数向量:\n'); disp(theta_hat); ``` 请注意,这只是一个简单的示例,可能需要根据您的具体需求进行修改和优化。希望对您有所帮助!
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值