大家好!这一篇要写的是稀疏矩阵的快速转置和快速相乘的操作,我们都知道,在稀疏矩阵中对零元素进行存储和算术运算没什么必要,所以我们要减小这两个操作的复杂度,比如矩阵相乘的传统的经典算法中:
for(i=1;i<=M.行;++i)
for(j=1;j<=N.列;++j)
{
Q[i][j]=0;
for(k=1;k<=M.列;++k)
Q[i][j]+=M[i][k]*N[k][j];
}
光是Q[i][j]+=M[i][k]*N[k][j];这一句就要M的行数*N的列数再乘以M的列数这样一个复杂度,所以我们有了以下的算法,在快速相乘中我会写比较多的注释方便大家理解,具体详细的说明请各位参考严的数据结构,下面有三个文件,一个是Matrix.h,一个是Matrix.c,最后一个是main.c
//Matrix.h
#ifndef MATRIX_H
#define MATRIX_H
#include <stdio.h>
#include <stdlib.h>
//#include <windows.h>
#define TRUE 1
#define FALSE 0
#define OK 1
#define ERROR 0
#define INFEASIBLE -1
#define OVERFLOW -2
#define MAXSIZE 7
#define MAXRC 7
#define ROW 4
#define COL 5
#define FORMAT "%3d"
#define TURNLINE printf("/n/n")
typedef int ElemType;
typedef int status;
typedef struct
{
int i,j;//行,列
ElemType e;//非零元素
}Triple;
typedef struct
{
Triple data[MAXSIZE+1];//非零元三元组表,data[0]不用
int mu,nu,tu;//矩阵的行数,列数及非零元个数
int rpos[MAXRC+1];//用于行逻辑链接的顺序表,即事先知道非零元在矩阵中每一行的位置,方便两个矩阵的快速相乘
}Matrix;
void test();
status createMatrix(Matrix *);//创建一个矩阵,这里留成空函数让大家自由发挥(因为这不是讨论的重点),下面自会有测试数据。
void printMatrix(Matrix M,int size);//打印矩阵
status findPostion(Matrix *);//查找矩阵中