結果
| 問題 |
No.8030 ミラー・ラビン素数判定法のテスト
|
| ユーザー |
Kiri8128
|
| 提出日時 | 2020-05-06 00:13:05 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 537 ms / 9,973 ms |
| コード長 | 715 bytes |
| コンパイル時間 | 200 ms |
| コンパイル使用メモリ | 82,048 KB |
| 実行使用メモリ | 77,312 KB |
| 最終ジャッジ日時 | 2024-11-16 23:24:48 |
| 合計ジャッジ時間 | 2,804 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 10 |
ソースコード
def isPrimeMR(n):
d = n - 1
d = d // (d & -d)
L = [2, 7, 61] if n < 1<<32 else [2, 3, 5, 7, 11, 13, 17] if n < 1<<48 else [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37]
if n in L: return 1
for a in L:
t = d
y = pow(a, t, n)
if y == 1: continue
while y != n - 1:
y = y * y % n
if y == 1 or t == n - 1: return 0
t <<= 1
return 1
P = [2, 3, 5, 7, 11, 13, 17, 19]
setP = set(P)
def isPrime(n):
if n in setP:
return 1
if n < 20:
return 0
for p in P:
if n % p == 0:
return 0
return isPrimeMR(n)
N = int(input())
for _ in range(N):
n = int(input())
print(n, isPrime(n))
Kiri8128