Introduction

　正整数は素因数分解によって、たとえば $18 = 2 \times 3^2$ のように、素数の正整数乗の積として一意に表すことができます。
　一方、指数に実数を許すと、この表現は一意ではなくなります。実際、$18 = 2^{1+\log_2 3} \times 3 = 2 \times 3 \times 5^{\log_5 3}$ のように、無数の表し方が存在します。
　このように、指数に実数、特に対数（logarithm）を用いた場合に、2つの異なる表現が等しいかどうかを判定するプログラムを、PythonのSymPyを使って簡単に作ってみました。しかし、意外にも正しく評価できないケースがあったため、今回この記事を書くに至りました。

Summary

以下では具体的な実例を扱いますが、全体像をつかんでいただくため、先にサマリーをまとめます。

指数と対数を組み合わせた等式が、正しく評価できない場合がありました
実例では前処理を工夫することで、正しく評価できる例を示しました
ただし、これらは実例に対するadhocな前処理方法であり、すべての等式に適用できるわけではありません
等式評価を汎用的に行える便利なライブラリがあれば、ぜひ知りたいです

Example

以下の$x,y$が等しいことを確かめよ

$x=2^{\log_6 18}$

$y=3^{(2+\log_3 2) \cdot \log_6 2}$

コードはGoogle colaboratoryで作り、以下のGithubに置きました
https://github.com/wagyan/math/blob/master/ExponentsAndLogarithms.ipynb

Codes

sympyを使っています。$\log$関数の引数は（真数、底）の順番であることに注意して下さい。

from sympy import log, simplify, N, Eq
x=2**(log(18,6))
y=3**((2+log(2,3))*log(2,6))

評価１

$x,y$と、等号評価「==」と、sympyの等号評価「equals」、「Eq」(simplify付き)の結果を表示します。

equalsは代数的に等しいかを評価し、Eqは2式が等しければTrueを返しますが基本的には（Falseであれば）等式そのものを返します。

print("x:: {}".format(x))
print("y:: {}".format(y))

print("x==y:: {}".format(x==y))
print("sympy.x.equals(y):: {}".format(x.equals(y)))
print("sympy.Eq(x,y).simplify():: {}".format(Eq(x,y).simplify()))

結果は、Falseです。

x:: 2**(log(3)/log(6) + 1)
y:: 3**((log(2)/log(3) + 2)*log(2)/log(6))
x==y:: False
sympy.x.equals(y):: False
sympy.Eq(x,y).simplify():: Eq(2**(log(3)/log(6) + 1), 3**(log(2)*log(18)/(log(3)*log(6))))

評価２（simplify）

次に$x,y$をsimplifyで整理してみます。

#数式を整理する
x=simplify(x)
y=simplify(y)

ここで$y$の式は少し整理されましたが、やはりFalseとなりました。なお、equalsでNoneと表示されているのは判定不能だそうです。

x:: 2**(log(3)/log(6) + 1)
y:: 3**(log(2)*log(18)/(log(3)*log(6)))
x==y:: False
sympy.x.equals(y):: None
sympy.Eq(x,y).simplify():: Eq(2**(log(3)/log(6) + 1), 3**(log(2)*log(18)/(log(3)*log(6))))

評価３（差分を取ってからsimplify）

$x,y$をそれぞれ整理するよりも、差分を取ってから式を整理すると、上手く0になることもあるので試してみました。

#xとyの差分をとって、0になるか確かめる
x=x-y
y=0

やはりここでもFalse、Noneとなってしまいます。

x:: -3**(log(2)*log(18)/(log(3)*log(6))) + 2**(log(18)/log(6))
y:: 0
x==y:: False
sympy.x.equals(y):: None
sympy.Eq(x,y).simplify():: Eq(-3**(log(2)*log(18)/(log(3)*log(6))) + 2**(log(3)/log(6) + 1), 0)

評価４（log_2 をとってから評価）

式が複雑過ぎて評価出来ない可能性があるので、少なくとも$x$の式を簡単にできるように$\log_2$を取ってみます。

#両辺のlog_2を取ると、途端に正解する
x=log(x,2)
y=log(y,2)

これでやっと、Trueと評価されました。ただし、「==」はFalseと評価されています。

x:: log(2**(log(3)/log(6) + 1))/log(2)
y:: log(3**((log(2)/log(3) + 2)*log(2)/log(6)))/log(2)
x==y:: False
sympy.x.equals(y):: True
sympy.Eq(x,y).simplify():: True

評価５（log_2 をとり、simplifyをしてから評価）

評価４に加えて、さらに両辺をsiplifyで整理してから評価します。

#両辺のlog_2を取ると、途端に正解する
x=log(x,2)
y=log(y,2)

#数式を整理する
x=simplify(x)
y=simplify(y)

ついに、$x,y$の式が一致しました。
これで「==」もTrueを返します。

x:: log(3)/log(6) + 1
y:: log(3)/log(6) + 1
x==y:: True
sympy.x.equals(y):: True
sympy.Eq(x,y).simplify():: True

評価６（数値的一致を見る）

sympyのNを用いて数値計算します。

ここでは、一番最初に定義した$x,y$の式をそのまま計算させています。

#数値計算する
x=N(x)
y=N(y)

計算方法や順序によっては誤差が出るので、本来は完全に0と一致するとは限りませんが、今回は完全一致したようです。

x:: 3.05918465698377
y:: 3.05918465698377
x==y:: True
sympy.x.equals(y):: True
sympy.Eq(x,y).simplify():: True

Conclusion

前述のSummaryのとおりです

それでは良きPythonライフを！