結果
問題 | No.6 使いものにならないハッシュ |
ユーザー |
|
提出日時 | 2024-02-28 00:35:59 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 101 ms / 5,000 ms |
コード長 | 2,101 bytes |
コンパイル時間 | 291 ms |
コンパイル使用メモリ | 82,488 KB |
実行使用メモリ | 79,488 KB |
最終ジャッジ日時 | 2024-09-29 12:10:53 |
合計ジャッジ時間 | 4,034 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 32 |
ソースコード
import sysimport mathimport bisectfrom heapq import heapify, heappop, heappushfrom collections import deque, defaultdict, Counterfrom functools import lru_cachefrom itertools import accumulate, combinations, permutations, productsys.setrecursionlimit(1000000)MOD = 10 ** 9 + 7MOD99 = 998244353input = lambda: sys.stdin.readline().strip()NI = lambda: int(input())NMI = lambda: map(int, input().split())NLI = lambda: list(NMI())SI = lambda: input()SMI = lambda: input().split()SLI = lambda: list(SMI())EI = lambda m: [NLI() for _ in range(m)]# 高速エラストテネス sieve[n]はnの最小の素因数def make_prime_table(n):sieve = list(range(n + 1))sieve[0] = -1sieve[1] = -1for i in range(4, n + 1, 2):sieve[i] = 2for i in range(3, int(n ** 0.5) + 1, 2):if sieve[i] != i:continuefor j in range(i * i, n + 1, i * 2):if sieve[j] == j:sieve[j] = ireturn sieve# Nの素因数分解を辞書で返す(単体)def prime_fact(n):root = int(n**0.5) + 1prime_dict = {}for i in range(2, root):cnt = 0while n % i == 0:cnt += 1n = n // iif cnt:prime_dict[i] = cntif n != 1:prime_dict[n] = 1return prime_dict# 約数列挙(単体)def divisors(x):res = set()for i in range(1, int(x**0.5) + 2):if x % i == 0:res.add(i)res.add(x//i)return resdef h(x):while x >= 10:x = sum(int(s) for s in str(x))return xdef main():K = NI()N = NI()prime_table = make_prime_table(N+1)primes = [p for i, p in enumerate(prime_table) if i == p and K <= p <= N]H = [h(p) for p in primes]ans = 0for g in range(9, 0, -1):ok = Falsefor i in range(len(H)+1-g):S = set(H[i:i+g])if len(S) == g:ans = max(ans, primes[i])ok = Trueif ok:breakprint(ans)if __name__ == "__main__":main()