[51nod1188]最大公约数之和 V2

Description

求

$$\sum_{i=1}^{n-1}\sum_{j=i+1}^{n}gcd(i,j)$$
n<=5*10^6,数据组数<=50000

Solution

一眼看到莫比乌斯反演（感觉中毒了）
但是会超时！！！
想想其他方法。
我们知道，如果要在区间[1..i]中询问gcd(i,x)=t的个数。那么也就相当于询问gcd(i/t,x/t)=1的个数。用一个欧拉函数就可以了。
那么能不能运用到这道题里呢？
我们发现，i一定小于j。枚举gcd那么我们就相当于求对于每一个j，[1..j-1]中有多少个数和它的gcd为那个数。
于是

$$Ans=\sum_{j=2}^{n}\sum_{t|j,t<j}\varphi({j\over t})*t$$
那么我们可以先枚举t，然后再上界里面与j有关的答案都刷出来。
时间复杂度和普通筛法一样。
然后求一遍前缀和，就可以O(1)回答了。

Code

#include<cmath>
#include<cstdio>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define N 5000000
using namespace std;
typedef long long ll;
ll a[N+5];
int p[N+5],phi[N+5],n,ty;
bool bz[N+5];
int main() {
    fo(i,2,N) {
        if (!bz[i]) p[++p[0]]=i,phi[i]=i-1;
        fo(j,1,p[0]) {
            int k=i*p[j];if (k>N) break;
            bz[k]=1;if (!(i%p[j])) {phi[k]=phi[i]*p[j];break;}
            phi[k]=phi[i]*(p[j]-1);
        }
    }phi[1]=1;
    fo(i,2,N) fo(j,1,N/i) a[i*j]+=phi[i]*j;
    fo(i,1,N) a[i]+=a[i-1];
    for(scanf("%d",&ty);ty;ty--) scanf("%d",&n),printf("%lld\n",a[n]);
}