【数据结构】数组和字符串（五）：特殊矩阵的压缩存储：稀疏矩阵——压缩稀疏行（CSR）

news2024/11/20 9:27:28

文章目录

  • 4.2.1 矩阵的数组表示
  • 4.2.2 特殊矩阵的压缩存储
    • a. 对角矩阵的压缩存储
    • b~c. 三角、对称矩阵的压缩存储
    • d. 稀疏矩阵的压缩存储——三元组表
    • e. 压缩稀疏行（Compressed Sparse Row，CSR）矩阵
      • 结构体
      • 创建CSR矩阵
      • 元素设置
      • 初始化
      • 打印矩阵
      • 销毁CSR矩阵
      • 主函数
      • 代码整合

4.2.1 矩阵的数组表示

【数据结构】数组和字符串（一）：矩阵的数组表示

4.2.2 特殊矩阵的压缩存储

  矩阵是以按行优先次序将所有矩阵元素存放在一个一维数组中。但是对于特殊矩阵，如对称矩阵、三角矩阵、对角矩阵和稀疏矩阵等, 如果用这种方式存储，会出现大量存储空间存放重复信息或零元素的情况，这样会造成很大的空间浪费。为节约存储空间和算法（程序）运行时间，通常会采用压缩存储的方法。

  • 对角矩阵：指除了主对角线以外的元素都为零的矩阵，即对 任意 i ≠ j (1≤ i , j ≤n)，都有M(i, j)=0。由于只有主对角线上有非零元素，只需存储主对角线上的元素即可。
  • 三角矩阵：指上三角或下三角的元素都为零的矩阵。同样地，只需存储其中一部分非零元素，可以节省存储空间。
  • 对称矩阵：指矩阵中的元素关于主对角线对称的矩阵。由于对称矩阵的非零元素有一定的规律，可以只存储其中一部分元素，从而减少存储空间。
  • 稀疏矩阵：指大部分元素为零的矩阵。传统的按行优先次序存储方法会浪费大量空间来存储零元素，因此采用压缩存储的方法更为合适。常见的压缩存储方法有：压缩稠密行（CSR）、压缩稠密列（CSC）、坐标列表（COO）等。

a. 对角矩阵的压缩存储

【数据结构】数组和字符串（二）：特殊矩阵的压缩存储：对角矩阵——一维数组

b~c. 三角、对称矩阵的压缩存储

【数据结构】数组和字符串（三）：特殊矩阵的压缩存储：三角矩阵、对称矩阵——一维数组

d. 稀疏矩阵的压缩存储——三元组表

【数据结构】数组和字符串（四）：特殊矩阵的压缩存储：稀疏矩阵——三元组表

e. 压缩稀疏行（Compressed Sparse Row，CSR）矩阵

  压缩稀疏行（Compressed Sparse Row，CSR）是一种常用的稀疏矩阵存储格式。CSR存储格式通过压缩非零元素的行指针和列索引，以及存储非零元素的值，来有效地表示稀疏矩阵。它包含以下几个关键组成部分：

  • row_ptr（行指针数组）：它是一个长度为rows + 1的数组，用于存储每一行在col_indices和elements数组中的起始索引位置。row_ptr[i]表示第i行的元素在col_indices和elements数组中的起始位置，而row_ptr[i+1] - row_ptr[i]表示第i行的非零元素个数。
  • col_indices（列索引数组）：它是一个长度为num_elements的数组，用于存储每个非零元素对应的列索引。col_indices[i]表示第i个非零元素所在的列索引。
  • elements（元素数组）：它是一个长度为num_elements的数组，用于存储每个非零元素的值。elements[i]表示第i个非零元素的值。

  CSR存储格式的主要优点是有效地压缩了稀疏矩阵的存储空间，只存储非零元素及其对应的行和列信息。此外，CSR格式还支持高效的稀疏矩阵向量乘法和稀疏矩阵乘法等操作。

结构体

typedef struct {
    int row;
    int col;
    int value;
} Element;

typedef struct {
    int rows;
    int cols;
    int num_elements;
    Element* elements;
    int* row_ptr;
    int* col_indices;
} CSRMatrix;

  • Element 结构体表示矩阵中的一个元素，包含三个成员变量：row（行索引）、col（列索引）和 value（元素值）。

  • CSRMatrix 结构体表示一个CSR矩阵，包含了矩阵的行数 rows、列数 cols、非零元素的个数 num_elements，以及三个指针成员变量 elements、row_ptr 和 col_indices。

创建CSR矩阵

CSRMatrix createCSRMatrix(int rows, int cols, int num_elements) {
    CSRMatrix matrix;
    matrix.rows = rows;
    matrix.cols = cols;
    matrix.num_elements = num_elements;
    matrix.elements = (Element*)malloc(num_elements * sizeof(Element));
    matrix.row_ptr = (int*)malloc((rows + 1) * sizeof(int));
    matrix.col_indices = (int*)malloc(num_elements * sizeof(int));
    return matrix;
}
  • createCSRMatrix 函数用于创建一个CSR矩阵。
  • 接受矩阵的行数、列数和非零元素的个数作为参数，并返回创建的CSR矩阵。
  • 在函数内部，通过动态内存分配分别为 elements、row_ptr 和 col_indices 分配内存空间，并将 row_ptr 数组的所有元素初始化为0，最后返回创建的矩阵。

元素设置

void setElement(CSRMatrix* matrix, int row, int col, int value) {
    if (row < 0 || row >= matrix->rows) {
        printf("Invalid row index.\n");
        return;
    }
    int index = matrix->row_ptr[row];
    matrix->elements[index].row = row;
    matrix->elements[index].col = col;
    matrix->elements[index].value = value;
    matrix->col_indices[index] = col;
    matrix->row_ptr[row]++;  // 递增索引值
}

  setElement 函数可用于设置（修改）CSR矩阵中某个位置的元素值。

  • 接受一个指向CSR矩阵的指针 matrix，以及要设置的元素的行索引、列索引和值作为参数。
  • 在函数内部，首先检查行索引是否有效，如果无效则打印错误信息并返回。
  • 然后，根据行索引找到对应行的起始位置，将元素的行索引、列索引和值分别赋给对应的矩阵元素，并更新 col_indices 数组和 row_ptr 数组中的值。

初始化

void initializeCSRMatrix(CSRMatrix* matrix, int* values, int* row_indices, int* col_indices, int num_elements) {
    for (int i = 0; i < num_elements; i++) {
        matrix->elements[i].value = values[i];
        matrix->elements[i].row = row_indices[i];
        matrix->elements[i].col = col_indices[i];
        matrix->col_indices[i] = col_indices[i];
        matrix->row_ptr[row_indices[i]]++;
    }

    int sum = 0;
    for (int i = 0; i <= matrix->rows; i++) {
        int temp = matrix->row_ptr[i];
        matrix->row_ptr[i] = sum;
        sum += temp;
    }
}

  initializeCSRMatrix 函数用于初始化CSR矩阵的数据。

  • 接受一个指向CSR矩阵的指针 matrix，以及包含非零元素的值、行索引和列索引的数组，以及非零元素的个数作为参数。
  • 通过遍历非零元素数组，将值、行索引和列索引分别赋给对应的矩阵元素，并更新 col_indices 数组和 row_ptr 数组中的值。row_ptr 数组的每个元素表示对应行的非零元素在 elements 数组中的起始位置，通过累加非零元素的个数来计算每行的结束位置。

打印矩阵

void printCSRMatrix(CSRMatrix matrix) {
    printf("CSR Matrix:\n");
    printf("Rows: %d, Cols: %d, Num Elements: %d\n", matrix.rows, matrix.cols, matrix.num_elements);
    printf("Values: ");
    for (int i = 0; i < matrix.num_elements; i++) {
        printf("%d ", matrix.elements[i].value);
    }
    printf("\n");
    printf("Row Pointer: ");
    for (int i = 0; i <= matrix.rows; i++) {
        printf("%d ", matrix.row_ptr[i]);
    }
    printf("\n");
    printf("Column Indices: ");
    for (int i = 0; i < matrix.num_elements; i++) {
        printf("%d ", matrix.col_indices[i]);
    }
    printf("\n");
}
void printMatrixForm(CSRMatrix matrix) {
    printf("Matrix Form:\n");
    for (int i = 0; i < matrix.rows; i++) {
        for (int j = 0; j < matrix.cols; j++) {
            int value = 0;
            for (int k = matrix.row_ptr[i]; k < matrix.row_ptr[i + 1]; k++) {
                if (matrix.elements[k].col == j) {
                    value = matrix.elements[k].value;
                    break;
                }
            }
            printf("%d ", value);
        }
        printf("\n");
    }
}
  • printCSRMatrix 函数用于打印CSR矩阵的信息：它接受一个CSR矩阵作为参数，并打印矩阵的行数、列数、非零元素的个数以及 elements、row_ptr 和 col_indices 数组的内容。
  • printMatrixForm 函数用于以矩阵形式打印CSR矩阵。它接受一个CSR矩阵作为参数，并按矩阵的行数和列数遍历矩阵元素，通过遍历 row_ptr 数组和 col_indices 数组来获取每个位置的元素值，并打印出矩阵的形式。

销毁CSR矩阵

void destroyCSRMatrix(CSRMatrix* matrix) {
    free(matrix->elements);
    free(matrix->row_ptr);
    free(matrix->col_indices);
    matrix->elements = NULL;
    matrix->row_ptr = NULL;
    matrix->col_indices = NULL;
}

主函数

int main() {
    int rows = 9;
    int cols = 3;
    int num_elements = 5;

    CSRMatrix matrix = createCSRMatrix(rows, cols, num_elements);

    int col_indices[] = {0, 0, 0, 0, 1};
    int row_indices[] = {3, 5, 7, 8, 7};
    int values[] =      {2, 1, 3, 1, 4};

    initializeCSRMatrix(&matrix, values, row_indices, col_indices, num_elements);

    printCSRMatrix(matrix);
    printMatrixForm(matrix);

    destroyCSRMatrix(&matrix);
    return 0;
}

在这里插入图片描述

代码整合

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

typedef struct {
    int row;
    int col;
    int value;
} Element;

typedef struct {
    int rows;
    int cols;
    int num_elements;
    Element* elements;
    int* row_ptr;
    int* col_indices;
} CSRMatrix;

CSRMatrix createCSRMatrix(int rows, int cols, int num_elements) {
    CSRMatrix matrix;
    matrix.rows = rows;
    matrix.cols = cols;
    matrix.num_elements = num_elements;
    matrix.elements = (Element*)malloc(num_elements * sizeof(Element));
    matrix.row_ptr = (int*)malloc((rows + 1) * sizeof(int));
    matrix.col_indices = (int*)malloc(num_elements * sizeof(int));

    memset(matrix.row_ptr, 0, (rows + 1) * sizeof(int));

    return matrix;
}
void initializeCSRMatrix(CSRMatrix* matrix, int* values, int* row_indices, int* col_indices, int num_elements) {
    for (int i = 0; i < num_elements; i++) {
        matrix->elements[i].value = values[i];
        matrix->elements[i].row = row_indices[i];
        matrix->elements[i].col = col_indices[i];
        matrix->col_indices[i] = col_indices[i];
        matrix->row_ptr[row_indices[i]]++;
    }


    int sum = 0;
    for (int i = 0; i <= matrix->rows; i++) {
        int temp = matrix->row_ptr[i];
        matrix->row_ptr[i] = sum;
        sum += temp;
    }
}

void setElement(CSRMatrix* matrix, int row, int col, int value) {
    if (row < 0 || row >= matrix->rows) {
        printf("Invalid row index.\n");
        return;
    }
    int index = matrix->row_ptr[row];
    matrix->elements[index].row = row;
    matrix->elements[index].col = col;
    matrix->elements[index].value = value;
    matrix->col_indices[index] = col;
    matrix->row_ptr[row]++;  // 递增索引值
}

void printCSRMatrix(CSRMatrix matrix) {
    printf("CSR Matrix:\n");
    printf("Rows: %d, Cols: %d, Num Elements: %d\n", matrix.rows, matrix.cols, matrix.num_elements);
    printf("Values: ");
    for (int i = 0; i < matrix.num_elements; i++) {
        printf("%d ", matrix.elements[i].value);
    }
    printf("\n");
    printf("Row Pointer: ");
    for (int i = 0; i <= matrix.rows; i++) {
        printf("%d ", matrix.row_ptr[i]);
    }
    printf("\n");
    printf("Column Indices: ");
    for (int i = 0; i < matrix.num_elements; i++) {
        printf("%d ", matrix.col_indices[i]);
    }
    printf("\n");
}
void printMatrixForm(CSRMatrix matrix) {
    printf("Matrix Form:\n");
    for (int i = 0; i < matrix.rows; i++) {
        for (int j = 0; j < matrix.cols; j++) {
            int value = 0;
            for (int k = matrix.row_ptr[i]; k < matrix.row_ptr[i + 1]; k++) {
                if (matrix.elements[k].col == j) {
                    value = matrix.elements[k].value;
                    break;
                }
            }
            printf("%d ", value);
        }
        printf("\n");
    }
}

void destroyCSRMatrix(CSRMatrix* matrix) {
    free(matrix->elements);
    free(matrix->row_ptr);
    free(matrix->col_indices);
    matrix->elements = NULL;
    matrix->row_ptr = NULL;
    matrix->col_indices = NULL;
}

int main() {
    int rows = 9;
    int cols = 3;
    int num_elements = 5;

    CSRMatrix matrix = createCSRMatrix(rows, cols, num_elements);

    int col_indices[] = {0, 0, 0, 0, 1};
    int row_indices[] = {3, 5, 7, 8, 7};
    int values[] =      {2, 1, 3, 1, 4};

    initializeCSRMatrix(&matrix, values, row_indices, col_indices, num_elements);

    printCSRMatrix(matrix);
    printMatrixForm(matrix);

    destroyCSRMatrix(&matrix);
    return 0;
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/1137867.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

位图算法经典剖析（未完）

位图算法经典剖析（未完）

第一题 只出现一次的数字|| 给你一个整数数组 nums &#xff0c;除某个元素仅出现 一次 外&#xff0c;其余每个元素都恰出现 三次 。请你找出并返回那个只出现了一次的元素。 力扣&#xff08;LeetCode&#xff09;官网 - 全球极客挚爱的技术成长平台 class Solution { publ…
阅读更多...
35二叉树-树的最小深度

35二叉树-树的最小深度

目录 LeetCode之路——111. 二叉树的最小深度 分析 解法一&#xff1a;广度优先查询 解法二&#xff1a;深度优先查询 LeetCode之路——111. 二叉树的最小深度 给定一个二叉树&#xff0c;找出其最小深度。 最小深度是从根节点到最近叶子节点的最短路径上的节点数量。 说…
阅读更多...
基于STM32与ESP8266的太空人WiFi天气时钟（代码开源）

基于STM32与ESP8266的太空人WiFi天气时钟（代码开源）

前言&#xff1a;本文为手把手教学ESP8266著名开源项目——太空人WiFi天气时钟&#xff0c;不同的是本次项目采用的是STM32作为MCU。两者开发过程中有因为各自芯片的特点&#xff08;时钟频率&#xff0c;内存大小等&#xff09;&#xff0c;导致开发程序大不相同&#xff0c;很…
阅读更多...
Redis 持久化配置，RDB和AOF方式配置说明

Redis 持久化配置，RDB和AOF方式配置说明

文章目录 一、概述二、RDB 持久化方式配置三、AOF 持久化方式配置 如果您对Redis的了解不够深入请关注本栏目&#xff0c;本栏目包括Redis安装&#xff0c;Redis配置文件说明&#xff0c;Redis命令和数据类型说明。 一、概述 Redis 持久化是指将 Redis 数据存储到磁盘上&#x…
阅读更多...
KVM实验之动态迁移

KVM实验之动态迁移

前言 一台NFS&#xff08;192.168.184.132&#xff09; 一台KVM-a&#xff08;192.168.184.133&#xff09; 一台KVM-b&#xff08;192.168.184.134&#xff09; NFS配置&#xff1a; [rootlocalhost ~]# setenforce 0 //关闭selinux [rootlocalhost ~]# service iptables st…
阅读更多...
python re 匹配所有字段名称相同的值

python re 匹配所有字段名称相同的值

import retext {"code": 200,"message": "success","traceId": "da0b668c-4d67-44bf-907f-c072fc63839a","data": {"list": [{"articleId": 121862102,"title": "python 目录…
阅读更多...
k8s 金丝雀发布与声明式管理

k8s 金丝雀发布与声明式管理

Deployment控制器支持自定义控制更新过程中的滚动节奏&#xff0c;如“暂停(pause)”或“继续(resume)”更新操作。比如等待第一批新的Pod资源创建完成后立即暂停更新过程&#xff0c;此时&#xff0c;仅存在一部分新版本的应用&#xff0c;主体部分还是旧的版本。然后&#xf…
阅读更多...
用图说话——流程图进阶

用图说话——流程图进阶

目录 一、基本流程图 二、时序流程图 一、基本流程图 经常阅读歪果仁绘制的流程图&#xff0c;感觉比较规范&#xff0c;自己在工作中也尝试用他们思维来绘图&#xff0c;这是一个小栗子&#xff1a; 二、时序流程图 在进行Detail设计过程中&#xff0c;一般的绘图软件显得…
阅读更多...
微信小程序笔记功能（富文本editor功能）开发

微信小程序笔记功能（富文本editor功能）开发

一、效果图展示 模拟器及pc 端效果图 手机端就不贴了 二、把官方示例copy到自己的项目传送门 三、改一改上传图片功能 insertImage() {const that thiswx.chooseImage({// count:1 一次选择图片的个数success: function (response) {// 多张图片上传response.tempFilePaths…
阅读更多...
基于opencv的selenium滑动验证码的实现

基于opencv的selenium滑动验证码的实现

这篇文章主要介绍了基于opencv的selenium滑动验证码的实现&#xff0c;文中通过示例代码介绍的非常详细&#xff0c;对大家的学习或者工作具有一定的参考学习价值&#xff0c;需要的朋友们下面随着小编来一起学习学习吧 基于selenium进行动作链 由于最近很多人聊到滑动验证码…
阅读更多...
Python通过pyecharts对爬虫房地产数据进行数据可视化分析（一）

Python通过pyecharts对爬虫房地产数据进行数据可视化分析（一）

一、背景 对Python通过代理使用多线程爬取安居客二手房数据&#xff08;二&#xff09;中爬取的房地产数据进行数据分析与可视化展示 我们爬取到的房产数据&#xff0c;主要是武汉二手房的房源信息&#xff0c;主要包括了待售房源的户型、面积、朝向、楼层、建筑年份、小区名称…
阅读更多...
图像去噪滤波算法汇总(Python)

图像去噪滤波算法汇总(Python)

前言 上篇文章&#xff1a;图像数据噪音种类以及Python生成对应噪音&#xff0c;汇总了常见的图片噪音以及噪音生成方法&#xff0c;主要用在数据增强上面&#xff0c;作为数据集填充的方式&#xff0c;可以避免模型过拟合。想要了解图像数据增强算法的可以去看本人所撰这篇文…
阅读更多...
IO,库-10.24.25

IO,库-10.24.25

库-10.24.25 一、概念 头文件&#xff1a;.h&#xff1a;函数声明&#xff0c;结构体定义&#xff0c;宏定义&#xff0c;外部引用&#xff0c;重定义&#xff0c;条件编译 #include <>:从系统路径&#xff08;/usr/include&#xff09;下查找 #include " "&am…
阅读更多...
如何将Mysql数据库的表导出并导入到另外的架构

如何将Mysql数据库的表导出并导入到另外的架构

如何将Mysql数据库的表导出并导入到另外的架构 准备一、解决方法1.右键->导出->用mysqldump导出2.注意路径一般为&#xff1a;C:/Program Files/MySQL/MySQL Server 8.0/bin/mysqldump.exe和导出的sql文件位置3.右键->SQL脚本->运行SQL脚本4.找到SQL脚本并点击确定…
阅读更多...
[moeCTF 2023] REV

[moeCTF 2023] REV

逆向这东西&#xff0c;不太好说。 base64 这是个pyc文件&#xff08;python编译后的字节码文件&#xff09;&#xff0c;这东西可以直接用各种方法反编译。也可以不弄&#xff0c;必竟这应该签到级别的。用notepad打开&#xff0c;可以看到base64的编译和两个码表。显然猜是…
阅读更多...
mac系统u盘启动盘制作教程，更新至macOS Sonoma 14

mac系统u盘启动盘制作教程，更新至macOS Sonoma 14

mac系统怎么制作装系统的u盘,如果您要在多台电脑上安装 macOS&#xff0c;而又不想每次都下载安装器&#xff0c;这时可引导安装器就会很有用。一起来看苹果电脑u盘启动盘制作教程吧。 Macos系统安装包合集包揽macos 10.15&#xff0c;macos 11和苹果最新系统等多个版本 1、A…
阅读更多...
Windows 和 Linux 这2个系统在进行编程实现的时候的一些区别：

Windows 和 Linux 这2个系统在进行编程实现的时候的一些区别：

很惭愧&#xff0c;学了很多年才意识到&#xff0c;噢&#xff0c;原来这两个系统实现一些功能的时候会使用到不同的库&#xff0c;使用不同的函数。 那么&#xff0c;也会延伸出一些问题&#xff1a; 比如&#xff0c;如何实现版本的迁移。一个在Linux上运行的代码如何可以比…
阅读更多...
Instant-NGP中的多分辨率哈希编码

Instant-NGP中的多分辨率哈希编码

Instant-NGP的出现&#xff0c;无疑给神经表达领域带来了新的生命力。可认为是NeRF诞生以来的关键里程碑了。首次让我们看到了秒级的重建、毫秒级的渲染的NeRF工作。 作为如此顶到爆的工作&#xff0c;Instant-NGP毫无疑问斩获SIGGRAPH 2022的最佳论文。虽然只是五篇最佳论文之…
阅读更多...
微信公众号推送封面图制作：专业技巧大揭秘

微信公众号推送封面图制作：专业技巧大揭秘

在微信公众号推送中&#xff0c;一个吸引人的封面设计是吸引读者点击的重要因素。本文将指导你如何使用免费在线海报制作工具&#xff0c;如乔拓云&#xff0c;来制作一个适合节日的海报模板&#xff0c;并编辑文本、图片以及调整字体、颜色、布局等&#xff0c;最后导出并保存…
阅读更多...
【Linux前篇 】VMWare虚拟机安装与环境配置及远程连接 —— windows版

【Linux前篇 】VMWare虚拟机安装与环境配置及远程连接 —— windows版

目录 一、操作系统 1.1 什么是操作系统 1.2 常见操作系统 1.3 个人版本和服务器版本的区别 1.4 Linux的各个版本 二、VMWare 虚拟机安装配置流程 2.1 安装 2.2 配置虚拟网络编辑器 三、安装配置 Windows Server 2012 R2 3.1 创建虚拟机 3.2 安装 Windows Server 2012…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023