T = int(input())
P = 10 ** 9 + 7

def f(x,y):
    return (x + y + x * y)%P
for _ in range(T):
    N = int(input())
    A = list(map(int,input().split()))
    ans = 0
    for a in A:
        ans = f(ans,a)
    print(ans)