三元组存储结构需定义为按行主序、列次序升序排列的非零元素序列，如vector<Triple>或vector<tuple<int,int,double>>，否则双指针合并会漏项或重复；相加时须严格比较(row,col)字典序，相同位置值相加后需用容差判断是否保留（如abs(sum) > 1e-12），并预先reserve容量以避免频繁realloc。
三元组存储结构怎么定义才方便相加
稀疏矩阵用三元组（行号、列号、值）存储时，vector<tuple<int, int, double>> 或自定义结构体都行，但必须保证：所有非零元按 row 主序、col 次序升序排列。否则后续双指针合并会出错——不是结果错，而是漏项或重复插入。
常见错误是直接 push_back 未排序的元素，或者从文件读入后没调用 sort。正确做法是插入后统一排序，或在插入时二分查找定位（小规模数据直接 sort 更稳妥）：
sort(triples.begin(), triples.end(), [](const auto& a, const auto& b) {
if (get<0>(a) != get<0>(b)) return get<0>(a) < get<0>(b);
return get<1>(a) < get<1>(b);
});
双指针合并时如何处理相同位置的元素
两个已排序三元组向量 A 和 B 相加，核心是模拟归并过程，但需合并相同 (i,j) 位置的值。不能简单“谁小谁先走”，必须显式比对行列坐标：
- 当
A[k]和B[l]的row和col都相等 → 把值相加，若结果非零才存入结果；为零则跳过 - 当
A[k]坐标字典序小于B[l]→ 取A[k]，k++ - 反之取
B[l]，l++
注意：C++ 中 tuple<int,int,double> 默认按字典序比较，可直接用 < 判断大小关系，但必须确保行列顺序一致（即 make_tuple(row, col, val)）。
立即学习“C++免费学习笔记（深入）”；
为什么不能用 map<pair<int,int>, double> 替代三元组 vector
虽然 map<pair<int,int>, double> 插入自动去重+排序，看似省事，但实际性能差很多：
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
- 每次插入是
O(log n)，而预排序 + 合并是O(n + m) - 内存碎片多，缓存不友好，尤其在千级非零元以上时，实测慢 2–5 倍
-
map无法保证遍历顺序严格按行主序（除非用std::map且 key 是pair，但仍有常数开销）
更关键的是：三元组快速相加算法依赖顺序遍历和双指针跳跃，map 迭代器不支持随机跳转，没法高效实现“跳过整行”这类优化。
结果中零元素怎么安全剔除
相加后某位置值为 0 是常见情况（比如 5.0 + (-5.0)），必须剔除，否则破坏稀疏性。但不能边遍历边 erase —— 会破坏迭代器有效性。
推荐做法是：先用双指针生成带零的结果 vector，再用 remove_if + erase 两步清除：
auto it = remove_if(res.begin(), res.end(), [](const auto& t) { return abs(get<2>(t)) < 1e-12; });
res.erase(it, res.end());
注意浮点比较要用容差（1e-12），别用 == 0.0；整数场景可直接判等。这个步骤不能省，否则后续所有基于三元组的操作（如转 CSR、矩阵乘）都会出错。
真正麻烦的是边界情况：两个空矩阵相加、一个全零结果、大量抵消导致结果比输入还稀疏——这些都要靠上述剔除逻辑兜底，而不是靠“假设不会发生”。
