結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー simansiman
提出日時 2020-01-14 18:38:00
言語 Ruby
(3.3.0)
結果
AC  
実行時間 4,956 ms / 9,973 ms
コード長 852 bytes
コンパイル時間 36 ms
コンパイル使用メモリ 7,296 KB
実行使用メモリ 19,200 KB
最終ジャッジ日時 2024-11-16 23:21:23
合計ジャッジ時間 11,793 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 115 ms
18,688 KB
testcase_01 AC 115 ms
18,560 KB
testcase_02 AC 116 ms
18,688 KB
testcase_03 AC 111 ms
18,944 KB
testcase_04 AC 2,430 ms
19,072 KB
testcase_05 AC 2,283 ms
19,200 KB
testcase_06 AC 320 ms
19,072 KB
testcase_07 AC 321 ms
19,200 KB
testcase_08 AC 331 ms
19,072 KB
testcase_09 AC 4,956 ms
19,200 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Syntax OK

ソースコード

diff # 

require 'openssl'

def modular_exponentiation(base, n, mod)
  res = 1

  while n > 0
    if n[0] == 1
      res = (res * base) % mod
    end

    base = (base ** 2) % mod
    n >>= 1
  end

  res
end

def witness(a, n)
  v = n - 1
  t = 0

  while v % 2 == 0
    t += 1
    v /= 2
  end

  u = n / (2 ** t)

  x = modular_exponentiation(a, u, n)

  t.times do
    bx = x
    x = (x ** 2) % n

    if x == 1 && bx != 1 && bx != n - 1
      return true
    end
  end

  return true if x != 1

  false
end

def miller_rabin(n, s = 10)
  return false if n <= 1
  return true if n == 2
  return false if n.even?

  s.times do
    a = rand(n - 1) + 1

    return false if witness(a, n)
  end

  true
end

N = gets.to_i
N.times do
  x = gets.to_i
  bn = OpenSSL::BN.new(x)

  if bn.prime?
    puts [x, 1].join(' ')
  else
    puts [x, 0].join(' ')
  end
end
0