#! ruby -Ks
hj=Time.now
def sh(a)
if a==1 then
return 0
end
b=a**(1/2.0)
i=0
while x[i]<=b
if a%x[i]==0 then
return 0
end
i+=1
end
return 1
end
x=[2]
print " ", x[0]," "
x[1]=3
print " ", x[1]," "
cn=2
for i in 5.step(10000,2)
if i%3!=0 && sh(i,x)==1 then
if i<10
print " ", i," "
else
if i<100 then
print " ", i," "
else
if i<1000 then
print " ", i," "
else
print i," "
end
end
end
x[cn]=i
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
end
ow=Time.now
print "\n"
print "見つかった素数の個数は",cn,"個です。\n"
print "素数探索時間は",ow-hj,"秒です。\n"
参考ダウンロード添付ファイル
実行結果
バージョン１の場合
バージョン２の場合
バージョン３の場合
バージョン４の場合
バージョン５の場合
今回
残念ながら、バージョンアップとは認定できません。
今回の試みを改善するために、
奇数を３つに分類します。
３，９，１５，２１，２７，３３，・・・
５，１１，１７，２３，２９，３５，・・・
７，１３，１９，２５，３１，３７，・・・
分類基準は何であるか分かりますか。
答は３０行下。
答え
３，９，１５，２１，２７，３３，・・・
５，１１，１７，２３，２９，３５，・・・
７，１３，１９，２５，３１，３７，・・・
第１グループ−−−３で割ったときに余りが０（すなわち３の倍数）
第２グループ−−−３で割ったときの余りが２
第３グループ−−−３で割ったときの余りが１
検索対象は、第２グループと第３グループに実で良いですね。
この点に注目すると、
for i in 5.step(10000,2)
を
for i in 5.step(10000,6)
とすることが出来ます。
改良しましょう。
結構難しいですよ。
第６話へ 第８話へ
