大整数乘法中的分治思想(TOOM-COOK的一种使用方法)

2023-11-24 13:30

本文主要是介绍大整数乘法中的分治思想(TOOM-COOK的一种使用方法),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

算法分析与设计学习中,接触到一道大整数乘法问题,分享出来,原题目如下:

算法分析在用分治法求两个n位大整数u和v的乘积时,将u和v都分割为长度为n/3的3段。证明可以用5次n/3位整数的乘法求得uv的值。按此思想设计大整数乘积的分治方法,并分析算法的计算复杂性。

先参考一道较为简单的题目:设有两个n位二进制数X,Y,求它们的乘积XY。
分析:按照一般算法,根据小学数学乘法规律,两个数中每位数都需要相应做乘法,则需要的时间复杂度是O(n^2)。

此时,我们考虑将X,Y分为高位、低位,即
X = | A| B | , Y = | C| D|
A,B,C,D均为n/2位,求XY的问题可以转换为:

	X * Y = (A * 2^(n/2) +B )  *(C * 2^(n/2) +D )= AC * 2^n + (AD + BC) * 2^(n/2)+BD			 

则,可以看出利用分治法后,进行了4次n/2位的乘法运算。但是并没有涉及算法性能优化,时间复杂度依然是O(n^2)。

这个算法进行性能提升可以使用如下方法:
先计算

U = (A + B)(C + D), V = AC, W = BD则 Z = XY = V * 2^n +(U - V- W ) * 2^(n/2) + W

上面过程中,由于只使用了U、V、W涉及的3次乘法,比没有优化的少了一种,时间复杂度就降低到了在这里插入图片描述

具体 过程可以参考 Karatsuba算法

Karatsuba算法伪代码实现如下

procedure karatsuba(num1, num2)if (num1 < 10) or (num2 < 10)return num1*num2/* calculates the size of the numbers */m = max(size_base10(num1), size_base10(num2))m2 = m/2/* split the digit sequences about the middle */high1, low1 = split_at(num1, m2)high2, low2 = split_at(num2, m2)/* 3 calls made to numbers approximately half the size */z0 = karatsuba(low1,low2)z1 = karatsuba((low1+high1),(low2+high2))z2 = karatsuba(high1,high2)return (z2*10^(2*m2))+((z1-z2-z0)*10^(m2))+(z0)

具体例子可以参照这个图片
Karatsuba算法

那参照上述步骤,在解决最开始提到的5次n/3位整数乘法中也可以设法将两个大整数U、V分割为

U = |A| B| C
V =|D| E| F

然后利用(A +B +C )(D +E +F)的方法,在分别细分合并其中几项,来简化复杂度,以求达到5次乘法的要求。

PS,这个方法是一开始我的想法,实现后发现只能简化到6次乘法实现,所以最后寻找别的算法。

最后发现有TOOM-COOK方法来做这个工作,而题目涉及的分为3等份的过程就是TOOM-COOK中当n = 3 的特例,也称为TOOM3。

那先来看看TOOM-COOK的一般实现,即不规定分为多少份,设分为m份。

设有两个大整数U,V,利用分治思维将U、V分为如下部分

U  = |U-(m-1)|……|U2 |U1 |U0
V   =|V-(m-1)|……|V2 |V1 |V0
设X = 10^(n/m)

则可以将u和v及其乘积w=uv表示为
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

将U,V和W都看作关于变量X的多项式,可以得到
在这里插入图片描述

此时我们可以取2m-1个不同的数x1,x2,…,x2m-1代入上多项式,可得W(xi) 与 W0,W1,W2……W(2m-2)的关系,转换为矩阵为

在这里插入图片描述
那设B为红框部分的矩阵可以推得
在这里插入图片描述

此时,取x1、 x2、 x3、 x4…. 为不同的数,再结合以下两个式子,
在这里插入图片描述
可以得出W0 、W1、 W2、 W3、 W4与U0、 U1、U2 ,V0、 V1、 V2 的关系
最后移位相加得到最终结果,移位相加可以参考如下的图,因为设的是10进制的数,所以每次移位就是10^(n/m)
在这里插入图片描述

那其实说到这里,相信有人已经明白了TOOM-COOK算法的整体核心思想,整体实现流程可以简单概括为
在这里插入图片描述

那回归到本题,在具体的3等分n位大整数U、V问题中,可以详细描述为以下细节
首先,分治划分U、V

U = |U2 |U1 |U0
V =|V2 |V1 |V0

同样的,设X = 10^(n/3),那U、V及其乘积W = UV可以表示为
在这里插入图片描述
分别取X为5个不同的数,即

X1 = 0,X2 = 1,X3 = -1, X4 = 2,X5 =-2

代入多项式中
在这里插入图片描述
第一个结论,是W(Xi)与U0,U1,U2,U3,V0,V1,V2的关系,暂且将5种X取值下的W(Xi)标记为a,b,c,d,e
在这里插入图片描述
第二个结论,是W(Xi)即a,b,c,d,e与W = UV中各项参数W0,W1,W2,W3,W4的,其中可以看到矩阵B(TOOM-COOK一般方法中提到)
在这里插入图片描述
至此,建立了W0,W1,W2,W3,W4与U1、U2 ,V0、 V1、 V2 的关系,经过运算后,可以求出两个式子

在这里插入图片描述
在这里插入图片描述
再带入W=UV=W0+W1X+W2X2 +W3X3+W4X4 中,可以求出W
到了这一步,相信大家都已经对这道题有了更深的理解,下面我们来看一道具体的应用。

123456*987654

按照程序实现流程,划分UV后,分别求得a,b,c,d,e以及W0,W1,W2,W3,W4如下
在这里插入图片描述
在这里插入图片描述

UV=W0+W1X+W2X2+W3X3+W4X4即W0, W1 , W2, W3, W4移位相加可得W结果
在这里插入图片描述
对比运算,可知结果无误
在这里插入图片描述

最后再说一下时间复杂度,a,b,c,d,e所涉及的一共5次乘法,加减不计,所以得到的递归方程如下:
在这里插入图片描述
求解得复杂度为在这里插入图片描述

TOOM-COOK实现思路和算法可以参考大数乘法问题

/*** Java8中的 Toom-Cook multiplication 3路乘法*/
private static BigInteger multiplyToomCook3(BigInteger a, BigInteger b) {int alen = a.mag.length;int blen = b.mag.length;int largest = Math.max(alen, blen);// k is the size (in ints) of the lower-order slices.int k = (largest+2)/3;   // Equal to ceil(largest/3)// r is the size (in ints) of the highest-order slice.int r = largest - 2*k;// Obtain slices of the numbers. a2 and b2 are the most significant// bits of the numbers a and b, and a0 and b0 the least significant.BigInteger a0, a1, a2, b0, b1, b2;a2 = a.getToomSlice(k, r, 0, largest);a1 = a.getToomSlice(k, r, 1, largest);a0 = a.getToomSlice(k, r, 2, largest);b2 = b.getToomSlice(k, r, 0, largest);b1 = b.getToomSlice(k, r, 1, largest);b0 = b.getToomSlice(k, r, 2, largest);BigInteger v0, v1, v2, vm1, vinf, t1, t2, tm1, da1, db1;v0 = a0.multiply(b0);da1 = a2.add(a0);db1 = b2.add(b0);vm1 = da1.subtract(a1).multiply(db1.subtract(b1));da1 = da1.add(a1);db1 = db1.add(b1);v1 = da1.multiply(db1);v2 = da1.add(a2).shiftLeft(1).subtract(a0).multiply(db1.add(b2).shiftLeft(1).subtract(b0));vinf = a2.multiply(b2);// The algorithm requires two divisions by 2 and one by 3.// All divisions are known to be exact, that is, they do not produce// remainders, and all results are positive.  The divisions by 2 are// implemented as right shifts which are relatively efficient, leaving// only an exact division by 3, which is done by a specialized// linear-time algorithm.t2 = v2.subtract(vm1).exactDivideBy3();tm1 = v1.subtract(vm1).shiftRight(1);t1 = v1.subtract(v0);t2 = t2.subtract(t1).shiftRight(1);t1 = t1.subtract(tm1).subtract(vinf);t2 = t2.subtract(vinf.shiftLeft(1));tm1 = tm1.subtract(t2);// Number of bits to shift left.int ss = k*32;BigInteger result = vinf.shiftLeft(ss).add(t2).shiftLeft(ss).add(t1).shiftLeft(ss).add(tm1).shiftLeft(ss).add(v0);if (a.signum != b.signum) {return result.negate();} else {return result;}
}

参考资料
大数乘法问题
Karatsuba
两个n位大整数相乘算法
算法分析设计习题2-5
TOOM-COOK算法
TOOM-COOK
TOOM-3

这篇关于大整数乘法中的分治思想(TOOM-COOK的一种使用方法)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/422021

相关文章

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python实现无痛修改第三方库源码的方法详解

《Python实现无痛修改第三方库源码的方法详解》很多时候,我们下载的第三方库是不会有需求不满足的情况,但也有极少的情况,第三方库没有兼顾到需求,本文将介绍几个修改源码的操作,大家可以根据需求进行选择... 目录需求不符合模拟示例 1. 修改源文件2. 继承修改3. 猴子补丁4. 追踪局部变量需求不符合很

java中使用POI生成Excel并导出过程

《java中使用POI生成Excel并导出过程》:本文主要介绍java中使用POI生成Excel并导出过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录需求说明及实现方式需求完成通用代码版本1版本2结果展示type参数为atype参数为b总结注:本文章中代码均为

mysql出现ERROR 2003 (HY000): Can‘t connect to MySQL server on ‘localhost‘ (10061)的解决方法

《mysql出现ERROR2003(HY000):Can‘tconnecttoMySQLserveron‘localhost‘(10061)的解决方法》本文主要介绍了mysql出现... 目录前言:第一步:第二步:第三步:总结:前言:当你想通过命令窗口想打开mysql时候发现提http://www.cpp

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

Spring Boot3虚拟线程的使用步骤详解

《SpringBoot3虚拟线程的使用步骤详解》虚拟线程是Java19中引入的一个新特性,旨在通过简化线程管理来提升应用程序的并发性能,:本文主要介绍SpringBoot3虚拟线程的使用步骤,... 目录问题根源分析解决方案验证验证实验实验1:未启用keep-alive实验2:启用keep-alive扩展建

MySQL INSERT语句实现当记录不存在时插入的几种方法

《MySQLINSERT语句实现当记录不存在时插入的几种方法》MySQL的INSERT语句是用于向数据库表中插入新记录的关键命令,下面:本文主要介绍MySQLINSERT语句实现当记录不存在时... 目录使用 INSERT IGNORE使用 ON DUPLICATE KEY UPDATE使用 REPLACE

使用Java实现通用树形结构构建工具类

《使用Java实现通用树形结构构建工具类》这篇文章主要为大家详细介绍了如何使用Java实现通用树形结构构建工具类,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录完整代码一、设计思想与核心功能二、核心实现原理1. 数据结构准备阶段2. 循环依赖检测算法3. 树形结构构建4. 搜索子

GORM中Model和Table的区别及使用

《GORM中Model和Table的区别及使用》Model和Table是两种与数据库表交互的核心方法,但它们的用途和行为存在著差异,本文主要介绍了GORM中Model和Table的区别及使用,具有一... 目录1. Model 的作用与特点1.1 核心用途1.2 行为特点1.3 示例China编程代码2. Tab