#include using namespace atcoder; using mint = modint1000000007; // using mint = modint998244353; // using mint = modint;//mint::set_mod(MOD); const long long MOD = 1000000007; // const long long MOD = 998244353; #include #define rep(i, a, b) for (ll i = (ll)(a); i < (ll)(b); i++) #define repeq(i, a, b) for (ll i = (ll)(a); i <= (ll)(b); i++) #define repreq(i, a, b) for (ll i = (ll)(a); i >= (ll)(b); i--) #define each(a, b) for (auto &(a) : (b)) #define endl '\n' // fflush(stdout); #define cYes cout << "Yes" << endl #define cNo cout << "No" << endl #define sortr(v) sort(v, greater<>()) #define pb push_back #define mp make_pair #define mt make_tuple #define tget(a, b) get(a) #define FI first #define SE second #define ALL(v) (v).begin(), (v).end() #define INFLL 3000000000000000100LL #define INF 1000000100 #define PI acos(-1.0L) #define TAU (PI * 2.0L) using namespace std; typedef long long ll; typedef pair Pll; typedef tuple Tlll; typedef vector Vi; typedef vector VVi; typedef vector Vl; typedef vector VVl; typedef vector VVVl; typedef vector VTlll; typedef vector Vm; typedef vector VVm; typedef vector Vs; typedef vector Vd; typedef vector Vc; typedef vector Vb; typedef vector VPll; typedef priority_queue PQl; typedef priority_queue, greater> PQlr; /* print */ template ostream &operator<<(ostream &os, const vector &V) { int N = V.size(); if (N == 0) { os << endl; return os; } rep(i, 0, N - 1) { os << V[i] << ' '; } os << V[N - 1] << endl; return os; } template ostream &operator<<(ostream &os, const vector> &V) { int N = V.size(); rep(i, 0, N) os << V[i]; return os; } template ostream &operator<<(ostream &os, pair const &P) { os << P.FI << ' ' << P.SE; return os; } ostream &operator<<(ostream &os, mint const &M) { os << M.val(); return os; } /* useful */ template void Vin(vector &v) { int n = v.size(); rep(i, 0, n) cin >> v[i]; } template int SMALLER(vector &a, T x) { return lower_bound(a.begin(), a.end(), x) - a.begin(); } template int orSMALLER(vector &a, T x) { return upper_bound(a.begin(), a.end(), x) - a.begin(); } template int BIGGER(vector &a, T x) { return a.size() - orSMALLER(a, x); } template int orBIGGER(vector &a, T x) { return a.size() - SMALLER(a, x); } template int COUNT(vector &a, T x) { return upper_bound(ALL(a), x) - lower_bound(ALL(a), x); } template bool chmax(T &a, T b) { if (a < b) { a = b; return 1; } return 0; } template bool chmin(T &a, T b) { if (a > b) { a = b; return 1; } return 0; } template void press(T &v) { v.erase(unique(ALL(v)), v.end()); } template vector zip(vector b) { pair p[b.size() + 10]; int a = b.size(); vector l(a); for (int i = 0; i < a; i++) p[i] = mp(b[i], i); sort(p, p + a); int w = 0; for (int i = 0; i < a; i++) { if (i && p[i].first != p[i - 1].first) w++; l[p[i].second] = w; } return l; } template vector vis(vector &v) { vector S(v.size() + 1); rep(i, 1, S.size()) S[i] += v[i - 1] + S[i - 1]; return S; } ll dem(ll a, ll b) { return ((a + b - 1) / (b)); } ll dtoll(double d, int g) { return round(d * pow(10, g)); } const double EPS = 1e-10; void init() { cin.tie(0); cout.tie(0); ios::sync_with_stdio(0); cout << fixed << setprecision(12); } // Comm + K -> Comm + 0 // Comm + K -> Comm + C = add// // Comm + K -> Comm + U = del// /********************************** START **********************************/ mint st_add(mint a, mint b) { return a + b; } mint st_zero() { return 0; } struct RAQ { int n; segtree st; RAQ(int input, mint init = 0) { n = input; st = segtree(vector(n + 1, init)); } void set(int i, mint k) { apply(i, i + 1, k - get(i)); } void apply(int l, int r, mint val) { st.set(l, st.get(l) + val); st.set(r, st.get(r) - val); } mint get(int i) { return st.prod(0, i + 1); } void out() { for (int i = 0; i < n; i++) cout << get(i) << " \n"[i == n - 1]; } void err() { for (int i = 0; i < n; i++) cerr << get(i) << " \n"[i == n - 1]; } mint operator[](int i) { return get(i); } }; void sol() { int n; cin >> n; Vl a(n); Vin(a); // cerr << zip(a) << endl; mint ans = 0; auto f = [&]() { map ma; mint pre = 0; RAQ st(n), str(n); rep(i, 0, n) { ma[a[i]].pb(i); } for (auto p : ma) { auto key = p.first; auto val = p.second; rep(i, 0, val.size()) { mint l, r, s; l = st[val[i]]; // s = st[n - 1]; // r = s - l; // r /= pow_mod(2, val[i] + 1, MOD); mint lr, rr, sr; lr = str[val[i]]; sr = str[n - 1]; rr = sr - lr; pre += l * rr; // /cerr << lr << ' ' << r << endl; } rep(i, 0, val.size()) { mint ad = pow_mod(2, val[i], MOD); st.apply(val[i], n, ad); mint adr = pow_mod(2, n - 1 - val[i], MOD); str.apply(val[i], n, adr); } } return pre; }; ans += f(); rep(i, 0, n) a[i] *= -1; ans += f(); cout << ans << endl; } int main() { init(); int q = 1; // cin >> q; while (q--) sol(); return 0; }