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

#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,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); 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分别记录三个矩阵的元素的存储位置 g = m = t = 1; /从data1开始 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-datag.e=sum; G-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; /记录运算完成后，第三方矩阵非零元的个数 /*结束*/ 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) /行下标相等时 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+; 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 = 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; /记录矩阵的各行的第一个非零元的位置 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; 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; 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 = 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; else x = T-tu+1; for(q = T-rpostrow; qdataq.j; tmpgcol += M-dataw.e * T-dataq.e; /对应非零元累加 for(gcol = 1; 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 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; 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); printf