結果
問題 | No.308 素数は通れません |
ユーザー |
|
提出日時 | 2016-06-29 09:48:29 |
言語 | Python3 (3.13.1 + numpy 2.2.1 + scipy 1.14.1) |
結果 |
AC
|
実行時間 | 39 ms / 1,000 ms |
コード長 | 676 bytes |
コンパイル時間 | 137 ms |
コンパイル使用メモリ | 12,544 KB |
実行使用メモリ | 11,520 KB |
最終ジャッジ日時 | 2024-11-27 18:46:41 |
合計ジャッジ時間 | 6,040 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 107 |
ソースコード
import random import queue def is_prime(n): if n==1: return False if n==2: return True if n%2==0: return False q = (n-1)>>1 while q&1 == 0: q >>= 1 # print("n={}, q={}".format(n, q)) for i in range(30): a = random.randint(1, n-1) a_pow = pow(a, q, n) if a_pow == 1: continue j = 0 while a_pow != n-1 and (1<<(j+1))*q != n-1: a_pow *= a_pow a_pow %= n j+=1 if (1<<(j+1))*q==n-1 and a_pow != n-1: return False return True magic = [0, 1, 0, 0, 3, 0, 5, 0, 7, 7, 7, 0, 11, 0, 13, 7, 7, 0, 8, 0, 19, 19, 7, 0, 23, 23, 8, 8, 8, 0] N = int(input()) if N < 30: print(magic[N]) else: if N%8 == 1 and is_prime(N-8): print(14) else: print(8)