高手帮忙做下“两个稀疏矩阵相乘Q=M*N过程求精”

Question

高手帮忙做下“两个稀疏矩阵相乘Q=M*N过程求精”

Answer 1

#include
#include
#define OK 1
#define ERROR 0
#define MAXSIZE 100 //最多非0元素的个数
#define MAXR 50 //rpos所能处理的最大行数
#define MAXC 50 //系数矩阵相乘时,保留临时列结果的数组temp[MAXC]
typedef struct NODE{ //定义稀疏矩阵结点
int i;
int j;
int data;
} Node;
typedef struct MATRIX{ //定义稀疏矩阵（可以快速访问）
int mu,nu,tu;
Node matrix[MAXSIZE+1];
int rpos[MAXR+1];
} Matrix;
int CreatSMatrix( Matrix* M ); //创建一个矩阵（由用户输入原始矩阵,转化为稀疏矩阵方式储存）
int Print( Matrix M ); //打印一个稀疏矩阵
int Mul_SMatrix( Matrix M,Matrix N,Matrix *Q); //两个稀疏矩阵相乘
main(){
Matrix A1,A2,A3; //定义矩阵
CreatSMatrix( &A1 );
CreatSMatrix( &A2 );
if( A1.nu==A2.mu ){ //判断能否相乘
Mul_SMatrix( A1,A2,&A3 );
printf("两矩阵相乘得：\n"); Print(A3);
}
system("pause");
return 0;
}
//稀疏矩阵相乘
int Mul_SMatrix( Matrix M,Matrix N,Matrix *Q)
{
int i,Mj;
int arow,Mlim,Nlim,Mcol,Nrow;
int ctemp[MAXC];
Q->tu=0;
//初始化Q
Q->mu=M.mu; Q->nu=M.nu;
if(M.tu*N.tu!=0){
//非零矩阵
for(arow=1; arowtu+1; //给Q->rpos[]数组赋值
Mlim = arowtu].i = arow;
Q->matrix[Q->tu].j = i;
Q->matrix[Q->tu].data = ctemp[i];
}
}
}
}
return 1;
}
//构建稀疏矩阵
int CreatSMatrix( Matrix* M ){
int temp,i,j;
printf("输入矩阵的行列数：");
scanf("%d%d",&M->mu,&M->nu);
M->tu=0;
printf("按行序输入矩阵：\n");
for( i=1; imu; i++){
M->rpos[i]=M->tu+1; //每计算完一行,给rpos[]赋值
for( j=1; jnu; j++){
scanf("%d",&temp );
if( temp ){ //非0值保存
M->matrix[M->tu+1].i= i;
M->matrix[M->tu+1].j= j;
M->matrix[M->tu+1].data=temp;
M->tu++;
}
}
}
return OK;
}
//打印稀疏矩阵
int Print( Matrix M){
int i;
if(M.tu==0){
printf("空矩阵\n\n");
return ERROR;
}
printf("i\tj\tdata\n");
for( i=1; i