如何找到一系列n * 0.1(在Python中)最接近反对数(10base)的每个整数

时间:2013-06-02 16:42:57

标签: python algorithm logarithm

在我的算法中,我想评估当前整数n是否恰好最接近0.1的任何正数的反对数,所以0.1,1.0,1.1,7.9,21.5等。

反对数? AntiLog(x)= 10 x

我找到了一个在线计算器,可以在这里找到antilog:http://ncalculators.com/number-conversion/anti-log-logarithm-calculator.htm但我无处可找到任何用Python或任何其他编程语言执行此操作的示例。

如果我不能使用Python来查找一系列数字的反对数,我将不得不求助于在我的程序中存储一个反对数值列表,为了性能,我可能会考虑到,但不管怎样它会很棒找出如何在代码中执行此操作。

更新 使用第一个答案的代码我能够做到这一点,这段代码完全展示了我想要做的所有事情:

#!/usr/bin/python
import sys
import math

for x in xrange(1, 1000000):
    target = round(math.log10(x),1)
    int_antilog = int(10**target+0.5) #fast round() alternative
    if (x == int_antilog):
        print 'do something at call ' + str(x)

2 个答案:

答案 0 :(得分:5)

Python有power operator执行你想要的取幂:

def antilog(x):
    return 10 ** x

示例:

>>> antilog(0.1)
1.2589254117941673
>>> antilog(3)
1000
>>> inputs = [0.1, 1.0, 1.1, 7.9, 21.5]
>>> outputs = [antilog(x) for x in inputs]
>>> print outputs
[1.2589254117941673, 10.0, 12.589254117941675, 79432823.47242822, 3.1622776601683794e+21]

以你喜欢的方式回合并转换为整数,你会很高兴。

答案 1 :(得分:1)

如果k是10的最接近整数,则是0.1的倍数,那么k ≈ 10ᵐʳ,其中m是整数r=0.1,所以ln k ≈ m r ln 10和{{1} }。如果m ≈ ln k / (r ln 10) = log₁₀(k)/r,那么10的幂与最接近k的0.1的倍数相比,是你的最大值u = log₁₀(k)/r倍或u的上限的10倍,如以下代码。

r

这是一个结果示例,通过ipython解释器:

import math

def aloghilo(k):
    u = 10*math.log10(k)
    return (10**(math.floor(u)/10), 10**(math.ceil(u)/10))

def isantilog(k):
   (lo,hi) = aloghilo(k)
   return k==round(lo) or k==round(hi)