#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int MAXN=1000+5,MOD=100000007;
int cnt[MAXN],_2[MAXN],gcd[MAXN];

int main(){
#ifdef bhiaibogf
    freopen("in.in","r",stdin);
#endif
    int t;
    _2[0]=1;
    for(int i=1;i<MAXN;i++)
        _2[i]=2*_2[i-1]%MOD;
    scanf("%d",&t);
    while(t--){
        int n,_max=0;
        scanf("%d",&n);
        memset(cnt,0,sizeof cnt);
        for(int i=0;i<n;i++){
            int ls;
            scanf("%d",&ls);
            _max=max(_max,ls);
            for(int j=1;j<=ls;j++) if(ls%j==0) cnt[j]++;
        }
        int ans=0;
        for(int i=_max;i>0;i--){
            gcd[i]=_2[cnt[i]]-1;
            for(int j=2;i*j<=_max;j++) (gcd[i]-=gcd[i*j])%=MOD;
            (ans+=1ll*gcd[i]*i%MOD)%MOD;
        }
        printf("%d\n",(ans+MOD)%MOD);
    }
    return 0;
}