結果
| 問題 |
No.8030 ミラー・ラビン素数判定法のテスト
|
| ユーザー |
nagiss
|
| 提出日時 | 2019-09-21 02:25:19 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
WA
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 821 bytes |
| コンパイル時間 | 377 ms |
| コンパイル使用メモリ | 12,672 KB |
| 実行使用メモリ | 11,392 KB |
| 最終ジャッジ日時 | 2024-11-18 17:04:56 |
| 合計ジャッジ時間 | 5,497 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 9 WA * 1 |
ソースコード
def miller_rabin(n):
# 確率的素数判定(ミラーラビン素数判定法)
# 素数なら確実に True を返す、合成数なら確率的に False を返す
# True が返ったなら恐らく素数で、False が返ったなら確実に合成数である
primes = [2, 325, 9375, 28178, 450775, 9780504, 1795265022]
if n == 2: return True
if n <= 1 or n & 1 == 0: return False
d = m1 = n-1
r = (d & -d).bit_length() - 1
d >>= r
for a in primes:
if a >= n: return True
x = pow(a, d, n)
if x == 1 or x == m1: continue
for _ in range(r-1):
x = x * x % n
if x == m1: break
else: return False
return True
N = int(input())
X = [int(input()) for _ in range(N)]
for x in X:
print(x, int(miller_rabin(x)))
nagiss