mod = 10 ** 9 + 7
t = int(input())
for _ in range(t):
    n = int(input())
    a = list(map(int, input().split()))
    a.sort()
    a2 = []
    for i in a:
        if i == 0:
            continue
        a2.append(i)
    ans = 0
    n = len(a2)
    for i in range(n):
        for j in range(i + 1, n):
            ans += (a2[i] + a2[j]) + (a2[i] * a2[j])
            ans %= mod
    print(ans)