C++标准库无直方图函数，需手动用std::vector分桶；对自定义边界桶（如{0.0,0.5,1.2,2.0,5.0}），应以std::lower_bound二分查找首个≥value的边界，索引减1得桶号，并确保value∈[bound[0], bound.back())。
用 std::vector 手动分桶，别依赖 std::histogram
C++ 标准库至今（C++23）**没有内置直方图函数**，std::histogram 是提案（P2301R3），尚未进入标准。想统计浮点数在桶区间的分布，必须自己实现分桶逻辑。核心是：给定数据、左闭右开（或统一约定）的桶边界，对每个数查它落在哪个索引位置。
常见错误是直接用 value / bin_width 取整——这仅适用于从 0 开始、等宽、连续无缺口的桶。实际中桶常为自定义边界（如 {0.0, 0.5, 1.2, 2.0, 5.0}），必须用二分查找。
- 推荐用
std::lower_bound查找第一个 ≥value的边界位置，该位置减 1 即为所属桶索引（需检查是否越界） - 桶定义为左闭右开区间：
[bound[i], bound[i+1])，这样能避免相邻桶边界值归属歧义 - 所有输入值必须满足
value ≥ bound[0]且value ，否则属于“溢出”，需单独计数或丢弃
处理边界和溢出：用哨兵值 + 预检查
如果桶边界是 std::vector<double> bounds = {0.0, 1.0, 2.0, 3.0}</double>，那只有 3 个有效桶（索引 0~2），对应区间 [0,1)、[1,2)、[2,3)。任何 ≥3.0 或 的值都不属于这些桶。
不加检查直接调 std::lower_bound 可能返回 bounds.end()，解引用会崩溃；用 std::distance 得到索引后若等于 bounds.size()，说明值 ≥ 最大边界。
立即学习“C++免费学习笔记（深入）”；
- 初始化直方图向量为
std::vector<size_t>(bounds.size() - 1)</size_t>，长度 = 桶数 - 对每个
val，先判断val = bounds.back()，真则归入underflow/overflow计数器 - 否则用
auto it = std::lower_bound(bounds.begin(), bounds.end(), val)，桶索引为std::distance(bounds.begin(), it) - 1
性能关键：预排序桶边界 + 用 std::lower_bound
桶边界必须严格递增，否则 std::lower_bound 行为未定义。如果边界来自用户输入或计算，务必先调 std::sort 去重并排序（但注意：去重会合并桶，通常不应自动去重，应报错或告警）。
每次查找复杂度是 O(log N)，N 是桶数。若桶数极少（如 ≤10），线性扫描反而更快；但超过 20 桶，二分优势明显。别用 std::find_if 遍历——尤其当数据量大时，整体复杂度会从 O(M log N) 退化到 O(M×N)（M 是数据点数）。
- 边界只须排序一次，可缓存
std::vector<double></double>并复用 - 若桶等宽且从某起点开始，可用公式
index = static_cast<int>((val - min_val) / bin_width)</int>，但必须确保val不越界，且用std::floor处理负数（避免截断误差） - 浮点比较本身有精度风险，但
std::lower_bound在有序容器上基于比较，只要边界和值都用相同类型（如全 <code>double），就足够可靠
完整最小可运行示例
#include <vector>
#include <algorithm>
#include <iostream>
<p>std::vector<size_t> histogram(const std::vector<double>& data,
const std::vector<double>& bounds) {
if (bounds.size() < 2) return {};
std::vector<size_t> hist(bounds.size() - 1, 0);
size_t under = 0, over = 0;</p><pre class='brush:php;toolbar:false;'>for (double val : data) {
if (val < bounds[0]) { ++under; continue; }
if (val >= bounds.back()) { ++over; continue; }
auto it = std::lower_bound(bounds.begin(), bounds.end(), val);
int idx = std::distance(bounds.begin(), it) - 1;
if (idx >= 0 && idx < (int)hist.size()) {
++hist[idx];
}
}
// 可选：返回 under/over 作为额外信息
return hist;}
// 使用：
// auto bounds = std::vector
真正容易被忽略的是：桶边界必须由使用者保证有序且无重复；std::lower_bound 返回的是“插入位置”，不是“匹配位置”——哪怕 val 恰好等于某个边界点，它也会落入左侧桶（因区间为左闭右开）。这点和直觉可能相反，但恰恰是避免边界值二义性的设计前提。
