#include using namespace std; const long long INF = (3LL << 59); class segtree { private: int sz; vector val; public: segtree() : sz(0), val(vector()) {} segtree(int n) : sz(n >= 3 ? 1 << (32 - __builtin_clz(n)) : n), val(vector(sz * 2, -INF)) {} void update(int pos, long long x) { pos += sz; val[pos] = x; while (pos > 1) { pos /= 2; val[pos] = max(val[pos * 2], val[pos * 2 + 1]); } } long long rangemax(int l, int r) const { l += sz; r += sz; long long ans = -INF; while (l != r) { if (l & 1) ans = max(ans, val[l++]); if (r & 1) ans = max(ans, val[--r]); l /= 2; r /= 2; } return ans; } }; int main() { // step #1. input cin.tie(0); ios::sync_with_stdio(false); int N; cin >> N; vector A(N); for (int i = 0; i < N; i++) { cin >> A[i]; } // step #2. calculate sum long long suma = 0; for (int i = 0; i < N; i++) { suma += A[i]; } // step #3. dynamic programming vector dp(N, -INF); segtree seg(N); dp[N - 1] = 0; seg.update(N - 1, dp[N - 1]); for (int i = N - 2; i >= 0; i--) { int posl = i + 1; for (int j = A[i]; j <= A[N - 1]; j += A[i]) { int posr = lower_bound(A.begin(), A.end(), j + A[i]) - A.begin(); long long res = seg.rangemax(posl, posr); dp[i] = max(dp[i], res - j); posl = posr; } seg.update(i, dp[i]); } // step #4. output long long ans = suma + dp[0]; cout << ans << endl; return 0; }