稀疏矩阵的创建与相加.docVIP

//本程序实现稀疏矩阵的创建与相加可在VS2008环境中运行 #includestdio.h #includestdlib.h #define MAXSIZE 100//非零元素个数最大值为 #define OK 1 typedef int ElemType;//元素值的类型定义为int型 typedef int Status; typedef struct{ int i,j;//该非零元素的行下标和列下标 ElemType e;//存储该非零元素的值 }Triple; typedef struct{ Triple data[MAXSIZE+1];//非零元三元组表，data[0]未用 int mu,nu,tu;//矩阵的行数、列数和非零元个数 }TSMatrix; Status CreatSMatrix(TSMatrix M) {//创建稀疏矩阵M int i,m,n; ElemType e; printf(请输入矩阵的行数，列数和非零元素个数(可用空格隔开):); scanf(%d,M.mu);scanf(%d,M.nu);scanf(%d,M.tu); if(M.tuMAXSIZE) {printf(非零元素个数太多请重新输入\n);exit(0);} M.data[0].i=0;//为下面比较顺序作准备 for(i=1;i=M.tu;i++) { printf(请按行序输入第%d个非零元素所在行、列、元素的值:,i); scanf(%d,m);scanf(%d,n);scanf(%d,e); if(m1||mM.mu||n1||nM.nu)//行或列超出范围 {printf(行或列超出范围请重新输入\n);exit(0);} if(mM.data[i-1].i||m==M.data[i-1].in=M.data[i-1].j)//行或列的顺序有错 {printf(行或列的顺序有错请重新输入\n);exit(0);} M.data[i].i=m; M.data[i].j=n; M.data[i].e=e; } return OK; } Status comp(int c1,int c2) {//比较行列的大小AddSMatrix函数将用到 if(c1c2) return -1; if(c1==c2) return 0; return 1; } Status AddSMatrix(TSMatrix M,TSMatrix N,TSMatrix Q) {//求稀疏矩阵的和Q=M+N int m=1,n=1,q=0;//分别为矩阵M、N、Q的编号 if(M.mu!=N.mu||M.nu!=N.nu)//M、N两稀疏矩阵行或列数不同 {printf(不满足矩阵相加的条件请重新输入\n);exit(0);} Q.mu=M.mu;Q.nu=M.nu;//矩阵Q的行数和列数与矩阵M(矩阵N)相同 while(m=M.tun=N.tu) {//矩阵M和N的元素都没有处理完 switch(comp(M.data[m].i,N.data[n].i)){ case -1: Q.data[++q]=M.data[m++];//将矩阵M的当前元素值赋给矩阵Q break; case 0://M、N矩阵当前行元素相等继续比较 switch(comp(M.data[m].j,N.data[n].j)){ case -1: Q.data[++q]=M.data[m++]; break; case 0://M、N矩阵当前非零元素的行列均相等 Q.data[++q]=M.data[m++]; Q.data[q].e+=N.data[n++].e;//矩阵M、N当前元素求和并赋值给矩阵Q if(Q.data[q].e==0)//元素值为不存入压缩矩阵 q--; break; case 1: Q.data[++q]=N.data[n++];//将矩阵N当前元素赋值给矩阵Q } break; case 1: Q.data[++q]=N.data[n++]; } } while(m=M.tu)//矩阵M元素全部处理完毕 Q.data[++q]=M.data[m++]; while(n=N.tu)//矩阵N元素全部处理完毕 Q.data[++q]=N.data[n++]; Q.tu=q;//矩阵Q的非零元素个数 if(qMAXSIZE)//非零元素个