#if __INCLUDE_LEVEL__ == 0 #include __BASE_FILE__ void Solve() { int n; IN(n); vector a(n); IN(a); int A = ranges::max(a); linear_sieve::init(A); vector> edges; for (int i : Rep(0, n)) { auto f = linear_sieve::factor(a[i]); for (int l = 0; l < Sz(f);) { int e = 1; while (l + e < Sz(f) && f[l] == f[l + e]) { ++e; } int p = f[l]; edges.emplace_back(p * e, i, n + p); l += e; } } ranges::sort(edges); atcoder::dsu dsu(n + A + 1); int64_t ans = 0; for (auto [w, i, j] : Rev(edges)) { if (dsu.same(i, j)) { ans += w; } else { dsu.merge(i, j); } } OUT(ans); } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); Solve(); } #elif __INCLUDE_LEVEL__ == 1 #include #include namespace linear_sieve { std::vector primes, lpf; void init(int n) { if (n < int(std::size(lpf))) return; if (n < 2 * int(std::size(lpf))) n = 2 * std::size(lpf); lpf.resize(n + 1, -1); for (int d = 2; d <= n; ++d) { if (lpf[d] == -1) lpf[d] = d, primes.push_back(d); for (int p : primes) { if (p * d > n or p > lpf[d]) break; lpf[p * d] = p; } } } std::vector factor(int n) { __glibcxx_assert(n >= 1); std::vector res; for (init(n); n > 1; n /= res.back()) res.push_back(lpf[n]); return res; } } // namespace linear_sieve template concept MyRange = std::ranges::range && !std::convertible_to; template concept MyTuple = std::__is_tuple_like::value && !MyRange; namespace std { istream& operator>>(istream& is, MyRange auto&& r) { for (auto&& e : r) is >> e; return is; } istream& operator>>(istream& is, MyTuple auto&& t) { apply([&](auto&... xs) { (is >> ... >> xs); }, t); return is; } ostream& operator<<(ostream& os, MyRange auto&& r) { auto sep = ""; for (auto&& e : r) os << exchange(sep, " ") << e; return os; } ostream& operator<<(ostream& os, MyTuple auto&& t) { auto sep = ""; apply([&](auto&... xs) { ((os << exchange(sep, " ") << xs), ...); }, t); return os; } } // namespace std using namespace std; #define Rev views::reverse #define Rep(...) [](int l, int r) { return views::iota(min(l, r), r); }(__VA_ARGS__) #define Sz(r) int(size(r)) #define IN(...) (cin >> forward_as_tuple(__VA_ARGS__)) #define OUT(...) (cout << forward_as_tuple(__VA_ARGS__) << '\n') #endif // __INCLUDE_LEVEL__ == 1