我想找到一个数字的所有基数和指数。
示例:
Number = 64
2^6=64
4^3=64
8^2=64
64^1=64
Number = 1845.28125
4.5^5=1845.28125
Number = 19683
3^9=19683
27^3=19683
19683^1=19683
我现在要做的是制作一个'数字'的整数。只看到多次计算的结果就会给出正确的结果:
basehits, expohits = [], []
if eval(Number) > 1000:
to = 1000 #base max 1000 in order to avoid too many calculations
else:
to = int(eval(Number))
for n in range(1,to):
for s in range(1,31): #just try with exponents from 1 to 30
calcres = pow(n,s)
if calcres == eval(Number):
basehits.append(n)
expohits.append(s)
elif calcres > eval(Number):
break
问题是,这从未找到浮动数字,例如1845.28125
(见上文)
当只知道结果时,是否有更好的方法来找到指数和基数?
答案 0 :(得分:3)
您的问题需要更多限制,但这里有一些帮助:
>>> from math import log
>>> help(log)
Help on built-in function log in module math:
log(...)
log(x[, base])
Return the logarithm of x to the given base.
If the base not specified, returns the natural logarithm (base e) of x.
>>> for base in range(2, 10):
... exp = log(64, base)
... print('%s ^ %s = %s' % (base, exp, base ** exp))
...
2 ^ 6.0 = 64.0
3 ^ 3.785578521428744 = 63.99999999999994
4 ^ 3.0 = 64.0
5 ^ 2.5840593484403582 = 63.99999999999999
6 ^ 2.3211168434072493 = 63.99999999999998
7 ^ 2.1372431226481328 = 63.999999999999964
8 ^ 2.0 = 64.0
9 ^ 1.892789260714372 = 63.99999999999994
答案 1 :(得分:1)
怎么样
import math
num=64
for i in range(2,int(math.sqrt(num))+1):
if math.log(num,i).is_integer():
print i,int(math.log(num,i))
输出是:
2 6
4 3
8 2
当然,你总是可以添加:
print num,1
获取
64,1
如果要添加分数,在点后面加n个十进制数字,可以使用:
from __future__ import division
import math
num=1845.28125
decimal_digits=1
ans=3
x=1
while(ans>=2):
ans=num**(1/x)
if (ans*10**decimal_digits).is_integer():
print ans,x
x+=1
其中decimal_digits
表示点后面的位数。
对于这个例子,答案是
4.5 5
,
如果您将num
更改为39.0625
而将decimal_digits
更改为2,则输出将为:
2.5 4
6.25 2
答案 2 :(得分:1)
对于整数,您可以查看数字的prime factors。一旦你知道64是2**6
,就可以很容易地列出你想要的所有结果。
现在,您期望哪个数字至少有两个不同的素因子?例如:应该写成3*5
,3**1 * 5**1
还是15**1
?
目前尚不清楚如何为Floats定义问题。
4.5
有什么特别之处?
如果计算1845.28125**(1.0/5)
,Python将返回4.5
,但对于其他输入数字,结果可能会偏离1e-16。
import math
def find_possible_bases(num, min_base = 1.9, max_decimals = 9, max_diff = 1e-15):
max_exponent = int(math.ceil(math.log(num,min_base)))
for exp in range(1,max_exponent):
base = round(num**(1.0/exp),max_decimals)
diff = abs(base**exp-num)
if diff < max_diff:
print('%.10g ** %d = %.10g' % (base, exp, base ** exp))
find_possible_bases(64)
# 64 ** 1 = 64
# 8 ** 2 = 64
# 4 ** 3 = 64
# 2 ** 6 = 64
find_possible_bases(19683)
# 19683 ** 1 = 19683
# 27 ** 3 = 19683
# 3 ** 9 = 19683
find_possible_bases(1845.28125)
# 1845.28 ** 1 = 1845.28
# 4.5 ** 5 = 1845.28
find_possible_bases(15)
# 15 ** 1 = 15
迭代可能的指数,并计算基数。它将其舍入为9位小数,并检查错误变为什么。如果它足够小,它会显示结果。您可以使用参数并找到最适合您的问题。 作为奖励,它也适用于整数(例如64和15)。
使用Rational numbers可能会更好。