"素数階段"をPARI/GPで描いてみる。
NHKスペシャルでリーマン予想に関する番組をやっていたそうで、そこで素数階段というものをやっていたそうだ(こういう表現からわかるように、見てません)。素数階段は自然数を横軸に、縦軸にその自然数までにある素数の数を取ったもの(らしい)。PARI/GPは素晴らしいことに、素数を求めるルーチンを書く必要すらなくて、たとえば素数の和なんかはforprime文で簡単にかけるし、prime(n)でn番目の素数を求められたりする。ということで素数階段を求めてみた。プログラムは一行(もっと簡単にできると思うけど、めんどくさい)。
for (i=2,1000,c=0;forprime(x=2,i,c=c+1;);write("prime.txt",i,"\t", c););
で、1000までの素数階段は、
160000とかまでは(時間かかるよ)
こんな感じ。素数定理からのn/log(n)と比較してます。
なんて簡単にできるんだろう。夢のよう。
« 新田辺付近のNTTのアンテナ | トップページ | Pepsi NEXのCMは懐かしい曲ばかり(Part2) »
「学問・資格」カテゴリの記事
- Interface2025年8月号Pythonで体験!はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10/10、RCR wireless newsのウェビナー2件(6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
- 高周波・RFニュース 2025年7月2日 5G Americasが6Gに向けセンシングと通信ホワイトペーパー発行、KYOCERA AVXが3dBハイブリッドカプラ発表、TDKが車載薄膜インダクタ発表、Nordicが1次電池向けPMIC発表、ローデ・シュワルツの6GとAI/ML解説記事(2025.07.02)
- 高周波・RFニュース 2025年6月30日 QualcommがAIを用いた6Rxアンテナ解説、Next G Allianceと日本のXGMFが5G,6Gで協力、5G Americasが25Q1で5G加入者増加と発表、TechInsigtsがHuawei Pura 80 Pro+分解、Qorvoが5-7GHzのWi-Fi 7 FEM発表(2025.06.30)
- Gemini CLIが使えるようになっていたので早速VSCodeのターミナルから使って、JavaScriptで連立一次方程式を計算するコードを書いてもらった。普通にガウスの消去法で計算するhtmlを作ってくれた。(2025.06.27)
コメント