做了一个矩阵压缩的东西,然后简单记录一下
矩阵输入
矩阵显式(压缩之后的和压缩之前的)
矩阵转置
矩阵乘法
只有矩阵乘法需要稍微说一下,由于矩阵要求按照列序为主存储,所以矩阵相乘的时候可以将左边的矩阵转置,然后直接左边的一列和右边的一列相乘然后存储。
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
// 三元组实现矩阵的压缩以及运算
#define OK 1
#define ERROR 0
typedef struct Element{
int row,col;
int val;
}element;
typedef struct Matrix {
// 矩阵
int row,col,num; // 矩阵的行数 、列数、 非零元素
element *base;
}matrix;
void MatShow(matrix *mat) {
int i = 0;
int ptr = 0,row = 0,col = 0;
{
// 显式压缩矩阵
for(int i=0; i<mat->num; i++) {
printf("%3d %3d %3d\n",mat->base[i].row,mat->base[i].col,mat->base[i].val);
}
}
{
// 显示原矩阵
printf("\nInit Matrix:\n");
printf("%d %d %d\n",mat->row,mat->col,mat->num);
int *p = (int*)malloc(sizeof(int) * (mat->row * mat->col));
//memset(p,0,(mat->row * mat->col));
for(int i=0; i<mat->row * mat->col; i++) p[i] = 0;
for(int i=0; i<mat->num; i++) {
int r = mat->base[i].row,c = mat->base[i].col,val = mat->base[i].val;
p[r * mat->col + c] = val;
}
for(int i=0; i<mat->row; i++) {
for(int j=0; j<mat->col; j++) {
printf("%d ",p[i * mat->col + j])