#pragma GCC optimize("Ofast")
#define _USE_MATH_DEFINES
#include "bits/stdc++.h"

using namespace std;

using u8 = uint8_t;
using u16 = uint16_t;
using u32 = uint32_t;
using u64 = uint64_t;
using i8 = int8_t;
using i16 = int16_t;
using i32 = int32_t;
using i64 = int64_t;

constexpr char newl = '\n';
constexpr double eps = 1e-10;

#define FOR(i,a,b) for (int i = (a); i < (b); i++)
#define F0R(i,b) FOR(i,0,b)
#define RFO(i,a,b) for (int i = ((b)-1); i >=(a); i--)
#define RF0(i,b) RFO(i,0,b)
//#define fi first
//#define se second
#define debug(x) cout << #x << ": " << x << '\n';
#define rng(a) a.begin(),a.end()
#define rrng(a) a.rbegin(),a.rend()

template<typename T1, typename T2> inline void chmin(T1& a, T2 b) { if (a > b) a = b; }
template<typename T1, typename T2> inline void chmax(T1& a, T2 b) { if (a < b) a = b; }
template<class T> void Print(vector<T> v) {
  F0R(i, v.size()) {
    cout << v[i] << ' ';
  }
  cout << newl;
}

#if 1

template<int MOD> struct Modint {
  long long val;
  constexpr Modint(long long v = 0) noexcept : val(v% MOD) {
    if (val < 0) val += MOD;
  }
  constexpr int getmod() { return MOD; }
  constexpr Modint operator - () const noexcept {
    return val ? MOD - val : 0;
  }
  constexpr Modint operator + (const Modint& r) const noexcept { return Modint(*this) += r; }
  constexpr Modint operator - (const Modint& r) const noexcept { return Modint(*this) -= r; }
  constexpr Modint operator * (const Modint& r) const noexcept { return Modint(*this) *= r; }
  constexpr Modint operator / (const Modint& r) const noexcept { return Modint(*this) /= r; }
  constexpr Modint& operator += (const Modint& r) noexcept {
    val += r.val;
    if (val >= MOD) val -= MOD;
    return *this;
  }
  constexpr Modint& operator -= (const Modint& r) noexcept {
    val -= r.val;
    if (val < 0) val += MOD;
    return *this;
  }
  constexpr Modint& operator *= (const Modint& r) noexcept {
    val = val * r.val % MOD;
    return *this;
  }
  constexpr Modint& operator /= (const Modint& r) noexcept {
    long long a = r.val, b = MOD, u = 1, v = 0;
    while (b) {
      long long t = a / b;
      a -= t * b; swap(a, b);
      u -= t * v; swap(u, v);
    }
    val = val * u % MOD;
    if (val < 0) val += MOD;
    return *this;
  }
  constexpr bool operator == (const Modint& r) const noexcept {
    return this->val == r.val;
  }
  constexpr bool operator != (const Modint& r) const noexcept {
    return this->val != r.val;
  }
  friend constexpr ostream& operator << (ostream& os, const Modint<MOD>& x) noexcept {
    return os << x.val;
  }
  friend constexpr Modint<MOD> modpow(const Modint<MOD>& a, long long n) noexcept {
    if (n == 0) return 1;
    auto t = modpow(a, n / 2);
    t = t * t;
    if (n & 1) t = t * a;
    return t;
  }

  friend constexpr Modint<MOD> Ncr(Modint<MOD> a, int n) {
    Modint<MOD> p = 1, q = 1;

    while (n > 0) {
      p *= a;
      a -= 1;
      q *= n--;
    }

    return p / q;
  }
};
template<int MOD>
using mint = Modint<MOD>;


constexpr int MOD = 1e9 + 7; // 1e9 + 7

// INSERT ABOVE HERE
signed main() {
  cin.tie(0);
  ios_base::sync_with_stdio(false);

  int N;
  cin >> N;

  vector<int> cs(N);
  F0R(i, N) {
    cin >> cs[i];
  }
  sort(rng(cs));

  mint<MOD> result = 0;

  // 数列の違う部分を変更する際、コストが少ない部分から順に変更していくことで、
  // 和が最小になる。

  // それぞれのCiにつき、払われるコストを全てのパターンで計算していけばいい

  // とある Ci が違う数列の組は 2 * 2 ^ (N-1) * 2 ^ (N-1) = 2 ^ (1 + N*2 - 2) 通り
  // 数列をソートし、降順にC1, C2 ... CN とする

  // すると C1 は明らかに最後の順番で使われるので、この部分のコストは C1 * 1 * 2^1 * 2 ^ (N-1) * 2 ^ (N-1)
  // C1 * A1 とする

  // C2 は、C1が違うときのみ最後から2番目(D=2), それ以外のとき最後(D=1)に使われる
  // D=2 のときは 2^2 * 2 ^ (N-2) * 2 ^ (N-2) 通り
  // A1 にその数を足せば A2 が得られるので　C2 * A2

  // D=a のときは 2^a * 2 ^ (N-a) * 2 ^ (N-a) = 2 ^ (N * 2 - a)

  // 解けないので解説を読む
  // pdfでも分からないので配信を見る

  // S xor T = X とすると、XとTがあればSは一意に定まる
  // 違いだけ考えればいいので、T=0000... に対する答えを出して、それ * 2^N すればいい
  // 言い換えると、全Xに対するコストの総和を取って、それ * 2^N

  // 数列Cは昇順にする
  
  // X = 00 Cost=0
  // X = 01 Cost=C2
  // X = 10 Cost=C1
  // X = 11 Cost = C1 * 2 + C2

  // X = ?????
  //       Ci
  //       ??? にある 1の数を全て足したのが、Ciの係数
  //     ?? の部分の1の数は関係ない

  // i の部分の1の合計: 2 ^ (N-i)
  // i より右の1の合計: 2 ^ (N-i-1) * (N-i)

  // 左の部分: 2^(i-1) 通り

  vector<mint<MOD>> two(N + 1);
  two[0] = 1;
  F0R(i, N) {
    two[i + 1] = two[i] * 2;
  }

  F0R(i, N) {
    mint<MOD> a = i < N - 1 ? two[N-i-2] : 0;
    result += two[i] * (two[N-i-1] + (a * (N-i-1))) * cs[i];
  }
  result *= two[N];

  cout << result.val;
}
#endif


Measured memory_bytes: 2560