結果

問題 No.3030 ミラー・ラビン素数判定法のテスト
ユーザー xpicxpic
提出日時 2024-11-29 17:01:00
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 323 ms / 9,973 ms
コード長 3,366 bytes
コンパイル時間 218 ms
コンパイル使用メモリ 82,372 KB
実行使用メモリ 78,152 KB
最終ジャッジ日時 2024-11-29 17:01:02
合計ジャッジ時間 2,477 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
55,060 KB
testcase_01 AC 40 ms
55,768 KB
testcase_02 AC 39 ms
56,300 KB
testcase_03 AC 39 ms
55,084 KB
testcase_04 AC 219 ms
78,068 KB
testcase_05 AC 223 ms
78,084 KB
testcase_06 AC 159 ms
78,112 KB
testcase_07 AC 153 ms
78,152 KB
testcase_08 AC 156 ms
78,048 KB
testcase_09 AC 323 ms
77,752 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import os
import sys
from io import BytesIO, IOBase
from random import randint

BUFSIZE = 1 << 15


class FastIO(IOBase):
    newlines = 0

    def __init__(self, file):
        self._fd = file.fileno()
        self.buffer = BytesIO()
        self.writable = "x" in file.mode or "r" not in file.mode
        self.write = self.buffer.write if self.writable else None

    def read(self):
        while True:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            if not b:
                break
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines = 0
        return self.buffer.read()

    def readline(self):
        while self.newlines == 0:
            b = os.read(self._fd, max(os.fstat(self._fd).st_size, BUFSIZE))
            self.newlines = b.count(b"\n") + (not b)
            ptr = self.buffer.tell()
            self.buffer.seek(0, 2), self.buffer.write(b), self.buffer.seek(ptr)
        self.newlines -= 1
        return self.buffer.readline()

    def flush(self):
        if self.writable:
            os.write(self._fd, self.buffer.getvalue())
            self.buffer.truncate(0), self.buffer.seek(0)


class IOWrapper(IOBase):
    def __init__(self, file):
        self.buffer = FastIO(file)
        self.flush = self.buffer.flush
        self.writable = self.buffer.writable
        self.write = lambda s: self.buffer.write(s.encode("ascii"))
        self.read = lambda: self.buffer.read().decode("ascii")
        self.readline = lambda: self.buffer.readline().decode("ascii")


def output(*args, **kwargs):
    sep, file = kwargs.pop("sep", " "), kwargs.pop("file", sys.stdout)
    at_start = True
    for x in args:
        if not at_start:
            file.write(sep)
        file.write(str(x))
        at_start = False
    file.write(kwargs.pop("end", "\n"))
    if kwargs.pop("flush", False):
        file.flush()

sys.stdin, sys.stdout = IOWrapper(sys.stdin), IOWrapper(sys.stdout)
input = lambda: sys.stdin.readline().rstrip("\r\n")

def miller_rabin(n, check):
    d, s = n - 1, 0
    while d % 2 == 0:
        d >>= 1
        s += 1
    for a in check:
        if n <= a:
            return True
        a = pow(a, d, n)
        if a == 1:
            continue
        r = 1
        while a != n - 1:
            if r == s:
                return False
            a = a * a % n
            r += 1
    return True


def is_prime32(n):
    return miller_rabin(n, [2, 7, 61])


def is_prime64(n):
    return miller_rabin(n, [2, 325, 9375, 28178, 450775, 9780504, 1795265022])


def is_prime(n):
    if n <= 1:
        return False
    if n <= 3:
        return True
    if n & 1 == 0:
        return False
    if n < 4759123141:
        return is_prime32(n)
    if n < 18446744073709551615:
        return is_prime64(n)
    d = (n - 1) >> 1
    while d & 1 == 0:
        d >>= 1
    for _ in range(100):
        a = randint(1, n - 1)
        t = d
        y = pow(a, t, n)
        while t != n - 1 and y != 1 and y != n - 1:
            y = (y * y) % n
            t <<= 1
        if y != n - 1 and t & 1 == 0:
            return False
    return True


def main():
    for _ in range(int(input())):
        a = int(input())
        output(a, end=' ')
        output(int(is_prime(a)))


if __name__ == '__main__':
    main()
0