結果
問題 | 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 randomimport queuedef is_prime(n):if n==1: return Falseif n==2: return Trueif n%2==0: return Falseq = (n-1)>>1while 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:continuej = 0while a_pow != n-1 and (1<<(j+1))*q != n-1:a_pow *= a_powa_pow %= nj+=1if (1<<(j+1))*q==n-1 and a_pow != n-1:return Falsereturn Truemagic = [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)