« 大戸屋で大戸屋風チキン南蛮 チキン1.5倍をいただく。ものすごいタルタル+チキンの量でかなりサクサクチキンだった。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2023年12月19日 WRC-23で新周波数帯は1GHz以下,3.5GHz,6GHz,6Gは7~8.5GHz帯を検討。IEEE Microwave Magazineでミリ波の軌道角運動量、Microwave Journalでミリ波3Dスキャナ、QSPICE受賞、GSA今年のまとめ、など。 »

2023年12月18日 (月)

NHK笑わない数学でバーチ・スウィンナートン=ダイアー予想(BSD予想)をやっていたので当時EDSACIIで計算したことをPythonで計算(その2) 現在わかっている最大ランク28(以上)とランク20の楕円曲線で計算をしてみる。確かに傾きが同じだ!

前回は番組で紹介していたランク0,1,2,3の楕円曲線で計算してみた。

NHK笑わない数学でバーチ・スウィンナートン=ダイアー予想(BSD予想)をやっていたので当時EDSACIIで計算したことをPython+SymPy+Numbaで計算。番組で取り上げられていた4つの楕円曲線がランクごとに傾きが予想通りになった! 

ところで今、楕円曲線のランクっていくつまで分かっているんだろう?と思って調べると

History of elliptic curves rank records 

Rank of an elliptic curve

最大はランク28以上のこれ。

y2 + xy + y = x3 − x2 − 20067762415575526585033208209338542750930230312178956502x + 34481611795030556467032985690390720374855944359319180361266008296291939448732243429

はっきりわかっていて大きいのはランク20らしい。

y2 + xy + y = x3  x2 -

244537673336319601463803487168961769270757573821859853707x + 961710182053183034546222979258806817743270682028964434238957830989898438151121499931

めちゃくちゃ係数がでかいが、Pythonな任意桁数の整数が扱えるから楽勝…と思ったらなぜかNumbaがエラーを出す。

Numbaなら10000までの素数で1分半くらいで1条件計算できるのに…仕方ないのでNumba外して計算。2条件で5時間以上かかった…

Bsdpython4

結果はこんな感じ。前回やったランク0,1,2,3も一緒に描いてみた。確かに傾きは予想通り!

Bsdpython3

 

« 大戸屋で大戸屋風チキン南蛮 チキン1.5倍をいただく。ものすごいタルタル+チキンの量でかなりサクサクチキンだった。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2023年12月19日 WRC-23で新周波数帯は1GHz以下,3.5GHz,6GHz,6Gは7~8.5GHz帯を検討。IEEE Microwave Magazineでミリ波の軌道角運動量、Microwave Journalでミリ波3Dスキャナ、QSPICE受賞、GSA今年のまとめ、など。 »

パソコン・インターネット」カテゴリの記事

学問・資格」カテゴリの記事

日記・コラム・つぶやき」カテゴリの記事

映画・テレビ」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« 大戸屋で大戸屋風チキン南蛮 チキン1.5倍をいただく。ものすごいタルタル+チキンの量でかなりサクサクチキンだった。 | トップページ | 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2023年12月19日 WRC-23で新周波数帯は1GHz以下,3.5GHz,6GHz,6Gは7~8.5GHz帯を検討。IEEE Microwave Magazineでミリ波の軌道角運動量、Microwave Journalでミリ波3Dスキャナ、QSPICE受賞、GSA今年のまとめ、など。 »

最近の記事

最近のコメント

2024年12月
1 2 3 4 5 6 7
8 9 10 11 12 13 14
15 16 17 18 19 20 21
22 23 24 25 26 27 28
29 30 31        
フォト
無料ブログはココログ