このワークシートは[url=https://www.geogebra.org/m/twxxx3yq]Math by Code[/url]の一部です。[br][br]楕円曲線をさらにさぐろう[br][br]今回は、有限体上での楕円曲線の解[br]です。
とりあえず実験する。[br][br]前回、楕円曲線は実平面から複素平面の格子のぷつぷつに移し替えると、[br]挙動が単純になることがわかりました。[br]今回は、楕円曲線のもとの数体を実数から有限体にぷつぷつの点に移し替えて調べてみましょう。[br][br][b][size=150]<有限体の決まりは簡単>[/size][/b][br]有限体というのは剰余類の体バージョンです。[br]群の基本となったZ3は3で割った余りによる分類でした。[br]割る数を素数にすることで、たし算だけでなく、[u]0を除けば[/u]かけ算でも群になります。[br]おまけにアーベル群で分配法則も使えるもの[br]最強のグループである体になります。[br]体といえば、有理数Q,実数R,複素数Cがあります。[br]これと同じしくみの有限バージョンです。[br]要するに、有限だから力技で調べることもできるということですね。[br]線形代数の知識もいらないし、[br]とても扱いやすいはずです。[br][br]正式な記法はたとえば、次のようになりますね。[br][b][color=#0000ff][size=150]Fp=(Z/(p),+,×) 素数pを法とする剰余類の体[br][/size][/color][/b]ゼロを除く[u]かけ算の群[/u]として使うときは、[b]F[sup]×[/sup][sub]p[/sub][/b][br]とかくこともありますね。[br][br][b][size=150]<有限体上の楕円曲線の解の個数>[br][/size][/b][br]・素数pを法として合同な整数が作る[b]体Fp上の楕円曲線Eの解の個数をNp[/b]と名付けよう。[br]たとえば、[br]曲線E:y^2=x^3-xのN2とは、[br]方程式E:y^2=x^3-xの解(x,y)を無限遠点O∞と(0,0),(0,1),(1,0),(1,1)の2^2+1=5点から数えたものですね。[br]曲線E:y^2=x^3-xのN3とは、[br]Eの解をO∞と(x,y) for x in range(3) for y in range(3)の3^2+1=10個から数えたものですね。[br][br]・無限遠点O∞はn等分点でもO+O=Oとなることから、フリーパスでE[n]の仲間に入りました。[br]理由は違いますが、楕円曲線Eは上下対称で上も下も∞につきぬけてます。この点をO∞としたわけです。[br]調べる前から「[b]O∞は楕円曲線の解の1つとなる[/b]」と言えます。[br]だから、解の候補の1つであり、チェックしなくても最後に調べた個数に1をたすことで[br][b]正しいNp[/b]が求められるということですね。[br][br][b][size=150][color=#9900ff]課題:geogebraのZip、Keepif, Lengthなどを使って、Np(pが素数)を求めよう。[br][/color][/size][/b][br]たとえば、次のようにするとスライダーで調べられます。[br][br]タイトル「有限体Fp上の楕円曲線の解となる点数Np」[br]Ps={2,3,5,7,11,13,17,19,23,29,31,37,41,43,47}[br]idx=slider(1,15,1) #どの素数を選ぶかのインデックス[br]n=Element(Ps,idx) #選んだ素数p[br]Rs= Sequence(k,k,0,n-1) #pを法とする剰余リスト [br]PPs=flatten(Zip(Zip((p,q), p, Rs),q,Rs)) #テスト前の(x,y) for x in range(p) for y in range(p) グレーでうすく[br]passed=KeepIf(y(A)^2==x(A)^3-x(A),A,PPs) "緑のまま[br]Np=Length(passed)[b]+1[/b] [br]text1="N" +n + "=" + Np +"" ♯サイズを大にして表示する。[br][br]実際にやてみると[br]N2=3 [br]N3=3[br]N5=3[br]N7=3[br]エラーは起きませんがなんか変ですね。[br]x(A)がおかしいわけではありませんでした。[br]Remが多項式か数値かで2通りの動きしかできないことが原因です。[br]でも多項式の剰余が欲しいわけではないので、きれいにRemを使って代数的にプログラムするのではなく、泥臭く、数値リストを追いかけることにしましょう。[br]PPs=flatten(Zip(Zip((p,q), p, Rs),q,Rs)) #これがベースになります。n^2個のリストです。[br]PPsx=Zip(x(k),k,PPs)[br]PPsy=Zip(y(k),k,PPs)[br]Es=Zip(E(p,q),p,PPsx,q,PPsy) #E値です。[br]Esp=Zip(Mod(k,n),k,Es) #E値のpによる剰余です。[br]Passed=KeepIf(k==0,k,Esp) ♯解0だけのリストです。[br]NP=Length(Passed)[b]+1[/b] #無限遠点O∞の1個を入れるのをお忘れなく。 [br][b]text1="N" +n + "=" + NP +"" #NPの部分だけ修正します。[/b][br][br]変顔をマークしよう。[br][br]N2 = 3[br]N3 = 4[br][b]N5 = 8[/b][br]N7 = 8[br]N11= 12 [br][b]N13= 8[br]N17= 16[br][/b]N19= 20[br]N23= 24[br][b]N29= 40[br][/b]N31= 32[br][b]N37= 40[br]N41= 32[br][/b]N43= 44[br]N47= 48[br][br]素数pが2、または、4で割って3余る数のときは[b]Np=p+1ピッタになる[/b]。[br]でも、[br]素数pが5,13,17,29,37,41のように4で割って1余る型だとズレができる。[br]このズレをap=(p+1)ーNpと定義しよう。[br]apに何か決まりとかはあるのだろうか。[br]これは見たらすぐわかるね。[br][b]素数pが4で割って1余る型ならどうなる?[br]ズレap=(p+1)-Np=2qになることをガウスが見つけていた。[br][br]しかも、pの平方数の和分解とのつながりもつけて。[br](qは負の数でもよい。ただし、p=q^2+r^2と平方数の和に分解したときに、q+r-1が4の倍数になる)[br][br][size=150]<振り返り>[br][/size][/b]ガウスの定理をたしかめよう。[br]5+1-8=2(-1)、5=(-1)^2+2^2から奇数q=-1,r=2 とするとq+r-1=0は4の倍数[br]13+1-8=2(3)、13=3^2+2^2から、奇数q=3,r=2 とするとq+r-1=4は4の倍数[br]17+1-16=2*1、 17=1^2+4^2から、奇数q=1,r=4とするとq+r-1=4は4の倍数[br]............[br]キリがないですね。[br]証明はしないが、面白い結果だ。[br]ガウスの定理というと、発散の積分で有名だけど、数論にもあるんだね。[br][br]面白いのはガウスの定理だけではない。[br]今は曲線E:y^2=x^3-xでの実験でしたが、[br]楕円曲線を変更してもp+1になったりちょっとずれたりと似たような結果になるようです。[br]詳しい研究は、ガウスの100年後の「[b]ハッセの定理[/b]」があります。[br][br]「どんな楕円曲線の[b][color=#0000ff]解の個数のズレapの絶対値は2√p[/color][/b]に必ず収まる」[br][br][b][size=200][color=#0000ff]|a[sub]p[/sub]|≦2√p[br][/color][/size][/b][br]この定理からわかるように、解の個数[b]Np[/b]はおよその範囲は決まるのですが、[br]いくつになるかのカンタンなアルゴリズムがなく、[br]予測はできるけれど、正確な個数を当てるには泥臭く数えるしかない。[br]しかも、特定の点を何倍もした位置を逆算するのは、スーパーコンピュータでも何百年もかかるという[br]ことがわかっています。[br]この、「有限体上の楕円曲線が持つ複雑さ」が[br]スマホでLINEを送ったり、ネットショッピングをしたりするときの[br]通信を守る「[b]楕円曲線暗号(ECC)[/b]」につながるらしいのです。[br][br]暗号は、素因数分解の難しさを利用した[b]公開鍵暗号[/b]だけではないんですね。[br]
さて、次に進むとしたら、やはり複素平面に広げることですよね。[br][br]でも、その前に足場を固めましょう。[br]日々、このようなことに触れている人からしたら、超基本レベルかもしれませんが、[br]たまにしかやらない人やふらっとやりだした人からすると、このページの内容自体に[br]違和感を感じるかもしれないからです。[br][br]次に進むために現在地の足場を固めよう。[br][br][b][size=150]<4で割って1余る素数が平方数の和に分解できること>[/size][/b][br]さきほどのガウスの定理で、4で割って1余る素数が平方数の和に分解できることを[br]さりげなく使ってます。[br]この土台を一応確認しますね。[br]整数を虚数単位iを使って拡大した環、ガウスの整数環Z[i]の世界、[br]a+bi (for a,b in Z for i is imaginary unit) では[br]4で割って1余る素数は、この世界では、分解できてしまいます。[br]5 = (2+i)(2-i),[br]13 = (3+2i)(3-2i),[br]17 = (4+i)(4-i) [br]これらを計算すると(a+bi)(a-bi)=a^2+b^2と平方数の和で表されるのでしたね。[br]つまり、2+i,2-i,3+2i,3-2i,4+i,4-iが素数のようなもの、分解できないもの(素の元、素元)です。[br]ということは、4で割って3余る素数の方は、この世界でも分解できなもの、素元だと推定できますね。[br]たとえば、素数3,7,...は平方数0,1,4,9の和では表せませんよね。[br]だから、3=(a+bi)(a-bi)とは分解できないということです。[br][br]わたしたちが素数として1まとめにしている数たちが、[br]ガウス整数の世界では4を法とした剰余で水と油のように、きっちりと[br]2分されているということがわかりますね。[br]くわしくはこちら、[url=https://www.geogebra.org/m/twxxx3yq#material/pfhtmgn7]https://www.geogebra.org/m/twxxx3yq#material/pfhtmgn7[br][br][/url][b][size=150]<そもそも楕円曲線をしらべることが何につながるのか>[/size][/b][br]楕円曲線という名前自体が楕円積分とからむという理由で背負った名前らしいけれども、[br]そもそも楕円曲線は何につながるのか?[br]これがわからないと、モチベーションが上がりませんよね。[br]ズバリ、楕円曲線の解は数論の方程式の解とつながるのです。[br]さっきのガウスの手法は整数Zを拡大して環を作るというやり方で代数学の手法だった。[br][br][b][size=150]フェルマーの大定理 [br]x^n + y^n=z^n( for x,y,z, n in Z ,n≧3)となるx,y,zはない。[/size][/b] (1630年代)[br]ピタゴラス数のオンパレードみたいなことは起こりませんよ。[br]でも、私フェルマーはたぶん無限降下法で証明できた。でも、余白が少ないから書きませんよ。[br][br]この式もガウス流の素元分解することができれば、[br]証明できるのではないかと考えたクンマーという数学者がいた。[br]しかし、この代数的な試みはあと少しのところまでいったが証明まではできなかった。[br]1のn乗根をgeogebra流にζn=(1;2π/n)とかくと、x^nがζnを利用した数の積に直せる。[br]このζnを使って整数を拡大したZ[ζn]の素元を探求すれば証明できるかもしれない、[br]というアイデアからクンマーは『素イデアル論』を樹立しました。[br][br]しかし、このチャレンジのおかげで、p進数、オイラーのζ関数を関連づける景色が生まれ、[br]それが楕円曲線へとつながり、最終的には、1994年9月のワイルズWilesによるフェルマーの大定理の証明につながったのです。[br][br][color=#0000ff][b][size=150]チャレンジは無駄ではなかった。[br][/size][/b][/color][br]ワイルズの証明の超概略は、斎藤毅さんによる整理によると、次の流れにあるといいます。[br][b][size=150]フェルマー方程式「x^n + y^n=z^n( for x,y,z, n in Z ,n≧3)」[br]の自明解はx、y、zの少なくとも1つは0だ。[br][/size][/b]もし フェルマーの方程式が自明解以外をもつ(反例)とするならば、[br]それに対応する楕円曲線が定義できる。[br]楕円曲線が定義できると、それにつながる保型形式があるはず。[br]ところが、そのような保型形式は存在しないので、自明解はない。[br]という背理法で説明できるのです。[br]もうすこし、歴史に寄せて語るとこんな感じです。[br][br][*][b][color=#0000ff]1.フェルマーの反例から「フライの楕円曲線」を作る。[br][/color][/b][/*][*][b][color=#0000ff]2.谷山・志村予想「すべての楕円曲線は保型形式から作れる」をワイルズが証明した。[br][/color][/b][/*][*][b][color=#0000ff]3.1と2から、フライの楕円曲線に対応する保型形式が「存在するはず」となる。[br][/color][/b][/*][*][b][color=#0000ff]4.しかし、「リベットの定理」により、そのような性質を持つ保型形式は「存在し得ない(矛盾)」。[/color][/b][/*][*][br][/*]楕円曲線から保型形式に至る部分が超絶難所の連続なので、そこはスルーします。[br][br]ただ、楕円曲線の解自体の分析周りをすることで、登山せずに[br]ドローンから見た登山者のようすが多少でも感じられるかもしれない。[br]そんな希望を持てたらいいですね。[br][br]そのために、今回も含めた4回連続で楕円曲線を扱ってきたのです。[br]代数的数論(群論)、曲線の解析、位相、線形代数[br]ここからは、数学の総力戦です。[br][br][color=#0000ff][b][size=150]でも、チャレンジは無駄になりません。[br]明確化や新しいつながりを見つけられるでしょう。[br][/size][/b][/color]