結果

問題 No.308 素数は通れません
ユーザー iwashi31iwashi31
提出日時 2015-12-01 22:46:16
言語 Python2
(2.7.18)
結果
AC  
実行時間 15 ms / 1,000 ms
コード長 1,496 bytes
コンパイル時間 1,144 ms
コンパイル使用メモリ 7,040 KB
実行使用メモリ 7,168 KB
最終ジャッジ日時 2024-09-26 12:12:03
合計ジャッジ時間 4,867 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 107
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from random import seed, randrange

def modpow(base, power, mod):
    result = 1
    while power > 0:
        if power & 1 == 1:
            result = (result * base) % mod
        base = (base * base) % mod
        power >>= 1
    return result

def isPrime(n):
    if n % 2 == 0:
        return False
    
    def MR_core(d, a, s = 1):
        if s == 1:
            return MR_core(d >> 1, a, 2) == n
        elif s == 2:
            if d & 1:
                m = MR_core(d >> 1, a, 3)
                s = m * m % n * a % n
                return n if s == n - 1 or s == 1 else s
            else:
                m = MR_core(d >> 1, a, 2)
                if m == n:
                    return n
                else:
                    s = m * m % n
                    return n if s == n - 1 else s
        elif s == 3:
            if d == 0:
                return 1
            else:
                m = MR_core(d >> 1, a, 3)
                if d & 1:
                    return m * m % n * a % n
                else:
                    return m * m % n
    
    k = 3
    return all(MR_core(n - 1, randrange(2, n)) for i in range(k))

n = input()
umekomi = [ 1, 1, 2, 3, 3, 5, 3, 7, 7, 7, \
            5,11,11,13, 7, 7, 7, 8, 8,19, \
           19, 7, 7,23,23, 8, 8, 8, 7, 8]
if n <= 30:
    print umekomi[n-1]
    sys.exit()

m = 8
while True:
    if isPrime(m+1) == False:
        if n%m != 1 or isPrime(n-m) == False:
            print m
            sys.exit()
    m += 2
0