mod = 10**9 + 7
n = int(input())
x = list(map(int,input().split()))
r = sum(x)
ans = 0
for i in range(n):
	ans += x[i] * (r-x[i]) % mod
	ans %= mod

ans *= pow(r, mod-2, mod)
ans %= mod
ans += 1
ans %= mod

print(ans)