数据结构--稀疏矩阵的快速转置及快速相乘操作

      大家好！这一篇要写的是稀疏矩阵的快速转置和快速相乘的操作，我们都知道，在稀疏矩阵中对零元素进行存储和算术运算没什么必要，所以我们要减小这两个操作的复杂度，比如矩阵相乘的传统的经典算法中：

 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 *);//查找矩阵中