bzoj2818: Gcd

求一遍欧拉函数前缀和，直接就可以O(n)解决。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
bool check[10000005];int prime[3000005],fai[10000005],tot;long long s[10000005];
int main()
{
    int n;
    cin>>n;
    for(int i=2;i<=n;i++)
    {
        if(!check[i])
        {
            prime[tot++]=i;
            fai[i]=i-1;
        }
        for(int j=0;j<tot;j++)
        {
            if(i*prime[j]>n)
            break;
            check[i*prime[j]]=1;
            if(i%prime[j])
            {
                fai[i*prime[j]]=fai[i]*(prime[j]-1);
            }
            else
            {
                fai[i*prime[j]]=fai[i]*prime[j];
                break;
            }
        }
    }
    fai[1]=1;
    for(int i=1;i<=n;i++)
    {
        s[i]=s[i-1]+fai[i];
    }
    long long ans=0;
    for(int i=0;i<tot;i++)
    {
        ans+=(s[n/prime[i]]<<1)-1;
    }
    cout<<ans;
}