如何快速查询文件中存储的矩阵索引对应的数值？

方案1：内存预构建查询索引（绝大多数学术场景首选，实现简单性能最高） #

你的需求本质是COO格式稀疏矩阵的点查询，只要三元组总规模没有超过内存上限（常规PC可轻松应对千万级三元组），优先采用本方案：

• 第一步：程序初始化阶段一次性读取所有三元组数据，根据矩阵维度大小选择对应存储结构：
  • 若矩阵最大行列索引≤10^4（总元素数≤1e8）：直接初始化二维数组，默认值全为0，读取三元组时给对应位置赋值，查询时直接取arr[c][d]即可，查询耗时为常数级。
  • 若矩阵维度极大、非零元素占比极低：使用哈希映射存储，将行列索引拼接为唯一键，非零值作为值存入。以C++为例，可以用std::unordered_map，参考实现逻辑：

#include <unordered_map>
#include <fstream>
#include <string>

// 定义哈希键，将两个int拼接为long long
inline uint64_t get_key(int row, int col) {
    return (static_cast<uint64_t>(row) << 32) | static_cast<uint32_t>(col);
}

std::unordered_map<uint64_t, int> mat_map;

// 预加载阶段执行一次
void load_mat(const std::string& file_path) {
    std::ifstream fin(file_path);
    int r, c, v;
    while (fin >> r >> c >> v) {
        // 可选：如果0值不需要存储可以加判断 if(v != 0) 再插入，节省空间
        mat_map[get_key(r, c)] = v;
    }
    fin.close();
}

// 查询阶段直接调用
int get_val(int c, int d) {
    auto it = mat_map.find(get_key(c, d));
    return it == mat_map.end() ? 0 : it->second;
}

• 第二步：在你的嵌套循环中直接调用get_val(c,d)即可，不需要再操作文件指针，查询平均复杂度为O(1)，相比每次遍历文件性能提升至少1000倍以上。

方案2：外存二分查询（仅适用于三元组规模超内存上限的极端场景） #

如果三元组总量超过10亿条、内存无法完全加载，可以提前对文件做预处理：

1. 离线将三元组文件按「行索引升序、列索引升序」排序，生成新的有序三元组文件
2. 查询时采用二分查找定位对应行的起始位置，再在该行内二分查找对应列的值，查询复杂度为O(logM)，M为三元组总条数，远低于全文件遍历的O(M)复杂度。

内容的提问来源于stack exchange，提问作者Safayet Hossain Sobuj