#include int ri() { int n; scanf("%d", &n); return n; } int main() { int n = ri(); std::map a; std::vector > sum; for (int i = 0; i < n; i++) a[ri()]++; std::pair cur = {0, 0}; for (auto &i : a) { sum.push_back(std::tuple (i.first, cur.first, cur.second)); cur.first += i.second; cur.second += (int64_t) i.first * i.second; } int cache[801]; for (int i = 0; i <= 800; i++) cache[i] = -1; int64_t ans = 0; for (auto &i : a) { int64_t res = 0; if (i.first <= 800) { if (cache[i.first] == -1) { cache[i.first] = 0; for (auto &j : a) cache[i.first] += (j.first % i.first) * j.second; } res = cache[i.first]; } else { for (int j = 0; j < 200000; j += i.first) { auto low = std::lower_bound(sum.begin(), sum.end(), std::tuple(j, 0, 0)); auto high = std::lower_bound(sum.begin(), sum.end(), std::tuple(j + i.first, 0, 0)); int64_t xsum = (high == sum.end() ? cur.second : std::get<2>(*high)) - (low == sum.end() ? cur.second : std::get<2>(*low)); int num = (high == sum.end() ? cur.first : std::get<1>(*high)) - (low == sum.end() ? cur.first : std::get<1>(*low)); // std::cerr << xsum << " " << num * i.first * (j / i.first) << std::endl; res += xsum; res -= (int64_t) num * i.first * (j / i.first); } } res *= i.second; ans += res; } std::cout << ans << std::endl; return 0; }