稀疏矩阵乘法.docVIP

稀疏矩阵乘法 #include stdio.h #include stdlib.h //错误（1） //宏变量定义:根据算法增加的语句 #define MAXSIZE 100 //定义Elemtype 和Status类型 #define MAXRC 10 typedef int ElemType; typedef int Status; //定义数据结构参见书上P98 int num[MAXSIZE]; int row; typedef struct { int i,j; ElemType e; }Triple; typedef struct { Triple data[MAXSIZE+1]; int rpos[MAXRC+1]; int mu,nu,tu; } RLSMatrix; void InitRLSMatrix(RLSMatrix M); void OutputRLSMatrix(RLSMatrix M);//按照三元组输出稀疏矩阵 void OutputRLSMatrix1(RLSMatrix M);//按照矩阵格式输出稀疏矩阵 Status MultSMaxtrix(RLSMatrix M,RLSMatrix T,RLSMatrix Q); void main() { RLSMatrix M,N,Q; printf(输入矩阵M（三元组形式）\n); InitRLSMatrix(M); printf(输入矩阵N（三元组形式）\n); InitRLSMatrix(N); MultSMaxtrix(M,N,Q); printf(矩阵Q（三元组形式）\n); OutputRLSMatrix(Q); printf(矩阵Q（矩阵格式）\n); OutputRLSMatrix1(Q); } void InitRLSMatrix(RLSMatrix M) { int m,n,t; int k; int i,j,e; printf(请输入稀疏矩阵的行数\n); scanf(%d,m); printf(请输入稀疏矩阵的列数\n); scanf(%d,n); printf(请输入稀疏矩阵的元素个数\n); scanf(%d,t); M.mu=m; M.nu=n; M.tu=t; printf(请依次输入稀疏矩阵的各个元素，输入格式：\n); printf(行号 列号 元素值\n); for(k=1;k=t;k++) { scanf(%d%d%d,i,j,e); M.data[k].i=i; M.data[k].j=j; M.data[k].e=e; } //计算Rpos for(row=1;row=M.mu;row++) num[row]=0; for(k=1;k=M.tu;k++) ++num[M.data[k].i]; M.rpos[1]=1; for(row=2;row=M.mu;row++) M.rpos[row]=M.rpos[row-1]+num[row-1]; //for(row=1;row=M.mu;row++) printf(%d,M.rpos[row]); } void OutputRLSMatrix(RLSMatrix M) { int k; printf(依次输出稀疏矩阵的各个元素\n); printf(行号 列号 元素值\n) ; for(k=1;kM.tu;k++) printf(%4d%4d%4d\n,M.data[k].i,M.data[k].j,M.data[k].e); } Status MultSMaxtrix(RLSMatrix M,RLSMatrix N,RLSMatrix Q) { int tp,t,ctemp[MAXSIZE];//注意这些变量定义 int arow,bcol,brow,ccol; int p,q; if(M.nu!=N.mu) return 0; Q.mu=M.mu ; Q.nu=N.nu; Q.tu=0; if(M.tu*N.tu !=0) { for(arow=1;arow=M.mu;arow++) { for(bcol=1;bcol=N.nu;bcol++) //将算法中的ctemp[]=0修改为可以执行的代码 ctemp[bcol]=0; // 当前行各元素累加器清零 Q.rpos[arow]=Q.tu+1; if(arowM.mu) tp=M.rpos[arow+1]; else tp=M.