結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー 👑 KazunKazun
提出日時 2020-09-18 02:51:23
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 946 bytes
コンパイル時間 269 ms
コンパイル使用メモリ 82,176 KB
実行使用メモリ 79,104 KB
最終ジャッジ日時 2024-04-29 14:19:10
合計ジャッジ時間 2,370 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 45 ms
53,760 KB
testcase_01 AC 44 ms
53,888 KB
testcase_02 AC 46 ms
53,632 KB
testcase_03 AC 46 ms
54,016 KB
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 AC 199 ms
79,104 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

#Miller-Rabinの素数判定法
def Miller_Rabin_Primality_Test(N,Times=20):
    """Miller-Rabinによる整数Nの素数判定を行う.

    N:整数
    ※:Trueは正確にはProbably Trueである(Falseは確定False).
    """
    from random import randint as ri

    if N==2:
        return True

    if N==1 or N&1==0:
        return False

    x=N-1
    k=0
    while x&1==0:
        k+=1
        x>>=1
    q=(N-1)//(1<<k)

    for _ in range(Times):
        m=ri(2,N-1)
        if pow(m,q,N)==1:
            return True

        y=pow(m,q,N)

        for i in range(k):
            if (y+1)%N==0:
                return True
            y*=y
            y%=N

        return False
#================================================
N=int(input())
Y=[0]*N
for i in range(N):
    x=int(input())
    if Miller_Rabin_Primality_Test(x):
        Y[i]="{} {}".format(x,1)
    else:
        Y[i]="{} {}".format(x,0)

print("\n".join(map(str,Y)))
0