DSB第4部分—新.pptVIP

row col value 0 0 16 0 4 91 1 1 12 2 1 3 2 6 15 3 0 22 3 2 -8 5 0 -16 方法二： 对数组 A (矩阵Mm ? n的三元组） 扫描 n 遍，每遍扫描t次。第 i 遍扫描，找到A中列号为i的元素（即M中第i列中的非0元素），将这些元素转置后依次放入B中，B即为转置矩阵。 row col value 0 0 16 0 3 22 0 5 -16 1 1 12 1 2 3 2 3 -8 4 0 91 6 2 15 0 1 2 3 4 5 6 7 0 1 2 3 4 5 6 7 时间复杂度为O(n*t), t是A的长度。 程序段 k=0; for(i=0;in;i++) for(j=0;jt;j++) if(A[j].col==i) { B[k].col=A[j].row; B[k].row=A[j].col; B[k++].value=A[j].value; } 分析方法二可以发现，它只用了1个指针k，用来指示转置后元素在B中的存放位置，因此在B中必须按照顺序存放，导致要反复扫描A，效率不高。 考虑：若扫描三元组A的任何元素时，都能够知道其在转置后的三元组B中的存放位置，就可以不对A反复扫描了。 方法三: 快速转置算法 快速转置算法通过增加适量的存储空间，达到节省时间的目的。 快速转置算法使用n个指针k[n] (n是矩阵的列数)，指向A中每一列的第一个非零元素在T中的存放位置。 row col value 0 0 16 0 3 22 0 5 -16 1 1 12 1 2 3 2 3 -8 4 0 91 6 2 15 0 1 2 3 4 5 6 7 k[0]=0 k[3]=5 k[5]=7 k[1]=2 k[2]=3 k[3]=6 k[0]=1 k[2]=4 row col value 0 0 16 0 4 91 1 1 12 2 1 3 2 6 15 3 0 22 3 2 -8 5 0 -16 0 1 2 3 4 5 6 7 显然k[0]=0。计算k[i]时，可先考虑计算每一列中 非零元素的个数num[i],得到如下计算公式： k[0]=0 k[i]=k[i-1]+num[i-1] (1≤i≤n) col 0 1 2 3 4 5 num[col] 2 1 2 2 0 1 k[col] 0 2 3 5 7 7 表4-1 num数组和k数组 快速转置算法基本步骤： (1)计算num (2)计算k (3)快速转置 row col value 0 0 16 0 3 22 0 5 -16 1 1 12 1 2 3 2 3 -8 4 0 91 6 2 15 0 1 2 3 4 5 6 7 num[0]++ - 1 num[3]++ - 1 num[5]++ - 1 num[1]++ - 1 num[2]++ - 1 num[3]++ - 2 num[0]++ - 2 num[2]++ - 2 col 0 1 2 3 4 5 num[col] 2 1 2 2 0 1 k[col] 0 2 3 5 7 7 表4-1 num数组和k数组 num[i]初始值为0。 计算num[i]时也只需 对三元组扫描一遍。 快速转置算法 程序4.6 稀疏矩阵的快速转置 template class T void SeqTripleT::Transpose(SeqTripleT B)const { int *num=new int[n]; int *k=new int[n]; B.m=n;B.n=m;B.t=t; if(t0) { for (int i=0; in; i++) num[i]=0;//初始化num for (i=0; it; i++) num[trip[i].col]++;//计算num k[0]=0; f