結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー JashinchanJashinchan
提出日時 2022-09-07 18:21:17
言語 Ruby
(3.3.0)
結果
AC  
実行時間 3,547 ms / 9,973 ms
コード長 1,048 bytes
コンパイル時間 39 ms
コンパイル使用メモリ 7,552 KB
実行使用メモリ 12,416 KB
最終ジャッジ日時 2024-05-02 10:40:14
合計ジャッジ時間 11,453 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 92 ms
12,160 KB
testcase_01 AC 93 ms
12,160 KB
testcase_02 AC 92 ms
12,032 KB
testcase_03 AC 94 ms
12,160 KB
testcase_04 AC 2,014 ms
12,288 KB
testcase_05 AC 1,951 ms
12,416 KB
testcase_06 AC 794 ms
12,288 KB
testcase_07 AC 780 ms
12,288 KB
testcase_08 AC 783 ms
12,288 KB
testcase_09 AC 3,547 ms
12,288 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
Syntax OK

ソースコード

diff # 

def pow_mod(a, k, n)
    ret = 1
    while k > 0
        if k.odd?
            ret = (ret * a) % n
        end
        a = (a ** 2) % n
        k >>= 1
    end
    ret
end

def miller_rabin(n, bases)
    d = n - 1
    s = 0
    while d.even?
        d = d >> 1
        s = s + 1
    end
    for a in bases do
        if n <= a
            return true
        end
        a = pow_mod(a, d, n)
        if a == 1 then
        else
            r = 1
            while a != n - 1 do
                if r == s then
                    return false
                end
                a = a * a % n
                r = r + 1
            end
        end
    end
    return true
end

def is_prime(n)
    return false if n < 2
    return true if n < 4
    return false if n.even?
    return miller_rabin(n, Array[2, 7, 61]) if n < 4759123141
    return miller_rabin(n, Array[2, 325, 9375, 28178, 450775, 9780504, 1795265022])
end


N = gets.to_i
N.times do
  x = gets.to_i
  if is_prime(x)
    puts [x, 1].join(' ')
  else
    puts [x, 0].join(' ')
  end
end
0