求和公式降维变形;
欧拉函数 利用欧拉反演替换 gcd;分配性质中的 ai 和 aj,依次同步,对 i 枚举,对 j 使用前缀计数统计;
#include <bits/stdc++.h>
#define int long long
using namespace std;
const int max_n=8e4+5,max_m=1e6;
int cnt,n;
int a[max_n],prime[max_m],phi[max_m],vis[max_m],sum[max_m];
void init()
{
phi[1]=1;
for (int i=2;i<=max_m;++i){
if(!vis[i])prime[++cnt]=i,phi[i]=i-1;
for (int j=1;j<=cnt&&i*prime[j]<=max_m;++j){
vis[i*prime[j]]=1;
if(i%prime[j]==0){phi[i*prime[j]]=phi[i]*prime[j];break;}
phi[i*prime[j]]=phi[i]*phi[prime[j]];
}
}
}
int fuck(int x)
{
int res=0;
for(int i=1;i*i<=x;i++){
if(!(x%i)){
res+=phi[i]*(sum[i]++);
if(i*i!=x)res+=phi[x/i]*(sum[x/i]++);
}
}
return res;
}
void clear()
{
for(int i=1;i<=max_m;i++)sum[i]=0;
}
void solve()
{
cin>>n;
for(int i=1;i<=n;i++)cin>>a[i];
sort(a+1,a+1+n);
int ans=0;
for(int i=1;i<=n;i++){
ans+=fuck(a[i])*(n-i);
}
cout<<ans<<endl;
clear();
}
signed main()
{
ios::sync_with_stdio(false);
cin.tie(nullptr);
init();
int _t;cin>>_t;
while(_t--){
solve();
}
return 0;
}