矩阵乘法的一个最佳算法.pdfVIP

矩阵乘法 的一个最佳算法 蒋 昌俊 吴哲挥 (山东矿亚学院 ，豪安) 关键词 算法、矩阵乘法、运算次数、阶 一 、 引 言 矩阵乘法是线性代数中常见的问题之一，许多数值计算 问题都包含着矩阵乘法的计算．因 此，降低矩阵乘法算法的时间复杂度 问题，多年来一直引起算法研究者们的高度重视 ． 1969年，Strassen提出了一个时间复杂度为 0(一 t)的矩阵乘法算法 ，第一次突破了 0()的界限，被誉为 “在代数 复杂性理论中最激动人心的结果 ．以后，又出现了一系列祈 的矩阵乘法算法 ． 据了解，迄今最好的算法是 1981年 由 Pan提出的，其时间复杂度为 O( )”． 一 一 口 占 矩阵乘法的算法最少需要多少时间量呢? Baase指出， 对于两个 阶方阵¨的乘法来说， 阶的乘除法运算是必须的．但他认为，达到这个下界的算法大概是不存在的 ． ， ， ～ 本文给出矩阵乘法的一个算法．该算法适用于有理数矩阵．对于两个 阶方阵相乘，即 (，m，一)问题 ，只需要 阶的运算次数．对于 行m列矩阵和m行 l列矩 阵的乘法 ，即 (一，m， f)问题 ，运算次数的阶为 0(m(一+ f))． 这样 ，本文提 出的算法其运算次数的阶 已达到 了理 论下 界 ． 二、非负整数矩阵乘法的算法 已知两个非负整数矩阵 求矩阵 c一 [fj，】_，， 使得 一 ∑ ^btj．i一1，2，·一，；，一1，2，… ，． 々= 1 对于这个 问题，可以通过下面的算法步骤来实现． 第 1步：计算 口一 max{口j^1i∈{1，2，·一，}，{∈{1，2，·一；m}1， b— max{b^j1{∈{1，2，·-．，m)，“ {1，2，·--，¨， 本文 198g年6月22日收到． 1988年9月19日收到任改稿． ’国家 自然科 学基盎资助项 目 第 4期 科 学 通 报 = ，” 5 卜 1． Y— nmabx一 十 1． 第 2步 ：对 i一 1，2，·一，m，用 Hornet方法计算 口^ 一 ∑ ， … f 6t一∑ b,iY ． 第3步：计算 一∑ 。6I． 第 {步：做除法．傲次求商 d，d 一，d㈦ ，r‘和余数 r】，… ，… 使满 足： f— dIy+ rI， O ≤ rL y dI-I— diY+ ，，， 0≤ ，r Y ，一 2， 3， ·-·， 一 2 以一 一 r,y+ ，f 0≤ {一l y． 令