結果

問題 No.308 素数は通れません
ユーザー rpy3cpprpy3cpp
提出日時 2015-12-05 10:48:58
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 64 ms / 1,000 ms
コード長 998 bytes
コンパイル時間 222 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 11,264 KB
最終ジャッジ日時 2024-11-27 18:40:18
合計ジャッジ時間 6,185 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 107
権限があれば一括ダウンロードができます

ソースコード

diff # 

import random

def MillerRabin(q, k = 1000):
    '''
    素数判定の ミラー・ラビン テスト 以下のコードをもとに作成。
    http://d.hatena.ne.jp/pashango_p/20090704/1246692091
    '''
    if q == 2:
        return True
    if q == 1 or not (q & 1):
        return False
    d = (q - 1) >> 1
    while not (d & 1):
        d >>= 1
    for i in range(k):
        a = random.randint(1, q - 1)
        t = d
        y = pow(a, t, q)
        while t != q - 1 and y != 1 and y != q - 1: 
            y = pow(y, 2, q)
            t <<= 1
        if y != q - 1 and not (t & 1):
            return False
    return True

def solve(N):
    dic = {4:3, 6:5, 8:7, 9:7, 10:7, 12:11, 14:13, 15:7, 16:7, 18:8, 20:19, 21:19, 22:7, 24:23, 25:23, 26:8}
    if N in dic:
        return dic[N]
    if N == 49:
        return 14
    if N < 57:
        return 8
    if N % 8 == 1:
        M = N - 8
        if MillerRabin(M):
            return 14
    return 8

N = int(input())
print(solve(N))
0