三元组相加不能暴力合并，因时间复杂度退化为O(nnz_A×nnz_B)；需先按行主序排序再双指针归并，复杂度降至O(nnz_A+nnz_B)，并跳过结果为零项；转置宜用计数排序替代比较排序，再局部排序保证列内行序；应抽象坐标比较器与映射器实现逻辑复用；须严控空矩阵、越界、浮点精度等边界条件。
三元组相加为什么不能直接按行主序暴力合并？
因为稀疏矩阵的非零元位置高度不规则，暴力遍历两组三元组（vector<tuple<int,int,double>>）逐个比对行列索引，时间复杂度会退化到 O(nnz_A * nnz_B)，实际中常卡在 10⁵ 级非零元时就超时。真正高效的相加必须利用行列索引的局部有序性——哪怕原始三元组未排序，也应先按行优先（再列）排序，然后用双指针归并，复杂度降为 O(nnz_A + nnz_B)。
实操建议：
立即学习“C++免费学习笔记（深入）”；
- 输入三元组务必先调用
std::sort，排序谓词为[&](const auto& a, const auto& b) { return get<0>(a) < get<0>(b) || (get<0>(a) == get<0>(b) && get<1>(a) < get<1>(b)); } - 归并时用两个游标
i、j，比较(row_i, col_i)与(row_j, col_j)：相等则值相加（注意浮点精度容差，比如fabs(val_i + val_j) > 1e-12才保留）；前者小则推i，后者小则推j - 别忘了跳过结果为零的项——相加后值被抵消掉的三元组必须丢弃，否则破坏稀疏性
转置时为什么不能只交换行列索引再重排序？
能交换，但重排序代价高。标准三元组转置本质是“按列主序重排”，而列主序 ≡ 按 (col, row) 排序。如果原始三元组已是行主序（(row, col)），直接 std::sort 一次没问题；但若后续还要频繁访问（如多次转置或乘法），这种每次 O(nnz log nnz) 的开销不可接受。
实操建议：
立即学习“C++免费学习笔记（深入）”；
- 用计数排序替代比较排序：先扫描一遍统计每列非零元个数 → 得到各列在转置结果中的起始偏移
col_ptr[col]→ 再扫一遍原三元组，根据col直接写入目标位置。总时间O(nnz + ncol)，无比较、稳定、缓存友好 - 转置结果仍需保证列内按行升序（即转置后是列主序），所以计数排序后要再对每列内部按
row排序——但这时每列数据已连续存放，可用std::sort局部排序，或更优：原始输入若按行主序，转置时用“行号作为键、列号作为值”建哈希桶，再对每个桶内row排序 - 避免复制：用
std::vector<std::vector<std::pair<int, double>>>存每列的(row, val)，转置后直接构造新三元组，减少中间内存分配
如何让相加和转置共享同一套索引管理逻辑？
硬编码两套排序/归并逻辑会导致维护困难，尤其当需要支持不同存储格式（COO、CSR、CSC）时。核心是抽象出“坐标比较器”和“索引映射器”两个组件。
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
实操建议：
立即学习“C++免费学习笔记（深入）”；
- 定义统一结构体
Triplet，含row、col、val成员，而非tuple，便于扩展字段（如是否已归一化） - 用模板参数传入比较策略，例如
struct RowMajorCmp { bool operator()(const Triplet& a, const Triplet& b) const { return a.row < b.row || (a.row == b.row && a.col < b.col); } };，相加用它，转置则用ColMajorCmp - 转置本身可视为一种“坐标变换操作”，封装为
transform_coords(const vector<Triplet>& in, vector<Triplet>& out, function<void(int&, int&)> f)，其中f是[&](int& r, int& c) { swap(r, c); }—— 这样未来支持旋转、镜像等变换只需换 lambda
实际部署时最容易被忽略的边界条件
不是算法逻辑错，而是工程细节崩盘：空矩阵、全零结果、行列号越界、浮点溢出、内存对齐导致的 vector realloc 失败。
实操建议：
立即学习“C++免费学习笔记（深入）”；
- 输入前检查
row >= 0 && row < rows && col >= 0 && col < cols，否则抛std::out_of_range或静默截断（根据场景选） - 相加后若结果为空，返回空
vector<Triplet>，不要留一个(0,0,0.0)占位——下游解析器可能误判维度 - 转置时若原始矩阵列数未知（仅给三元组），必须从所有
col中取max_element+ 1，不能假设列数等于行数 - 用
reserve()预分配结果容器：相加最多nnz_A + nnz_B项，转置数量不变，避免多次扩容拷贝
稀疏矩阵运算的性能瓶颈往往不在算法复杂度，而在内存访问模式和分支预测失败。把排序换成计数、把 tuple 换成 struct、把 inline lambda 拆成命名函数——这些微小调整，在百万级非零元场景下，差异就是秒级与毫秒级的区别。
