4

巡回対称性による1~4次方程式の統一的導出

169
1
$$$$

巡回対称性による1~4次方程式の統一的導出

はじめに

 四半世紀ほど前、外回り営業で車を運転中に夜空の星座のようなイメージが頭に浮かんできました。星々の配置関係を対称性をもって並び替えたらどうなるだろうという思いつきから、気づけば$n$乗根や$n$次方程式の解という発想に広がり、もしかしたらその巡回対称性から逆算すれば、統一的な手法で解の公式を導出できるのでは?などとおぼろげに考えたのです。

 そのアイデアから少なくとも四次方程式までは実際に導出できることを確認するに至り、その過程を再現したものが本稿となります。数学を専門的に学んだことがないもので、当時はカルダノもフェラーリもよく知らず(調べたのはここ10年の間くらい)、せいぜい二次方程式の解の公式を知っていた程度でした。ですので、私個人にとっては斬新なアイデアだと思う一方で、私が思いつく程度のことですから、もしかしたらとっくに先人が通ってきた道だったとしても何ら不思議はないでしょう。先行研究があるかどうかは存じませんが、新規性があろうとなかろうと独力で導出できたという感動体験を記録として残しておきたいと考え、この記事を書き上げました。

 というわけで、巡回対称性という視点から、一次方程式から四次方程式までを共通の枠組みで導出する過程を示します。

$1$$n$次方程式の解を表す関数を予想

$1$$n$次方程式の一般形は
$$\sum_{j=0}^na_jx^{n-j}=0\quad(a_j\ne0)$$
と表すことができ、複素数範囲では高々$n$個の解に集約されることが知られています。それらの解は
この方程式に対して互いに優劣なく完全な対称性を持つと思われますので、$n$次方程式の解のうち
インデックス$m_{(\operatorname{mod}n)}$における一般解を表す関数$x(m_{(\operatorname{mod}n)})$を次のように予想してみました。
$$\begin{align} &\mathrm{P}=e^{2\pi i}(=1)\text{として}\\ &x(m_{(\operatorname{mod}n)})=\sum_{k=0}^{n-1}s_k\cdotp\left(\mathrm{P}^\frac kn\right)^m\quad\text{(予想)}\\ \end{align}$$
$s_k$$a_j$を用いて『全ての$a_j$が実数だった場合$s_k$も実数となるような式
(この記事内では「実式」と呼ぶことにします)で表せる』という予想です。
先の方程式はモニックとなるように式全体の倍率を変えても解に影響しないため、
$$\begin{align} &a_0=1\text{として}\\ &\sum_{j=0}^na_jx^{n-j}=\prod_{m=0}^{n-1}\left(x-\sum_{k=0}^{n-1}s_k\mathrm{P}^\frac{km}n\right)=0\quad\text{(予想)} \end{align}$$
のように表現可能だろうと考えました。

少なくとも$n\lt5$においては有限回の四則演算&冪根範囲で解の公式を導出できることが知られていますので、
この予想を元に$1$次から$4$次までの各$s_k$を求められるかどうか、検証してみましょう。
$ $

$n=1$の場合($1$次方程式)

$$\begin{align} &x+a_1\\ =&\prod_{m=0}^0\left(x-s_0\mathrm{P}^\frac{0m}1\right)\\ =&\left(x-(s_0)\right)\\ =&x~\underbrace{-s_0}_{a_1}\\ =&0\\ \end{align}$$
$$\begin{align} s_0&=-a_1\\ &=\left(-\frac{a_1}1\right)\ \end{align}$$
$$\begin{align} &x(m_{(\operatorname{mod}1)})=\left(-\frac{a_1}1\right)\mathrm{P}^\frac{0m}1\\ \end{align}$$
$ $
$x+b=0$ の解
$\begin{align} x&=\left(-\frac b1\right)\\ &=-b\\ \end{align}$
$ $
$ax+b=0$ の解
$\begin{align} x &=-\left(\frac ba\right)\\ &=-\frac ba\\ \end{align}$
$ $

$n=2$の場合(2次方程式)

$$\begin{align} &x^2+a_1x+a_2\\ =&\prod_{m=0}^1\left(x-(s_0\mathrm{P}^\frac{0m}2+s_1\mathrm{P}^\frac{1m}2)\right)\\ =&\left(x-(s_0+s_1\mathrm{P}^\frac12)\right)\left(x-(s_0+s_1\mathrm{P}^\frac22)\right)\\ =&(x-(s_0+s_1))(x-(s_0-s_1))\\ =&x^2~\underbrace{-2s_0}_{a_1}x~\underbrace{+(s_0^2-s_1^2)}_{a_2}\\ =&0\\ \end{align}$$
$$\begin{align} s_0&=\left(-\frac{a_1}2\right)\\ s_1^2&=s_0^2-a_2\\ &=\left(-\frac{a_1}2\right)^2+2\left(-\frac{a_2}2\right)\\ s_1&=\left[\left(-\frac{a_1}2\right)^2+2\left(-\frac{a_2}2\right)\right]^\frac12\quad\text{(実式)}\\ \end{align}$$
$$\begin{align} &x(m_{(\operatorname{mod}2)})\\ &=\left(-\frac{a_1}2\right)\mathrm{P}^\frac{0m}2+\left[\left(-\frac{a_1}2\right)^2+2\left(-\frac{a_2}2\right)\right]^\frac12\mathrm{P}^\frac{1m}2\\ &=\left(-\frac{a_1}2\right)+\left[\left(-\frac{a_1}2\right)^2+2\left(-\frac{a_2}2\right)\right]^\frac12\mathrm{P}^\frac{m}2\\ \end{align}$$
$ $
$x^2+bx+c=0$ の解
$$\begin{align} x &=\left(-\frac b2\right)\pm\sqrt{\left(-\frac b2\right)^2+2\left(-\frac c2\right)}\\ &=\frac{-b\pm\sqrt{b^2-4c}}2\\ \end{align}$$
$ $
$ax^2+bx+c=0$ の解
$$\begin{align} x &=\frac{-\left(\frac ba\right)\pm\sqrt{\left(\frac ba\right)^2-4\left(\frac{c}a\right)}}{2}\\ &=\frac{-b\pm\sqrt{b^2-4ac}}{2a}\\ \end{align}$$
$ $

$n=3$の場合(3次方程式)

$$\begin{align} &x^3+a_1x^2+a_2x+a_3\\ =&\prod_{m=0}^2\left(x-(s_0\mathrm{P}^\frac{0m}3+s_1\mathrm{P}^\frac{1m}3+s_2\mathrm{P}^\frac{2m}3)\right)\\ =&\left(x-(s_0+s_1+s_2)\right)\left(x-(s_0+s_1\mathrm{P}^\frac13+s_2\mathrm{P}^\frac23)\right)\left(x-(s_0+s_1\mathrm{P}^\frac23+s_2\mathrm{P}^\frac43)\right)\\ =&x^3~\underbrace{-(3s_0)}_{a_1}x^2~~\underbrace{+3(s_0^2-s_1s_2)}_{a_2}x~\underbrace{-(s_0^3+s_1^3+s_2^3-3s_0s_1s_2)}_{a_3}\\ =&0\\ \end{align}$$
$s_1$$s_2$は可換で、値を入れ替えても(式の並び替えで)結果的に同一の式となることを確認しておきます。
$$\begin{align} s_0&=\left(-\frac{a_1}3\right)\\ s_1s_2&=s_0^2-\frac{a_2}3\\ &=\left(-\frac{a_1}3\right)^2+\left(-\frac{a_2}3\right)\\ s_1^3+s_2^3&=3s_0s_1s_2-s_0^3-a_3\\ &=3\left(-\frac{a_1}3\right)\left(\left(-\frac{a_1}3\right)^2+\left(-\frac{a_2}3\right)\right)-\left(-\frac{a_1}3\right)^3+3\left(-\frac{a_3}3\right)\\ &=2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\\ \end{align}$$
$(S^3-s_1^3)(S^3-s_2^3)=(S^3)^2~\underbrace{-(s_1^3+s_2^3)}_{a'_1}~S^3\underbrace{+s_1^3s_2^3}_{a'_2}=0$ より
$$\begin{align} S^3(m'_{(\operatorname{mod}2)}) &=\left(-\frac{a'_1}2\right)+\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\mathrm{P}^\frac{m'}2\\ \end{align}$$
$\mathrm{P}^\frac{m'}2$は、$m'$が偶数なら$\mathrm{P}^\frac{m'}2=+1$$m'$が奇数なら$\mathrm{P}^\frac{m'}2=-1$です。$s_1^3$$s_2^3$$S^3(m'_{(\operatorname{mod}2)})$における
 $2$次共役解であるため$s_k^3$と表したときの$k$を利用してどちらかをどちらかに対応させます。$s_1$$s_2$は可換であり
 どちらに対応させても順不同で結果は同じであるため、ここでは便宜上$m'=k-1$を採用することにします。
$$\begin{align} s_{k(\in\{1,2\})}^3 &=\left(-\frac{a'_1}2\right)+\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\mathrm{P}^\frac{k-1}2\\ s_{k(\in\{1,2\})} &=\left[\left(-\frac{a'_1}2\right)+\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\mathrm{P}^\frac{k-1}2\right]^\frac13\quad\text{(実式)}\\ \end{align}$$
$$\begin{cases} a_1'&=-(s_1^3+s_2^3)\\ &=-\left(2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\right)\\ a_2'&=(s_1s_2)^3\\ &=\left(\left(-\frac{a_1}3\right)^2+\left(-\frac{a_2}3\right)\right)^3\\ \end{cases}$$
$$\begin{cases} s_0&=\left(-\frac{a_1}3\right)\\ s_1&=\left[\left(-\frac{a'_1}2\right)+\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\right]^\frac13\\ s_2&=\left[\left(-\frac{a'_1}2\right)-\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\right]^\frac13\\ \end{cases}$$
$$\begin{align} x(m_{(\mathrm{mod}3)})&=\left(-\frac{a_1}3\right)\mathrm{P}^\frac{0m}3\\ &+\left[\left(-\frac{a'_1}2\right)+\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{1m}3\\ &+\left[\left(-\frac{a'_1}2\right)-\left[\left(-\frac{a'_1}2\right)^2+2\left(-\frac{a'_2}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{2m}3\\ &=\left(-\frac{a_1}3\right)\\ &+\left[\left(-\frac{-\left(2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a_1}3\right)^2+\left(-\frac{a_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{+m}3\\ &+\left[\left(-\frac{-\left(2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{a_1}3\right)^3+3\left(\left(-\frac{a_1}3\right)\left(-\frac{a_2}3\right)+\left(-\frac{a_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a_1}3\right)^2+\left(-\frac{a_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{-m}3\\ \end{align}$$
$ $
$x^3+bx^2+cx+d=0$ の解($1$$3$乗根とその共役を$\omega$$\overline{\omega}$で表します)
$$\begin{align} x&=\left(-\frac{b}3\right)\\ &+\omega\sqrt[3]{\left(-\frac{-\left(2\left(-\frac{b}3\right)^3+3\left(\left(-\frac{b}3\right)\left(-\frac{c}3\right)+\left(-\frac d3\right)\right)\right)}2\right)+\sqrt{\left(-\frac{-\left(2\left(-\frac{b}3\right)^3+3\left(\left(-\frac{b}3\right)\left(-\frac{c}3\right)+\left(-\frac d3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{b}3\right)^2+\left(-\frac c3\right)\right)^3}2\right)}}\\ &+\overline{\omega}\sqrt[3]{\left(-\frac{-\left(2\left(-\frac{b}3\right)^3+3\left(\left(-\frac{b}3\right)\left(-\frac{c}3\right)+\left(-\frac d3\right)\right)\right)}2\right)-\sqrt{\left(-\frac{-\left(2\left(-\frac{b}3\right)^3+3\left(\left(-\frac{b}3\right)\left(-\frac{c}3\right)+\left(-\frac d3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{b}3\right)^2+\left(-\frac c3\right)\right)^3}2\right)}}\\ &=-\frac{b}3\\ &+\frac{\omega\sqrt[3]{ \left(-b^3+\frac{9(bc-3d)}2\right)+\sqrt{\left(-b^3+\frac{9(bc-3d)}2\right)^2-\left(b^2-3c\right)^3}}}3\\ &+\frac{\overline{\omega}\sqrt[3]{ \left(-b^3+\frac{9(bc-3d)}2\right)-\sqrt{\left(-b^3+\frac{9(bc-3d)}2\right)^2-\left(b^2-3c\right)^3}}}3\\ &=-\frac{b}3\\ &+\frac{\omega\sqrt[3]{ \frac{\left(9(bc-3d)-2b^3\right)+\sqrt{\left(9(bc-3d)-2b^3\right)^2-4\left(b^2-3c\right)^3}}2}}3\\ &+\frac{\overline{\omega}\sqrt[3]{ \frac{\left(9(bc-3d)-2b^3\right)-\sqrt{\left(9(bc-3d)-2b^3\right)^2-4\left(b^2-3c\right)^3}}2}}3\\ \end{align}$$
$ $
$ax^3+bx^2+cx+d=0$ の解($1$$3$乗根とその共役を$\omega$$\overline{\omega}$で表します)
$$\begin{align} x&=-\frac{b}{3a}\\ &+\frac{\omega\sqrt[3]{ \left(-b^3+\frac{9a(bc-3ad)}2\right)+\sqrt{\left(-b^3+\frac{9a(bc-3ad)}2\right)^2-\left(b^2-3ac\right)^3}}}{3a}\\ &+\frac{\overline{\omega}\sqrt[3]{ \left(-b^3+\frac{9a(bc-3ad)}2\right)-\sqrt{\left(-b^3+\frac{9a(bc-3ad)}2\right)^2-\left(b^2-3ac\right)^3}}}{3a}\\ &=-\frac{b}{3a}\\ &+\frac{\omega\sqrt[3]{ \frac{\left(9a(bc-3ad)-2b^3\right)+\sqrt{\left(9a(bc-3ad)-2b^3\right)^2-4\left(b^2-3ac\right)^3}}2}}{3a}\\ &+\frac{\overline{\omega}\sqrt[3]{ \frac{\left(9a(bc-3ad)-2b^3\right)-\sqrt{\left(9a(bc-3ad)-2b^3\right)^2-4\left(b^2-3ac\right)^3}}2}}{3a}\\ \end{align}$$
$ $

$n=4$の場合(4次方程式)

$$\begin{align} &x^4+a_1x^3+a_2x^2+a_3x+a_4\\ =&\prod_{m=1}^4\left(x-(s_0\mathrm{P}^\frac{0m}4+s_1\mathrm{P}^\frac{1m}4+s_2\mathrm{P}^\frac{2m}4+s_3\mathrm{P}^\frac{3m}4)\right)\\ =&(x-(s_0+s_1+s_2+s_3))(x-(s_0+s_1\mathrm{P}^\frac14+s_2\mathrm{P}^\frac24+s_3\mathrm{P}^\frac34))(x-(s_0+s_1\mathrm{P}^\frac24+s_2\mathrm{P}^\frac44+s_3\mathrm{P}^\frac64))(x-(s_0+s_1\mathrm{P}^\frac34+s_2\mathrm{P}^\frac64+s_3\mathrm{P}^\frac94))\\ =&(x-(s_0+s_1+s_2+s_3))(x-(s_0+s_1i-s_2-s_3i))(x-(s_0-s_1+s_2-s_3))(x-(s_0-s_1i-s_2+s_3i))\\ =&x^4\\ &\underbrace{-\left(4s_0\right)}_{a_1}x^3\\ &\underbrace{+2\left(3s_0^2-s_2^2-2s_1s_3\right)}_{a_2}x^2\\ &\underbrace{-4\left(s_0^3-s_0s_2^2+s_1^2s_2+s_2s_3^2-2s_0s_1s_3\right)}_{a_3}x\\ &\underbrace{+\left(s_0^4-s_1^4+s_2^4-s_3^4+4(s_0s_1^2s_2-s_0^2s_1s_3+s_0s_2s_3^2-s_1s_2^2s_3)-2(s_0^2s_2^2-s_1^2s_3^2)\right)}_{a_4}\\ =&0\\ \end{align}$$
$s_1$$s_3$は可換で、値を入れ替えても(式の並び替えで)結果的に同一の式となることを確認しておきます。
$$\begin{align} s_0&=\left(-\frac{a_1}4\right)\\ s_2^2+2s_1s_3&=3s_0^2-\frac{a_2}2\\ (s_1^2+s_3^2)s_2 &=-s_0^3+s_0(s_2^2+2s_1s_3)-\frac{a_3}4\\ a_4&=s_0^4-s_1^4+s_2^4-s_3^4+4(s_0s_1^2s_2-s_0^2s_1s_3+s_0s_2s_3^2-s_1s_2^2s_3)-2(s_0^2s_2^2-s_1^2s_3^2)\\ &=s_0^4-s_1^4+s_2^4-s_3^4+4s_0s_1^2s_2-4s_0^2s_1s_3+4s_0s_2s_3^2-4s_1s_2^2s_3-2s_0^2s_2^2+2s_1^2s_3^2\\ &=\textcolor{#f37}{4s_2^4-4s_2^4-8s_1s_2^2s_3}+s_0^4+s_2^4\textcolor{#f37}{+4s_1^2s_3^2+4s_1s_2^2s_3}+4s_0s_1^2s_2+4s_0s_2s_3^2-2s_0^2s_2^2-4s_0^2s_1s_3-s_1^4-s_3^4\textcolor{#f37}{-2s_1^2s_3^2}\\ &=\textcolor{#f37}{4s_2^4-4s_2^2(s_2^2+2s_1s_3)}+s_0^4+(s_2^2\textcolor{#f37}{+2s_1s_3})^2+4s_0s_2(s_1^2+s_3^2)-2s_0^2(s_2^2+2s_1s_3)-(s_1^2+s_3^2)^2\\ a_4\cdotp\frac{s_2^2}4 &=s_2^6-(s_2^2+2s_1s_3)s_2^4+\frac{s_0^4-2(s_2^2+2s_1s_3)s_0^2+4(s_1^2+s_3^2)s_2s_0+(s_2^2+2s_1s_3)^2}4s_2^2-\frac{(s_1^2+s_3^2)^2s_2^2}4\\ 0&=s_2^6-(s_2^2+2s_1s_3)s_2^4+\frac{s_0^4-2(s_2^2+2s_1s_3)s_0^2+4(s_1^2+s_3^2)s_2s_0+(s_2^2+2s_1s_3)^2-a_4}4s_2^2-\frac{(s_1^2+s_3^2)^2s_2^2}4\\ \end{align}$$
$(s_2^2)^3~\underbrace{-(s_2^2+2s_1s_3)}_{a'_1}~(s_2^2)^2\underbrace{+\frac{s_0^4-2(s_2^2+2s_1s_3)s_0^2+4(s_1^2+s_3^2)s_2s_0+(s_2^2+2s_1s_3)^2-a_4}4}_{a'_2}(s_2^2)\underbrace{-\frac{(s_1^2+s_3^2)^2s_2^2}4}_{a'_3}=0$ より
$$\begin{align} (s_2^2)(m'_{(\operatorname{mod}3)})&=\left(-\frac{a'_1}3\right)\\&+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{+m'}3\\ &+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\mathrm{P}^\frac{-m'}3 \\ \end{align}$$
$m'$は任意の整数ですが、$s_2$が実式であることを担保するためここでは$m'\equiv0\pmod3$としておきます。
$$\begin{align} s_2^2 &=\left(-\frac{a'_1}3\right)\\ &+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\\ &+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\\ s_2&=\left[ \left(-\frac{a'_1}3\right) +\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13 +\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13 \right]^\frac12\quad\text{(実式)} \end{align}$$
$$\begin{align} s_1s_3&=\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)-s_2^2}2\\ \end{align}$$
$$\begin{cases} s_2\neq0~\text{の場合}\\ \quad s_1^2+s_3^2=\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2}\\ s_2=0~\text{の場合}\\ \begin{align} a_4&=s_0^4-s_1^4-s_3^4-4s_0^2s_1s_3+2s_1^2s_3^2\\ &=s_0^4-(s_1^4+s_3^4)-4s_0^2(s_1s_3)+2(s_1s_3)^2\\ &=s_0^4-((s_1^2+s_3^2)^2-2(s_1s_3)^2)-4s_0^2(s_1s_3)+2(s_1s_3)^2\\ &=s_0^4-(s_1^2+s_3^2)^2-4s_0^2(s_1s_3)+4(s_1s_3)^2\\ (s_1^2+s_3^2)^2 &=s_0^4-4s_0^2(s_1s_3)+4(s_1s_3)^2-a_4\\ s_1^2+s_3^2 &=\left[s_0^4-4s_0^2\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)}2\right)+4\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)}2\right)^2-a_4\right]^\frac12\quad\text{(実式)}\\ \end{align} \end{cases}$$

$(S^2-s_1^2)(S^2-s_3^2)=(S^2)^2-\left(s_1^2+s_3^2\right)S^2+s_1^2s_3^2$ より
$$\begin{align} S^2(m''_{(\operatorname{mod}2)})&=\left(-\frac{-(s_1^2+s_3^2)}2\right)+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2-(s_1s_3)^2\right]^\frac12\mathrm{P}^\frac{m''}2\\ \end{align}$$
$\mathrm{P}^\frac{m''}2$は、$m''$が偶数なら$\mathrm{P}^\frac{m''}2=+1$$m''$が奇数なら$\mathrm{P}^\frac{m''}2=-1$です。$s_1^2$$s_3^2$$S^2(m''_{(\operatorname{mod}2)})$における
 $2$次共役解であるため$s_k^2$と表したときの$k$を利用してどちらかをどちらかに対応させます。$s_1$$s_3$は可換であり
 どちらに対応させても順不同で結果は同じであるため、ここでは便宜上$m''=\frac{k-1}2$を採用することにします。
$$\begin{align} s_{k(\in\{1,3\})}^2&=\left(-\frac{-(s_1^2+s_3^2)}2\right)+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\mathrm{P}^\frac{\frac{k-1}2}2\\ s_{k(\in\{1,3\})}&=\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\mathrm{P}^\frac{k-1}4\right]^\frac12\quad\text{(実式)}\\ \end{align}$$
$$\begin{cases} a'_1&=-(s_2^2+2s_1s_3)\\ &=-\left(3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)\right)\\ a'_2&=\frac{s_0^4-2\left(s_2^2+2s_1s_3\right)s_0^2+4\left(s_1^2+s_3^2\right)s_2s_0+\left(s_2^2+2s_1s_3\right)^2-a_4}4\\ &=\frac{\left(-\frac{a_1}4\right)^4-2\left(3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)\right)\left(-\frac{a_1}4\right)^2+4\left(2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)\right)\left(-\frac{a_1}4\right)+\left(3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)\right)^2-a_4}4\\ a'_3&=-\frac{\left(s_1^2+s_3^2\right)^2s_2^2}4\\ &=-\frac{\left(2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)\right)^2}4\\ \end{cases}$$
$$\begin{cases} s_0&=\left(-\frac{a_1}4\right)\\ s_1&=\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\right]^\frac12\\ s_2&=\left[ \left(-\frac{a'_1}3\right) +\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13 +\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{a'_1}3\right)^3+3\left(\left(-\frac{a'_1}3\right)\left(-\frac{a'_2}3\right)+\left(-\frac{a'_3}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{a'_1}3\right)^2+\left(-\frac{a'_2}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13 \right]^\frac12\\ s_3&=\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)-\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\right]^\frac12\\ \end{cases}$$
$$\begin{cases} s_2\neq0~\text{の場合}\\ \quad s_1^2+s_3^2=\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2}\\ s_2=0~\text{の場合}\\ \begin{align} \quad s_1^2+s_3^2 &=\left[s_0^4-4s_0^2(s_1s_3)+4(s_1s_3)^2-a_4\right]^\frac12\quad\text{(実式)}\\ &=\left[\left(-\frac{a_1}4\right)^4-4\left(-\frac{a_1}4\right)^2\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)}2\right)+4\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)}2\right)^2-a_4\right]^\frac12 \end{align} \end{cases}$$
$$\begin{align} x(m_{(\mathrm{mod}4)})&=\left(-\frac{a_1}4\right)\mathrm{P}^\frac{0m}4\\ &+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\right]^\frac12\mathrm{P}^\frac{1m}4\\ &+s_2\mathrm{P}^\frac{2m}4\\ &+\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)-\left[\left(-\frac{-(s_1^2+s_3^2)}2\right)^2+2\left(-\frac{(s_1s_3)^2}2\right)\right]^\frac12\right]^\frac12\mathrm{P}^\frac{3m}4\\ &=\left(-\frac{a_1}4\right)\\ &+i\left[\left(-\frac{-\left(\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2} \right)}2\right)+\left[\left(-\frac{-\left(\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2} \right)}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)-s_2^2}2\right)^2}2\right)\right]^\frac12\right]^\frac12\\ &-s_2\mathrm{P}^\frac{2m}4\\ &-i\left[\left(-\frac{-\left(\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2} \right)}2\right)-\left[\left(-\frac{-\left(\frac{2\left(-\frac{a_1}4\right)^3+2\left(-\frac{a_1}4\right)\left(-\frac{a_2}4\right)+\left(-\frac{a_3}4\right)}{s_2} \right)}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{a_1}4\right)^2+2\left(-\frac{a_2}4\right)-s_2^2}2\right)^2}2\right)\right]^\frac12\right]^\frac12\\ \end{align}$$
$ $
$x^4+bx^3+cx^2+dx+e=0$ の解
$$\begin{cases} B&=-\left(3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)\right)\\ C&=\frac{\left(-\frac{b}4\right)^4-2\left(3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)\right)\left(-\frac{b}4\right)^2+4\left(2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)\right)\left(-\frac{b}4\right)+\left(3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)\right)^2-e}4\\ D&=-\frac{\left(2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)\right)^2}4\\ \alpha&=\left(-\frac{B}3\right) +\left[\left(-\frac{-\left(2\left(-\frac{B}3\right)^3+3\left(\left(-\frac{B}3\right)\left(-\frac{C}3\right)+\left(-\frac{D}3\right)\right)\right)}2\right)+\left[\left(-\frac{-\left(2\left(-\frac{B}3\right)^3+3\left(\left(-\frac{B}3\right)\left(-\frac{C}3\right)+\left(-\frac{D}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{B}3\right)^2+\left(-\frac{C}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13 +\left[\left(-\frac{-\left(2\left(-\frac{B}3\right)^3+3\left(\left(-\frac{B}3\right)\left(-\frac{C}3\right)+\left(-\frac{D}3\right)\right)\right)}2\right)-\left[\left(-\frac{-\left(2\left(-\frac{B}3\right)^3+3\left(\left(-\frac{B}3\right)\left(-\frac{C}3\right)+\left(-\frac{D}3\right)\right)\right)}2\right)^2+2\left(-\frac{\left(\left(-\frac{B}3\right)^2+\left(-\frac{C}3\right)\right)^3}2\right)\right]^\frac12\right]^\frac13\\ \end{cases}$$
[$\alpha\neq0$ の場合]
$$\begin{align} x&=\left(-\frac{b}4\right)\\ &+i\sqrt{\left(-\frac{-\left(\frac{2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)}{\sqrt\alpha}\right)}2\right)+\sqrt{\left(-\frac{-\left(\frac{2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)}{\sqrt\alpha}\right)}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)-\alpha}2\right)^2}2\right)}}\\ &-\sqrt{\alpha}\\ &-i\sqrt{\left(-\frac{-\left(\frac{2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)}{\sqrt\alpha}\right)}2\right)-\sqrt{\left(-\frac{-\left(\frac{2\left(-\frac{b}4\right)^3+2\left(-\frac{b}4\right)\left(-\frac{c}4\right)+\left(-\frac{d}4\right)}{\sqrt\alpha}\right)}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)-\alpha}2\right)^2}2\right)}}\\ &=\cdots\quad(\text{導出は完了し、以降は係数表示への代入・整理のみとなるため省略}) \end{align}$$
[$\alpha=0$ の場合]
$$\begin{align} x&=\left(-\frac{b}4\right)\\ &+i\sqrt{\left(-\frac{\sqrt{\left(-\frac{b}4\right)^4-4\left(-\frac{b}4\right)^2\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)+4\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2-e}}2\right)+\sqrt{\left(-\frac{\sqrt{\left(-\frac{b}4\right)^4-4\left(-\frac{b}4\right)^2\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)+4\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2-e}}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2}2\right)}}\\ &-i\sqrt{\left(-\frac{\sqrt{\left(-\frac{b}4\right)^4-4\left(-\frac{b}4\right)^2\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)+4\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2-e}}2\right)-\sqrt{\left(-\frac{\sqrt{\left(-\frac{b}4\right)^4-4\left(-\frac{b}4\right)^2\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)+4\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2-e}}2\right)^2+2\left(-\frac{\left(\frac{3\left(-\frac{b}4\right)^2+2\left(-\frac{c}4\right)}2\right)^2}2\right)}}\\ &=\cdots\quad(\text{導出は完了し、以降は係数表示への代入・整理のみとなるため省略}) \end{align}$$

投稿日:11日前
更新日:4時間前
数学の力で現場を変える アルゴリズムエンジニア募集 - Mathlog served by OptHub

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

https://mathlog.info/articles/323         数学を愛する会 副会長 CCO / ガラパゴ数学 開拓者 / 猫舌・甘党・薄味派

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中