#! 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]," "
x[2]=5
print " ", x[2]," "
cn=3
for i in 7.step(10000,30)
if 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
if i+4<10000 && sh(i+4,x)==1 then
if i+4<10
print " ", i+4," "
else
if i<100 then
print " ", i+4," "
else
if i<1000 then
print " ", i+4," "
else
print i+4," "
end
end
end
x[cn]=i+4
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
if i+6<10000 && sh(i+6,x)==1 then
if i+6<10
print " ", i+6," "
else
if i<100 then
print " ", i+6," "
else
if i<1000 then
print " ", i+6," "
else
print i+6," "
end
end
end
x[cn]=i+6
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
if i+10<10000 && sh(i+10,x)==1 then
if i+10<10
print " ", i+10," "
else
if i<100 then
print " ", i+10," "
else
if i<1000 then
print " ", i+10," "
else
print i+10," "
end
end
end
x[cn]=i+10
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
if i+12<10000 && sh(i+12,x)==1 then
if i+12<10
print " ", i+12," "
else
if i<100 then
print " ", i+12," "
else
if i<1000 then
print " ", i+12," "
else
print i+12," "
end
end
end
x[cn]=i+12
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
if i+16<10000 && sh(i+16,x)==1 then
if i+16<10
print " ", i+16," "
else
if i<100 then
print " ", i+16," "
else
if i<1000 then
print " ", i+16," "
else
print i+16," "
end
end
end
x[cn]=i+16
cn+=1
if cn>1 && cn%20==0 then
print "\n"
end
end
if i+22<10000 && sh(i+22,x)==1 then
if i+22<10
print " ", i+22," "
else
if i<100 then
print " ", i+22," "
else
if i<1000 then
print " ", i+22," "
else
print i+22," "
end
end
end
x[cn]=i+22
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"
参考ダウンロード添付ファイル
実行結果
バージョン１の場合
バージョン２の場合
バージョン３の場合
バージョン４の場合
バージョン５の場合
バージョン６の場合
今回
返ってスピードが落ちてしまいました。
バージョンアップ認定ならず！
検索対象を外すのは、偶数と３の倍数しか効果がないようです。
私には、残念ながらこれ以上のアイデアはありません。
これで素数探索の第10講は終わりとします。
本講義は、第１部を終了して第２部へと進みます。
第２部第11講では、
偶数版魔方陣の法則による４次魔方陣と６次魔方陣の作成に挑戦します。
第８話へ 第11講第１話へ
