O - GCD - Extreme (II)(欧拉函数在gcd里的小运用)

O

emmm没思路。。。看题解。。。暴力会T飞的哈哈哈哈
gcd(i,j)=d得gcd(i/d,j/d)=1,就转化成phi(j/d)*d。故对于每个j：gcd(i,j)，该j对整个题目上的式子的贡献就是sigma(phi(j/d)*d),d是j的因子！！！
你以为到了这里对我这种智障就完了吗！！！
我竟然要外层循环for(i,1,N)一个一个的枚举出来每个的sigma。。。每个i会有好多不是i的因子的无效值。
为什么不是外层是for(d,1,N) 对于每个因子都把它对不同j的贡献算到j里面去。这样枚举每次都是有用的。
最有弄一个前缀和就好了。。

void init()
{
    for(int d=1;d<N-5;++d)//枚举因子d
        for(int j=d*2;j<=N-5;j+=d)//d是j的因子
            f[j]+=d*phi[j/d];//f[j]+=该因子d对j的贡献.
    for(int i=2;i<=N-5;++i)//前缀和
        ans[i]=ans[i-1]+f[i];
}

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#define en '\n'
typedef long long ll;
const int N=4e6+5;
using namespace std;
int phi[N],prime[N],tot;
bool is[N];
void getphi()
{
    tot=0,phi[1]=1;
    for(int i=2;i<=N-5;++i)
    {
        if(!is[i]) prime[++tot]=i,phi[i]=i-1;
        for(int j=1;j<=tot;++j)
        {
            if(i*prime[j]>N-5) break;
            is[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
}
ll f[N+5],ans[N+5];
void init()
{
    for(int d=1;d<N-5;++d)
        for(int j=d*2;j<=N-5;j+=d)
            f[j]+=d*phi[j/d];
    for(int i=2;i<=N-5;++i)
        ans[i]=ans[i-1]+f[i];
}
int main()
{
    getphi();init();
    int n;
    while(~scanf("%d",&n),n)
        printf("%lld\n",ans[n]);
}