Question

这是我的代码，我正在尝试使用math.log来查找int的日志。我尝试了解决不准确日志的大部分方法，但似乎没有任何代码可以使用我的代码：

num = 810333333333333333253432343224234
print(num)
getcontext().prec = 100 # set the precision  (no. of digits)
getcontext().rounding = ROUND_DOWN  # this will effectively truncate
logA = Decimal(math.log(num,2))

print(logA)
#find num using component of logs
seedfrac = Decimal(logA - int(logA))
nonfrac = 2**int(logA)
frac = Decimal(2**seedfrac)
print(Decimal(frac*nonfrac))
#find num directly 
print(Decimal(2**logA ))

输出以检查值：

810333333333333333253432343224234
109.320214523928456173962331376969814300537109375
810333333333340688230991571062453.5576735355138989553226940926052307835750320132600744995177525695124
810333333333340688230991571062453.5576735355138989553226940926052307835750320132600744995177525695128

如果您知道任何解决方法，请回复，谢谢。

Answer 1

math.log仅适用于浮点数，因此当您执行Decimal(math.log(num,2))时，您只需要将num转换为浮点数然后采用其浮点精度调用Decimal日志中。

首先将您的数字转换为Decimal实例，保留精度，然后使用其对数方法之一（下面插入换行符以提高可读性）：

In [22]: Decimal(num)
Out[22]: Decimal('810333333333333333253432343224234')

In [23]: Decimal(num).ln()
Out[23]: Decimal('75.77499847546938418941086142055648421904630259496362044
                  157156139031463040418047508186428539214239394')

In [24]: Decimal(num).ln().exp()
Out[24]: Decimal('810333333333333333253432343224234.0000000000000000000000
                  000000000000000000000000000000000000000000020')

In [25]: Decimal(num).ln() / Decimal(2).ln()
Out[25]: Decimal('109.32021452392844307936383214097698765413874994582696830
                  23528366628242511675596034347551786781907708')

In [26]: 2**(Decimal(num).ln() / Decimal(2).ln())
Out[26]: Decimal('810333333333333333253432343224233.999999999999999999999
                  9999999999999999999999999999999999999999999548')

python math.log在较大的int处输出不准确？

1 个答案: