数据结构实验报告-特殊矩阵和稀疏矩阵.doc

精品文档 实验五 特殊矩阵和稀疏矩阵【实验目的】1、掌握数组的结构类型（静态的内存空间配置）；通过数组的引用下标转换成该数据在内存中的地址；2、掌握对称矩阵的压缩存储表示；3、掌握稀疏矩阵的压缩存储-三元组表表示，以及稀疏矩阵的转置算法。【实验学时】2学时【实验预习】回答以下问题：1、什么是对称矩阵？写出对称矩阵压缩存储sak与aij之间的对应关系。2、什么是稀疏矩阵？稀疏矩阵的三元组表表示。【实验内容和要求】1、编写程序exp5_1.c，将对称矩阵进行压缩存储。（1）对称矩阵数组元素Aij转换成为以行为主的一维数组sak，请描述k与ij的关系。（注意C程序中，i，j，k均从0开始）（2）调试程序与运行。对称矩阵存储下三角部分即i=j。对称矩阵为 3,9,1,4,7 9,5,2,5,8 1,2,5,2,4 4,5,2,1,7 7,8,4,7,9参考程序如下：#include#define N 5int main() int upperNN= 3,9,1,4,7, 9,5,2,5,8, 1,2,5,2,4, 4,5,2,1,7, 7,8,4,7,9 ; /*对称矩阵*/ int rowMajor15; /*存储转换数据后以行为主的数组*/ int Index; /*数组的索引值*/ int i,j; printf(Two dimensional upper triangular array:n); for (i=0; iN; i+) /*输出对称矩阵*/ for(j=0; jN; j+) printf(%3d,upperij); printf(n); for(i=0; iN; i+) /*进行压缩存储*/ for(j=0; j=j_) /*下三角元素进行存储*/ Index=_i*(i+1)/2+j_; /*ij与index的转换*/ rowMajorIndex=upperij; printf(nRow Major one dimensional array:n); for(i=0; i15; i+) /*输出转换后的一维数组*/ printf(%3d, rowMajori); printf(n); return 1;2、完成程序exp5_2.c，实现稀疏矩阵的三元组表存储及稀疏矩阵的转置。调试并给出结果： 补充完整程序，运行稀疏矩阵的一般转置算法； 完成稀疏矩阵的快速转置算法，并修改主函数的转置调用算法，验证快速转置算法的正确性。exp5_2.c部分代码如下：#include#define MAXSIZE 20 /*非零元素个数最大值*/typedef int ElemType;typedef struct int i,j; ElemType e;Triple;typedef struct Triple dataMAXSIZE+1; /*三元组表，data0不用*/ int mu,nu,tu;/*矩阵的行数、列数、非零元个数*/TSMatrix;void TransposeSMatrix(TSMatrix *T,TSMatrix *M);/*一般转置算法*/void FastTransposeSMatrix(TSMatrix *M,TSMatrix *T);/*快速转置算法*/int main() /int i,j,k,q,col,p; int i,j,k; int temp67=0,12,9,0,0,0,0, /*稀疏矩阵*/ 0,0,0,0,0,0,0, -3,0,0,0,0,14,0, 0,0,24,0,0,0,0, 0,18,0,0,0,0,0, 15,0,0,-7,0,0,0, ; TSMatrix T,M; M.mu=6; M.nu=7; M.tu=0; k=1; for (i=0;i M.mu;i+) /*转换为稀疏矩阵的三元组表示*/ for (j=0;j M.nu;j+) if (tempij!=0) M.datak.i=i+1; M.datak.j=j+1; M.datak.e=tempij; k+; M.tu=k-1; FastTransposeSMatrix(&M,&T); /*调用转置算法进行转置*/ /*输出转置结果*/ printf(稀疏矩阵：n); for (i=0;i M.mu;i+) /*转换为稀疏矩阵的三元组表示*/ for (j=0;j M.nu;j+) printf(%3d,tempij); printf(n); printf(转置前M三元组表：nmutnuttun); printf(%dt%dt%dn,M.mu,M.nu,M.tu); printf(nitjten); for (i=1;i=M.tu;i+) printf(%dt%dt%dn,M.datai.i,M.datai.j,M.datai.e); printf(转置后T三元组表：nmutnuttun); printf(%dt%dt%dn,T.mu,T.nu,T.tu); printf(nitjten); for (i=1;imu=M-nu; T-nu=M-mu; T-tu=M-tu; if (T-tu) q=1; for (col=1;colnu;+col) for (p=1;ptu;+p) if (M-datap.j=col) T-dataq.i=M-datap.j; T-dataq.j=M-datap.i; T-dataq.e=M-datap.e; +q; /*稀疏矩阵的快速转置算法*/void FastTransposeSMatrix(TSMatrix *M,TSMatrix *T) int t,q,col,p,numMAXSIZE,cpotMAXSIZE; T-mu=M-nu; T-nu=M-mu; T-tu=M-tu; if (T-tu) /*快速转置过程的实现，请补充代码*/ for (col=1;colnu;+col) numcol=0; for(t=1;ttu;+t) +numM-datat.j; cpot1=1; for(col=2;colnu;+col) cpotcol=cpotcol-1+numcol-1; for(p=1;ptu;+p) col=M-datap.j; q=cpotcol;