数据结构(C++)(第05章.pptVIP

每一列的第一个非零元的位置为 pot[0]=0 第0列第一个非零元位置 pot[1]=pot[0]+pot[1]=2 第1列第一个非零元位置 pot[2]=pot[1]+pot[2]=4 第2列第一个非零元位置 pot[3]=pot[2]+pot[3]=6 第3列第一个非零元位置 pot[4]=pot[3]+pot[4]=7 第4列第一个非零元位置 pot[5]=pot[4]+pot[5]=7 第5列第一个非零元位置 pot[6]=pot[5]+pot[6]=8 第6列第一个非零元位置 则M稀疏矩阵的转置矩阵N的三元组表很容易写出（见图5-8），算法描述如下： void sparmatrix::fastrans(sparmatrix a,sparmatrix b) { int pot[100],col,ano,bno; b.rows=a.cols; b.cols=a.rows; b.terms=a.terms; if (b.terms0) { for(col=0;cola.cols;col++) pot[col]=0; for(int t=0;ta.terms;t++) //求出每一列的非零元个数 { col=a.data[t].j; pot[col+1]=pot[col+1]+1; } pot[0]=0; for(col=1;cola.cols;col++) //求出每一列的第一个非零元在转置后的位置 pot[col]=pot[col-1]+pot[col]; for( ano=0;anoa.terms;ano++) //转置 { col=a.data[ano].j; bno=pot[col]; b.data[bno].j=a.data[ano].i; b.data[bno].i=a.data[ano].j; b.data[bno].v=a.data[ano].v; pot[col]=pot[col]+1; }} for(int i=0;ia.terms;i++) coutb.data[i].i b.data[i].j b.data[i].vendl; //输出转置后的三元组 } void main() { sparmatrix a,b; cina.rowsa.colsa.terms; //输入稀疏矩阵的行、列数及非零元的个数 for( int i=0;ia.terms;i++) cina.data[i].ia.data[i].ja.data[i].v; //输入转置前的三元组 for(i=0;ia.terms;i++) couta.data[i].i a.data[i].j a.data[i].vendl; //输出转置前的三元组 a.fastrans( a,b); //调用快速转置算法 } 该算法比按列转置多用了辅助向量空间pot，但它的时间为四个单循环，故总的时间复杂度为O(a.cols+a.terms)，比按列转置算法效率要高。 2．稀疏矩阵的相加运算 当稀疏矩阵用三元组表进行相加时，有可能出现非零元素的位置变动，这时候，不宜采用三元组表作存储结构，而应该采用十字链表较方便(为描述方便，假设稀疏矩阵的行列下标从1开始而非0)。 （1）十字链表的建立 下面分两步讨论十字链表的建立算法： 第一步，建立表头的循环链表： 依次输入矩阵的行、列数和非零元素个数：m,n和t。由于行、列链表共享一组表头结点，因此，表头结点的个数应该是矩阵中行、列数中较大的一个。假设用s 表示个数，即s=max（m，n）。依次建立总表头结点（由hm指针指向）和s个行、列表头结点，并使用next域使s+1个头结点组成一个循环链表，总表头结点的行、列域分别为稀疏矩阵的行、列数目，s个表头结点的行列域分别为0。并且开始时，每一个行、列链表均是一个空的循环链表，即s个行、列表头结点中的行、列指针域rptr和cptr均指向头结点本身。 第二步，生成表中结点： 依次输入t个非零元素的三元组（i,j,v），生成一个结点，并将它插入到第i行链表和第j列链表中的正确位置上，使第i个行链表和第j个列链表变成一个非空的循环链表。 算法描述如下： class linknode { public: int i,j; union vnext { int v; linknode *next; } k; linknode *rptr,*cptr; linknode *creatlindmat(); }; linknode *linknode::creatlindmat