uva 11426 GCD - Extreme (II) 欧拉函数

https://vjudge.net/problem/UVA%20-%2011426

题目大意：如上图所示，就是让你求那个式子的结果。

思路：不妨引入一个新函数$f$，且有$f(n)={\sum }_{i=1}^{n-1}gcd(n,i)$，那么对于任意给定的数$n$，答案就是${\sum }_{i=2}^{n}f(i)$。现在我们思考一下如何求$f(n)$，考虑枚举$n$的因子，假设某个因子为$i$，那么只要求出$gcd(n,j)=i,1<=j<=n$的个数就可以统计这个因子的贡献了，由$gcd(n,j)=i$可得$gcd(n/i,j/i)=1$，因此只要求出$gcd(n/i,j/i)=1,1<=j<=n$的个数就可以了，显然这个值等于$\phi (n/i)$，那么这道题就解决了。

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;

const int maxn=4e6+5;

int prime[maxn];
bool vis[maxn];
ll phi[maxn],f[maxn],s[maxn];
int k=-1,n;

void euler()
{
    for(int i=2;i<=n;i++)
    {
        if(!vis[i])
        {
            prime[++k]=i;
            phi[i]=i-1;
        }
        for(int j=0;j<=k&&i*prime[j]<=n;j++)
        {
            vis[prime[j]*i]=1;
            if(i%prime[j]==0)
            {
                phi[prime[j]*i]=phi[i]*prime[j];
                break;
            }
            else
                phi[prime[j]*i]=phi[i]*phi[prime[j]];
        }
    }
}

int main()
{
    n=4e6;
    int nn=n<<1;
    euler();
    for(int i=1;i<=nn;i++)//枚举因子
        for(int j=i<<1;j<=n;j+=i)
            f[j]+=i*phi[j/i];
    for(int i=2;i<=n;i++)
        s[i]=s[i-1]+f[i];
    while(~scanf("%d",&n)&&n)
        printf("%lld\n",s[n]);
    return 0;
}