Problem:
Score: 400 points

Problem Statement

Quickly after finishing the tutorial of the online game ATChat, you have decided to visit a particular place with N-1 players who happen to be there. These N players, including you, are numbered 1 through N, and the friendliness of Player i is A_i.

The N players will arrive at the place one by one in some order. To make sure nobody gets lost, you have set the following rule: players who have already arrived there should form a circle, and a player who has just arrived there should cut into the circle somewhere.

When each player, except the first one to arrive, arrives at the place, the player gets comfort equal to the smaller of the friendliness of the clockwise adjacent player and that of the counter-clockwise adjacent player. The first player to arrive there gets the comfort of 0.

What is the maximum total comfort the N players can get by optimally choosing the order of arrivals and the positions in the circle to cut into?

Constraints

All values in input are integers.

2 \leq N \leq 2 \times 10^5

1 \leq A_i \leq 10^9

Input

Input is given from Standard Input in the following format:

N
A_1 A_2 \dots A_N

Output

Print the maximum total comfort the N players can get.

Sample Input 1

4
2 2 1 3

Sample Output 1

7

By arriving at the place in the order Player 4, 2, 1, 3, and cutting into the circle as shown in the figure, they can get the total comfort of 7.

They cannot get the total comfort greater than 7, so the answer is 7.

Sample Input 2

7
1 1 1 1 1 1 1

Sample Output 2

6

Solution in C++ (GCC 9.2.1) (judge verdict: Runtime Error)

#include <bits/stdc++.h>

using namespace std;
using ll = long long;
using Graph = vector<vector<int>>;
using WGraph = vector<vector<pair<int, ll>>>;

template<class T>inline bool chmax(T &a, const T &b) { if (b > a) { a = b; return true; } return false; }
template<class T>inline bool chmin(T &a, const T &b) { if (b < a) { a = b; return true; } return false; }

constexpr int dx[4] = {-1, 0, 0, 1}, dy[4] = {0, -1, 1, 0};
constexpr int INF = 1e9;
constexpr int MOD = 1e9 + 7;
constexpr long long LINF = 1e18;
constexpr double EPS = 1e-10;
constexpr double PI = M_PI;

template<class RandomAccessIterator>
auto medianOfMedians(RandomAccessIterator first, RandomAccessIterator last, int k) {
    auto _first = *first;
    using T = decltype(_first);

    int len = last - first;
    if (len == 1 || k == 0) return *first;

    int ms = (len + 4) / 5;
    std::vector<T> medians(ms);
    for (int i=0; i<ms; ++i) {
        auto iter = first + 5*i;
        if (5*i + 4 < len) {
            for (int j=0; j<3; ++j) {
                int minInd = j;
                for (int k=j+1; k<5; ++k) {
                    if (*(iter+k) < *(iter+minInd)) minInd = k;
                }
                swap(*(iter+j), *(iter+minInd));
            }
            medians[i] = *(iter+2);
        } else {
            medians[i] = *iter;
        }
    }

    T pivot = medianOfMedians(medians.begin(), medians.end(), ms/2);

    for (int i=0; i<len; ++i) {
        if (*(first+i) == pivot) {
            swap(*(first+i), *(first+len-1));
            break;
        }
    }

    int store = 0;
    for (int i=0; i<len; ++i) {
        if (*(first+i) < pivot) {
            swap(*(first+i), *(first+store));
            ++store;
        }
    }
    swap(*(first+store), *(first+len-1));

    if (store == 0 || store == k) {
        return pivot;
    } else if (k < store) {
        return medianOfMedians(first, first+k, k);
    } else {
        return medianOfMedians(first+k+1, last, k-store-1);
    }
}

void solve() {
    int N;
    cin >> N;
    vector<int> a(N);
    for (int i=0; i<N; ++i) cin >> a[i];

    int maxi = -INF;
    for (int x : a) chmax(maxi, x);
    int medi = medianOfMedians(a.begin(), a.end(), N/2);
    ll res = 0;
    int cnt = -1;
    for (int x : a) {
        if (medi < x) {
            res += 2 * x;
            cnt += 2;
        }
    }
    res -= maxi;
    res += (ll)medi * (N - 1 - cnt);
    cout << res << '\n';
}

int main() {
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout << fixed << setprecision(20);

    solve();

    return 0;
}