- iOS/Objective-C 元类和类别
- objective-c - -1001 错误,当 NSURLSession 通过 httpproxy 和/etc/hosts
- java - 使用网络类获取 url 地址
- ios - 推送通知中不播放声音
今天,我来求助于一道小数学题。我提出了一个程序,该程序计算具有任意长度数字的方程式。 https://bitbucket.org/Dermenslof/bistro在优化方面,我正在为两个操作数超过 10,000 位的运算集成乘法算法。这个算法是 toom-cook 第 3 部分。以下是原则:我们要进行如下操作:123 456 789 * 987 654 321我们将把这两个数字分成三部分:
------------------- ------------------
| 123 | 456 | 789 | | x2 | x1 | x0 |
------------------- = ------------------
| 987 | 654 | 321 | | y2 | y1 | y0 |
------------------- ------------------
我们被告知:
the result P(X) = A(X) * B(X)
A(X) = x2 * X² + x1 * X + x0
B(X) = y2 * X² + y1 * X + y0
P(X) = p4 * X⁴ + p3 * X³ + p2 * X² + p1 * X¹ + p0
由此我们可以得到:
X = -2 A = 369 B = 2961 A * B = AB(2) = 1092609 = (4 * x2 + 2 * x1 + x0) * (4 * y2 + 2 * y1 + y0)
X = -1 A = 456 B = 654 A * B = AB(-1) = 298224 = (x2 - x1 + x0) * (y2 - y1 + y0)
X = 0 A = 789 B = 321 A * B = AB(0) = 253269 = x0 * y0 (gives us directly p0)
X = 1 A = 1368 B = 1962 A * B = AB(1) = 2684016 = (x2 + 2 * x1 + x0) * (y2 + 2 * y2 + y0)
X = 2 A = 2193 B = 5577 A * B = AB(2) = 12230361 = (4 * x2 + 2 * x1 + x0) * (4 * y2 + 2 * y1 + y0)
X = inf A = 123 B = 987 B * B = AB(inf) = 121401 = x2 * y2 (gives us directly p4)
我们推断:
AB(-2) = 16 * p4 - 8 * p3 + 4 * p2 + 2 * p1 + p0
AB(-1) = p4 - p3 + p2 - p1 + p0
AB(0) = p0
AB(1) = p4 + p3 + p2 + p1 + p0
AB(2) = 16 * p4 + 8 * p3 + 4 * p2 + 2 * p1 + p0
我们现在必须隔离 p3、p2 和 p1
p3 = (2 * AB(-1) + AB(2) - 18p4 - 6p2 - 3p0) / 6
p2 = (AB(1-) + AB(1) - 2 * p4 - 2 * p0) / 2
p1 = AB(-1) - p4 - p3 - p2 - p0
我的结果是:
121 401 000 000 000 000
+ 000 530 514 000 000 000
+ 000 001 116 450 000 000
+ 000 000 000 662 382 000
+ 000 000 000 000 253 269
-------------------------
121 932 631 112 635 269
一切正常。现在让我烦恼的是:
p3 = (2 * AB(-1) + AB(2) - 18 * p4 - 6 * p2 - 3 * p0) / 6
由于技术原因,我们发现除法的计算速度非常慢。除以 2 或除以基数的幂不是问题。但是将一个可能非常大的数字除以 6 是一个巨大的数字。所以你必须找到 X 的不同值,或者这个多项式的另一种形式以获得更合适的方程。
我得到了一些信息,显然以下方程式更合适:( https://stuff.mit.edu/afs/athena/astaff/source/src-9.0/third/gmp/doc/multiplication )
AB(0) = p0
16 * AB(1/2) = p4 + 2 * p3 + 4 * p2 + 8 * p1 + 16 * p0
AB(1) = p4 + p3 + p2 + p1 + p0
AB(2) = 16 * p4 + 8 * p3 + 4 * p2 + 2 * p1 + p0
AB(inf) = p4
这就是我请求你帮助的地方,我无法从这些等式中分离出 p3、p2 和 p1预先感谢那些对这个问题感兴趣的人
最佳答案
我无法从这些[方程式中分离出] p3、p2 和 p1
我开始编辑你的问题以修复奇数符号错误(不是唯一的一个;参见AB(-2)),继续对齐“ps”(发现结果与您链接的旧 GMP 文档中的表格非常相似):
AB(0) = p0
16AB(1/2) = p4 + 2p3 + 4p2 + 8p1 + 16p0
AB(1) = p4 + p3 + p2 + p1 + p0
AB(2) = 16p4 + 8p3 + 4p2 + 2p1 + p0
AB(inf) = p4
“显然,我们可以”减去第一行和最后一行的适当倍数得到
h = 2p3 + 4p2 + 8p1
u = p3 + p2 + p1
d = 8p3 + 4p2 + 2p1
类似地,10u -h-d 给出
6p2
,2h-12u+ d
6p1
,并且对称地 h-12u+2d
6p3
.
这仍然需要至少除以 6(或三,如果你认为/2 “免费”) - 据我所知,它是 state of the art截至 17 年底。
关于algorithm - toom-cook 第 3 部分分辨率多项式,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/27470561/
我在一本书(Interview Question)中读到这个问题,想在这里详细讨论这个问题。请点亮它。 问题如下:- 隐私和匿名化 马萨诸塞州集团保险委员会早在 1990 年代中期就有一个绝妙的主意
我最近接受了一次面试,面试官给了我一些伪代码并提出了相关问题。不幸的是,由于准备不足,我无法回答他的问题。由于时间关系,我无法向他请教该问题的解决方案。如果有人可以指导我并帮助我理解问题,以便我可以改
这是我的代码 public int getDist(Node root, int value) { if (root == null && value !=0) return
就效率而言,Strassen 算法应该停止递归并应用乘法的最佳交叉点是多少? 我知道这与具体的实现和硬件密切相关,但对于一般情况应该有某种指南或某人的一些实验结果。 在网上搜索了一下,问了一些他们认为
我想学习一些关于分布式算法的知识,所以我正在寻找任何书籍推荐。我对理论书籍更感兴趣,因为实现只是个人喜好问题(我可能会使用 erlang(或 c#))。但另一方面,我不想对算法进行原始的数学分析。只是
我想知道你们中有多少人实现了计算机科学的“ classical algorithms ”,例如 Dijkstra's algorithm或现实世界中的数据结构(例如二叉搜索树),而不是学术项目? 当有
我正在解决旧编程竞赛中的一些示例问题。在这个问题中,我们得到了我们有多少调酒师以及他们知道哪些食谱的信息。制作每杯鸡尾酒需要 1 分钟,我们需要使用所有调酒师计算是否可以在 5 分钟内完成订单。 解决
关闭。这个问题是opinion-based .它目前不接受答案。 想要改进这个问题? 更新问题,以便 editing this post 可以用事实和引用来回答它. 关闭 8 年前。 Improve
我开始学习 Nodejs,但我被困在中间的某个地方。我从 npm 安装了一个新库,它是 express -jwt ,它在运行后显示某种错误。附上代码和错误日志,请帮助我! const jwt = re
我有一个证书,其中签名算法显示“sha256rsa”,但指纹算法显示“sha1”。我的证书 SHA1/SHA2 的标识是什么? 谢谢! 最佳答案 TL;TR:签名和指纹是完全不同的东西。对于证书的强度
我目前在我的大学学习数据结构类(class),并且在之前的类(class)中做过一些算法分析,但这是我在之前的类(class)中遇到的最困难的部分。我们现在将在我的数据结构类(class)中学习算法分
有一个由 N 个 1x1 方格组成的区域,并且该区域的所有部分都是相连的(没有任何方格无法到达的方格)。 下面是一些面积的例子。 我想在这个区域中选择一些方块,并且两个相邻的方块不能一起选择(对角接触
我有一些多边形形状的点列表,我想将其包含在我页面上的 Google map 中。 我已经从原始数据中删除了尽可能多的不必要的多边形,现在我剩下大约 12 个,但它们非常详细以至于导致了问题。现在我的文
我目前正在实现 Marching Squares用于计算等高线曲线,我对此处提到的位移位的使用有疑问 Compose the 4 bits at the corners of the cell to
我正在尝试针对给定算法的约束满足问题实现此递归回溯函数: function BACKTRACKING-SEARCH(csp) returns solution/failure return R
是否有包含反函数的库? 作为项目的一部分,我目前正在研究测向算法。我正在使用巴特利特相关性。在 Bartlett 相关性中,我需要将已经是 3 次矩阵乘法(包括 Hermitian 转置)的分子除以作
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 这个问题似乎与 help center 中定义的范围内的编程无关。 . 关闭 8 年前。 Improve
问题的链接是UVA - 1394 : And There Was One . 朴素的算法是扫描整个数组并在每次迭代中标记第 k 个元素并在最后停止:这需要 O(n^2) 时间。 我搜索了一种替代算法并
COM 中创建 GUID 的函数 (CoCreateGUID) 使用“分散唯一性算法”,但我的问题是,它是什么? 谁能解释一下? 最佳答案 一种生成 ID 的方法,该 ID 具有一定的唯一性保证,而不
在做一个项目时我遇到了这个问题,我将在这个问题的实际领域之外重新措辞(我想我可以谈论烟花的口径和形状,但这会使理解更加复杂).我正在寻找一种(可能是近似的)算法来解决它。 我有 n 个不同大小的容器,
我是一名优秀的程序员,十分优秀!