稀疏矩阵的压缩存储及运算.docxVIP

稀疏矩阵的压缩存储及运算 一、 实验内容 实现稀疏矩阵的压缩存储方法以及在特定存储方法下的基本运算 二、 实验母的 掌握数组的应用，包括稀疏矩阵、特殊矩阵的压缩存储方法。 矩阵的基本运算的实现， 包括 矩阵相加、转置、乘法等。 三、 问题描述 用行逻辑链接顺序表和十字链表分别实现稀疏矩阵的压缩存储 编程实现矩阵的转置运算和乘法运算 (运用行逻辑链接顺序表或十字链表作为存储结构) 四、 问题的实现 稀疏矩阵的抽象数据类型定义： ADT SpareseMatrix( 数据对象：D ={.日 |i =1,2，,m; j =1,2, , n; ai产ElemSet m和n分别称为矩阵的行数和列数} 数据关系： R ={ Row, Col} Row={ai,j,ai,j*|1 ^*m,1 勺 n—1} Col ={ai,j,ai*j〉|1i 4m-1,1 j ￡n} 基本操作： CreateSMatrix(M); 操作结果: 创建稀疏矩阵 M。 PrintSMatrix(M) ； 初始条件： 稀疏矩阵 M 存在。 操作结果: 输出稀疏矩阵 M。 AddSMatrix(M ,N , Q); 初始条件： 稀疏矩阵 M 和N的行数和列数对应相等。 操作结果: 求稀疏矩阵的和 Q=M+N 。 MultSMatrix(M ,N, Q); 初始条件： 稀疏矩阵 M 的列数等于N的行数。 操作结果：求稀疏矩阵乘积 Q=M*N 。 TransposeSMatrix(M,T); 初始条件: 稀疏矩阵 M存在。 操作结果: 求稀疏矩阵 M的转置矩阵 }ADT SpareseMatrix 五、主要源程序代码 #includeiostream using namespacestd; #define MAXSIZE 100; int num[100]; typedef struct OLNode( int i,j; int e; struct OLNode *right,*down; }OLNode,*OLink; typedef struct( OLink *rhead,*chead; int mu,nu,tu; }CrossList; 〃十字链表结构体定义 int CreatSMatrix_OL(CrossList M){ int i,j,e; OLink q; OLink p; cout请输入稀疏矩阵的行数，列数，非零元素的个数： ; cinM.mu;cinM.nu;cinM.tu; M.rhead=(OLink *)malloc((M.mu+1)* sizeof(OLNode)); M.chead=(OLink *)malloc((M.nu+1)* sizeof(OLNode)); for(i=1;i=M.mu;i++) M.rhead[i]=NULL; for(i=1;i=M.nu;i++) M.chead[i]=NULL; cout请输入稀疏矩阵，若为空，则退出 ”endl; cini;cinj;cine; while (i!=0){ p=(OLink)malloc( sizeof(OLNode)); p-i=i;p-j=j;p-e=e; if (M.rhead[i]==NULL || M.rhead[i]-jj) {p-right=M.rhead[i];M.rhead[i]=p;} else{ q=M.rhead[i]; while (q-right q-right-jj) q=q-right; p-right=q-right; q-right=p;} if (M.chead[j]==NULL || M.chead[j]-ii) {p-down=M.chead[j];M.chead[j]=p;} else{ q=M.chead[j]; while (q-down q-down-ii) q=q-down; p-down=q-down; q-down=p; } cini;cinj;cine; } return 1; } //创建十字链表 void TurnSMatrix_OL(CrossList M)( int a,b; OLink p,q; for(a1;a=M.mu;a++) ( q=p=M.rhead[a]; while (q)( b=p-i; p-i=p-j; p-j=b; q=p-right; p-right=p-down; p-down=q; } } } 〃十字链表实现稀疏矩阵转置 int AddSMatrix_OL(CrossList *A,CrossList *B){ OLNode *pa,*pb,*p,*pre,*cp[100]; int i,j,t; t=A-tu+B-tu; for(j=1;j=A-nu;j++) cp[j