import sys
input=sys.stdin.readline

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

K=0;Mod=10**9+7;two_inv=pow(2,Mod-2,Mod)
A_max=max(A)

X=[0]*(A_max+1)
for a in A:
    X[a]+=1

for a in A:
    x1=x2=0
    for b in range(a):
        x1+=X[b]
        x2+=X[b]*X[b]
    K+=x1*x1-x2
    K%=Mod

K=(K*two_inv)%Mod
print(K)