#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>
#include <queue>
using namespace std;
typedef long long ll;
const int N=1005;
const int MOD=1E8+7;
int d[N][N],f[N][N],n,T,ans,x,w;

void updata(int &u,int w)
{
    u=(u+w)%MOD;
}

int main()
{
    for (int i=0;i<N;i++)
        for (int j=0;j<N;j++) d[i][j]=__gcd(i,j);
    cin>>T;
    while (T--)
    {
        cin>>n;
        memset(f,0,sizeof(f));
        f[0][0]=1; ans=0;
        for (int i=1;i<=n;i++)
        {
            scanf("%d",&x);
            for (int j=0;j<N;j++) if (w=f[i-1][j])
            {
                updata(f[i][j],w);
                updata(f[i][d[j][x]],w);
            }
        }
        for (int j=1;j<N;j++) updata(ans,(ll)f[n][j]*j%MOD);
        cout<<ans<<endl;
    }
}