D. Small GCD

Problem - 1900D - Codeforces  

求和公式降维变形；

欧拉函数 利用欧拉反演替换 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;
}