#include #define rep(i, n) for(int(i) = 0; (i) < (n); (i)++) #define FOR(i, m, n) for(int(i) = (m); (i) < (n); (i)++) #define ALL(v) (v).begin(), (v).end() #define LLA(v) (v).rbegin(), (v).rend() #define PB push_back #define MP(a, b) make_pair((a), (b)) using namespace std; template inline vector make_vec(size_t a, T val) { return vector(a, val); } template inline auto make_vec(size_t a, Ts... ts) { return vector(a, make_vec(ts...)); } template inline T read() { T t; cin >> t; return t; } template inline vector readv(size_t sz) { vector ret(sz); rep(i, sz) cin >> ret[i]; return ret; } template inline tuple reads() { return {read()...}; } template struct edge { int to; T cost; edge(int t, T c) : to(t), cost(c) {} }; using ll = long long; using pii = pair; using pll = pair; using Graph = vector>; template using WGraph = vector>>; const int INF = 1 << 30; const ll LINF = 1LL << 60; const int MOD = 1e9 + 7; // 参考:https://qiita.com/drken/items/1b7e6e459c24a83bb7fd // 1-indexed BIT // add:log(N),sum:log(N),get_k-th:log(N) // https://codeforces.com/contest/1354/problem/D struct BIT { vector dat; BIT(int n) { dat.resize(n + 1); for(int i = 0; i < (int)dat.size(); i++) dat[i] = 0; } // a is 1-indexed void add(int a, ll x) { for(int i = a; i < (int)dat.size(); i += i & -i) dat[i] = dat[i] + x; } // [1,a], a is 1-indexed ll sum(int a) { ll res = 0; for(int i = a; i > 0; i -= i & -i) { res = res + dat[i]; } return res; } // [a,b), a and b are 1-indexed ll sum(int a, int b) { return sum(b - 1) - sum(a - 1); } // k-th number (k is 1-indexed) // lower_bound int get(ll k) { if(k <= 0) return 0; int res = 0; int N = 1; while(N < (int)dat.size()) N *= 2; for(int i = N / 2; i > 0; i /= 2) { if(res + i < (int)dat.size() && dat[res + i] < k) { k = k - dat[res + i]; res = res + i; } } return res + 1; } }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int N = read(); auto A = readv(N); BIT bit(200010); ll sum = 0; rep(i, N) { bit.add(A[i], 1); sum += A[i]; } map table; ll sub = 0; for(int i = 0; i < N; i++) { if(table.count(A[i])) { sub += A[i] * table[A[i]]; } else { ll modsum = 0; for(int j = 1; A[i] * j <= 200000; j++) { modsum += bit.sum(j * A[i], min(200001LL, (j + 1) * A[i])) * j; } sub += A[i] * modsum; table[A[i]] = modsum; } } ll res = sum * N - sub; cout << res << endl; }