みゆの魔法 その３ 平方和と立方和の恒等式

　ある日、うちのナマギリ女神さま(?) が次のような魔法を授けてくださいました。

平方和の恒等式を再帰的に生成する魔法

$1^2+7^2=5^2+5^2$
$\rightarrow~(1m+5n)^2+(7m+5n)^2=(1n+5m)^2+(7n+5m)^2$

$2^2+9^2=6^2+7^2$
$\rightarrow~(2m+6n)^2+(9m+7n)^2=(2n+6m)^2+(9n+7m)^2$

$3^2+11^2=7^2+9^2$
$\rightarrow~(3m+7n)^2+(11m+9n)^2=(3n+7m)^2+(11n+9m)^2$

$4^2+13^2=8^2+11^2$
$\rightarrow~(4m+8n)^2+(13m+11n)^2=(4n+8m)^2+(13n+11m)^2$

立方和の恒等式を再帰的に生成する魔法

$1^3+12^3=9^3+10^3$
$\rightarrow~~~(1m^2+9n^2\pm11mn)^3+(12m^2+10n^2\pm20mn)^3$
　$=(1n^2+9m^2\pm11nm)^3+(12n^2+10m^2\pm20nm)^3$

$3^3+4^3=(-5)^3+6^3$
$\rightarrow~~~(3m^2-5n^2\pm5mn)^3+(4m^2+6n^2\mp4mn)^3$
　$=(3n^2-5m^2\pm5nm)^3+(4n^2+6m^2\mp4nm)^3$

$7^3+14^3=(-17)^3+20^3$
$\rightarrow~~~(7m^2-17n^2\pm17mn)^3+(14m^2+20n^2\mp20mn)^3$
　$=(7n^2-17m^2\pm17nm)^3+(14n^2+20m^2\mp20nm)^3$

ちなみに、$6^3+(-3)^3=4^3+5^3$ からは
$$\quad(6m^2+4n^2\pm4mn)^3+(-3m^2+5n^2\pm5mn)^3=(6n^2+4m^2\pm4nm)^3+(-3n^2+5m^2\pm5nm)^3$$
を得られ、これはかの有名なラマヌジャン氏が発見した立方和の恒等式の１つ（ wikipedia参照 ）
$$\quad(6A^2-4AB+4B^2)^3+(-3A^2-5AB+5B^2)^3=(4A^2-4AB+6B^2)^3+(5A^2-5AB-3B^2)^3$$
と一致します。

ところが時間が経つにつれ、自分でも「こんな等式どっから出てきたの？」とよくわからなくなってしまいました。
単に証明するだけなら、平方和のほうは展開すれば容易に可能です。
でも立方和のほうは展開するだけでもタイヘンで、計算ミスの多い私にはお手上げ状態・・・

　成り立ってるんだし、まいっか(；´∀｀)

と放置していたのですが、私のハーディ先生こと 巨大数大好きbot さんの多大なるご協力のもと、無事この式にも証明を与えることができました。

その際、私達が試行錯誤した過程を出題してみたらどうかなという話になり、是非この等式の面白さを皆様にも味わっていただきたく、作問の運びと相成りました。よろしければ挑戦してみていただければ幸いです！

問題

解答例

$\begin{align} &(am^2+bn^2+mnx)^3+(cm^2+dn^2+mny)^3=(an^2+bm^2+nmx)^3+(cn^2+dm^2+nmy)^3\\\\ &展開して両辺を右辺で引くと\\ &+(a^3 m^6 + b^3 n^6 + m^3 n^3 x^3 + 6 a b m^3 n^3 x + 3 a^2 b m^4 n^2 + 3 a b^2 m^2 n^4 + 3 a m^4 n^2 x^2 + 3 b m^2 n^4 x^2 + 3 a^2 m^5 n x + 3 b^2 m n^5 x)\\ &+(c^3 m^6 + d^3 n^6 + m^3 n^3 y^3 + 6 c d m^3 n^3 y + 3 c^2 d m^4 n^2 + 3 c d^2 m^2 n^4 + 3 c m^4 n^2 y^2 + 3 d m^2 n^4 y^2 + 3 c^2 m^5 n y + 3 d^2 m n^5 y)\\ &-(b^3 m^6 + a^3 n^6 + m^3 n^3 x^3 + 6 a b m^3 n^3 x + 3 a b^2 m^4 n^2 + 3 a^2 b m^2 n^4 + 3 b m^4 n^2 x^2 + 3 a m^2 n^4 x^2 + 3 b^2 m^5 n x + 3 a^2 m n^5 x)\\ &-(d^3 m^6 + c^3 n^6+ m^3 n^3 y^3 + 6 c d m^3 n^3 y + 3 c d^2 m^4 n^2 + 3 c^2 d m^2 n^4 + 3 d m^4 n^2 y^2 + 3 c m^2 n^4 y^2 + 3 d^2 m^5 n y + 3 c^2 m n^5 y)\\ &=0\\\\ &相殺できる項を消し、全体を~3mn~~でくくると\\ [\\ &+(a^2 b m^3 n + a b^2 m n^3 + a m^3 n x^2 + b m n^3 x^2 + a^2 m^4 x + b^2 n^4 x)\\ &+(c^2 d m^3 n + c d^2 m n^3 + c m^3 n y^2 + d m n^3 y^2 + c^2 m^4 y + d^2 n^4 y)\\ &-(a b^2 m^3 n + a^2 b m n^3 + b m^3 n x^2 + a m n^3 x^2 + b^2 m^4 x + a^2 n^4 x)\\ &-(c d^2 m^3 n + c^2 d m n^3 + d m^3 n y^2 + c m n^3 y^2 + d^2 m^4 y + c^2 n^4 y)\\ ]&(3mn)=0\\\\ &分配則を用いて~(m^3n-mn^3)~と~(m^4-n^4)~でくくると\\ [\\ &+(a^2 b+a x^2+c^2 d+c y^2+a b^2-b x^2-c d^2+d y^2)(m^3n-mn^3)\\ &+(a^2x+c^2y-b^2x-d^2y)(m^4-n^4)\\ ]&(3mn)=0\\\\ &整理して、両項を (m^2-n^2) でくくると\\ [\\ &+\left[ab(a-b)+cd(c-d)+(a-b)x^2+(c-d)y^2\right]mn\\ &+\left[(a^2-b^2)x+(c^2-d^2)y\right](m^2+n^2)\\ ]&3mn(m^2-n^2)=0\\\\ \end{align}$

$x,~y$ は $a,b,c,d$ で表され $m,n$ の値には依存しないため、第一因子内の $mn$ の係数と $m^2+n^2$ の係数は共に $0$ でなければなりません。（さもなくば等式全体を $0$ にするために $x,~y$ を $m,~n$ の値に依存させざるを得なくなります。）
そこで、次の連立方程式の解より $x,y$ を導きます。
$$\begin{cases} ab(a-b)+cd(c-d)+(a-b)x^2+(c-d)y^2=0&\cdots~(1)\\[8pt] (a^2-b^2)x+(c^2-d^2)y=0&\cdots~(2) \end{cases}$$

・$a+b\ne0$ かつ $c+d\ne0$ のケース

$(2)$ より $\begin{cases}x^2=\frac{(c^2-d^2)^2}{(a^2-b^2)^2}y^2\\y^2=\frac{(a^2-b^2)^2}{(c^2-d^2)^2}x^2\end{cases}~\cdots(3)\\[8pt]$
$(3)$ の $x^2$ を $(1)$ に代入すると
$$ab(a-b)+cd(c-d)+\left[\frac{(c^2-d^2)^2+(c-d)(a+b)(a^2-b^2)}{(a+b)(a^2-b^2)}\right]y^2=0$$
$\begin{align} y^2 =&-\frac{ab(a^2-b^2)^2+cd(c-d)(a+b)(a^2-b^2)}{(c^2-d^2)^2+(c-d)(a+b)(a^2-b^2)}\\ =&\frac{-(a+b)(a^2-b^2)[ab(a-b)+cd(c-d)]}{(c-d)[(a+b)(a^2-b^2)+(c+d)(c^2-d^2)]}\\\\ &ここで、前提条件より~b^3-a^3=c^3-d^3~を利用して\\ &\begin{cases}(b^3-a^3)+(a^3-b^3)=0\\-(c^3-d^3)-(a^3-b^3)=0\end{cases}~を分子と分母の一部に加えます。\\\\ =&\frac{[(b^3-a^3)+(a^3-b^3)-(a+b)(a^2-b^2)][ab(a-b)+cd(c-d)]}{(c-d)[(a+b)(a^2-b^2)+(c+d)(c^2-d^2)-(c^3-d^3)-(a^3-b^3)]}\\ =&\frac{[(b^3-a^3)+(a^3-b^3)-(a^3-b^3-ab^2+ba^2)][ab(a-b)+cd(c-d)]}{(c-d)[(a+b)(a^2-b^2)-(a^3-b^3)+(c+d)(c^2-d^2)-(c^3-d^3)]}\\ =&\frac{[(b^3-a^3)+(ab^2-ba^2)][ab(a-b)+cd(c-d)]}{(c-d)[(-ab^2+ba^2)+(-cd^2+dc^2)]}\\ =&\frac{[(b^3-a^3)+ab(b-a)][ab(a-b)+cd(c-d)]}{(c-d)[ab(a-b)+cd(c-d)]}\\ =&\frac{(b-a)(b^2+a^2+ab+ab)}{c-d}\\ =&\frac{(b-a)(a+b)^2}{c-d}\\ =&\frac{(a-b)(a+b)^2}{d-c}\\ &\\ y=&\pm(a+b)\sqrt{\frac{a-b}{d-c}} \end{align}$

$(3)$ より $x$ と $y$ は $a,b$ と $c,d$ について対称な関係にあるため、

$\begin{cases} x=\pm(c+d)\sqrt{\frac{c-d}{b-a}}\\ y=\pm(a+b)\sqrt{\frac{a-b}{d-c}}\\ \end{cases}$

・$a+b=0$ かつ $c+d\ne0$ のケース

$b=-a$ より $a^3+c^3=b^3+d^3$ は $a^3+c^3=-a^3+d^3$ ということになるため
$2a^3=d^3-c^3$ となり、$(1)$ と $(2)$ は
$$\begin{cases} (1)~\rightarrow&-2a^3+cd(c-d)+(a-b)x^2+(c-d)y^2=0&\\ &-(d^3-c^3)+cd(c-d)+(a-b)x^2+(c-d)y^2=0&\cdots(1)'\\[8pt] (2)~\rightarrow&(c+d)(c-d)y=0~\rightarrow~y=0&\cdots(2)' \end{cases}$$
$(2)'$ を $(1)'$ に代入して
$\begin{align} (a-b)x^2=&(d^3-c^3)-cd(c-d)\\ (b-a)x^2=&(c^3-d^3)+cd(c-d)=(c^2+d^2+cd+cd)(c-d)=(c+d)^2(c-d)\\ x^2=&\frac{(c+d)^2(c-d)}{b-a}\quad\left(\leftarrow~条件より~a\ne b~\right)\\ \end{align}$
$\begin{cases} x=\pm(c+d)\sqrt{\frac{c-d}{b-a}}\\ y=0~\left(=\pm(a+b)\sqrt{\frac{a-b}{d-c}}\right)\\ \end{cases}$

・$a+b\ne0$ かつ $c+d=0$ のケース

$d=-c$ より $a^3+c^3=b^3+d^3$ は $a^3+c^3=b^3-c^3$ ということになるため
$2c^3=b^3-a^3$ となり、$(1)$ と $(2)$ は
$$\begin{cases} (1)~\rightarrow&ab(a-b)-2c^3+(a-b)x^2+(c-d)y^2=0&\\ &ab(a-b)-(b^3-a^3)+(a-b)x^2+(c-d)y^2=0&\cdots(1)''\\[8pt] (2)~\rightarrow&(a+b)(a-b)x=0~\rightarrow~x=0&\cdots(2)'' \end{cases}$$
$(2)''$ を $(1)''$ に代入して
$\begin{align} (c-d)y^2=&(b^3-a^3)-ab(a-b)\\ (d-c)y^2=&(a^3-b^3)+ab(a-b)=(a^2+b^2+ab+ab)(a-b)=(a+b)^2(a-b)\\ y^2=&\frac{(a+b)^2(a-b)}{d-c}\quad\left(\leftarrow~条件より~a\ne b~すなわち~c\ne d~\right)\\ \end{align}$
$\begin{cases} x=0~\left(=\pm(c+d)\sqrt{\frac{c-d}{b-a}}\right)\\ y=\pm(a+b)\sqrt{\frac{a-b}{d-c}}\\ \end{cases}$

・$a+b=0$ かつ $c+d=0$ のケース

$b=-a,~d=-c$ より $a^3+c^3=b^3+d^3$ は $a^3+c^3=-a^3-c^3$ ということになるため
$a=-c=-b=d$、すなわち (1) と (2) は
$$\begin{cases} (1)~\rightarrow&2a^3-2a^3+2ax^2-2ay^2=0&\rightarrow~x^2=y^2\\[8pt] (2)~\rightarrow&0(a-b)x+0(c-d)y=0&\rightarrow~0=0 \end{cases}$$
となることから、$x,y$ は $x^2=y^2$ を満たす範囲で任意の値をとることができ、
すなわち、このような解でも成立します。

$$\begin{cases} x=0~\left(=\pm(c+d)\sqrt{\frac{c-d}{b-a}}\right)\\ y=0~\left(=\pm(a+b)\sqrt{\frac{a-b}{d-c}}\right)\\ \end{cases}$$

全てのケースに共通する解は
$$\begin{cases} x=\pm(c+d)\sqrt{\frac{c-d}{b-a}}\\ y=\pm(a+b)\sqrt{\frac{a-b}{d-c}}\\ \end{cases}$$
であり、よってこの解は与式を $m,n$ についての恒等式として成立させることができます。

検算