// author: hotman78 // date: 2023/09/29-22:23:39 // --- begin raw code ----------------- // #include"cpplib/util/template.hpp" // #include // using namespace atcoder; // int op(int a, int b) { return min(a, b); } // int e() { return 1<<30; } // int main(){ // lint n; // cin>>n; // vec a(n); // rep(i,n)cin>>a[i]; // vector v(n,1<<30); // segtree seg(v); // seg.set(n-1,-a[n-1]); // lint ans=SUM(a); // rrep(i,n-1){ // int e=a[i]; // int val=1<<30; // rep(j,e,1000001,e){ // lint l=j; // lint r=j+e; // auto p=lower_bound(all(a),l)-a.begin(); // chmax(p,i+1); // auto q=lower_bound(all(a),r)-a.begin(); // auto res=seg.prod(p,q); // if(p using namespace std; #line 1 "cpplib/util/ioutil.hpp" // template // std::ostream& output(std::ostream& out,const Head& head,const Args&... args){ // out>>head; // return output(head,args...); // } // template // std::ostream& output(std::ostream& out,const Head& head){ // out>>head; // return out; // } template std::ostream &operator<<(std::ostream &out, std::pair v) { out << "(" << v.first << "," << v.second << ")"; return out; } // template // ostream& operator<<(ostream& out,std::tuplev){ // std::apply(output,v); // return out; // } #line 11 "cpplib/util/template.hpp" struct __INIT__ { __INIT__() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(15); } } __INIT__; typedef long long lint; constexpr long long INF = 1LL << 60; constexpr int IINF = 1 << 30; constexpr double EPS = 1e-10; #ifndef REACTIVE #define endl '\n'; #endif typedef vector vec; typedef vector> mat; typedef vector>> mat3; typedef vector svec; typedef vector> smat; template using V = vector; template using VV = V>; template inline void output(T t) { bool f = 0; for (auto i : t) { cout << (f ? " " : "") << i; f = 1; } cout << endl; } template inline void output2(T t) { for (auto i : t) output(i); } template inline void debug(T t) { bool f = 0; for (auto i : t) { cerr << (f ? " " : "") << i; f = 1; } cerr << endl; } template inline void debug2(T t) { for (auto i : t) debug(i); } #define loop(n) for (long long _ = 0; _ < (long long)(n); ++_) #define _overload4(_1, _2, _3, _4, name, ...) name #define __rep(i, a) repi(i, 0, a, 1) #define _rep(i, a, b) repi(i, a, b, 1) #define repi(i, a, b, c) \ for (long long i = (long long)(a); i < (long long)(b); i += c) #define rep(...) _overload4(__VA_ARGS__, repi, _rep, __rep)(__VA_ARGS__) #define _overload3_rev(_1, _2, _3, name, ...) name #define _rep_rev(i, a) repi_rev(i, 0, a) #define repi_rev(i, a, b) \ for (long long i = (long long)(b)-1; i >= (long long)(a); --i) #define rrep(...) _overload3_rev(__VA_ARGS__, repi_rev, _rep_rev)(__VA_ARGS__) #define all(n) begin(n), end(n) template bool chmin(T &s, const E &t) { bool res = s > t; s = min(s, t); return res; } template bool chmax(T &s, const E &t) { bool res = s < t; s = max(s, t); return res; } const vector dx = {1, 0, -1, 0, 1, 1, -1, -1}; const vector dy = {0, 1, 0, -1, 1, -1, 1, -1}; #define SUM(v) accumulate(all(v), 0LL) #if __cplusplus >= 201703L template auto make_vector(T x, int arg, Args... args) { if constexpr (sizeof...(args) == 0) return vector(arg, x); else return vector(arg, make_vector(x, args...)); } #endif #define extrep(v, ...) for (auto v : __MAKE_MAT__({__VA_ARGS__})) #define bit(n, a) ((n >> a) & 1) vector> __MAKE_MAT__(vector v) { if (v.empty()) return vector>(1, vector()); long long n = v.back(); v.pop_back(); vector> ret; vector> tmp = __MAKE_MAT__(v); for (auto e : tmp) for (long long i = 0; i < n; ++i) { ret.push_back(e); ret.back().push_back(i); } return ret; } using graph = vector>; template using graph_w = vector>>; #if __cplusplus >= 201703L constexpr inline long long powll(long long a, long long b) { long long res = 1; while (b--) res *= a; return res; } #endif template pair &operator+=(pair &s, const pair &t) { s.first += t.first; s.second += t.second; return s; } template pair &operator-=(pair &s, const pair &t) { s.first -= t.first; s.second -= t.second; return s; } template pair operator+(const pair &s, const pair &t) { auto res = s; return res += t; } template pair operator-(const pair &s, const pair &t) { auto res = s; return res -= t; } #define BEGIN_STACK_EXTEND(size) \ void *stack_extend_memory_ = malloc(size); \ void *stack_extend_origin_memory_; \ char *stack_extend_dummy_memory_ = (char *)alloca( \ (1 + (int)(((long long)stack_extend_memory_) & 127)) * 16); \ *stack_extend_dummy_memory_ = 0; \ asm volatile("mov %%rsp, %%rbx\nmov %%rax, %%rsp" \ : "=b"(stack_extend_origin_memory_) \ : "a"((char *)stack_extend_memory_ + (size)-1024)); #define END_STACK_EXTEND \ asm volatile("mov %%rax, %%rsp" ::"a"(stack_extend_origin_memory_)); \ free(stack_extend_memory_); int floor_pow(int n) { return n ? 31 - __builtin_clz(n) : 0; } #line 2 "main.cpp" #include #include #include #include #ifdef _MSC_VER #include #endif #if __cplusplus >= 202002L #include #endif namespace atcoder { namespace internal { #if __cplusplus >= 202002L using std::bit_ceil; #else unsigned int bit_ceil(unsigned int n) { unsigned int x = 1; while (x < (unsigned int)(n)) x *= 2; return x; } #endif int countr_zero(unsigned int n) { #ifdef _MSC_VER unsigned long index; _BitScanForward(&index, n); return index; #else return __builtin_ctz(n); #endif } constexpr int countr_zero_constexpr(unsigned int n) { int x = 0; while (!(n & (1 << x))) x++; return x; } } // namespace internal } // namespace atcoder namespace atcoder { #if __cplusplus >= 201703L template struct segtree { static_assert(std::is_convertible_v>, "op must work as S(S, S)"); static_assert(std::is_convertible_v>, "e must work as S()"); #else template struct segtree { #endif public: segtree() : segtree(0) {} explicit segtree(int n) : segtree(std::vector(n, e())) {} explicit segtree(const std::vector &v) : _n(int(v.size())) { size = (int)internal::bit_ceil((unsigned int)(_n)); log = internal::countr_zero((unsigned int)size); d = std::vector(2 * size, e()); for (int i = 0; i < _n; i++) d[size + i] = v[i]; for (int i = size - 1; i >= 1; i--) { update(i); } } void set(int p, S x) { assert(0 <= p && p < _n); p += size; d[p] = x; for (int i = 1; i <= log; i++) update(p >> i); } S get(int p) const { assert(0 <= p && p < _n); return d[p + size]; } S prod(int l, int r) const { assert(0 <= l && l <= r && r <= _n); S sml = e(), smr = e(); l += size; r += size; while (l < r) { if (l & 1) sml = op(sml, d[l++]); if (r & 1) smr = op(d[--r], smr); l >>= 1; r >>= 1; } return op(sml, smr); } S all_prod() const { return d[1]; } template int max_right(int l) const { return max_right(l, [](S x) { return f(x); }); } template int max_right(int l, F f) const { assert(0 <= l && l <= _n); assert(f(e())); if (l == _n) return _n; l += size; S sm = e(); do { while (l % 2 == 0) l >>= 1; if (!f(op(sm, d[l]))) { while (l < size) { l = (2 * l); if (f(op(sm, d[l]))) { sm = op(sm, d[l]); l++; } } return l - size; } sm = op(sm, d[l]); l++; } while ((l & -l) != l); return _n; } template int min_left(int r) const { return min_left(r, [](S x) { return f(x); }); } template int min_left(int r, F f) const { assert(0 <= r && r <= _n); assert(f(e())); if (r == 0) return 0; r += size; S sm = e(); do { r--; while (r > 1 && (r % 2)) r >>= 1; if (!f(op(d[r], sm))) { while (r < size) { r = (2 * r + 1); if (f(op(d[r], sm))) { sm = op(d[r], sm); r--; } } return r + 1 - size; } sm = op(d[r], sm); } while ((r & -r) != r); return 0; } private: int _n, size, log; std::vector d; void update(int k) { d[k] = op(d[2 * k], d[2 * k + 1]); } }; } // namespace atcoder using namespace atcoder; int op(int a, int b) { return min(a, b); } int e() { return 1 << 30; } int main() { lint n; cin >> n; vec a(n); rep(i, n) cin >> a[i]; vector v(n, 1 << 30); segtree seg(v); seg.set(n - 1, -a[n - 1]); lint ans = SUM(a); rrep(i, n - 1) { int e = a[i]; int val = 1 << 30; rep(j, e, 1000001, e) { lint l = j; lint r = j + e; auto p = lower_bound(all(a), l) - a.begin(); chmax(p, i + 1); auto q = lower_bound(all(a), r) - a.begin(); auto res = seg.prod(p, q); if (p < q) { // cerr<