結果
| 問題 | 
                            No.2039 Copy and Avoid
                             | 
                    
| コンテスト | |
| ユーザー | 
                             tktk_snsn
                         | 
                    
| 提出日時 | 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)
            
            
            
        
            
tktk_snsn