#include using namespace std; #define FOR(i,m,n) for(int i=(m);i<(n);++i) #define REP(i,n) FOR(i,0,n) using ll = long long; constexpr int INF = 0x3f3f3f3f; constexpr long long LINF = 0x3f3f3f3f3f3f3f3fLL; constexpr double EPS = 1e-8; constexpr int MOD = 998244353; // constexpr int MOD = 1000000007; constexpr int DY4[]{1, 0, -1, 0}, DX4[]{0, -1, 0, 1}; constexpr int DY8[]{1, 1, 0, -1, -1, -1, 0, 1}; constexpr int DX8[]{0, -1, -1, -1, 0, 1, 1, 1}; template inline bool chmax(T& a, U b) { return a < b ? (a = b, true) : false; } template inline bool chmin(T& a, U b) { return a > b ? (a = b, true) : false; } struct IOSetup { IOSetup() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << fixed << setprecision(20); } } iosetup; template requires requires { typename T::Monoid; {T::id()} -> std::same_as; {T::merge(std::declval(), std::declval())} -> std::same_as; } struct SegmentTree { using Monoid = typename T::Monoid; explicit SegmentTree(const int n) : SegmentTree(std::vector(n, T::id())) {} explicit SegmentTree(const std::vector& a) : n(a.size()), p2(std::bit_ceil(a.size())) { dat.assign(p2 << 1, T::id()); std::copy(a.begin(), a.end(), dat.begin() + p2); for (int i = p2 - 1; i > 0; --i) { dat[i] = T::merge(dat[i << 1], dat[(i << 1) + 1]); } } void set(int idx, const Monoid val) { idx += p2; dat[idx] = val; while (idx >>= 1) dat[idx] = T::merge(dat[idx << 1], dat[(idx << 1) + 1]); } Monoid get(int left, int right) const { Monoid res_l = T::id(), res_r = T::id(); for (left += p2, right += p2; left < right; left >>= 1, right >>= 1) { if (left & 1) res_l = T::merge(res_l, dat[left++]); if (right & 1) res_r = T::merge(dat[--right], res_r); } return T::merge(res_l, res_r); } Monoid operator[](const int idx) const { return dat[idx + p2]; } template int find_right(int left, const G g) const { if (left >= n) [[unlikely]] return n; Monoid val = T::id(); left += p2; do { while (!(left & 1)) left >>= 1; Monoid nxt = T::merge(val, dat[left]); if (!g(nxt)) { while (left < p2) { left <<= 1; nxt = T::merge(val, dat[left]); if (g(nxt)) { val = nxt; ++left; } } return left - p2; } val = nxt; ++left; } while (!std::has_single_bit(static_cast(left))); return n; } template int find_left(int right, const G g) const { if (right <= 0) [[unlikely]] return -1; Monoid val = T::id(); right += p2; do { --right; while (right > 1 && (right & 1)) right >>= 1; Monoid nxt = T::merge(dat[right], val); if (!g(nxt)) { while (right < p2) { right = (right << 1) + 1; nxt = T::merge(dat[right], val); if (g(nxt)) { val = nxt; --right; } } return right - p2; } val = nxt; } while (!std::has_single_bit(static_cast(right))); return -1; } private: const int n, p2; std::vector dat; }; namespace monoid { template struct RangeMinimumQuery { using Monoid = T; static constexpr Monoid id() { return std::numeric_limits::max(); } static Monoid merge(const Monoid& a, const Monoid& b) { return std::min(a, b); } }; template struct RangeMaximumQuery { using Monoid = T; static constexpr Monoid id() { return std::numeric_limits::lowest(); } static Monoid merge(const Monoid& a, const Monoid& b) { return std::max(a, b); } }; template struct RangeSumQuery { using Monoid = T; static constexpr Monoid id() { return 0; } static Monoid merge(const Monoid& a, const Monoid& b) { return a + b; } }; } // namespace monoid int main() { int n; cin >> n; vector a(n); REP(i, n) cin >> a[i]; SegmentTree> dp(a.back() + 1); dp.set(a.back(), 0); for (int i = n - 2; i >= 0; --i) { ll dp_ai = LINF; for (int d = 1; ; ++d) { const ll tmp = dp.get(a[i] * d, min(a[i] * (d + 1), a.back() + 1)); if (tmp != monoid::RangeMinimumQuery::id()) { chmin(dp_ai, tmp + 1LL * d * a[i]); } if (a[i] * (d + 1) > a.back()) break; } dp.set(a[i], dp_ai); } // REP(i, a.back() + 1) cout << (dp[i] == monoid::RangeMinimumQuery::id() ? -1 : dp[i]) << " \n"[i == a.back()]; cout << reduce(a.begin(), a.end(), 0LL) - dp[a.front()] << '\n'; return 0; }