前两天突发奇想,写一个乘法的实现,但不用乘号*。并测试一下性能如何。因此就有了下面的代码:(本文主要目的是为了玩递归和位移,因此仅限自然数)

首先,标准乘法:

 int commonMultiplication(int a, int b) => a * b;

第二,从数学的角度,乘法其实就是加法,只是加法的简写而已,因此 a * b 可以理解为 b 个 a 相加;故得出用加法代替的乘法。为了减少加法的次数,取 a, b 两数的最小值进行循环:

 int plusMultiplication(int a, int b) {
var t = a > b ? b : a, r = 0;
if (t == a) {
a = b;
b = t;
} while (b-- > 0) r += a;
return r;
}

最后,位移乘法。计算机进行位移操作的速度是非常快的,因此,如果能把乘法换算成位移操作的话,速度会不会更快?但是根据计算机二进制的特点,每次向左位移,相当于是乘了2的n次方,怎么办?我们看看例子:

  31 * 18

怎么用位移来实现呢?首先我们取小的数:18,进行拆解。拆解的算法是,任何一个自然数,在计算机中都可以用二进制表示,那么,任何一个自然数,也可以拆解为2的幂的和。比如 18 的二进制表示 为  10010,拆解成十进制,即为 18 = 16 + 2;那么上面可以拆解为:

  31 * 18 = 31 * (16 + 2) = 31 * 16 + 31 * 2

此时即可换成位移的算法:

  31 * 18 = 31 * (16 + 2) = 31 * 16 + 31 * 2 = 31 * 2^4 + 31 * 2^1 = (31 << 4) + (31 << 1)

是不是转成了位移?这里要注意对乘数 18 的处理,位移的关键是找到乘数 b 如何拆解为 2 的幂,这里要递归拆解,以18为例,第一次的关键是找到16这个2的幂,然后用18 - 16,结果等于 2,继续找2等于多少2的幂,直到差小于2为止。同时,考虑到a*b,在实际计算中即使小的乘数也很大的情况,采用类似二分搜索的方式,给定阶数,从大的阶数向小的阶数递归。阶数的设定可以参考计算机的位数,即32位还是64位。因此,计算一个数可以拆分出来的最大2的幂的代码如下:

 int _bitShift(int d, int order) {
if (order == 0) return 0;
var i = 0;
for (var t = d; (t >>= order) > 0; d = t, t = d) i += order;
return i + _bitShift(d, order >> 1);
}

位移乘法实现的代码

int bitMultiplication(int a, int b, int order) {
var t = a > b ? b : a, r = 0;
if (t == a) {
a = b;
b = t;
} for (var i = _bitShift(b, order); i > 0; i = _bitShift(b, order)) {
r += a << i;
b -= 1 << i;
} return b == 1 ? r + a : r;
}

经过采用随机大数进行多次循环测试,结果显示,用加法实现的乘法速度最慢,最快的仍然是 普通乘法,即用乘号的乘法;位移乘法居中。看来底层的计算机指令还是对*乘法做了优化。

最后说明一下,本文纯属原创,没有使用dart 自带math包,尤其是在求一个数可以拆出来的最大的2的幂的数字时。否则直接用math包中的log函数即可得出结果,但这样就起不到自己练习位运算和递归的目的了。

递归算法之不用乘号的乘法——用位移实现乘法(dart语言实现)的更多相关文章

  1. 矩阵乘法 codevs 1287 矩阵乘法

    1287 矩阵乘法  时间限制: 1 s  空间限制: 128000 KB  题目等级 : 黄金 Gold       题目描述 Description 小明最近在为线性代数而头疼,线性代数确实很抽象 ...

  2. 补码一位乘法(Booth算法,C语言实现)

    补码一位乘法 首先了解下什么是补码? 补码概念的理解,需要先从“模”的概念开始. 我们可以把模理解为一个容器的容量.当超出这个 容量时,会自动溢出.如:我们最常见到的时钟,其容量 是 12,过了 12 ...

  3. 大整数算法[11] Karatsuba乘法

    ★ 引子         前面两篇介绍了 Comba 乘法,最后提到当输入的规模很大时,所需的计算时间会急剧增长,因为 Comba 乘法的时间复杂度仍然是 O(n^2).想要打破乘法中 O(n^2) ...

  4. 大整数算法[09] Comba乘法(原理)

    ★ 引子          原本打算一篇文章讲完,后来发现篇幅会很大,所以拆成两部分,先讲原理,再讲实现.实现的话相对复杂,要用到内联汇编,要考虑不同平台等等. 在大整数计算中,乘法是非常重要的,因为 ...

  5. 【老鸟学算法】大整数乘法——算法思想及java实现

    算法课有这么一节,专门介绍分治法的,上机实验课就是要代码实现大整数乘法.想当年比较混,没做出来,颇感遗憾,今天就把这债还了吧! 大整数乘法,就是乘法的两个乘数比较大,最后结果超过了整型甚至长整型的最大 ...

  6. Algorithm --> 矩阵链乘法

    动态规划--矩阵链乘法 1.矩阵乘法       Note:只有当矩阵A的列数与矩阵B的行数相等时A×B才有意义.一个m×r的矩阵A左乘一个r×n的矩阵B,会得到一个m×n的矩阵C. #include ...

  7. [转]大整数算法[11] Karatsuba乘法

    ★ 引子         前面两篇介绍了 Comba 乘法,最后提到当输入的规模很大时,所需的计算时间会急剧增长,因为 Comba 乘法的时间复杂度仍然是 O(n^2).想要打破乘法中 O(n^2) ...

  8. python_九九乘法表

    # 九九乘法表 print(" 九九乘法表") for table_x in range(1,10): for table_y in range(1,table_x +1): pr ...

  9. For循环练习之99乘法表和转义字符

    之前说了for循环的概念以及常用到的操作,那么我们接下来做几个巩固练习: 1.打印99乘法表: 99乘法表的形式: 1*1 = 1 1*2 = 2 2*2 = 4 1*3 = 3 2*3 = 6 3* ...

随机推荐

  1. SQL语句 运算符

    6.2 运算符   6.2.1 算术运算符 加 / 减 / 乘 / 除 6.2.2 连接运算符 是用来连接字符串的.跟java中的 + 是一致的. select 'abc' || ' bcd ' as ...

  2. scikit-plot

    安装说明 安装Scikit-plot非常简单,直接用命令: pip install scikit-plot 即可完成安装. 仓库地址: https://github.com/reiinakano/sc ...

  3. C#任务调度——LimitedConcurrencyLevelTaskScheduler

    这是参考大佬分享的代码写的有问题请提出指正,谢谢. using Serilog; using System; using System.Collections.Generic; using Syste ...

  4. priority_queue与multiset

    感觉c++最有用的stl是bitset.堆(优先队列)和平衡树,其他的都可以手打 这里主要讲一下堆和平衡树的基本用法和区别所在 priority_queue 堆/优先队列 定义: priority_q ...

  5. guava中Multimap、Multiset使用

    guava中的Multimap接口 Multimap和java.util.Map接口没有任何继承关系.同Map一样,也是放键值对,但是Multimap的值是一个集合.同样支持泛型,假如键值对的key的 ...

  6. linux 硬件时间 系统时间

    建了一个虚拟机,发现每次重启后系统时钟总是跟现有时间相差10多个小时 用date -s进行修正之后,再次重启又出现该问题.于是怀疑跟硬件时钟有关 用hwclock看了下,发现的确硬件时钟的时间存在差异 ...

  7. BZOJ 4423: [AMPPZ2013]Bytehattan 平面图转对偶图 + 并查集

    Description 比特哈顿镇有n*n个格点,形成了一个网格图.一开始整张图是完整的.有k次操作,每次会删掉图中的一条边(u,v),你需要回答在删除这条边之后u和v是否仍然连通. Input 第一 ...

  8. [Noip模拟题]宠物之战senso

    Description 众所周知,moreD的宠物已经被moreD奴役得体无完肤.这只宠物实在忍无可忍,把自己每天走魔法树的经历告诉了 自己的宠物.同时他还说明了自己爬树是多么地慢,以至于moreD每 ...

  9. Android中StatFs获取系统/sdcard存储(剩余空间)大小

    1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 3 ...

  10. (49)LINUX应用编程和网络编程之四 Linux进程全解

    补充: 1.  C程序的执行过程: C编译器调用链接器,链接器设置可执行程序文件的启动起始地址(启动例程),启动例程获得内核传递来的 命令行参数和环境变量值,为调用main函数做准备.[实际上该启动例 ...