mod=10**9+7
def solve():
  N=int(input())
  A=list(map(int,input().split()))
  P=A[0]
  for i in range(1,N):
    P=(P+A[i]+P*A[i])%mod
  print(P)

for i in range(int(input())):
  solve()