python - 如何组合指数？ (x**a)**b => x**(a*b)?

Question

如何在sympy中简化方程中的指数

from sympy import symbols 
a,b,c,d,e,f=symbols('abcdef')
j=(a**b**5)**(b**10)
print j 
(a**(b**5))**(b**10) #ans even after using expand simplify 
# desired output 
a**(b**15)

如果不能使用 sympy，我应该在 python 中导入哪个模块？

编辑 即使我将“b”定义为实数，所有其他符号也是如此

b=symbols('b',real=True) 没有得到简化的指数 只有当指数是常量时它才会化简

a=symbols('a',real=True)
b=symbols('b',real=True)
(a**5)**10
 a**50  #simplifies only if exp are numbers
(a**b**5)**b**10


(a**(b**5))**b**10  #no simplification

Answer 1

(x^m)ⁿ = x^mn 为真 only if m, n are real .

>>> import math
>>> x = math.e
>>> m = 2j*math.pi
>>> (x**m)**m      # (e^(2πi))^(2πi) = 1^(2πi) = 1
(1.0000000000000016+0j)
>>> x**(m*m)       # e^(2πi×2πi) = e^(-4π²) ≠ 1
(7.157165835186074e-18-0j)

据我所知，sympy supports complex numbers ，所以我认为不应该进行这种简化，除非你能证明 b 是真实的。

---

编辑:如果 x 不是正数，它也为假。

>>> x = -2
>>> m = 2
>>> n = 0.5
>>> (x**m)**n
2.0
>>> x**(m*n)
-2.0

编辑(由 gnibbler 编写):这是应用了 Kenny 限制的原始示例

>>> from sympy import symbols 
>>> a,b=symbols('ab', real=True, positive=True)
>>> j=(a**b**5)**(b**10)
>>> print j
a**(b**15)