#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
map<int,ll>m;
int n;
void inc(int x, int idx){
for(int i=1;i<=sqrt(x);i++){
if(x%i==0)
m[i] += (n-idx);
if(i*i != x)
m[x/i] += (n-idx);
}
}
int main() {
int t;
scanf("%d",&t);
while(t--){
scanf("%d",&n);
int a[n];
for(int i=0;i<n;i++)
scanf("%d",&a[i]);
ll ans = 0;
for(int i=n-1;i>=0;i--){
ans += 1ll*(i+1)*m[a[i]];
inc(a[i],i);
}
cout<<ans<<endl;
m.clear();
}
return 0;
}