結果

問題 No.308 素数は通れません
ユーザー cielciel
提出日時 2015-12-02 03:03:15
言語 Ruby
(3.3.0)
結果
AC  
実行時間 100 ms / 1,000 ms
コード長 1,557 bytes
コンパイル時間 343 ms
コンパイル使用メモリ 11,336 KB
実行使用メモリ 15,576 KB
最終ジャッジ日時 2023-08-18 13:15:36
合計ジャッジ時間 14,163 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 100 ms
15,380 KB
testcase_01 AC 95 ms
15,400 KB
testcase_02 AC 94 ms
15,344 KB
testcase_03 AC 94 ms
15,296 KB
testcase_04 AC 93 ms
15,224 KB
testcase_05 AC 93 ms
15,232 KB
testcase_06 AC 93 ms
15,344 KB
testcase_07 AC 95 ms
15,372 KB
testcase_08 AC 93 ms
15,128 KB
testcase_09 AC 93 ms
15,232 KB
testcase_10 AC 93 ms
15,220 KB
testcase_11 AC 93 ms
15,360 KB
testcase_12 AC 95 ms
15,228 KB
testcase_13 AC 95 ms
15,228 KB
testcase_14 AC 94 ms
15,352 KB
testcase_15 AC 93 ms
15,228 KB
testcase_16 AC 94 ms
15,256 KB
testcase_17 AC 96 ms
15,372 KB
testcase_18 AC 96 ms
15,216 KB
testcase_19 AC 94 ms
15,228 KB
testcase_20 AC 93 ms
15,456 KB
testcase_21 AC 94 ms
15,116 KB
testcase_22 AC 94 ms
15,300 KB
testcase_23 AC 94 ms
15,224 KB
testcase_24 AC 94 ms
15,352 KB
testcase_25 AC 93 ms
15,324 KB
testcase_26 AC 95 ms
15,396 KB
testcase_27 AC 94 ms
15,372 KB
testcase_28 AC 96 ms
15,108 KB
testcase_29 AC 94 ms
15,128 KB
testcase_30 AC 94 ms
15,216 KB
testcase_31 AC 94 ms
15,308 KB
testcase_32 AC 93 ms
15,224 KB
testcase_33 AC 94 ms
15,208 KB
testcase_34 AC 95 ms
15,340 KB
testcase_35 AC 97 ms
15,260 KB
testcase_36 AC 95 ms
15,352 KB
testcase_37 AC 94 ms
15,112 KB
testcase_38 AC 100 ms
15,304 KB
testcase_39 AC 94 ms
15,300 KB
testcase_40 AC 95 ms
15,436 KB
testcase_41 AC 93 ms
15,372 KB
testcase_42 AC 93 ms
15,380 KB
testcase_43 AC 94 ms
15,224 KB
testcase_44 AC 93 ms
15,132 KB
testcase_45 AC 93 ms
15,344 KB
testcase_46 AC 94 ms
15,224 KB
testcase_47 AC 95 ms
15,312 KB
testcase_48 AC 94 ms
15,220 KB
testcase_49 AC 95 ms
15,392 KB
testcase_50 AC 96 ms
15,128 KB
testcase_51 AC 95 ms
15,300 KB
testcase_52 AC 94 ms
15,388 KB
testcase_53 AC 94 ms
15,348 KB
testcase_54 AC 94 ms
15,204 KB
testcase_55 AC 93 ms
15,128 KB
testcase_56 AC 93 ms
15,488 KB
testcase_57 AC 94 ms
15,336 KB
testcase_58 AC 95 ms
15,216 KB
testcase_59 AC 94 ms
15,384 KB
testcase_60 AC 96 ms
15,220 KB
testcase_61 AC 95 ms
15,456 KB
testcase_62 AC 94 ms
15,312 KB
testcase_63 AC 98 ms
15,396 KB
testcase_64 AC 96 ms
15,216 KB
testcase_65 AC 95 ms
15,208 KB
testcase_66 AC 94 ms
15,388 KB
testcase_67 AC 95 ms
15,140 KB
testcase_68 AC 93 ms
15,252 KB
testcase_69 AC 94 ms
15,212 KB
testcase_70 AC 95 ms
15,224 KB
testcase_71 AC 94 ms
15,112 KB
testcase_72 AC 96 ms
15,336 KB
testcase_73 AC 93 ms
15,356 KB
testcase_74 AC 94 ms
15,356 KB
testcase_75 AC 96 ms
15,112 KB
testcase_76 AC 93 ms
15,348 KB
testcase_77 AC 92 ms
15,236 KB
testcase_78 AC 92 ms
15,360 KB
testcase_79 AC 94 ms
15,376 KB
testcase_80 AC 93 ms
15,368 KB
testcase_81 AC 94 ms
15,392 KB
testcase_82 AC 94 ms
15,268 KB
testcase_83 AC 94 ms
15,260 KB
testcase_84 AC 94 ms
15,576 KB
testcase_85 AC 95 ms
15,120 KB
testcase_86 AC 95 ms
15,500 KB
testcase_87 AC 96 ms
15,560 KB
testcase_88 AC 93 ms
15,224 KB
testcase_89 AC 94 ms
15,224 KB
testcase_90 AC 96 ms
15,508 KB
testcase_91 AC 96 ms
15,208 KB
testcase_92 AC 96 ms
15,376 KB
testcase_93 AC 99 ms
15,344 KB
testcase_94 AC 94 ms
15,348 KB
testcase_95 AC 93 ms
15,320 KB
testcase_96 AC 94 ms
15,252 KB
testcase_97 AC 95 ms
15,244 KB
testcase_98 AC 93 ms
15,268 KB
testcase_99 AC 93 ms
15,204 KB
testcase_100 AC 93 ms
15,132 KB
testcase_101 AC 99 ms
15,452 KB
testcase_102 AC 93 ms
15,272 KB
testcase_103 AC 94 ms
15,212 KB
testcase_104 AC 94 ms
15,248 KB
testcase_105 AC 94 ms
15,128 KB
testcase_106 AC 95 ms
15,224 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Syntax OK

ソースコード

diff # 

#!/usr/bin/ruby
require 'prime'

# https://github.com/ChapterHouse/prime_miller_rabin
# Copyright (c) 2013 Frank Hall
# MIT License

class Prime::MillerRabin < Prime::PseudoPrimeGenerator
  def succ()
    self.last_prime = next_prime(last_prime || 1)
  end

  def rewind()
    self.last_prime = nil
  end

  def prime?(x)
    miller_rabin(x)
  end

  private

  attr_accessor :last_prime

  def likely_prime?(a, n)
    d = n - 1
    s = 0
    while d % 2 == 0 do
      d >>= 1
      s += 1
    end

    b = 1
    while d > 0
      u = d % 2
      t = d / 2
      b = (b * a) % n if u == 1
      a = a**2 % n
      d = t
    end

    if b == 1
      true
    else
      s.times do |i|
        return true if b == n - 1
        b = (b * b) % n
      end
      (b == n - 1)
    end
  end

  def miller_rabin(n)
    if n.abs < 2
      false
    else
      likely_prime = true
      # 26 Yields a probability of prime at 99.99999999999998% so lets kick it up a notch.
      27.times do |i|
        begin
          a = rand(n)
        end while a == 0
        likely_prime = likely_prime?(a, n)
        break unless likely_prime
      end
      likely_prime
    end
  end

  def next_prime(x)
    if x < 2
      2
    elsif x < 3
      3
    elsif x < 5
      5
    else
      x += (x.even? ? 1 : (x % 10 == 3 ? 4 : 2 ))
      x += (x % 10 == 3 ? 4 : 2 ) until x.prime?
      x
    end
  end
end

N=gets.to_i
if N>25
	p N%8==1&&Prime::MillerRabin.new.prime?(N-8) ? 14 : 8
else
	h={4=>3,6=>5,12=>11,14=>13,18=>8,20=>19,21=>19,24=>23,25=>23};h.default=7
	p h[N]
end
0