数据结构 系数矩阵的表示和转置.docVIP

实验四 稀疏矩阵的表示和转置 实验人： 学号: 时间： 分析与探讨 附录：源代码 源代码列在附录中，要求程序风格清晰易理解，有充分的注释。有意义的注释行不少于30%。 #includestdio.h #includestdlib.h #define P 4 #define N 5 #define MAXSIZE 100 typedef int ElementType; typedef struct { int row,col; ElementType e; }Triple; typedef struct { Triple data[MAXSIZE+1]; int m,n,len; }TSMatrix; void InitMatrix(int A[P+1][N+1])//初始化矩阵 { int i,j; printf(请输入一个4行5列的矩阵：\n); for(i=1;i=P;i++) { for(j=1;j=N;j++) { scanf(%d,A[i][j]); } } for(i=1;iP+1;i++) { for(j=1;jN+1;j++) { printf(%d\t,A[i][j]); } printf(\n); } printf(\n); } TSMatrix MatrixToTriple(int A[P+1][N+1],TSMatrix M)//将一般的矩阵转换成三元组形式的 { int r,c,k=0; for(r=1;r=P;r++) { for(c=1;c=N;c++) { if(A[r][c]!=0) { k++; M.data[k].row=r; M.data[k].col=c; M.data[k].e=A[r][c];//一个小小的错误，将A[r][c]写成A[c][r]，导致调了5个多小时 } //调程序一定要认真啊(⊙o⊙)！ } } M.m=P; M.n=N; M.len=k; return M; } TSMatrix FastTransposeTSMatrix(TSMatrix A,TSMatrix B) {//稀疏矩阵的快速转置算法 int lie,t,p,q; int num[MAXSIZE],position[MAXSIZE]; B.len=A.len; B.m=A.n; B.n=A.m; if(B.len) { for(lie=1;lie=A.n;lie++) { num[lie]=0; } for(t=1;t=A.len;t++) { num[A.data[t].col]++; } position[1]=1; for(lie=2;lie=A.n;lie++) { position[lie]=position[lie-1]+num[lie-1]; } for(p=1;p=A.len;p++) { lie=A.data[p].col; q=position[lie]; B.data[q].row=A.data[p].col; B.data[q].col=A.data[p].row; B.data[q].e=A.data[p].e; position[lie]++; } } return B; } void printMatrix(TSMatrix M) { int B[N+1][P+1]={0}; int i,j,k=1; for(i=1;iN+1;i++) { for(j=1;jP+1;j++) { if(i==M.data[k].row j==M.data[k].col) { B[i][j]=M.data[k].e; k++; } } } printf(稀疏矩阵转换成一般矩阵后，再进行打印：\n); for(i=1;iN+1;i++) { for(j=1;jP+1;j++) { printf(%4d\t,B[i][j]); } printf(\n); } } void main() { int A[P+1][N+1]={0}; TSMatrix M={{0},0,0,0},T={{0},0,0,0}; InitMatrix(A); M=MatrixToTriple(A,M); T=FastTransposeTSMatrix(M,T);