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

计算机科学与工程学院实验说明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。数组numcol用于统计矩阵A中第col列的非零元素个数，数组cpotcol用于表示矩阵A第col列中第一个非零元素在b-data中的恰当的存储位置。cpotcol值可按下列公式求得：cpot1=1；cpotcol=cpotcol-1+numcol-1； 2coln#define SMAX 16typedef int datatype;typedef struct threeint i,j; datatype v;node;typedef struct threelist int m,n,t; node dataSMAX+1;spmat;spmat *fast_transmat(spmat *a) int ano,bno; int col; spmat *b; int numSMAX; int cpotSMAX; b=(spmat *)malloc(sizeof(spmat); b-m=a-n; b-n=a-m; b-t=a-t;if (b-t0) for (col=1;coln;col+) numcol=0; for (ano=1;anot;ano+) col=a-dataano.j; numcol=numcol+1; cpot1=1; for (col=2;coln;col+) cpotcol=cpotcol-1+numcol-1; for (ano=1;anot;ano+) col=a-dataano.j; bno=cpotcol; b-databno.i= a-dataano.j; b-databno.j= a-dataano.i;b-databno.v= a-dataano.v;cpotcol=cpotcol+1; return b;给定稀疏矩阵A如下：0 -2 0 0 0 13 0 3 0 00 5 0 0 90 6 0 -7 00 4 0 0 0用快速转置算法实现对该矩阵进行转置并输出。int main() spmat *A=(spmat *)malloc(sizeof(spmat); spmat *B=(spmat *)malloc(sizeof(spmat); coutA-m;coutendl; coutA-n;coutendl; coutA-t;coutendl; for(int a=1;at ;+a) coutA-dataa.i; coutA-dataa.j; coutA-dataa.v; coutendl; coutafter transmat :n;B=fast_transmat(A);for(int k=1