N = int(input()) Ms = [int(input()) for i in range(N)] dp = [float('inf')] * (1 << N) dp[0] = 0 discount = [0] * (1 << N) full = (1 << N) - 1 for m in range(N): msb = 1 << m for mask in range(1 << m): discount[msb + mask] = (discount[mask] + Ms[m]) % 1000 frontiers = [0] for m in range(N): new_frontiers = [] for f in frontiers: dpf = dp[f] dsc = discount[f] for idx in range(N): bit = 1 << idx if bit & f: continue nf = bit | f if dp[nf] == float('inf'): new_frontiers.append(nf) dp[nf] = min(dp[nf], dpf + max(0, Ms[idx] - dsc)) frontiers = new_frontiers print(dp[full])