import sys
input = sys.stdin.readline

mod=10**9+7

N=int(input())
A=list(map(int,input().split()))

def solve(A):
    A=list(A)

    S=sorted(set(A))
    D={S[i]:i+1 for i in range(len(S))}

    for i in range(N):
        A[i]=D[A[i]]

    # BIT(BIT-indexed tree)
    LEN=len(S)+1

    BIT=[0]*(LEN+1) # 1-indexedなtree. 配列BITの長さはLEN+1にしていることに注意。

    def update(v,w): # index vにwを加える
        while v<=LEN:
            BIT[v]+=w
            v+=(v&(-v)) # v&(-v)で、最も下の立っているビット. 自分を含む大きなノードへ. たとえばv=3→v=4

    def getvalue(v): # [1,v]の区間の和を求める
        ANS=0
        while v!=0:
            ANS+=BIT[v]
            v-=(v&(-v)) # 自分より小さい自分の和を構成するノードへ. たとえばv=14→v=12へ
        return ANS


    BIT2=[0]*(LEN+1) # 1-indexedなtree. 配列BITの長さはLEN+1にしていることに注意。

    def update2(v,w): # index vにwを加える
        while v<=LEN:
            BIT2[v]+=w
            v+=(v&(-v)) # v&(-v)で、最も下の立っているビット. 自分を含む大きなノードへ. たとえばv=3→v=4

    def getvalue2(v): # [1,v]の区間の和を求める
        ANS=0
        while v!=0:
            ANS+=BIT2[v]
            v-=(v&(-v)) # 自分より小さい自分の和を構成するノードへ. たとえばv=14→v=12へ
        return ANS

    POW2=[1]
    for i in range(10**5+10):
        POW2.append(POW2[-1]*2%mod)


    for i in range(N):
        a=A[i]
        update(a,POW2[i])

    ANS=0

    for i in range(N-1,-1,-1):
        a=A[i]
        ANS+=getvalue(a-1)*getvalue2(a-1)
        ANS%=mod

        update(a,-POW2[i])
        update2(a,POW2[N-1-i])

    return ANS


B=[10**10-a for a in A]
print(solve(tuple(A))+solve(tuple(B)))