实验报告稀疏矩阵程序.docVIP

实验说明 1) 首先将稀疏矩阵表示成三元组形式---用数组A表示； 2） 接着利用稀疏矩阵的快速转置算法进行转换，保存在数组B中； 3） 最后输出经过转置后的稀疏矩阵形式。 实验详细说明 快速转置算法的基本思想：按照A的三元组a-data的顺序进行转置，并且将转置后的三元组一次置入b-data中恰当的位置上。这就需要预先确定矩阵A中各列第一个非零元素在b-data中的正确位置。注意到A中第0列的第一个非零元素一定存储在B中下标为0的位置上，该列中其它非零元素应存放在B中后面连续的位置上，那么第1列的第一个非零元素在B中的位置便等于第0列的第一个非零元素在B中的位置加上第0列的非零元素的个数，以此类推。 为了实现上述算法，需要设置两个一维数组num和cpot。数组num[col]用于统计矩阵A中第col列的非零元素个数，数组cpot[col]用于表示矩阵A第col列中第一个非零元素在b-data中的恰当的存储位置。cpot[col]值可按下列公式求得： #define SMAX 16 typedef int datatype; typedef struct three { int i,j; datatype v; }node; typedef struct threelist { int m,n,t; node data[SMAX+1]; }spmat; spmat *fast_transmat(spmat *a) { int ano,bno; int col; spmat *b; int num[SMAX]; int cpot[SMAX]; b=(spmat *)malloc(sizeof(spmat)); b-m=a-n; b-n=a-m; b-t=a-t; if (b-t0) { for (col=1;col=a-n;col++) num[col]=0; for (ano=1;ano=a-t;ano++) { col=a-data[ano].j; num[col]=num[col]+1; } cpot[1]=1; for (col=2;col=a-n;col++) cpot[col]=cpot[col-1]+num[col-1]; for (ano=1;ano=a-t;ano++) { col=a-data[ano].j; bno=cpot[col]; b-data[bno].i= a-data[ano].j; b-data[bno].j= a-data[ano].i; b-data[bno].v= a-data[ano].v; cpot[col]=cpot[col]+1; } } return b; } 给定稀疏矩阵A如下： 用快速转置算法实现对该矩阵进行转置并输出。 int main() { spmat *A=(spmat *)malloc(sizeof(spmat)); spmat *B=(spmat *)malloc(sizeof(spmat)); coutrow m=; cinA-m;coutendl; coutcol n=; cinA-n;coutendl; coutdatatype t=; cinA-t;coutendl; for(int a=1;a=A-t ;++a) { coutA enter i , j, v:; cinA-data[a].i; cout\t; cinA-data[a].j; cout\t; cinA-data[a].v; coutendl; } coutafter transmat :\n; B=fast_transmat(A); for(int k=1;k=B-t;++k){ couti=B-data [k].i\t j=B-data [k].j\t v=B-data[k] .v\n; } } 1 实 验 总 结（结果和心得体会） 计算机科学与工程学院 5 《数据结构》实验报告 0 -2 0 0 0 13 0 3 0 0 0 5 0 0 9 0 6 0 -7 0