import sys

sys.setrecursionlimit(10 ** 7)
input = sys.stdin.readline
f_inf = float('inf')
mod = 10 ** 9 + 7


class BIT:
    """ 1-indexで使用しないと無限ループするので注意 """

    def __init__(self, n):
        self.n = n
        self.bit = [0] * (n + 1)

    def update(self, idx, x):
        """ idx番目にxをO(logN)で加算する """
        while idx <= self.n:
            self.bit[idx] += x
            idx += idx & (-idx)

    def query(self, idx):
        """ 1 ~ idx番目までの区間和をO(logN)で返す """
        res = 0
        while idx > 0:
            res += self.bit[idx]
            idx -= idx & (-idx)
        return res

    def sec_sum(self, left, right):
        """ left ~ right番目までの区間和をO(logN)で返す """
        return self.query(right) - self.query(left - 1)

    def lower_bound(self, w):
        """ 区間和がw以上になる最小のindexをO(logN)で返す """
        if w <= 0:
            return 0
        x = 0
        k = 1 << self.n.bit_length()
        while k:
            if x + k <= self.n and self.bit[x + k] < w:
                w -= self.bit[x + k]
                x += k
            k //= 2
        return x + 1

    def debug(self):
        """ bitの状態を出力する """
        print(*[self.sec_sum(i, i) for i in range(1, self.n + 1)])


def resolve():
    n = int(input())
    A = list(map(int, input().split()))

    B = []
    for idx, a in enumerate(A):
        B.append([a, idx])
    B.sort()

    A_comp = [0] * n
    prev = -1
    i = 0
    for a, idx in B:
        if prev != a:
            i += 1
        A_comp[idx] = i
        prev = a

    bit_cnt = BIT(n)
    bit_total = BIT(n)
    cnt = [[0, 0] for _ in range(n)]
    total = [[0, 0] for _ in range(n)]
    for i in range(n):
        cnt[i][0] = bit_cnt.sec_sum(A_comp[i] + 1, n)
        total[i][0] = bit_total.sec_sum(A_comp[i] + 1, n) % mod
        bit_cnt.update(A_comp[i], 1)
        bit_total.update(A_comp[i], A[i])

    bit_cnt = BIT(n)
    bit_total = BIT(n)
    for i in reversed(range(n)):
        cnt[i][1] = bit_cnt.query(A_comp[i] - 1)
        total[i][1] = bit_total.query(A_comp[i] - 1) % mod
        bit_cnt.update(A_comp[i], 1)
        bit_total.update(A_comp[i], A[i])

    res = 0
    for i in range(n):
        left_cnt, right_cnt = cnt[i]
        left_tot, right_tot = total[i]
        res += left_cnt * right_tot % mod + right_cnt * left_tot % mod + A[i] * left_cnt % mod * right_cnt % mod
        res %= mod
    print(res)


if __name__ == '__main__':
    resolve()