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

1、#include#include #define MAXSIZE 12500/ 三元组结构typedef struct int i,j; int e;Triple;/矩阵结构/矩阵行下标和列下标/值typedef struct Triple 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,

2、Matrix* G);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 不用printf(“第%d个非零元的行号：“,i);/以下为数据初始化scanf(%d,&M-datai.i);

3、printf(“第%d个非零元的列号：,i); scanf(%d,&M-datai.j); printf(”第d个非零元的元素值：,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;的行、列数。/因为加减运算必须维数相等，所以M、T 行、列数相等，初始化第三方矩阵G-nu = M-nu;G-tu = 0; int g,m,t; g = m = t = 1;int k;int sum

4、;/g、 m、 t 分别记录三个矩阵的元素的存储位置 /从 data1 开始/和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-datag.e=sum;G-

5、tu+;g+;m+;t+;else / 矩阵相加为零时，不做存储操作m+;t+;结束 */ 行下标和下标不等的元素直接插入到 G 的相应位置开始 *while(M-datam.i = k) / 插入 M 剩余的元素 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) / 插入 T 剩余的元素G-datag.i = k;G-datag.j = T-datat.j;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;插入 M 剩余的元素G-datag.j = M-datam.j;G-datag.e = M-datam.e

8、;G-tu+; m+; g+;while(T-datat.i = k) / G-datag.i = k;插入 T 的元素G-datag.j = T-datat.j;G-datag.e = 0 -G-tu+;t+;g+;/G-tu=g;T-datat.e;/记录运算完成后，第三方矩阵非零元的个数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

9、numMAXSIZE;/ 记录矩阵的各行的第一个非零元的位置int m,t,g;m = t = g = 1;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

10、; numk+;M-rpos1 = 1;/第一行第一个非零元的位置在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 中各行非零元的个数for(trow = 1; trow mu; trow+)/清零numtrow = 0

11、;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);for(trow = 1;trow mu; trow+)printf(“第%d行的第一个非零元在data中的序号是：(d),trow,T-rpostrow);printf(n);结束 *for(mrow

12、 = 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)x = T-rpostrow+1;elsex = T-tu+1;for(q = T-rpostrow; qdataq.j;tmpgcol += M-dataw.e * T-dataq.e;/对应非零元累加for(gcol = 1

13、; gcol nu;gcol+)/压缩存储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)int k;printf( 非零元素数 %dn,M-tu); for(k = 1;k tu;k+) printf(%d,%d,%d)n,M-datak.i,M-datak.j,M-datak.e);/7、主函数int mai

14、n(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;printf( 请输入第一个矩阵的行数和列数： ); scanf(%d %d,&(M-mu),&(M-nu);printf( 非零元的个数： ); scanf(%d,&(M-tu);Init(M); printf(n);printf( 请输入第二个矩阵的行数和列数： ); scanf(%d %d,&(T-mu),&(T-nu); print