稀疏矩阵的表示和转置.doc

实验2 稀疏矩阵的表示和转置 实验人：杜国胜 学号:X 时间：11.9 实验目的 掌握稀疏矩阵的三元组顺序表存储结构 掌握稀疏矩阵的转置算法。 实验内容 采用三元组表存储表示，求稀疏矩阵M的转置矩阵T。 （算法5.1） 实验步骤： 构建稀疏矩阵M。 求稀疏矩阵M的转置矩阵T。 输出稀疏矩阵M和稀疏矩阵T。 算法说明 首先要创建稀疏矩阵和三元组顺序表，定义mu,mu,tu分别表示矩阵的行列数和非零元个数。在进行稀疏矩阵的转置时要做到1.将矩阵的行列值相互交换2.将每个三元组的I,j相互调换3.重排三元组之间的次序 测试结果 分析讨论 在此次程序中转置的方法称为快速转置，在转置前，应先求的M的每一列中非零元的个数，进而求得每一列的第一个非零元的位置 附录：源代码 #includestdio.h #define MAXSIZE 100 typedef struct { int i,j; int e; }Triple; typedef struct { Triple data[MAXSIZE+1]; int mu,nu,tu; }TSMatrix; //创建稀疏矩阵M void CreateSMatrix (TSMatrix *M) { int i,m,n,e,k; printf(输入矩阵M的行数、列数、非零元的个数（中间用逗号隔开）：); scanf(%d,%d,%d,(*M).mu,(*M).nu,(*M).tu); (*M).data[0].i=0; printf(\n); for(i=1;i=(*M).tu;i++) { do { printf(输入第%d个非零元素所在的行(1~%d)列(1~%d)值以及该数值：,i,(*M).mu,(*M).nu); scanf(%d,%d,%d,m,n,e); k=0; if(m1||m(*M).mu||n1||n(*M).nu) k=1; if(m(*M).data[i-1].i||m==(*M).data[i-1].in(*M).data[i-1].j) k=1; }while(k); (*M).data[i].i=m; (*M).data[i].j=n; (*M).data[i].e=e; } printf(\n); } //输出稀疏矩阵M void PrintSMatrix(TSMatrix M) { int i; printf(**************************************\n); for(i=1;i=M.tu;i++) printf(%2d%4d%8d\n,M.data[i].i,M.data[i].j,M.data[i].e); printf(**************************************\n); printf(\n); } //求稀疏矩阵M的转置矩阵T void TransposeSMatrix(TSMatrix M,TSMatrix *T) { int p,q,col; (*T).mu=M.nu; (*T).nu=M.mu; (*T).tu=M.tu; if((*T).tu) { q=1; for(col=1;col=M.nu;++col) for(p=1;p=M.tu;++p) if(M.data[p].j==col) { (*T).data[q].i=M.data[p].j; (*T).data[q].j=M.data[p].i; (*T).data[q].e=M.data[p].e; ++q; } } } void print(TSMatrix A) { int k=1,a,b; int M[MAXSIZE][MAXSIZE]; printf(非零元素所对应的位置：\n); printf(**************************************\n); for(a=0;aA.mu;a++) { for(b=0;bA.nu;b++) M[a][b]=0; } while(k=A.tu) { M[A.data[k].i-1][A.data[k].j-1]=A.data[k].e; k++; } for(a=0;aA.mu;a++) { printf( | ); for(b=0;bA.nu;b++) printf(%d ,M[a][b]); printf( | \n); } printf(**************************************\n); printf(\n); } //主函数 int main() { TSMatrix M,T; printf(创建矩阵M：