#include<bits/stdc++.h>
namespace {
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-function"
#include<atcoder/all>
#pragma GCC diagnostic pop
using namespace std;
using namespace atcoder;
#define rep(i,n) for(int i = 0; i < (int)(n); i++)
#define rrep(i,n) for(int i = (int)(n) - 1; i >= 0; i--)
#define all(x) begin(x), end(x)
#define rall(x) rbegin(x), rend(x)
template<class T> bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } else return false; }
template<class T> bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } else return false; }
using ll = long long;
using P = pair<int,int>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<ll>;
using VVL = vector<VL>;
using mint = modint1000000007;

mint op(mint x, mint y) { return x + y; }
mint e() { return mint(); }

template<class S>
struct value_compression : vector<S> {
  bool built = false;
  using VS = vector<S>;
  using VS::vector;
  value_compression(vector<S> v) : vector<S>(move(v)) {}
  void build() {
    sort(VS::begin(), VS::end());
    VS::erase(unique(VS::begin(), VS::end()), VS::end());
    built = true;
  }
  template<class T>
  void convert(T first, T last) { for (; first != last; ++first) *first = (*this)(*first); }
  int operator()(S x) {
    assert(built);
    return lower_bound(VS::begin(), VS::end(), x) - VS::begin();
  }
  void clear() { VS::clear(); built = false; }
};

} int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n;
  cin >> n;
  VI a(n);
  rep(i, n) cin >> a[i];
  value_compression vc = a;
  vc.build();
  vc.convert(all(a));
  int sz = vc.size();
  vector<mint> cf1(n), cf2(n), cb1(n), cb2(n);
  segtree<mint, op, e> seg(sz);
  mint p2 = 1;
  rep(i, n) {
    cf1[i] = seg.prod(a[i] + 1, sz);
    cf2[i] = seg.prod(0, a[i]);
    seg.set(a[i], seg.get(a[i]) + p2);
    p2 *= 2;
  }
  p2 = 1;
  rep(i, sz) seg.set(i, 0);
  rrep(i, n) {
    cb1[i] = seg.prod(a[i] + 1, sz);
    cb2[i] = seg.prod(0, a[i]);
    seg.set(a[i], seg.get(a[i]) + p2);
    p2 *= 2;
  }
  mint ans;
  rep(i, n) ans += cf1[i] * cb1[i] + cf2[i] * cb2[i];
  cout << ans.val() << '\n';
}