5次方程式の解の公式を求める

いまから5次方程式の解の公式を書き下してみます.途中で2次・3次・4次方程式を解く過程がありますが,それらの解の公式は既知のものとして扱います.5次方程式の解は一般に根号と四則演算によって書くことができないことが証明されているため,特殊関数を用いる必要があります.今回はその中でも楕円積分とヤコビのツェータ関数を使った表示を調べました.5次方程式の解の公式だけを見たい場合は一番最後の「まとめ」を御覧ください.


I. ニュートンの恒等式

方程式を解くための準備として,ニュートンの恒等式というものを見てみる.

■ニュートンの恒等式
x1,,xnn 個の変数とするとき, x1m++xnm=(1)m1mSm+i=1m1(1)m+i1Smi(x1i++xni) が成り立つ.ここで, Smm 番目の基本対称式を表すものとする.すなわち S1=x1++xnS2=x1x2+x1x3++xn1xnS3=x1x2x3+x1x2x4++xn2xn1xnSn=x1x2x3xn である.

以降 x1i++xni のような n 個の変数に対する和は Σ(xki) のように略記する(和を取るときの束縛変数は k とする). 以上のニュートンの恒等式を方程式の解に対して適用してみよう.いま n 次方程式 xn+an1xn1++a1x+a0=0 が与えられているとする.この方程式は複素数の範囲に重複度を込めて n 個の解を持ち, それを x1,,xn とする.方程式の係数は n 個の解の基本対象式と関連があり, S1=an1S2=an2Sm=(1)manmSn1=(1)n1a1Sn=(1)na0 のようになっている.したがってニュートンの恒等式により Σ(xkm)=(1)m1mSm+i=1m1(1)m+i1SmiΣ(xki)=manmi=1m1aniΣ(xki) が得られる.これは解の m 乗和と方程式の係数との間の関係を表す式になっている.5次方程式を簡単化していく過程で Σ(xk)Σ(yk3) などのようなものを大量に計算する必要があるため,ニュートンの恒等式を明示しておく必要があった.


II. 5次方程式の簡単化

5次方程式は四則演算と根号を用いた解の公式を持たないため,特殊関数による表示を与えるほかない.そのためにはまず5次方程式を 簡単な形にしていかなければならない. まず5次方程式 ax5+bx4+cx3+dx2+ex+f=0 は両辺を a0 で割ることにより,最高次係数は 1 の場合だけを考えればいいことになる.また方程式 x5+ax4+bx3+cx2+dx+e=0 に対して, x=ya/5 という変数変換をすることにより, y5+(b2a25)y3+(c3ab5+4a325)y2+(d2ac5+3a2b253a4125)y+(ead5+a2c25a3b125+4a53125)=0 のように4次の項を消すことができる.すなわち y5+py3+qy2+ry+s=0 という形の5次方程式を考えればよいことになる.実際にこの方程式を簡単な形に変形することを試みてみよう. まずは上の形から z5+Az2+Bz+C=0 という方程式へ変形する.その方法としては, y5+py3+qy2+ry+s=0 の5つの解を yk(k=1,,5)z5+Az2+Bz+C=0 の解を zk(k=1,,5) として,この2つの方程式の解の間に zk=yk2+αyk+β という関係があることを仮定することである.もし何かしらの方法によって α,β,A,B を決定することができれば,実際にそのような関係が存在することが分かるのである.そこで上のような関係が2つの方程式の解の間に存在したと仮定する.まず z5+Az2+Bz+C=0 の解に対して Σ(zk)=0Σ(zk2)=0Σ(zk3)=3AΣ(zk4)=4BΣ(zk5)=5C となることを計算により確かめることができる. いま zk=yk2+αyk+β であるから (1)Σ(yk2+αyk+β)=0(2)Σ((yk2+αyk+β)2)=0(3)Σ((yk2+αyk+β)3)=3A(4)Σ((yk2+αyk+β)4)=4B(5)Σ((yk2+αyk+β)5)=5C という式を得ることができる.式(1),(2)から 5β2p=04r6αq+2p24βp2α2p+5β2=0 という α,β に関する2次方程式が得られる.これにより α=3q2p±(3q2p)2+3p52rpβ=2p5 という解を得ることができる.すなわち y5+py3+qy2+ry+s=0z5+Az2+Bz+C=0 の 2つの方程式の解の間には zk=yk2+(3q2p±(3q2p)2+3p52rp)yk+2p5 という関係が存在することが分かった.また式(3)(4)(5)により, 3A=15αs+6pr12βr12α2r+3q2+15αpq18αβq3α3q2p3+6βp2+6α2p26β2p6α2βp+5β34B=8qs+28αps60αβs20α3s+4r2+28αqr8p2r+24βpr+36α2pr24β2r48α2βr4α4r8pq2+12βq2+18α2q228αp2q+60αβpq+20α3pq36αβ2q12α3βq+2p48βp312α2p3+12β2p2+24α2βp2+2α4p28β3p12α2β2p+5β45C=5s2+45αrs20pqs+40βqs+80α2qs45αp2s+140αβps+70α3ps150αβ2s100α3βs5α5s10pr2+20βr2+40α2r210q2r90αpqr+140αβqr+70α3qr+10p3r40βp2r80α2p2r+60β2pr+180α2βpr+30α4pr40β3r120α2β2r20α4βr15αq3+15p2q240βpq280α2pq2+30β2q2+90α2βq2+15α4q2+45αp3q140αβp2q70α3p2q+150αβ2pq+100α3βpq+5α5pq60αβ3q30α3β2q2p5+10βp4+20α2p420β2p360α2βp310α4p3+20β3p2+60α2β2p2+10α4βp210β4p20α2β3p+5β5 なので,これによって A,B,C を決定することができる.したがってより簡単化された方程式 z5+Az2+Bz+C=0 を解くことができれば,上記の解の関係 zk=yk2+αyk+β を経由して y を求めることができるようになる.したがって5次方程式は z5+Az2+Bz+C=0 という形にまで簡単化された. ここからさらに, w5+Pw+Q=0 という形にまで簡単化しようとするわけであるが,その方法の起点はいまとさほど変わらない.ただし,途中で出てくる方程式に少し技巧が必要となる.いま z5+Az2+Bz+C=0 の方程式の解を zk とし, w5+Pw+Q=0 の解を wk とおいたときに,それらの解の間に wk=zk4+ηzk3+θzk2+κz+μ という関係が存在すると仮定する.このとき Σ(wk)=0Σ(wk2)=0Σ(wk3)=0Σ(wk4)=4PΣ(wk5)=5Q という式を用いて (6)Σ(zk4+ηzk3+θzk2+κz+μ)=0(7)Σ((zk4+ηzk3+θzk2+κz+μ)2)=0(8)Σ((zk4+ηzk3+θzk2+κz+μ)3)=0(9)Σ((zk4+ηzk3+θzk2+κz+μ)4)=4P(10)Σ((zk4+ηzk3+θzk2+κz+μ)5)=5Q という5つの方程式を得ることができる.式(6)より 5μ3Aη4B=0μ=3Aη+4B5 となる.式(7)より (15Aθ20Bη25C)κ10Bθ225Cηθ(11)+15A2θ+3A2η2+23ABη+20AC+2B2=0 このとき,求めるべき変数が η,θ,κ,μ の4つに対して,それらを求めるための方程式が(6)(7)(8)の3つだけであることを考慮すると,解には1次元分の自由度が生まれることになる.そこでいま θ を上記の κ の項が 0 になるように定めることとする.すなわち 15Aθ20Bη25C=0 となるように θ を定めると θ=4Bη+5C3A となる.したがって(11)は 300ABCη2160B3η2+27A4η2+375AC2η400B2Cη+27A3Bη250BC245A3C+18A2B2=0 と変形され,この η の2次方程式を解くことによって η を求めることもできる.そして式(8)より 675A3κ33375A2Cηκ2+3600AB2ηκ2+4500ABCκ2+2025A4κ26000B2Cη2κ675A3Bη2κ15000BC2ηκ+4050A3Cηκ7200A2B2ηκ9375C3κ9675A2BCκ2025A5κ+225A2BCη3+320AB3η3+54A5η3+1125A2C2η2+3900AB2Cη2+960B4η2+756A4Bη2+9375ABC2η+2400B3Cη1485A4Cη+3843A3B2η+6250AC3+1500B2C2+4770A3BC+108A2B3+675A6=0 という κ の3次方程式が得られる. θ を上のように選んだのは,この κ の方程式が高次にならないようにするためであった.したがって η,θ,κ,μ の4つが方程式の係数により記述できたので, z5+Az2+Bz+C=0w5+Pw+Q=0 の解の間には,4次式の関係があることが分かった.また式(9)(10)により PQ も得ることができる.最終的な式は最後のまとめに記す.これにより,当初の5次方程式は w5+Pw+Q=0 という形にまで簡単化できることが判明した.この形はBring-Jerrard型と呼ばれる形である.ここで w=P4u とすることにより,方程式は P54u5+PP4u+Q=0P54u5P54u+Q=0u5u+QP54=0 という形にすることができる. R=Q/P54 とすれば,最終的に5次方程式は u5u+R=0 という形にまで持っていくことができる.


III. モジュラー関数による解の公式

さて, u5u+R=0 という形の5次方程式の解の公式を得ることができれば,最初の x5+ax4+bx3+cx2+dx+e=0 という方程式の解まで導出できることは分かった.問題は u5u+R=0 をどうやって解けばよいのかである.それには様々な方法が存在するが,ここではエルミートによるモジュラー関数を使った方法を紹介する. まずは,解を記述するために使うモジュラー関数を定義する.モジュラー関数は,一つの関数を指す用語ではなく,ある保型的性質を満たすような関数の総称である.今回は楕円関数から話を始める.まずパラメーター ττ=iK(1k2)K(k) とする.これは楕円関数論においてhalf-period ratioと呼ばれるものである.また k はモジュラスと呼ばれるパラメータである.さらに K(k)=0π211k2sin2φdφ はモジュラス k に対する完全楕円積分である.これを用いて 楕円モジュラー関数 φ(τ)=k4ψ(τ)=1k28 が定義される.ちなみにこの2つの関数 φ,ψ は,ヤコビのテータ関数 ϑ を使って φ(τ)=ϑ10(0;τ)ϑ00(0;τ)ψ(τ)=ϑ01(0;τ)ϑ00(0;τ) と表す事もできる.このとき関数のペア φ(τ),φ(5τ) は有理数体上代数的に従属であることが知られている.すなわち Φ5(φ(τ),φ(5τ))0 を満たすような有理数係数多項式 Φ5 が存在する.具体的にそれは Φ5(x,y):=x6y6+5x2y2(x2y2)+4xy(1x4y4) という形の6次多項式である.このとき Φ5(φ(τ),y) の6つの解は yj(τ)=φ(τ+16j5)(j=0,1,2,3,4)y5(τ)=φ(5τ) である.この6つの解を使って定義される J(τ)=(y5(τ)y0(τ))(y4(τ)y1(τ))(y3(τ)y2(τ))={φ(5τ)+φ(τ5)}{φ(τ+165)φ(τ+645)}{φ(τ+325)φ(τ+485)} という τ の関数は, J(τ)52000φ(τ)4ψ(τ)16J(τ)+16005φ(τ)3ψ(τ)16(1+φ(τ)8)=0 を任意の τ に対して満たすことが示されている.詳しくは参考文献「Briot-Bonquet, Théorie des fonctions elliptiques」のp.654-を参照していただきたい. 文献では J(τ) の各点の位数の勘定や, J(τ)yj(τ) に対する対称性などから u8β0J5+p=15u8βpJ5p(ap+bpu8+cpu16+)=0(u=φ(τ)) という u -級数の方程式を立て,そこから各係数 βp,ap,bp,cp, を決定していくことにより最終的な5次方程式を得ている. この5次方程式 T52000φ(τ)4ψ(τ)16T+16005φ(τ)3ψ(τ)16(1+φ(τ)8)=0 のその他の4つの解は J(τ+16),J(τ+32),J(τ+48),J(τ+64) である.いま τ を固定したときの一つの解を J とし, J=2534φ(τ)ψ(τ)4X というように変換を行うと,初等的な計算によって X5X+2(1+φ(τ)8)554φ(τ)2ψ(τ)4=0 となる.少し複雑になってきたが,つまり J(τ+16)(=0,1,2,3,4) という形の τ の関数は (J(τ+16)2534φ(τ)ψ(τ)4)5(J(τ+16k)2534φ(τ)ψ(τ)4)+2(1+φ(τ)8)554φ(τ)2ψ(τ)4=0 を任意の τ に対して満たしているということである. したがってまとめると,最初に与えられた方程式 u5u+R=0 に対して, R=2(1+φ(τ)8)554φ(τ)2ψ(τ)4 となるような τ の値を求めることができれば, u5u+R=0 の解は u=J(τ+16)2534φ(τ)ψ(τ)4(=0,1,2,3,4) というように求めることができることが分かった.したがって今度は R=2(1+φ(τ)8)554φ(τ)2ψ(τ)4 を満たすような τ を求めることになる.これを変形すると, 2(1+φ(τ)8)=R554φ(τ)2ψ(τ)44(1+φ(τ)8)2=R255φ(τ)4ψ(τ)8 これをモジュラス k を使って表すと 4(1+k2)2=R255k(1k2)4k4+R255k3+8k2R255k+4=0k4+(R5542)2k3+2k2(R5542)2k+1=0 となる.この方程式の解は三角関数の恒等式 tan4θ+4sin4θtan3θ+2tan2θ4sin4θtanθ+1=0 を用いることにより, sin4θ=4(2R554)2 を満たす θ に対して k=tanθ というように求めることができる.すなわち k=tan(14arcsin(16R255)) である.このモジュラス k を使って τ=iK(1k2)K(k) となる. したがって最終的に u5u+R=0 の解は u=J(τ+16)2534φ(τ)ψ(τ)4(=0,1,2,3,4)τ=iK(1k2)K(k),k=tan(14arcsin(16R255)) と求めることができる.



IV. まとめ

以上をまとめると,最初の x5+ax4+bx3+cx2+dx+e=0 の解の公式は以下のように与えられる.

x=ya5

y=α2+(α2)2β+z

α=3q2p±(3q2p)2+3p52rp

β=2p5

p=b2a25

q=c3ab5+4a325

r=d2ac5+3a2b253a4125

s=ead5+a2c25a3b125+4a53125

3A=3α(5pq5s)+6pr+(3β+3α2)(2p24r)+3q23(6αβ+α3)q2p32(3β2+3α2β)p+5β3

4B=4α(7ps+7qr7p2q)+8qs+(12αβ+4α3)(5pq5s)+4r2+(4β+6α2)(6pr+3q22p3)8p2r+(6β2+12α2β+α4)(2p24r)8pq23(12αβ2+4α3β)q+2p42(4β3+6α2β2)p+5β4

5C=5s2+5α(9rs9p2s18pqr3q3+9p3q)+(5β+10α2)(8qs+4r28p2r8pq2+2p4)+(20αβ+10α3)(7ps+7qr7p2q)20pqs+(30αβ2+20α3β+α5)(5pq5s)10pr2+(10β2+30α2β+5α4)(6pr+3q22p3)10q2r+10p3r+(10β3+30α2β2+5α4β)(2p24r)+15p2q23(20αβ3+10α3β2)q2p52(5β4+10α2β3)p+5β5

P54u=z4+ηz3+θz2+κz+μ

μ=3Aη+4B5

θ=4Bη+5C3A

0=d3κ3+d2κ2+d1κ+d0

d0=(225A2BC+320AB3+54A5)η3+(1125A2C2+3900AB2C+960B4+756A4B)η2+(9375ABC2+2400B3C1485A4C+3843A3B2)η+6250AC3+1500B2C2+4770A3BC+108A2B3+675A6

d1=(6000B2C675A3B)η2+(15000BC2+4050A3C7200A2B2)η9375C39675A2BC2025A5

d2=(3375A2C+3600AB2)η+4500ABC+2025A4

d3=675A3

η=e12e2+(e12e2)2e0e2

e0=250BC245A3C+18A2B2

e1=375AC2400B2C+27A3B

e2=300ABC160B3+27A4

4P=(8AC+4B2)(θ4+12μθ2+12ηκθ2+12η2μθ+12κ2θ+6μ2+24ηκμ+6η2κ2)+3A2(4μθ3+6κ2θ2+12μ2θ+24ηκμθ+6η2μ2+12κ2μ+4ηκ3)+7AB(4κθ3+12ημθ2+24κμθ+12ηκ2θ+12ημ2+12η2κμ+4κ3)+(9BC3A3)(4ηθ3+12κθ2+24ημθ+12η2κθ+12κμ+4η3μ+12ηκ2)+(5C210A2B)(4θ3+6η2θ2+12μθ+24ηκθ+12η2μ+6κ2+4η3κ)4B(6μ2θ2+12κ2μθ+4μ3+12ηκμ2+κ4)5C(12κμθ2+12ημ2θ+4κ3θ+12κμ2+12ηκ2μ)+(11A2C11AB2)(12ηθ2+12κθ+4η3θ+12ημ+12η2κ)+(24ABC4B3+3A4)(6θ2+12η2θ+4μ+12ηκ+η4)3A(12κμ2θ+4ημ3+4κ3μ)+(13AC213B2C+13A3B)(12ηθ+4κ+4η3)+(14BC2+14A3C+21A2B2)(4θ+6η2)+5μ4+4(5C3+45A2BC+15AB33A5)η+24A2C2+48AB2C+4B416A4B

5Q=(5C210A2B)(θ5+20μθ3+20ηκθ3+30η2μθ2+30κ2θ2+30μ2θ+120ηκμθ+30η2κ2θ+30η2μ2+30κ2μ+20η3κμ+20ηκ3)+(8AC+4B2)(5μθ4+10κ2θ3+30μ2θ2+60ηκμθ2+30η2μ2θ+60κ2μθ+20ηκ3θ+10μ3+60ηκμ2+30η2κ2μ+5κ4)+(9BC3A3)(5κθ4+20ημθ3+60κμθ2+30ηκ2θ2+60ημ2θ+60η2κμθ+20κ3θ+30κμ2+10η3μ2+60ηκ2μ+10η2κ3)+(11A2C11AB2)(5ηθ4+20κθ3+60ημθ2+30η2κθ2+60κμθ+20η3μθ+60ηκ2θ+30ημ2+60η2κμ+10κ3+10η3κ2)+(24ABC4B3+3A4)(5θ4+10η2θ3+30μθ2+60ηκθ2+60η2μθ+30κ2θ+20η3κθ+10μ2+60ηκμ+5η4μ+30η2κ2)+3A2(10μ2θ3+30κ2μθ2+20μ3θ+60ηκμ2θ+5κ4θ+10η2μ3+30κ2μ2+20ηκ3μ)+7AB(20κμθ3+30ημ2θ2+10κ3θ2+60κμ2θ+60ηκ2μθ+20ημ3+30η2κμ2+20κ3μ+5ηκ4)+(13AC213B2C+13A3B)(20ηθ3+30κθ2+10η3θ2+60ημθ+60η2κθ+20κμ+20η3μ+30ηκ2+5η4κ)+(14BC2+14A3C+21A2B2)(10θ3+30η2θ2+20μθ+60ηκθ+5η4θ+30η2μ+10κ2+20η3κ)4B(10μ3θ2+30κ2μ2θ+5μ4+20ηκμ3+5κ4μ)5C(30κμ2θ2+20ημ3θ+20κ3μθ+20κμ3+30ηκ2μ2+κ5)+(5C3+45A2BC+15AB33A5)(30ηθ2+20κθ+20η3θ+20ημ+30η2κ+η5)+(24A2C2+48AB2C+4B416A4B)(10θ2+30η2θ+5μ+20ηκ+5η4)3A(20κμ3θ+5ημ4+10κ3μ2)+(51ABC2+17B3C17A4C34A3B2)(20ηθ+5κ+10η3)+(18AC3+27B2C272A3BC36A2B3+3A6)(5θ+10η2)+5μ5+5(19BC338A3C2114A2B2C19AB4+19A5B)η+5C4120A2BC280AB3C+20A5C4B5+50A4B2

u=J(τ+16)253k(1k2)24(=0,1,2,3,4)

τ=iK(1k2)K(k)

k=tan(14arcsin(16R255))

R=QP54

J(τ)=(y5(τ)y0(τ))(y4(τ)y1(τ))(y3(τ)y2(τ))

yj(τ)=φ(τ+16j5)(j=0,1,2,3,4)

y5(τ)=φ(5τ)

φ(ω)=ϑ10(0;ω)ϑ00(0;ω)

ψ(ω)=ϑ01(0;ω)ϑ00(0;ω)

ϑab(z;ω)=n=eπi(n+a)2ω+2πi(n+a)(z+b)(a,bR)

K(k)=0π211k2sin2φdφ

ここで ϑ はヤコビのツェータ関数, K は第一種完全楕円積分という特殊関数である.5次方程式の解の公式はいざというときに必要になる可能性が僅かに存在するかもしれないので,記憶しておきましょう.



V. 参考文献

(1) V.S.Adamchik, David.J.Jeffrey, Polynomial Transformations of Tschirnhaus, Bring and Jerrard, ACM SIGSAM Bullentin, Vol37, No.3 (2003)
5次方程式をBring-Jerrard型にまで簡単化する方法が説明されています.

(2) Semjon Adlaj, Multiplication and division on elliptic curves, torsion points and roots of modular equations
レベル5の楕円モジュラー多項式について参考にしました.

(3) Briot, Bouquet, Théorie des fonctions elliptiques, p654-660
楕円モジュラー方程式から5次の標準形に持っていく部分の証明が載っています.

コメント