結果
問題 |
No.308 素数は通れません
|
ユーザー |
|
提出日時 | 2015-12-01 14:48:28 |
言語 | PyPy2 (7.3.15) |
結果 |
AC
|
実行時間 | 86 ms / 1,000 ms |
コード長 | 806 bytes |
コンパイル時間 | 2,860 ms |
コンパイル使用メモリ | 76,928 KB |
実行使用メモリ | 75,648 KB |
最終ジャッジ日時 | 2024-11-27 18:16:38 |
合計ジャッジ時間 | 12,752 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 107 |
ソースコード
def suspect(a, e, k, n): a = pow(a, e, n) if a == 1: return True for _ in xrange(k): if a == n - 1: return True a = a * a % n return False def isprime(n): if n <= 1 or (n > 2 and n % 2 == 0): return False e = n - 1 k = 0 while e % 2 == 0: e /= 2 k += 1 bases = [2, 3, 5, 7, 11, 13, 17, 19, 23] for i in bases: if n > i and not suspect(i, e, k, n): return False return True ans = [1,1,1,3,1,5,1,7,7,7,1,11,1,13,7,7,1,8,1,19,19,7,1,23,23,8,8,8,1,8,1,8,8,8,8,8,1,8,8,8,1,8,1,8,8,8,1,8,14,8,8,8,1,8,8,8,8,8,1,8,1,8,8,8,8,8,1,8,8,8,1,8,1,8,8,8,8,8,1,8,14,8,1,8,8,8,8,8,1,8,8,8,8,8,8,8,1,8,8] n = input() if n < 100: print ans[n - 1] else: if n % 8 != 1: print 8 elif not isprime(n - 8): print 8 elif n % 14 != 1: print 14 elif not isprime(n - 14): print 14 else: print 20