結果
問題 | No.2039 Copy and Avoid |
ユーザー |
![]() |
提出日時 | 2022-08-13 12:54:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 353 ms / 2,000 ms |
コード長 | 1,017 bytes |
コンパイル時間 | 231 ms |
コンパイル使用メモリ | 82,600 KB |
実行使用メモリ | 91,268 KB |
最終ジャッジ日時 | 2024-09-24 16:46:35 |
合計ジャッジ時間 | 4,834 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 31 |
ソースコード
import bisect import heapq inf = 2 * 10 ** 18 def gen_div(N): yield 1 for i in range(2, N+1): if i * i > N: break if N % i == 0: yield i if i * i != N: yield N // i yield N N, M, A, B = map(int, input().split()) C = sorted(map(int, input().split())) ans = inf divisors_of_N = sorted(gen_div(N)) K = len(divisors_of_N) dtoi = {d: i for i, d in enumerate(divisors_of_N)} dp = [[inf] * K for _ in range(K)] dp[0][0] = 0 for i in range(K): for j in range(i): dp[i][i] = min(dp[i][i], dp[i][j] + B) di = divisors_of_N[i] NG = [] for c in C: if c % di == 0: NG.append(c) for j in range(i+1, K): dj = divisors_of_N[j] if dj % di != 0: continue if bisect.bisect_right(NG, dj) - bisect.bisect_left(NG, di) > 0: break dp[j][i] = min(dp[j][i], dp[i][i] + (dj - di) // di * A) ans = min(dp[K-1]) if ans >= inf: ans = -1 print(ans)