方程式x^y=y^x

本稿では方程式x^y=y^xの代数的・幾何的な性質とその意味について数学的に解説する。

一般に、実数において冪乗は可換とはならないが、方程式x^y = y^xには(x,y)=(2,4)などの無数の実数解が存在する。

## 歴史
方程式x^y=y^xについて、少なくとも1728年6月29日にダニエル・ベルヌーイがクリスティアン・ゴールドバッハに宛てた手紙に綴られている。手紙の内容には、x\ne y の時、その自然数解は(2, 4)と(4, 2)のみであるが、その有理数解は(\tfrac{27}{8}, \tfrac{9}{4})や(\tfrac{64}{27},\tfrac{256}{81})など無数に存在するという記述がされている。

また、1729年1月31日のゴールドバッハの返信では、y=vxと代入することで得られる一般解について述べられていおり、同様の解がレオンハルト・オイラーによっても得られている。(後述)

1888年、オランダ人数学者のヨーハン・ファン・ヘンゲル は、二つの正整数r, nについて、r \geq 3のときr^{r+n} > (r+n)^rが常に成り立つことを示し、x^y=y^x(x\ne y)のすべての自然数解を求めるには（式の対称性により）x = 1,2の場合を考えれば十分であることを証明した。これにより、この方程式の自然数解は(2, 4)と(4, 2)のみであることが確定した。

その後もこの問題は多くの出版物で議論を呼んだ。1961年にニューヨーク市立大学シティ・カレッジのアルヴィン・ハウスナーは結果を代数体まで拡張することに成功した。

また、1960年にはウィリアム・ローウェル・プットナム数学競技大会の問題の一つとして出題されている。

## 正の実数解

### 解の陽函数表示
正の実数の自明な解の集合は\{(x,y) | x = y\}によって与えられる。非自明な解の集合はランベルトのW函数によって表される。具体的には、方程式をae^b = cの形に変形し、aとb を一致させ、ランベルトのW函数に関する公式a'e^{a'} = c' \Rightarrow a' = W(c')を適応することで得られる。

: \begin{align}
y^x &= x^y = \exp\left(y\ln x\right) & \\
y^x \exp\left(-y\ln x\right) &= 1 & \left(\mbox{multiply by } \exp\left(-y\ln x\right)\right) \\
y\exp\left(-y\frac{\ln x}{x}\right) &= 1 & \left(\mbox{raise by } 1/x\right) \\
-y\frac{\ln x}{x}\exp\left(-y\frac{\ln x}{x}\right) &= \frac{-\ln x}{x} & \left(\mbox{multiply by } \frac{-\ln x}{x}\right)
\end{align}

: \Rightarrow -y\frac{\ln x}{x} = W\left(\frac{-\ln x}{x}\right)
: \Rightarrow y = \frac{-x}{\ln x}\cdot W\left(\frac{-\ln x}{x}\right) = \exp\left(-W\left(\frac{-\ln x}{x}\right)\right)

最後のステップで使用した恒等式W(x)/x = \exp(-W(x))については、ランベルトのW函数の2つの函数に分割し、解を含むそれぞれの区間で恒等式を適用することで得ることができる。

: \begin{align}
W_0\left(\frac{-\ln x}{x}\right) &= -\ln x \quad&\text{for } &0 

* 0 :

: \Rightarrow \frac{-\ln x}{x} \ge 0
: \begin{align}\Rightarrow y &= \exp\left(-W_0\left(\frac{-\ln x}{x}\right)\right) \\
 &= \exp\left(-(-\ln x)\right) \\
 &= x \end{align}

* 1 :

: \Rightarrow \frac{-1}{e} 
: \Rightarrow y = \begin{cases}
\exp\left(-W_0\left(\frac{-\ln x}{x}\right)\right) = x \\
\exp\left(-W_{-1}\left(\frac{-\ln x}{x}\right)\right)
\end{cases}

* x = e:

: \Rightarrow \frac{-\ln x}{x} = \frac{-1}{e}
: \Rightarrow y = \begin{cases}
\exp\left(-W_0\left(\frac{-\ln x}{x}\right)\right) = x \\
\exp\left(-W_{-1}\left(\frac{-\ln x}{x}\right)\right) = x
\end{cases}

* x > e:

: \Rightarrow \frac{-1}{e} 
: \Rightarrow y = \begin{cases}
\exp\left(-W_0\left(\frac{-\ln x}{x}\right)\right) \\
\exp\left(-W_{-1}\left(\frac{-\ln x}{x}\right)\right) = x
\end{cases}

以上より、

が得られる。

### パラメータ解（陰函数表示）
非自明な解についてx \ne yと仮定しy = vx(v\ne 0)と置くことによって

: (vx)^x = x^{vx} = (x^v)^x

が得られ、両辺を\tfrac{1}{x}乗しxで割ることで、

: v = x^{v-1}

が得られる。これを使うことで、下記の正の実数における非自明な解のパラメータ解が得られる。

自明な解と併せて、正の実数での一般解は次のようになる。

(y=x) \cup \left(v^{1/(v-1)},v^{v/(v-1)}\right) \text{ for } v > 0, v \neq 1 .

この解に基づくと、自明な解 y=xのyに関する導関数dy/dxは1となり、非自明な解のyに関する導関数は(dy/dv)/(dx/dv)で求められ、次のように計算される。

: \frac{dy}{dx} = v^2\left(\frac{v-1-\ln v}{v-1-v\ln v}\right) \text{ for } v > 0, v \neq 1 .

また、v=2またはv=\tfrac{1}{2}を上記のパラメータ解に代入することで非自明な唯一の自然数解4^2=2^4が得られる。代数的数が解となる例は他にも(\sqrt 3, 3\sqrt 3)、( \sqrt[3]4, 4\sqrt[3]4)などがある。

上記のパラメータ解からは、x^y=y^xのグラフの幾何学的性質を導くことができる。例えば、先ほどの導関数の結果からは、v\neq 1を満たす正の実数v
に対してy=x^vは傾きがv^2となる点を持ち、それらはすべてx^y=y^xを通ることがわかる。

自明な解と非自明な解はv = 1のときに一致する。 非自明な解はv = 1で値を持たないが、v\to 1の極限を取ることで値を得られる。計算は、v = 1 + 1/nを代入してn\to\inftyに飛ばすことで簡単にできる。

: x = \lim_{v\to 1}v^{1/(v-1)} = \lim_{n\to\infty}\left(1+\frac 1n\right)^n = e.

よって、y = xのグラフとx^y = y^x(y \ne x)のグラフの交差する点は(x,y)=(e,e)であることがわかる。なお、ここでeはネイピア数である。

また、非自明な曲線は、x \to \inftyのときy = 1に漸近するため、漸近展開は次のようになる。

: y = 1 + \frac{\ln x}{x} + \frac{3}{2} \frac{(\ln x)^2}{x^2} + \cdots.

## その他の実数解
上記のパラメータ解のvに特定の実数を入れることで、xとyの一方が負であるような実数解も無数に得られる。たとえば、v=-2を代入することによってx=-\frac{1}{\sqrt[3]{2}}, y=\frac{2}{\sqrt[3]{2}}が得られ、xは負の値をとっている。同様に、自明な解y=xについてもx^xが実数となるならば、離散的ではあるがx=y=-1のような解が無数に得られる。

## 似たグラフ

### 
方程式\sqrt[x]y = \sqrt[y]xのグラフは、直線y=xと曲線が(x,y)=(\frac{1}{e},\frac{1}{e})で交差する。曲線は(0,1)と(1,0)を端点としている。

曲線部分のグラフの式は次のように陽函数表示される。

y=e^{W_0(x\ln(x))} \quad \mathrm{for} \quad 0

y=e^{W_{-1}(x\ln(x))} \quad \mathrm{for} \quad 1/e
### 
方程式\log_x(y) = \log_y(x)のグラフは、直線y=xと曲線が(x,y)=(1,1)で交差する。曲線は双曲線y=\frac{1}{x}の正部分と一致し、x \to \inftyのときy = 0に漸近する。

## 参考文献