O - GCD - Extreme (II)

题目链接https://vjudge.net/contest/233487#problem/O

思路：将GCD（n,x）为一个值得分类整理计算。

易知G[n]=G[n-1]+gcd(1,n)+gcd(2,n)+....gcd(n-1,n)

令gcd(1,n)+gcd(2,n)+....gcd(n-1,n)=A（n），对A（n）而言gcd（x,n） 1<=x<=n  若其值为d，则同除d ，gcd（x/d,n/d）=1。

这样的数的个数可以用欧拉函数计算出，故A（n）为 ∑phi(n/i)*i

 

---

Code：

#include<bits/stdc++.h>
#define LL long long
#define maxn 4000001
LL phi[maxn],gcd[maxn],res[maxn];
using namespace std;
void init()
{
    memset(phi,0,sizeof(phi));
    memset(gcd,0,sizeof(gcd));
    memset(res,0,sizeof(res));
    int i,j,index=0,x;
    phi[1]=1;
    for(i=2; i<=maxn; i++)
        if(!phi[i])
            for( j=i; j<=maxn; j+=i)
            {
                if(!phi[j])
                    phi[j]=j;
                phi[j]=phi[j]/i*(i-1);  //预处理欧拉函数值
            }
    for(i = 1; i<=maxn; i++)
        for(j = 2*i; j<=maxn; j+=i)
            gcd[j]+=phi[j/i]*i;         //求A（n）即gcd[i]

    for(i = 2; i<=maxn; i++)
        res[i]=res[i-1]+gcd[i];
}
int main()
{
    int n;
    init();
    while(scanf("%lld",&n)&&n)
        printf("%lld\n",res[n]);
}