結果
| 問題 |
No.8030 ミラー・ラビン素数判定法のテスト
|
| ユーザー |
nagiss
|
| 提出日時 | 2019-09-21 02:27:47 |
| 言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
| 結果 |
WA
(最新)
AC
(最初)
|
| 実行時間 | - |
| コード長 | 862 bytes |
| コンパイル時間 | 260 ms |
| コンパイル使用メモリ | 12,672 KB |
| 実行使用メモリ | 11,392 KB |
| 最終ジャッジ日時 | 2024-11-18 17:04:39 |
| 合計ジャッジ時間 | 5,546 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 9 WA * 1 |
ソースコード
def miller_rabin(n):
# 参考: http://tjkendev.github.io/procon-library/python/prime/probabilistic.html
# 確率的素数判定(ミラーラビン素数判定法)
# 素数なら確実に 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
d //= d & -d
for a in primes:
if a >= n: return True
t = d
y = pow(a, t, n)
while t!=m1 and y!=1 and y!=m1:
y = y * y % n
t <<= 1
if y!=m1 and t&1==0: 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