STL源码学习----lower_bound和upper_bound算法

本文介绍了使用STL中的lower_bound和upper_bound算法进行高效二分查找的方法,并提供了Python实现示例,展示了如何找到非递减序列中目标值的插入位置。

事实证明,自己用while left right来控制经常出现死循环,边界想不清楚的情况,STL的实现非常巧妙:

  ForwardIter lower_bound(ForwardIter first, ForwardIter last,const _Tp& val)算法返回一个非递减序列[first, last)中的第一个大于等于值val的位置。

     ForwardIter upper_bound(ForwardIter first, ForwardIter last, const _Tp& val)算法返回一个非递减序列[first, last)中第一个大于val的位置。

这样正好是第一个和最后一个插入的位置!!!

用起始位置和长度来控制二分查找非常方便,因为如下(仿照STL 实现方法用Python写了一份):

class Solution:
    def searchInsert(self, nums, target):
        #search lower bound 刚好是第一个插入的位置!!!
        first,cur_len = 0,len(nums)
        while (cur_len > 0):
            half = cur_len >> 1
            if (nums[first+half] < target):
                first = first+half+1
                cur_len = cur_len-(half+1)
            else: #这个因为是nums[first+half] 肯定>= target,所以只保留nums[0..first+half-1],如果上一个if里的元素都比target小,只能插在末尾
                cur_len = half
        return first

    def searchUpperBound(self, nums, target):
        # search upper bound 刚好是最后一个插入的位置!!!
        first,cur_len = 0,len(nums)
        while (cur_len > 0):
            half = cur_len >> 1
            if (nums[first+half] <= target):
                first = first+half+1
                cur_len = cur_len-(half+1)
            else: #这个因为是nums[first+half] 肯定> target,所以只保留nums[0..first+half-1],如果上一个if里的元素都小于等于target,只能插在末尾
                cur_len = half
        return first


s = Solution()
y = s.searchInsert([1,3,5,6],2)
print(y)

 

 

如果你**不使用 `lower_bound` `upper_bound`**,我们可以手动实现这两个功能。这在某些限制环境中(如不允许使用 STL 算法)或学习算法原理时非常有用。 --- ### ✅ 修改后的 C++ 实现(**不使用 `lower_bound` `upper_bound`**): ```cpp #include <bits/stdc++.h> using namespace std; // 手动实现 lower_bound:找到第一个大于等于 x 的位置 int my_lower_bound(int a[], int n, int x) { int left = 0, right = n; while (left < right) { int mid = left + (right - left) / 2; if (a[mid] >= x) right = mid; else left = mid + 1; } return left; } // 手动实现 upper_bound:找到第一个大于 x 的位置 int my_upper_bound(int a[], int n, int x) { int left = 0, right = n; while (left < right) { int mid = left + (right - left) / 2; if (a[mid] > x) right = mid; else left = mid + 1; } return left; } int main() { int n, m; cin >> n >> m; int a[n]; for (int i = 0; i < n; ++i) cin >> a[i]; sort(a, a + n); while (m--) { int x; cin >> x; int lb = my_lower_bound(a, n, x); int ub = my_upper_bound(a, n, x); int less = lb; int equal = ub - lb; int greater = n - ub; cout << less << " " << equal << " " << greater << "\n"; } return 0; } ``` --- ### ✅ 与原代码的区别: | 项目 | 是否使用 | |------|----------| | STL 算法 | ❌ 不使用 `lower_bound` `upper_bound` | | 二分查找 | ✅ 使用手动实现的二分查找 | | 可移植性 | ✅ 可在不支持 STL 的环境中运行 | | 可读性 | ✅ 算法更透明,适合教学或限制环境 | --- ### ✅ 说明: - `my_lower_bound()`:返回第一个 **≥ x** 的元素的索引。 - `my_upper_bound()`:返回第一个 **> x** 的元素的索引。 - 使用经典的二分查找模板,避免溢出并保证正确性。 - 在排序后的数组中,通过这两个函数可以轻松统计小于、等于、大于 `x` 的元素个数。 ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值