三元组稀疏矩阵加法不能直接按行遍历,因mata与matb非零元位置乱序,需先按(row,col)字典序排序再归并;转置需两次扫描统计列频与定位,而非简单交换行列;格式化输出应模拟稠密布局逐行填充;求和前须校验维度、内存及浮点零值。

三元组稀疏矩阵加法为什么不能直接按行遍历
因为两个三元组表 matA 和 matB 的非零元素位置是乱序存储的(即使各自按行优先排序,合并时仍需对齐行列索引),直接双指针逐行扫容易漏项或重复累加。正确做法是把两组三元组统一按 (row, col) 字典序排序后归并——这和归并两个有序链表逻辑一致,但必须先保证输入各自已按行主序、再按列主序排序。
常见错误现象:addTriplets() 结果中出现重复坐标、漏掉某处非零和、甚至越界写入 data[] 数组。
- 输入三元组表必须预先调用
sortTripletsByRowCol()排序(不是只按行号排) - 归并时用三个指针:
i指向matA,j指向matB,k指向结果matC - 当
matA[i].row == matB[j].row && matA[i].col == matB[j].col时才相加;否则取较小者插入 - 结果容量要预分配足够:最坏情况是所有非零元都不重叠,
matC.size = matA.size + matB.size
转置操作为何要两次扫描而非交换 row/col 字段
单纯交换每个三元组的 row 和 col 字段,只会得到按列主序排列的三元组表,而标准转置要求输出仍是行主序(即新矩阵的第 0 行对应原矩阵第 0 列所有非零元,且该行内按列号升序排列)。所以必须统计每列非零元个数,算出每行起始位置,再做一次填充。
典型错误:transpose() 后发现输出三元组顺序错乱,比如原矩阵第 2 列的两个非零元在转置后不在新矩阵第 2 行连续出现。
- 第一次扫描:统计原矩阵每列非零元数量,存入
numCols[col] - 计算
cpot[col](cumulative position):表示转置后第col行在结果中的起始下标 - 第二次扫描:对每个三元组
(r, c, v),查cpot[c]得到目标位置,填入(c, r, v),然后cpot[c]++ - 注意:
cpot数组长度应为原矩阵列数cols,而非行数
格式化输出时如何避免行列错位和对齐失真
稀疏矩阵打印不是简单循环输出三元组,而是模拟稠密矩阵布局:逐行扫描,遇到非零元就输出值,否则补空格或 0。但若行列范围大(如 1000×1000 矩阵只有 5 个非零元),全量遍历会严重低效。
用户常误以为“格式化”只是控制小数位数,实际核心是行列坐标映射与空白填充策略。
- 先提取所有唯一行号,排序后遍历每行
r - 对该行,收集所有
col满足triplet.row == r的三元组,再按col排序 - 用变量
lastCol = -1记录上一个输出列,中间空缺列补" "或"0 " - 推荐使用
std::set<int></int>存每行的列索引,避免手写二分查找;但若性能敏感,可用vector<pair double>></pair>+sort() - 不要用
printf("%6.2f", val)直接打,先转字符串再 pad,否则负数宽度错乱
求和运算容易忽略的边界与内存管理问题
多个稀疏矩阵连续相加(如 A + B + C)时,若复用同一块内存做中间结果,极易因容量不足导致越界写或静默截断。更隐蔽的问题是:不同矩阵行列维度不一致时,加法本身无定义,但三元组层面可能只报“结果为空”,不提示维度冲突。
实操中最常踩的坑不是算法错,而是没校验输入合法性。
- 每次加法前必须检查:
matA.rows == matB.rows && matA.cols == matB.cols - 动态分配的三元组数组(如
Trip* data)在多次加法中要delete[]前次结果,否则内存泄漏 - 若用
std::vector<trip></trip>,注意reserve()避免频繁 realloc;但shrink_to_fit()不保证释放内存 - 浮点数求和要考虑精度:
fabs(a + b) 应作零处理,否则转置或输出时带极小非零值 - 读取文件构造三元组时,行列索引是否从 0 开始?很多数据集从 1 开始,不减 1 会导致越界访问
真正麻烦的从来不是三元组归并或转置逻辑,而是维度校验、内存生命周期、以及浮点零值判定这些看似琐碎却决定结果可靠性的环节。
C++免费学习笔记(深入):立即使用
在学习笔记中,你将探索 C++ 的入门与实战技巧!











