////#include<bits/stdc++.h>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
#include <string>
#include <iostream>
#include <algorithm>
using namespace std;
#include <queue>
#include <stack>
#include <vector>
#include <bitset>
#include <deque>
#include <set>
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <map>
typedef long long LL;
const int INF = 0x3f3f3f;
const int MOD = 1e8+7;
const int MAXN = 400000;
const double eps =1e-6;
const double pi = acos(-1.0);
using namespace std;
int gcd(int a,int b)
{
    if(b==0) return a;
    else return gcd(b,a%b);
}
int a[1234],n;
LL dp[2][1010];
void Gao(LL &a,LL b)
{
    a=(a+b)%MOD;
}
int w[1010][1010];
int main()
{
    for(int i=1;i<=1000;i++)
        for(int j=1;j<=1000;j++)
            w[i][j]=gcd(i,j);
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%d",&a[i]);
        memset(dp,0,sizeof(dp));
        int now=0;
        for(int i=0;i<n;i++)
        {
            memset(dp[now],0,sizeof(dp[now]));
            for(int j=1;j<=1000;j++)
            {
                if(dp[now^1][j]){
                    Gao(dp[now][w[j][a[i]]],dp[now^1][j]);
                    Gao(dp[now][j],dp[now^1][j]);
                }
            }
            Gao(dp[now][a[i]],1LL);
            now^=1;
        }
        LL ans=0;
        for(int i=1;i<=1000;i++)
        {
            ans=(ans+dp[now^1][i]*i%MOD)%MOD;
        }
        printf("%I64d\n",ans);
    }
    return 0;
}