結果
問題 | No.308 素数は通れません |
ユーザー |
|
提出日時 | 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 |
ソースコード
import randomdef MillerRabin(q, k = 1000):'''素数判定の ミラー・ラビン テスト 以下のコードをもとに作成。http://d.hatena.ne.jp/pashango_p/20090704/1246692091'''if q == 2:return Trueif q == 1 or not (q & 1):return Falsed = (q - 1) >> 1while not (d & 1):d >>= 1for i in range(k):a = random.randint(1, q - 1)t = dy = pow(a, t, q)while t != q - 1 and y != 1 and y != q - 1:y = pow(y, 2, q)t <<= 1if y != q - 1 and not (t & 1):return Falsereturn Truedef 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 14if N < 57:return 8if N % 8 == 1:M = N - 8if MillerRabin(M):return 14return 8N = int(input())print(solve(N))