第4章补充矩阵乘法.pptVIP

* 第4章补充 稀疏矩阵乘法 矩阵相乘的数学公式 * 矩阵相乘的经典算法 两个矩阵相乘的经典算法若设C=A*B其中，A是m1*n1矩阵，B是m2*n2矩阵。当n1=m2时有： for (i=1;im1; ++i ) for ( j=1; j=n2; ++j){ C[i][j]=0; for(k=1; k=n1; ++k) C[i][j]+=A[i][k]*B[k][j]; } 此算法的时间复杂度是O(m1*n1*n2)。 矩阵相乘的经典算法的特点 * ci1 ci2 ci3 经典算法就是这样一项项加起来的：计算时要检索矩阵N的列元素。 当M和N是稀疏矩阵并用三元组表作存储结构时，就不能套用经典算法。 如何实现用二元组表作存储结构时的稀疏矩阵的乘法算法？关键是如何利用三元组和二元组检索行很方便的特点？ 经典算法中乘积项aikbkj计算次序 矩阵相乘的新算法 * 注意：对任意k=1,2,3,4，乘积项aikbkj 必是且只能是cij求和公式中的一项， 故可改变计算aikbkj的次序。 先取ai1，分别乘以N中的第1行的每个元素b1j(j=1,2,3)，并分别加到每个cij(j=1,2,3); 当M和N是稀疏矩阵并用三元组表作存储结构时，就不能套用经典算法。 ci1 ci2 ci3 为求得Q的第i行，按此次序计算乘积项aikbkj 矩阵相乘的新算法 * 注意：对任意k=1,2,3,4，乘积项aikbkj 必是且只能是cij求和公式中的一项， 故可改变计算aikbkj 的次序。 再取ai2，分别乘以N中的第2行的每个元素b2j(j=1,2,3)，并分别加到每个cij(j=1,2,3); 当M和N是稀疏矩阵并用三元组表作存储结构时，就不能套用经典算法。 ci1 ci2 ci3 为求得Q的第i行，按此次序计算乘积项aikbkj 矩阵相乘的新算法 * 注意：对任意k=1,2,3,4，乘积项aikbkj 必是且只能是cij求和公式中的一项， 故可改变计算aikbkj的次序。 再取ai3，分别乘以N中的第3行的每个元素b3j(j=1,2,3)，并分别加到每个cij(j=1,2,3); 如此，将M中第i行的每个元素都进行这样的计算，获得Q的第i行。 当M和N是稀疏矩阵并用三元组表作存储结构时，就不能套用经典算法。 ci1 ci2 ci3 为求得Q的第i行，按此次序计算乘积项aikbkj 矩阵相乘的新算法 * 如果矩阵M和N中有0元素，计算时，实际的项不需要计算。 当M和N是稀疏矩阵并用三元组表作存储结构时，就不能套用经典算法。 c11 c12 c13 三元组表表示的矩阵相乘 i j e 1 1 a11 1 2 a12 1 4 a14 2 1 a21 2 3 a23 3 2 a32 3 4 a34 A.data B.data C.data i j e 1 1 b11 1 2 b12 2 2 b22 2 3 b23 3 1 b31 3 3 b33 4 1 b41 i j e 1 1 c11 1 2 c12 1 3 c13 2 1 c21 2 2 c22 … … i j e 1 1 a11 1 2 a12 1 4 a14 2 1 a21 2 3 a23 3 2 a32 3 4 a34 A.data B.data ctemp i j e 1 1 b11 1 2 b12 2 2 b22 2 3 b23 3 1 b31 3 3 b33 4 1 b41 1 0 2 0 3 0 1 0 2 2 3 4 4 6 B.rowStart i j e C.data 计算C=A*B的第一行，首先将ctemp置0 三元组表表示的矩阵相乘 i j e 1 1 a11 1 2 a12 1 4 a14 2 1 a21 2 3 a23 3 2 a32 3 4 a34 ctemp i j e 1 1 b11 1 2 b12 2 2 b22 2 3 b23 3 1 b31 3 3 b33 4 1 a41 c11 0 c12 0 c13 0 1 0 2 2 3 4 4 6 i j e 计算Q=M*N的第一行，在将M第1行元素乘N相应元素并加到ctemp 三元组表表示的矩阵相乘 +a11b11 +a11b12 A.data B.data B.rowStart C.data i j e 1 1 a11 1 2 a12 1 4 a14 2 1 a21 2 3 a23 3 2 a32 3 4 a34 ctemp i j e 1 1 b11 1 2 b12 2 2 b22