#include<iostream>
#include<algorithm>

using namespace std;


int gcd(int a,int b)
{
    if(b==0) return a;
    return gcd(b,a%b);
}

int num[1024];
int dp[2][1024];
int gcd_[1024][1024];

int main()
{
    for(int i=1;i<=1010;i++)
    {
        for(int j=1;j<=1010;j++)
        {
            gcd_[i][j]=gcd(i,j);
        }
    }



    int t;
    cin>>t;
    while(t--)
    {
        int n;
        cin>>n;
        for(int i=0;i<n;i++)
        {
            scanf("%d",num+i);
        }

        int cur=0;
        memset(dp[cur],0,sizeof dp[cur]);

        for(int i=0;i<n;i++)
        {
            cur^=1;
            memset(dp[cur],0,sizeof dp[cur]);

            dp[cur][num[i]]=1;
            for(int j=1;j<=1000;j++)
            {
                if(dp[cur^1][j]!=0)
                {
                    dp[cur][gcd_[j][num[i]]]+=dp[cur^1][j];
                    dp[cur][gcd_[j][num[i]]]%=100000007;
                    dp[cur][j]+=dp[cur^1][j];
                    dp[cur][j]%=100000007;
                }
            }
        }
        long long ans=0;
        for(int i=1;i<=1000;i++)
        {
            ans+=1LL*dp[cur][i]*i%100000007;
            ans%=100000007;
        }
        cout<<ans<<endl;
    }
    return 0;
}