稀疏矩阵的操作(数据结构).doc

稀疏矩阵的操作(数据结构)稀疏矩阵的操作要求：1.稀疏矩阵采用三元组表示。2.两个相同行列数的稀疏矩阵A和B的相加矩阵C，并求出C.3.求A矩阵的转置矩阵D,输出D. 我恰好也做此题啊!你运气真好!下面的是C语言写的,不知是否可行?那矩阵的名字可能不符合你的要求,你自己改一下吧,呵呵 #include #define MAXSIZE 100 /* 非零元个数的最大值 */ typedef struct triple int i,j; /* 行下标,列下标 */ int e; /* 非零元素值 */ triple; typedef struct tsmatrix triple dataMAXSIZE+1; /* 非零元三元组表,data0未用 */ int mu,nu,tu; /* 矩阵的行数、列数和非零元个数 */* 各列第一个非零元的位置表rpos0未用 */ rlsmatrix;createsmatrix(rlsmatrix *M) /* 创建稀疏矩阵M */ int e,i,m,n; M-data0.i=0; /* 为以下比较顺序做准备 */ printf(请输入矩阵的行数，列数，和非零元素的个数:); scanf(%d,&M-mu);scanf(%d,&M-nu);scanf(%d,&M-tu); for(i=1;itu;i+) printf(请按行序顺序输入第%d个非零元素所在的行(1%d),列(1%d),元素值：,i,M-mu,M-nu); scanf(%d,&m);scanf(%d,&n);scanf(%d,&e); if(mM-mu|nM-nu) /*行或列超出范围 */ printf(行或列超出范围);getch();exit(); if(mdatai-1.i|m=M-datai-1.i&ndatai-1.j) /*行或列的顺序有错*/ printf(行或列的顺序有错);getch();exit(); M-datai.i=m; M-datai.j=n; M-datai.e=e; /* 求矩阵的快速转置 */void transposesmatrix(rlsmatrix M,rlsmatrix *T) /* cpos存放每列的第一个非零元素的地址，temp中间变量 */ int i,m,*cpos,*temp,k=0; T-mu=M.nu; T-nu=M.mu; T-tu=M.tu; cpos=(int *)malloc(M.mu*sizeof(int); if(cpos=NULL)exit(); temp=(int *)malloc(M.mu*sizeof(int); if(temp=NULL)exit(); /* 对cpos对初始化，初值为0 */ *(cpos+1)=0; for(i=1;i=M.nu;i+) for(m=1;m1) *(cpos+i)=*(temp+i-1)+1; free(temp); for(i=1;idata*(cpos+M.datai.j).i=M.datai.j; T-data*(cpos+M.datai.j).j=M.datai.i; T-data*(cpos+M.datai.j).e=M.datai.e; (*(cpos+M.datai.j)+; free(cpos); multsmatrix(rlsmatrix M,rlsmatrix N,rlsmatrix *T) int i,j,Qn=0; int *Qe; if(M.nu!=N.mu) printf(两矩阵无法相乘);getch();exit(); T-mu=M.mu; T-nu=N.nu; Qe=(int *)malloc(M.mu*N.nu*sizeof(int); /* Qe为矩阵Q的临时数组*/ for(i=1;i=M.mu*N.nu;i+) *(Qe+i)=0;/* 矩阵Q的第i行j列的元素值存于*(Qe+(M.datai.i-1)*N.nu+N.dataj.j)中，初值为0 */ for(i=1;i=M.tu;i+) /* 结果累加到Qe */ for(j=1;j=N.tu;j+) if(M.datai.j=N.dataj.i) *(Qe+(M.datai.i-1)*N.nu+N.dataj.j)+=M.datai.e*N.dataj.e; for(i=1;i=M.mu;i+) /*Qe矩阵中，因为M的每一行和N的每一列相乘都是T的一个元素，不管它是零或非零 */ for(j=1;jdataQn.e=*(Qe+(i-1)*N.nu+j); T-dataQn.i=i; T-dataQn.j=j; free(Qe); T-tu=Qn; return ; void printmatrix(rlsmatrix M)/输出 int i,m=1,n,k; printf(矩阵的简化模式是:n); for(i=1;i=M.tu;i+) printf(%d,%d,%dn,M.datai.i,M.datai.j,M.datai.e); printf(矩阵的行数是:%dn,M.mu); printf(矩阵的列数是:%dn,M.nu); printf(矩阵中非零元素个数是:%dn,M.tu); printf(矩阵的完整模式的:n); for(n=1;n=M.mu;n+) printf(|); for(k=1;k=M.nu;k+) if(M.datam.i=n&M.datam.j=k) printf(%3d,M.datam.e); m+; else printf( 0); printf( |n); return; void main() rlsmatrix M,N,T,K; printf(请为矩阵M赋值!n); createsmatrix(&M); printmatrix(M); printf(请为矩阵N赋值!n); createsmatrix(&N); printmatrix(N);getch(); printf(矩阵M*N得矩阵T!n); multsmat