查找PPT.pptVIP

第五章 数组和广义表 Status MlutSMatrix (RLSMatrix M, RLSMatrix N, RLSMatrix Q) { //求矩阵乘积Q = M×N，采用逻辑链接存储表示。 if (M.nu != N.mu) return ERROR; Q.mu = M.mu; Q.nu = M.nu; Q.tu = 0; //Q初始化 if (M.tu * N.tu !=0) { //Q是非零矩阵 for (arow = 1; arow = M.mu; ++ arow) { //处理M的每一行 ctemp[ ] = 0; //当前各元素累加器清零 Q.rpos[arow] = Q.tu + 1; if (arow M.mu) tp = M.rpos[arow + 1]; 两个稀疏矩阵相乘(Q = M×N)的过程可大致描述如下： Q初始化； if (Q是非零矩阵) { //逐行求积 for (arow = 1; arow = M.mu; ++ arow) { //处理M的每一行 ctemp[ ] = 0; //累加器清零 计算Q中第arow行的积并存入ctemp[ ]中； 将ctemp[ ]中非零元压缩存储到Q.data； } // for arow } // if 算法5.3是上述过程求精的结果。 算法5.3如下： else tp = M.tu + 1; for (p = M.rpos[arow]; p tp; ++p) { //为当前行中每一个非零元找到对应元在N中的行号 brow = M.data[p].j; if (brow N.mu) t = N.rpos[brow + 1]; else t = N.tu + 1; for (q = N.rpos[brow]; q t; ++q) { ccol = N.data[q].j; //乘积元素在Q中的列号 ctemp[ccol] + = M.data[p].e * N.data[q].e; } // for q } //求得Q中第crow( = arow)行的非零元 for (ccol = 1; ccol = Q.nu; ++ccol) //压缩存储该行非零元 if (ctemp[ccol]) { if (++Q.tu MAXSIZE) return ERROR; Q.data[Q.tu] = (arow, ccol, ctemp[ccol]); } // if } // for arow } // if return OK; } // MlutSMatrix 时间复杂度：O(M.mu×N.nu+M.tu×N.tu/N.mu)。 其中累加器ctemp初始化的时间复杂度为O(M.mu×N.nu)， 求Q的所有非零元的时间复杂度为O(M.tu×N.tu/N.mu)， 进行压缩存储的时间复杂度为O(M.mu×N.nu)。 ③十字链表 1．定义 十字链表：在链表中，每个非零元既是某个行链表中的一个结点， 又是某个列链表中的一个结点，这个矩阵构成了一个十字交叉的链表 结构。可用两个分别存储行链表的头指针和列链表的头指针的一维数 组表示之。 2．结点结构 row col val down right 其中，row，col和val分别表示该非零元所在的行、列和非零元的值。 right域，用以链接同一行中下一个非零元。down域，用以链接同一 列中下一个非零元。 3．图形表示 例如，矩阵M＝ 的十字链表如图5.6所示。 图5.6 稀疏矩阵M的十字链表 M.chead 3 1 2 M.rhead 1 1 3 1 4 5 2 2 -1 4．C语言描述 //-----------------------稀疏矩阵的十字链表存储表示----------------------------- typedef struct OLNode { int i, j; //该非零元的行和列下标