有一个很好的技巧可以计算给定浮点数的有理逼近(基于 Euclid 的 GCD 算法的一些属性)。我们可以使用它来确定“最佳”近似值是否为A/(2^a 5^b) 的形式,如果是,则浮点数终止(以 10 为基数),如果不是,它将有一些重复分量。棘手的一点是确定哪个近似值是正确的(由于浮点精度问题)。
以下是您获得近似有理表达式的方法。
首先迭代x = 1/x - floor(1/x) 跟踪int(x)
x = 0.12341234
1/x = 8.102917
x <= 1/x - 8 = 0.102917
1/x = 9.7165
x <= 1/x - 9 = 0.71265277
1/x = 1.3956
x < 1/x - 1 = 0.3956
...
接下来将 x 的 int 部分粘贴到该表的第一行,称它们为 k_i。
A_i = A_{i-2} + k_i * A_{i-1} 的值与 B_i 的值相同。
|| 8 | 9 | 1 | 2 | 1 | 1 | 8 | 1 | 1
A = 1 0 || 1 | 9 | 10 | 29 | 39 | 68 | 583 | 651 | 1234
B = 0 1 || 8 | 73 | 81 | 235 | 316 | 551 | 4724 | 5275 | 9999
然后有理近似是A_n/B_n。
1/8 = 0.12500000000000000 | e = 1.5e-3
9/73 = 0.12328767123287671 | e = 1.2e-4
10/81 = 0.12345679012345678 | e = 4.4e-5
29/235 = 0.12340425531914893 | e = 8.1e-6
39/316 = 0.12341772151898735 | e = 5.4e-6
68/551 = 0.12341197822141561 | e = 3.6e-7
583/4724 = 0.12341236240474174 | e = 2.2e-8
651/5275 = 0.12341232227488151 | e = 1.8e-8
1234/9999 = 0.12341234123412341 | e = 1.2e-9
因此,如果我们在 1234/9999 阶段确定我们的错误足够低,我们注意到 9999 不能写成 2^a 5^b 的形式,因此我们的十进制扩展正在重复。
请注意,虽然这似乎需要很多步骤,但如果我们使用
x = 1/x - round(1/x)(并改为跟踪 round(1/x))。在这种情况下,表格变为
8 10 -4 2 9 -2
1 0 1 10 -39 -68 -651 1234
0 1 8 81 -316 -551 -5275 9999
这会以更少的步骤为您提供先前结果的子集。
有趣的是,分数 A_i/B_i 总是使得 A_i 和 B_i 没有公因数,所以你不必担心抵消因数或类似的事情。
为了比较,让我们看看 x = 0.123 的展开式。我们得到的表格是:
8 8 -3 -5
1 0 1 8 -23 123
0 1 8 65 -187 1000
那么我们的近似序列是
1/8 = 0.125 e = 2.0e-3
8/65 = 0.12307.. e = 7.6e-5
23/187 = 0.12299.. e = 5.3e-6
123/1000 = 0.123 e = 0
我们看到 123/1000 正是我们想要的分数,因为 1000 = 10^3 = 2^3 5^3 我们的分数正在终止。
如果你真的想知道分数的重复部分是什么(什么数字和什么句点),你需要做一些额外的技巧。这涉及分解分母并找到具有所有这些因素(2 和 5 除外)的最小数字 (10^k-1),然后 k 将是您的周期。所以对于我们的顶级案例,我们发现 A = 9999 = 10^4-1 (因此 10^4-1 包含 A 的所有因素 - 我们在这里有点幸运......)所以重复部分的周期是 4 . 你可以找到关于这最后一部分的更多细节here。
此算法的最后一个重要方面是它不需要所有数字都将十进制扩展标记为重复。考虑 x = 0.34482,这有表格:
3 -10 -156
1 0 1 -10 .
0 1 3 -29 .
我们在第二个条目处得到一个非常准确的近似值并停在那里,得出的结论是我们的分数可能是 10/29(因为它在 1e-5 内使用)并且从上面链接中的表格中我们可以看出它的周期将是 28 位数字。这永远无法通过对数字的简短版本进行字符串搜索来确定,这需要至少 57 位数字才能知道。