数组的类型定义数组的顺序表示和实现矩阵的压.pptVIP

分析算法FastTransposeSMatrix的时间复杂度：时间复杂度为:O(M.nu+M.tu)for(col=1;col=M.nu;++col)……for(t=1;t=M.tu;++t)……for(col=2;col=M.nu;++col)……for(p=1;p=M.tu;++p)……5.3.3矩阵的压缩存储三元组顺序表又称有序的双下标法，它的特点是，非零元在表中按行序有序存储，因此便于进行依行顺序处理的矩阵运算。然而，若需随机存取某一行中的非零元，则需从头开始进行查找。行逻辑联接的顺序表5.3.3矩阵的压缩存储#defineMAXRC500typedefstruct{Tripledata[MAXSIZE+1];intrpos[MAXRC+1];intmu,nu,tu;}RLSMatrix;//行逻辑链接顺序表类型将上节快速转置矩阵的算法中创建的指示“行”信息的辅助数cpot固定在稀疏矩阵的存储结构中。5.3.3矩阵的压缩存储两个矩阵相乘的经典算法Q=M×NM是m1×n1,N是n1×n2,Q是m1×n2矩阵乘法的精典算法:for(i=1;i=m1;++i)for(j=1;j=n2;++j){Q[i][j]=0;for(k=1;k=n1;++k)Q[i][j]+=M[i][k]*N[k][j];}其时间复杂度为？O(m1×n2×n1)5.3.3矩阵的压缩存储稀疏矩阵相乘×＝5.3.3矩阵的压缩存储ije11314522-1312ije12221131-2324ije12621-1324row1234rpos[row]1245row1234rpos[row]13455.3.3矩阵的压缩存储StatusMultSMatrix(RLSMatrixM,RLSMatrixN,RLSMatrixQ){intarow,brow,p,q,t,ctemp[30],l,ccol,tp;if(M.nu!=N.mu)returnERROR;Q.mu=M.mu;Q.nu=N.nu;Q.tu=0;//Q初始化if(M.tu*N.tu!=0){//Q是非零矩阵5.3.3矩阵的压缩存储{for(arow=1;arow=M.mu;++arow){for(l=1;l=M.nu;++l)ctemp[l]=0;Q.rpos[arow]=Q.tu+1;if(arowM.mu)tp=M.rpos[arow+1];elsetp=M.tu+1;for(p=M.rpos[arow];ptp;++p){//对当前行中每一个非零元brow=M.data[p].j;//找到对应元在N中的行号if(browN.mu)t=N.rpos[brow+1];elset=N.tu+1;5.3.3矩阵的压缩存储for(q=N.rpos[brow];qt;++q){ccol=N.data[q].j;//乘积元素在Q中列号ctemp[ccol]+=M.data[p].e*N.data[q].e;}//forq}//求得Q中第crow(=arow)行的非零元5.3.3矩阵的压缩存储for(ccol=1;c