c++ 稀疏矩阵存储 c++如何实现csr或csc格式的稀疏矩阵-C++-PHP中文网

csr与csc的核心区别在于indptr的语义：csr中indptr[i]表示第i行首个非零元在values中的下标，适合按行访问；csc中indptr[j]表示第j列首个非零元下标，适合按列访问。

c++ 稀疏矩阵存储 c++如何实现csr或csc格式的稀疏矩阵

CSR 和 CSC 的核心区别在哪？别只看存储顺序

CSR（Compressed Sparse Row）和 CSC（Compressed Sparse Column）本质不是“行优先 vs 列优先”的简单类比，而是**访问模式决定结构**。CSR 适合按行遍历（比如矩阵-向量乘 y = A * x），CSC 适合按列遍历（比如 A.transpose() * x 或求解器中 LU 的列主元选主）。误用会导致缓存失效、循环嵌套反直觉、甚至逻辑错误——比如用 CSR 写一个按列求和的函数，性能可能差 5 倍以上。

二者共享三数组结构：values（非零值）、indices（对应行号或列号）、indptr（行/列起始偏移）。关键差异在 indptr 含义：
– CSR 中 indptr[i] 是第 i 行第一个非零元在 values 中的下标；
– CSC 中 indptr[j] 是第 j 列第一个非零元在 values 中的下标。

手写 CSR 类时，`indptr` 容易越界或漏填 1 个元素

indptr 长度必须是 rows + 1（CSR）或 cols + 1（CSC），最后一个元素恒为非零元总数 nnz。新手常犯两个错：
– 把 indptr 初始化成长度 rows，导致访问 indptr[rows] 时越界；
– 构造时只填了前 rows 个位置，忘了设 indptr[rows] = nnz，后续遍历某行会读到错误范围。

正确做法：先统计每行非零个数，做前缀和：

std::vector<int> indptr(rows + 1, 0);
for (int i = 0; i <p>注意：所有索引默认从 0 开始，别混用 1-based 的文献伪代码。</p><div class="aritcle_card flexRow artxards">
											<div class="artcardd flexRow">
												<a class="aritcle_card_img" rel="nofollow" href="/xiazai/skill5502" title="C++ Code Review Master"><img
														src="https://img.php.cn/upload/skill/000/000/081/179051228971575.jpg" alt="C++ Code Review Master" onerror="this.onerror='';this.src='/static/lhimages/moren/morentu.png'" ></a>
												<div class="aritcle_card_info flexColumn">
													<a rel="nofollow" href="/xiazai/skill5502" title="C++ Code Review Master" class="overflowclass">C++ Code Review Master</a>
													<p class="overflowclass">组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。</p>
												</div>
												<a rel="nofollow" href="/xiazai/skill5502" title="C++ Code Review Master" class="aritcle_card_btn flexRow flexcenter"><b></b><span>下载</span>
												</a>
											</div>
										</div>
<h3>插入新元素？别直接 push_back 到 <code>values</code> 和 <code>indices</code>
</h3>
<p>C++ 稀疏矩阵一旦用 CSR/CSC 格式构建完成，就**不支持高效随机插入**。因为插入会破坏 <code>indptr</code> 的连续性，还可能触发整行数据搬移。常见错误场景：<br>
– 读取 COO 格式（三元组）时边读边插，结果时间复杂度退化成 O(nnz²)；<br>
– 在已构造好的 CSR 上调用 <code>insert(i, j, val)</code>，内部反复 resize 向量。</p>
<p>正确路径只有两条：<br>
– 批量构建：先收集全部三元组，排序（按行→列对 CSR），再一次性构造三数组；<br>
– 用中间格式过渡：COO → 排序 → CSR。标准库没提供现成排序，得自己写：<br><code>std::sort(coo_triplets.begin(), coo_triplets.end(), [](const auto& a, const auto& b) { return a.row </code></p>
<h3>用 Eigen 或 Intel MKL？先确认你真需要它们</h3>
<p>Eigen 的 <code>SparseMatrix<double rowmajor></double></code> 默认就是 CSR，<code>ColMajor</code> 就是 CSC。但它封装太深：比如想直接访问 <code>indptr</code>，得调 <code>innerIndexPtr()</code>、<code>outerIndexPtr()</code>、<code>valuePtr()</code>，且这些指针只读；MKL 的 <code>sparse_matrix_t</code> 更是要求显式 convert，调试时看不到原始数组。</p>
<p>如果你只是做教学、小规模实验、或需深度定制（比如带符号压缩、混合精度索引），手写轻量 CSR 类更可控。示例骨架：</p>
<pre class="brush:php;toolbar:false;">struct CSRMatrix {
    std::vector<double> values;
    std::vector<int> indices;   // 行索引（CSR）或列索引（CSC）
    std::vector<int> indptr;     // 长度 rows+1（CSR）或 cols+1（CSC）
    int rows, cols, nnz;
};</int></int></double>

真正难的不是存数据，而是写对 matvec 和 transpose —— 这俩函数里 indptr 和 indices 的配合稍有偏差，结果就全错。建议先用小矩阵（3×4）手算一遍每步下标，再写代码验证。