Differences between revisions 5 and 6
Revision 5 as of 2018-10-10 15:05:03
Size: 3588
Comment:
Revision 6 as of 2018-10-10 15:11:45
Size: 3857
Comment:
Deletions are marked like this. Additions are marked like this.
Line 19: Line 19:
  *   * 与えられた整数以下の素数を列挙するのであれば,[[https://ja.wikipedia.org/wiki/エラトステネスの篩|エラトステネスの篩(ふるい)]]を使うのがよいでしょう
   * 後日,リストを取り上げてからやりましょう

基本問題(8)

  • 演習(8)

    基本課題(8a)

    【関数定義】

    試し割り法による素数判定

    BP_02_008a_judge_prime_number_by_trial_division.py

  • 解答は基本的に,下記 6 行です.
    • 2~√xまでの整数で順に割ってみて,割り切れたら,(1と)自分自身以外に約数があることを意味する.
      • 逆に言えば,割り切れる整数が,その範囲内に存在しなければ,素数であることを意味する
    • x**0.5は,平方根です.したがって,mathモジュールのsqrt()関数を使ってもかまいません
    • 与えられた整数以下の素数を列挙するのであれば,エラトステネスの篩(ふるい)を使うのがよいでしょう

      • 後日,リストを取り上げてからやりましょう

   1 def is_prime_number( x ):
   2     last = int( x ** 0.5 )
   3     for i in range( 2, last + 1 ):
   4         if is_divisor( i, x ):
   5             return( False )
   6     return( True )
  • ここで,is_divisor()関数は,単に,約数かどうか判定する以下の2 行の関数ですが,関数にするまでもないですね.

   1 def is_divisor( d, x ):
   2     return( x % d == 0 )
  • なので,もうすこし短く書けそうです.

   1 def is_prime_number( x ):
   2     for i in range( 2, int( x ** 0.5 + 1 ) ):
   3         if x % d == 0 :
   4             return( False )
   5     return( True )

   1 #!/usr/bin/env python
   2 # -*- coding: utf-8 -*-
   3 # ==============================================================================
   4 # * Copyright (c) 2018 IIJIMA, Tadashi
   5 # *       (IIJIMA Laboratory, Dept. of Science and Technology, Keio University).
   6 # ==============================================================================
   7 # ソフトウェア工学[02] 基本課題[02]-(008a)  BP_02_008a_judge_prime_number_by_trial_division.py
   8 # BP(Basic Problem) 02-008a: 試し割による素数判定
   9 #        2018-10-03 飯島 正 (iijima@ae.keio.ac.jp)
  10 # ==============================================================================
  11 # dはxの約数である
  12 def is_divisor( d, x ):
  13     return( x % d == 0 )
  14 # ==============================================================================
  15 # 2 ~ √xの間に,約数が存在しないなら,xは素数
  16 def is_prime_number( x ):
  17     last = int( x ** 0.5 )
  18     for i in range( 2, last + 1 ):
  19         if is_divisor( i, x ):
  20             return( False )
  21     return( True )
  22 # ==============================================================================
  23 # ===== 【メイン・プログラム】  =====
  24 # ----- オープニングメッセージ -----
  25 print( "試し割り法によって,素数か否かを判定する関数: " )
  26 print( "  50以下の素数列を表示するために使ってみる. " )
  27 print( "  指定の数以下の素数を列挙するには,「エラトステネスの篩」を使う方が効率がよいが," )
  28 print( "  その実装は,リストを取り上げてから,後日,行う" )
  29 
  30 # ----- 計算と実行結果の表示 -----
  31 for i in range( 2, 50 ):
  32     if is_prime_number( i ):
  33         print( i, end=", " )
  34 print()
  35 # ==============================================================================


SoftEng: Python/Prog/Practice/Basic/02/BP_008a (last edited 2018-10-28 14:44:43 by TadashiIijima)