第13講 素数探索
第5話 素数探索−素数で割っていくバージョン
前話解答例

Dim cn As Integer
Dim s(10000000) As Integer
Private Sub CommandButton1_Click()
  
  Call CommandButton2_Click
  Dim i As Integer, a As Integer, cns As Integer
  Dim hj As Single, ow As Single
  hj = Timer
  cn = 0
  For i = 1 To 10000
    If sh(i) = 1 Then
      a = cn Mod 10
      cns = Int(cn / 10)
      Cells(6 + cns, 1 + a) = i
      s(cn) = i
      cn = cn + 1
    End If
  Next
  ow = Timer
  Cells(4, 1) = "素数生成にかかった時間は"
  Cells(4, 5) = ow - hj
  Cells(4, 7) = "秒です。"
  Cells(5, 1) = "生成された素数は"
  Cells(5, 4) = cn
  Cells(5, 5) = "個です。"

End Sub
Function sh(a As Integer)

  If a = 1 Then
    sh = 0
    Exit Function
  End If
  If a = 2 Then
    sh = 1
    Exit Function
  End If
  If a Mod 2 = 0 Then
    sh = 0
    Exit Function
  End If
  Dim q As Integer, i As Integer
  q = Int(Sqr(a))
  For i = 0 To cn
    If s(i) > q Then
      sh = 1
      Exit Function
    End If
    If a Mod s(i) = 0 Then
      sh = 0
      Exit Function
    End If
  Next
  sh = 1
  
End Function

Private Sub CommandButton2_Click()

  Range("4:2000").Select
  Selection.ClearContents
  Cells(1, 1).Select
  
End Sub
実行例

参考ダウンロード添付ファイル
奇数で割っていくバージョン
as
は、素数を1229個生成するのに0.054688秒ですから、
少しですが確かにスピードアップしています。

それぞれをさらに、次のように変更して実験すると
奇数で割っていくバージョン
Private Sub CommandButton1_Click()
  
  Call CommandButton2_Click
  Dim i As Long, a As Long, cns As Long, cn As Long
  Dim hj As Single, ow As Single
  hj = Timer
  cn = 0
  For i = 1 To 10000000
    If sh(i) = 1 Then
      a = cn Mod 200
      cns = Int(cn / 200)
      Cells(6 + cns, 1 + a) = i
      cn = cn + 1
    End If
  Next
  ow = Timer
  Cells(4, 1) = "素数生成にかかった時間は"
  Cells(4, 5) = ow - hj
  Cells(4, 6) = "秒です。"
  Cells(5, 1) = "生成された素数は"
  Cells(5, 4) = cn
  Cells(5, 5) = "個です。"

End Sub
Function sh(a As Long)

  If a = 1 Then
    sh = 0
    Exit Function
  End If
  If a = 2 Then
    sh = 1
    Exit Function
  End If
  If a Mod 2 = 0 Then
    sh = 0
    Exit Function
  End If
  Dim o As Long, i As Long
  o = Int(Sqr(a))
  For i = 3 To o Step 2
    If a Mod i = 0 Then
      sh = 0
      Exit Function
    End If
  Next
  sh = 1
  
End Function

Private Sub CommandButton2_Click()

  Range("4:50000").Select
  Selection.ClearContents
  Cells(1, 1).Select
  
End Sub
参考ダウンロード添付ファイル

素数で割ってバージョン
Dim cn As Long
Dim s(10000000) As Long
Private Sub CommandButton1_Click()
  
  Call CommandButton2_Click
  Dim i As Long, a As Long, cns As Long
  Dim hj As Single, ow As Single
  hj = Timer
  cn = 0
  For i = 1 To 10000000
    If sh(i) = 1 Then
      a = cn Mod 200
      cns = Int(cn / 200)
      Cells(6 + cns, 1 + a) = i
      s(cn) = i
      cn = cn + 1
    End If
  Next
  ow = Timer
  Cells(4, 1) = "素数生成にかかった時間は"
  Cells(4, 5) = ow - hj
  Cells(4, 6) = "秒です。"
  Cells(5, 1) = "生成された素数は"
  Cells(5, 4) = cn
  Cells(5, 5) = "個です。"

End Sub
Function sh(a As Long)

  If a = 1 Then
    sh = 0
    Exit Function
  End If
  If a = 2 Then
    sh = 1
    Exit Function
  End If
  If a Mod 2 = 0 Then
    sh = 0
    Exit Function
  End If
  Dim q As Long, i As Long
  q = Int(Sqr(a))
  For i = 0 To cn
    If s(i) > q Then
      sh = 1
      Exit Function
    End If
    If a Mod s(i) = 0 Then
      sh = 0
      Exit Function
    End If
  Next
  sh = 1
  
End Function

Private Sub CommandButton2_Click()

  Range("4:50000").Select
  Selection.ClearContents
  Cells(1, 1).Select
  
End Sub
参考ダウンロード添付ファイル




とさらに差が広がりました。

さらに、改良しましょう。
2以外の偶数は、偶数ではありませんから、
2以外は、奇数のみを素数であるか、
判定するプログラムに変更して下さい。




第4話へ 第6話へ
004

eclipse c++ 入門
魔方陣 数独で学ぶ VBA 入門
数独のシンプルな解き方・簡単な解法の研究
vc++講義へ
excel 2013 2010 2007 vba入門へ
VB講義基礎へ
初心者のための世界で一番わかりやすいVisual C++入門基礎講座へ
初心者のための世界で一番わかりやすいVisual Basic入門基礎講座へ
専門用語なしの C言語 C++ 入門(Visual C++ 2010で学ぶ C言語 C++ 入門)
専門用語なしの excel vba マクロ 入門 2013 2010 2007 対応講義 第1部
eclipse java 入門へ
excel 2016 vba 入門へ
小学生からエンジニアまでのRuby入門へ
本サイトトップへ