#include #define REP(i, n) for (int i = 0; (i) < (int)(n); ++ (i)) #define REP3(i, m, n) for (int i = (m); (i) < (int)(n); ++ (i)) #define REP_R(i, n) for (int i = (int)(n) - 1; (i) >= 0; -- (i)) #define ALL(x) std::begin(x), std::end(x) using namespace std; template inline void chmin(T & a, U const & b) { a = min(a, b); } /** * @note O(\sqrt{n}) */ struct prepared_primes { int size; std::vector sieve; std::vector primes; prepared_primes(int size_) : size(size_) { sieve.resize(size); REP3 (p, 2, size) if (sieve[p] == 0) { primes.push_back(p); for (int k = p; k < size; k += p) { if (sieve[k] == 0) { sieve[k] = p; } } } } }; int64_t solve(int n, const vector & a) { // Gaussian elimination vector basis; for (int64_t a_i : a) { for (int k = 0; (a_i << k) < (1ull << 40); ++ k) { int64_t x = a_i << k; for (int64_t b : basis) { chmin(x, x ^ b); } if (x) { basis.push_back(x); } } } // list constructibles int64_t max_a = *max_element(ALL(a)); vector constructed(max_a + 1); constructed[0] = true; for (int64_t b : basis) { REP (x, max_a) if (constructed[x]) { if ((b ^ x) < constructed.size()) { constructed[b ^ x] = true; } } } // fast zeta transform prepared_primes primes(1e5 + 3); vector acc(max_a + 1); for (int64_t a_i : a) { acc[a_i] += a_i; } for (int64_t p : primes.primes) { REP_R (x, max_a / p + 1) { acc[x] += acc[x * p]; } } // answer int64_t sum = accumulate(ALL(a), 0ll); int64_t answer = INT64_MAX; REP3 (x, 1, max_a + 1) if (constructed[x]) { chmin(answer, sum - acc[x] + acc[x] / x); } return answer; } int main() { int n; cin >> n; vector a(n); REP (i, n) { cin >> a[i]; } cout << solve(n, a) << endl; return 0; }