三元组实现稀疏矩阵加减乘

1、#define MAXSIZE 12500 / 三元组结构typedef structint e; / Triple;/ 矩阵结构typedef structTriple dataMAXSIZE+1;int rposMAXSIZE+1; / 这是存放各行第一非零元在矩阵中的位置int mu,nu,tu; / 矩阵的行数、列数、非零元个数Matrix;void Init(Matrix* M);void Add(Matrix* M,Matrix* T,Matrix* G);void Jian(Matrix* M,Matrix* T,Matrix* G);int i,j; /矩阵行下标和列下标pri

2、ntf(第4 个非零元的列号：,i);void Cheng(Matrix* M,Matrix* T,Matrix* G);void Cheng(Matrix* M,Matrix* T,Matrix* G);void PrintMatrix(Matrix* M);/2 、初始化矩阵void Init(Matrix* M)int i;if(M-mu nu tu M-mu*M-nu)printf( 出错 !n); / 如果矩阵的行数、列数不符合要求，打印出错for(i=1;itu;i+) /data0 不用prin tf(第4 个非零元的行号：,i);/以下为数据初始化scanf(%d,&M

3、-datai.i);scanf(%d,&M-datai.j);printf(”第4 个非零元的元素值：,i);scanf(%d,&M-datai.e);printf(n);printf( 您创建的矩阵如下 :n);PrintMatrix(M);/3 、矩阵相加void Add(Matrix* M,Matrix* T,Matrix* G)G-mu = M-mu;/初始化第三方矩阵的行、列数。G-nu = M-nu;G-tu = 0;int g,m,t;/g、m、t 分别记录三个矩阵的元素的存储位置g = m = t = 1;/从 data1 开始因为加减运算必须维数相等，所以 M

4、、T 行、列数相等，printf(第4 个非零元的列号：,i);int k;int sum; / 和if(M-mu != T-mu | M-nu != T-nu)printf( 维数不等 !n);for(k = 1;k mu;k+)/ 行下标和下标相等的对应元素相加开始if(M-datam.i = k & T-datat.i = k) /if(M-datam.j = T-datat.j ) /sum=M-datam.e+T-datat.e; /行下标相等时列下标相等时矩阵向加if(sum != 0)/相加不为零时G-datag.i=k;G-datag.j=M-datam.j;G-dat

5、ag.e=sum;G-tu+;g+;m+;t+;else / 矩阵相加为零时，不做存储操作m+;t+;*/行下标和下标不等的元素直接插入到G 的相应位置while(M-datam.i = k) /G-datag.i = k;G-datag.j = M-datam.j;G-datag.e = M-datam.e;G-tu+;m+;g+;while(T-datat.i = k) /G-datag.i = k;G-datag.j = T-datat.j;开始插入 M 剩余的元素插入 T 剩余的元素G-datag.e = T-datat.e;G-tu+;t+;g+;/G-tu=g; /记录运算完成后，

6、第三方矩阵非零元的个数/*结束 */printf(n);printf( 相加后的矩阵如下 :n);PrintMatrix(G);/4 、矩阵相减void Jian(Matrix* M,Matrix* T,Matrix* G)G-mu = M-mu;G-nu = M-nu;G-tu = 0;int sum,k;int m,t,g;m = t = g = 1;if(M-mu != T-mu | M-nu != T-nu)printf( 维数不等 !n);for(k = 1;k mu;k+)/ 行下标和下标相等的对应元素相减if(M-datam.i = k & T-datat.i = k)

7、/if(M-datam.j = T-datat.j ) /开始行下标相等时列下标相等时sum=M-datam.e-T-datat.e; / 矩阵向减 if(sum != 0)/ 相加不为零时G-datag.i=k;G-datag.j=M-datam.j;G-datag.e=sum;G-tu+;g+;m+;t+;elsem+;t+;/行下标和下标不等的元素直接插入到G 的相应位置while(M-datam.i = k) /G-datag.i = k;G-datag.j = M-datam.j;G-datag.e = M-datam.e;G-tu+;m+;g+;while(T-datat.i =

8、k) /结束开始插入 M 剩余的元素插入 T 的元素G-datag.i = k;G-datag.j = T-datat.j;G-datag.e = 0 - T-datat.e;G-tu+;t+;g+;/G-tu=g; /*printf(n);printf( 相减后的矩阵如下 :n);PrintMatrix(G);/5 、矩阵相乘void Cheng(Matrix* M,Matrix* T,Matrix* G) int k,w,q,x,mrow,trow,gcol;记录运算完成后，第三方矩阵非零元的个数结束 */int tmpMAXSIZE; /累加器int numMAXSIZE; /记录矩阵的

9、各行的第一个非零元的位置int m,t,g; m = t = g w =1;if(M-nu != T-mu) /如果第一个矩阵的列数不等于第二个矩阵的行数printf( 这两个矩阵不能相乘 !n);/ 初始化G-mu = M-mu;G-nu = T-nu;G-tu = 0;if(M-tu * T-tu != 0)/ 记录各个矩阵的各行的第一个非零元的位置/*/记录矩阵 M 中各行非零元的个数for(mrow = 1; mrow mu; mrow+) /清零nummrow = 0;for(w = 1;w tu;w+)k = M-dataw.i;numk+;开始for(trow = 1; trow

10、 mu; trow+) /清零第一行第一个非零元的位置在 data 的第一位for(mrow = 2; mrow mu; mrow+) / 第 N 行第一个非零元的位置M-rposmrow = M-rposmrow-1 + nummrow-1;printf(n);printf( 矩阵 Mn);for(mrow = 1;mrow mu; mrow+)printf( 第 %d 行 的 第 一 个 非 零 元 在 data 中 的 序 号 是:(d),mrow,M-rposmrow);printf(n);/ 记录矩阵 T 中各行非零元的个数numtrow = 0;M-rpos1 = 1; /for(

11、trow = 1;trow mu; trow+)for(w = 1;w tu;w+)k = T-dataw.i;numk+;T-rpos1 = 1; / 第一行第一个非零元的位置在 data 的第一位for(trow = 2; trow mu; trow+) /第 N 行第一个非零元的位置T-rpostrow = T-rpostrow-1 + numtrow-1;printf(n);printf(n);printf(矩阵 Tn);printf(第 %d 行 的 第 一 个 非 零 元 在 data 中 的 序 号是:(d),trow,T-rpostrow);printf(n);I*for(mr

12、ow = 1; mrow mu; mrow+)for(k = 1; k nu;k+) /当前行各元素累加器清零tmpk = 0;G-rposmrow = G-tu + 1; /矩阵 G 的各行的第一个非零元在G-data 中的位置for(w = M-rposmrow;w rposmrow+1;w+)trow = M-dataw.j; /找到对应 T 的行号if(trow mu)结束for(gcol = 1; gcol nu;gcol+) /压缩存储x = T-rpostrow+1;elsex = T-tu+1;for(q = T-rpostrow; qdataq.j;tmpgcol += M-

13、dataw.e * T-dataq.e; / 对应非零元累加if(tmpgcol)G-dataG-tu+1.i = mrow;G-dataG-tu+1.j = gcol;G-dataG-tu+1.e = tmpgcol;G-tu+;printf(n);printf( 相乘后的矩阵如下 :n);PrintMatrix(G);/6 、矩阵打印void PrintMatrix(Matrix* M)printf(请输入第一个矩阵的行数和列数： );int k;printf( 非零元素数 %dn,M-tu);for(k = 1;k tu;k+)printf(%d,%d,%d)n,M-datak.i,M-

14、datak.j,M-datak.e);/7 、主函数int main(int argc, char* argv)Matrix* M =(Matrix*)malloc(sizeof(Matrix);Matrix* T =(Matrix*)malloc(sizeof(Matrix);Matrix* G =(Matrix*)malloc(sizeof(Matrix);int tag = 1;int n;scanf(%d %d,&(M-mu),&(M-nu);printf( 非零元的个数： );scanf(%d,&(M-tu);Init(M);printf(n);printf( 请输入第二个矩阵的行数和列数： );scanf(%d %d,&(T-mu),&(T-nu);printf( 非零元