bzoj2818: Gcd懵逼乌斯反演

由于是单组数据，强行不分块O(n)过

线性筛部分非常神奇，用了一个奇妙的推导（懒得写了）

#include <bits/stdc++.h>
using namespace std;
int n,m,mu[10000001],f[10000001],p[10000001];bool o[10000001];
int main()
{
    scanf("%d",&n);
    mu[1]=1;
    for(int i=2;i<=n;i++)
    {
        if(!o[i]) p[++m]=i,mu[i]=-1,f[i]=1;
        for(int j=1;j<=m;j++)
        if((long long)p[j]*i<=(long long)n)
        {
            o[p[j]*i]=1;
            if(i%p[j]==0) f[p[j]*i]=mu[i];
            else f[p[j]*i]=mu[i]-f[i];
            if(i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
        else break;
    }
    long long ret=0;
    for(int i=1;i<=n;i++)
        ret+=(long long)(n/i)*(n/i)*f[i];
    printf("%lld\n",ret);
    return 0;
}