数式処理ソフトによるガロア群の算出と、べき根を用いた厳密解の表現 その17

以前の記事のちょっとした補足を2つ述べます。ひとつは、各解を V で表す式の導出の仕方、もうひとつは g(x) を因数分解するときの、べき根でθ(x) を求める整合性のあるやり方についてです。

各解を V で表す式の導出・補足

以前、終結式の反復計算によって「V と同じ位置(先頭)に α を持つような Vk すべてを根に持つ多項式」を効率的に求める方法を説明しましたが、本 blog では(そして jurupapa さんの blog でも)それをどう使えば「αV で表す式」が得られるのか、具体的には説明していませんでした。以前の記事で

多項式の GCD 計算だけ(グレブナー基底を使ってもいいですね)で「αを根に持つ Q(V) 係数の 1 次多項式」を求めることができ

とだけ書いたやり方を説明します。

この計算法は三森明夫 “ガロア論文の古典的証明” で紹介されていたもので、おそらく、元々Galois が採用していた方法なのだと思います(考案者がGaloisなのか、それ以前から知られていたやり方なのかは私は知らない)。私が最初にこの辺りの話をいじくり回していた際、手計算で
α=13V2V2β=23V2+V+4γ=13V22
という式を求めた時も、実際は Lagrange 補間ではなくこの計算法を使っていました(Lagrange補間を使った方法はかなり式が複雑になり、3次方程式であっても手計算ではなかなかやる気が起きない式になります)。

この計算法ですが、同書の記述は私にとってはかなりわかりにくくて、最初は何が何だかさっぱり理解できませんでした。一連の考察がある程度進んでから再度読み返して、色々な記述が何を言おうとしているのかある程度推測できるようになって、初めてアルゴリズムが理解できました(jurupapa さんは最初から理解されていたようで、その理解力には脱帽です)。私にとって一番ハードルが高かったのは、同書では V という文字が多項式の不定元を表していることも α+2β+3γ のような「具体的な値」を表していることも両方あって、しかもすぐ近くでそのふたつの用法が断りなしに切り替えられていることがあったことでした。

そんなわけで、ここではその計算法を、私がわかりやすいと思える説明で述べてみます。やはり具体例を通じて説明することにして、方程式は例によって f(x)=x33x1=0 とし、V=α+2β+3γ とします。V の中の解 α, β, γ を対称群 S3 の元によって置換した値は
V1=α+2β+3γ=VV2=β+2γ+3αV3=γ+2α+3βV4=β+2α+3γV5=γ+2β+3αV6=α+2γ+3β
です。

まず、以前の記事でも説明した通り、Vk のうち「α を固定する置換」に対応する
V1=α+2β+3γV6=α+2γ+3β
を根に持つ多項式 Fα(x)=(xV1)(xV6) は終結式の反復計算で求められ、
Fα(x)=x2+3αx+3α23
でした。

Fα(x) の作り方から Fα(V)=0 だから、
V2+3αV+3α23=0
です。よって、t(x)=V2+3Vx+3x23 とおけば t(α)=0 だから、t(x)α を根に持つ多項式です。ここでさらに、β, γt(x) の根にならないことも言えます。その理由はと言うと…

Fα(x) を求める計算を、αβ を入れ替えて行ったとしましょう。するとそれは β+2α+3γ(=V2)β+2γ+3α(=V4) を根に持つ多項式を求めることになりますが、α, β, γ の対称性から、結果の式も αβ を入れ替えたもの、つまり Fα(x)αβ に置き換わった式になります。これはつまり
(xβ2α3γ)(xβ2γ3α)=x2+3βx+3β23
であることを意味します。この多項式の根は V2, V4 で、V=V1 は根ではありません。したがって右辺に x=V を代入しても 0 にならず、
V2+3βV+3β230
です。これより t(β)0 です。まったく同じようにして、t(γ)0 も言えます。

以上から、t(x) は「α は根に持つが、β, γ は根に持たない多項式」です。よって、f(x)t(x) の GCD は xα になります。t(x)=3x2+3Vx+V23 はその作り方から必然的に Q(V) 係数の多項式で、f(x)=x33x1Q 係数だから、f(x)t(x) で互除法によって GCD を求めれば、「α を根に持つ Q(V) 係数の 1 次多項式」が得られます。
gcd(x33x1,3x2+3Vx+V23)=(V23)x+3V+3(V39V9=0 を使って整理した)
よって αVQ 係数有理式で表す式が得られ、V39V9=0 を使ってそれを VQ 係数多項式に直せる、というわけです。
α=3V+3V23==13V2V2

βV で表す式を求めるには、「β を固定する置換に対応する Vk」を根に持つ多項式 Fβ(x) を求める所から始めて同様の手順を辿ればOKです(ここでは、終結式の反復計算は新規にやり直す必要がある。以前の計算結果を流用することはできない)。

θ(x) の求め方

以前は1次独立性に訴える議論をしていたので、1 のべき根が他数含まれるようになって「Galois群の部分群に対応する体の元」かどうかがよくわからなくなってくる状況で成立する話かどうかすっきりとは見やすくなくなっていましたが、そこをもっと簡潔にできることに気づきました。話としては「置換群としての Galois群で不変な解の多項式の値は、各解を VQ 係数多項式で表した式を代入して V について次数下げすれば V が自動的に消え去ってわかる」ということと同類の話に過ぎず、類似のことが解の「有理式」に対してもなりたつ、という話になります。「多項式」の話が多項式の係数の範囲に限定されなかったのと同様、「有理式」の方も係数の範囲にはまったくよりません。

今、
θ1(x)=h0(x)+ζh1(x)++ζp1hp1(x)
の係数は、解 α,β,K(ζ) 係数多項式として表されています。それらは下位の置換群 H ではどれも不変で、p 次巡回群 G/H の生成元に対応する置換 σ では一斉に ζ1 倍される、という関係にありました。なので、それらの比は σ で不変です。したがって、θ1(x) の最高次の係数を a とすると、θ1(x)a の係数は「G で不変な解の有理式」になっています。それらの値が求めたかったわけです。詳細は後述しますが、それには各解を V で表す式をその有理式に代入し、分母・分子それぞれを V について次数下げすると、余り難しい計算なしに求められます。特に、「分母の有理化」は必要ありません。

まず、分母、分子を次数下げして、g1(x) の次数未満にしたものをそれぞれ δ(x), ν(x) とします。ここで、g1(x)V の現在の最小多項式で、今因数分解したいと思っている式そのものであり、上位の置換群 G に対応する Vk すべてを根に持つ多項式です。さて、元の有理式が上位の置換群 G で不変であることは、

有理式 ν(x)δ(x) に、g1(x) の根を代入したら、どの根を代入しても分数としての値は同じ

ということを意味します。この「同じ」値が今求めたい値で、これを λ とおきましょう。λ は、K(ζ) の元になることだけが今わかっていて、具体的な値はまだ不明です。
ν(Vk)δ(Vk)=λ for all roots of g1(x)
すると、ν(x)λδ(x) は「g1(x) よりも低次なのに、g1(x) のどの根を代入しても 0 になる多項式」です。そのような多項式はゼロ多項式に限られます。よって多項式の等式として
ν(x)=λδ(x)
がなりたちます。すなわち有理式の等式として
ν(x)δ(x)=λ
がなりたつのです。つまり、元の解の有理式(θ(x)a の任意の係数)の分母・分子を V で表してそれぞれ次数下げした結果は、V の多項式として定数倍の関係にあり、その比 ν(V)δ(V) はただの多項式の割り算で求められます。これで、未知だった値 λ がわかるわけです(逆数 1δ(V)VK(ζ) 係数多項式として求める必要はない)。

さらに、このように

  • ν(x)δ(x) が多項式として定数倍の関係にある

ということから、次のふたつもなりたちます。

  • ν(x)δ(x) の次数が同じになる(※)
  • 両者の比 λ は、両者の最高次の係数の比に等しい

そうなっていることは、実際に ν(x)δ(x) を求める前から予めわかるので、実際の計算では ν(x)δ(x) は最高次の係数だけ求まれば十分です。つまり、次数下げで g1(x) による剰余を求める計算では、高次の係数から順に求まっていきますが、0 でない係数が見つかった時点で計算を打ち切ってよい、ということです。そうして、ν(x)δ(x) の最高次の係数(いずれも K(ζ) の元)の比を計算すれば λ が求まるわけです。ここでは「K(ζ) の元の逆数」の計算は必要ですが、より複雑な「K(ζ)[V] の元の逆数を K(ζ)[V] の元として求める」という計算は避けられる、というわけです。

※ より正確に言えば、ν(x) の方はゼロ多項式になる可能性はあるので、その場合のみ「同じ次数」ではなくなります。なので、「剰余は最高次係数だけ求めればいいから、0 でない係数が出た段階で計算を打ち切っていい」という話も、正確には δ(x) の方についてのみ言える話で、先にそちらを計算してから ν(x) の計算に入り、「δ(x) の最高次と同じ次数の係数が求まった時点で計算を打ち切っていい(もしそれが 0 だったら、それ以上計算しなくても ν(x) がゼロ多項式と確定する)」という話になります。(「δ(x) よりも高次の項で、ν(x)0 でない係数を持つことはない」ということも、実際に ν(x) を計算する前からわかる)

当然、今の計算法は θ2(x) 以降の計算にも利用できます。
θk(x)=h0(x)+(ζk)h1(x)++(ζk)p1hp1(x)
の各係数は解の K(ζ) 係数多項式の形で、いずれも H の置換で不変、σζk 倍されます。よって、引き続き aθ1(x) の最高次の係数とすれば θk(x)ak の各係数は G の置換で不変な解の有理式ですから、それぞれの値を先ほどと同じようにして計算することができます。


投稿日

カテゴリー:

,

投稿者:

タグ:

コメント

“数式処理ソフトによるガロア群の算出と、べき根を用いた厳密解の表現 その17” への2件のフィードバック

  1. みらいのアバター
    みらい

    井汲先生!YouTubeで先生の現職での動画を拝見しました。20年前に高校生の頃教わっていました。お元気ですか?全然老けないですね!!

    1. 井汲 景太のアバター
      井汲 景太

      ご覧になりましたか。ちょっと恥ずかしい(笑)。私の力不足から SEG からは離れることになって、今はああいうことをしています。当時はみらいさんのお力にはなれたでしょうか。数学がみらいさんの人生を豊かにすることにつながっていましたら幸いです。容姿は、幸い(?)あんまり変わらない体質みたいです。(白髪はぽつぽつ目につくようになってきました)

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

This site uses Akismet to reduce spam. Learn how your comment data is processed.