N = int(input())
A = [int(x) for x in input().split()]

mod = 10**9 + 7
A.sort()
c = []
prev = 0
for ai in A:
    if ai != prev:
        c += [1]
        prev = ai
    else:
        c[-1] += 1

res = 0
su_d = sum(c[2:])
for p in range(len(c) - 2):
    su = su_d
    for q in range(p + 1, len(c) - 1):
        res += c[p] * c[q] * su % mod
        res %= mod
        su -= c[q + 1]
    su_d -= c[p + 2]

print(res)