結果
問題 | No.12 限定された素数 |
ユーザー | rpy3cpp |
提出日時 | 2015-08-07 02:18:52 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 1,035 ms / 5,000 ms |
コード長 | 1,583 bytes |
コンパイル時間 | 197 ms |
コンパイル使用メモリ | 12,928 KB |
実行使用メモリ | 49,276 KB |
最終ジャッジ日時 | 2024-05-03 09:33:54 |
合計ジャッジ時間 | 14,920 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 385 ms
48,844 KB |
testcase_01 | AC | 437 ms
48,980 KB |
testcase_02 | AC | 400 ms
48,868 KB |
testcase_03 | AC | 1,035 ms
49,192 KB |
testcase_04 | AC | 421 ms
49,200 KB |
testcase_05 | AC | 491 ms
48,828 KB |
testcase_06 | AC | 649 ms
49,092 KB |
testcase_07 | AC | 697 ms
49,196 KB |
testcase_08 | AC | 460 ms
49,176 KB |
testcase_09 | AC | 459 ms
49,168 KB |
testcase_10 | AC | 444 ms
49,064 KB |
testcase_11 | AC | 845 ms
49,120 KB |
testcase_12 | AC | 688 ms
49,200 KB |
testcase_13 | AC | 583 ms
48,960 KB |
testcase_14 | AC | 458 ms
49,108 KB |
testcase_15 | AC | 583 ms
48,892 KB |
testcase_16 | AC | 733 ms
49,276 KB |
testcase_17 | AC | 398 ms
48,848 KB |
testcase_18 | AC | 379 ms
48,852 KB |
testcase_19 | AC | 383 ms
49,064 KB |
testcase_20 | AC | 401 ms
49,112 KB |
testcase_21 | AC | 437 ms
48,860 KB |
testcase_22 | AC | 397 ms
49,100 KB |
testcase_23 | AC | 394 ms
49,116 KB |
testcase_24 | AC | 392 ms
49,072 KB |
testcase_25 | AC | 483 ms
49,136 KB |
ソースコード
def primes2(limit): ''' returns a list of prime numbers upto limit. source: Rossetta code: Sieve of Eratosthenes http://rosettacode.org/wiki/Sieve_of_Eratosthenes#Odds-only_version_of_the_array_sieve_above ''' if limit < 2: return [] if limit < 3: return [2] lmtbf = (limit - 3) // 2 buf = [True] * (lmtbf + 1) for i in range((int(limit ** 0.5) - 3) // 2 + 1): if buf[i]: p = i + i + 3 s = p * (i + 1) + i buf[s::p] = [False] * ((lmtbf - s) // p + 1) return [2] + [i + i + 3 for i, v in enumerate(buf) if v] def read_data(): N = int(input()) A = list(map(int, input().split())) return N, A def solve(N, A): limit = 5 * 10**6 primes = primes2(limit) NGs = [True] * 10 for a in A: NGs[a] = False OKset = set(A) prev_p = 0 maxspan = -1 currentset = set() for p in primes: if is_safe(p, NGs): update_set(currentset, p) continue elif currentset == OKset: span = p - prev_p - 2 if maxspan < span: maxspan = span currentset = set() prev_p = p p = limit if currentset == OKset: span = p - prev_p - 1 if maxspan < span: maxspan = span return maxspan def is_safe(p, NGs): while p: p, r = divmod(p, 10) if NGs[r]: return False return True def update_set(currentset, p): while p: p, r = divmod(p, 10) currentset.add(r) N, A = read_data() print(solve(N, A))