#include // #include using namespace std; // using namespace atcoder; #define DEBUG(x) cerr<<#x<<": "< #define vl vector #define vii vector< vector > #define vll vector< vector > #define vs vector #define pii pair #define pis pair #define psi pair #define pll pair template pair operator+(const pair &s, const pair &t) { return pair(s.first + t.first, s.second + t.second); } template pair operator-(const pair &s, const pair &t) { return pair(s.first - t.first, s.second - t.second); } template ostream& operator<<(ostream& os, pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } #define X first #define Y second #define rep(i,n) for(int i=0;i<(int)(n);i++) #define rep1(i,n) for(int i=1;i<=(int)(n);i++) #define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--) #define rrep1(i,n) for(int i=(int)(n);i>0;i--) #define REP(i,a,b) for(int i=a;i bool chmax(T &a, const T &b) { if (a bool chmin(T &a, const T &b) { if (a>b) { a = b; return 1; } return 0; } #define UNIQUE(v) v.erase(std::unique(v.begin(), v.end()), v.end()); const ll inf = 1000000001; const ll INF = (ll)1e18 + 1; const long double pi = 3.1415926535897932384626433832795028841971L; #define Sp(p) cout<> n; vi a(n); rep (i, n) cin >> a[i]; sort(all(a)); vl sum(n + 1); rep (i, n) sum[i + 1] = sum[i] + a[i]; ll pre_value = 0; ll ans = 0; ll small_sum = 0; rep (i, n) { if (i > 0 and a[i] == a[i - 1]) { ans += pre_value; small_sum += a[i]; continue; } ll add = 0; add += small_sum; int left = i + 1; for (ll mul = 1; a[i] * mul <= 211111; mul++) { if (left == n) break; if (a[left] >= a[i] * (mul + 1)) continue; int ok = left, ng = n; while (ok + 1 < ng) { int mid = (ok + ng) / 2; if (a[mid] >= a[i] * (mul + 1)) ng = mid; else ok = mid; } add += (sum[ok + 1] - sum[left]) - (ok + 1 - left) * mul * a[i]; left = ok + 1; } ans += add; pre_value = add; small_sum += a[i]; } cout << ans << endl; }