bzoj 2818: Gcd

#include<cstdio>
#include<iostream>
#define M 10000009
#define ll long long
int tot,zhan[M],n,f[M];
ll sum[M],ans;
void ou()
{
    sum[1]=1;
    for(int i=2;i<=n;i++)
      {
        if(!f[i])
          {
            sum[i]=i-1;
            zhan[++tot]=i;
          }
        for(int j=1;j<=tot;j++)
          {
            if(zhan[j]*i>n)
              break;
            f[zhan[j]*i]=1;
            if(i%zhan[j])
              sum[i*zhan[j]]=sum[i]*sum[zhan[j]];
            else
              sum[i*zhan[j]]=sum[i]*zhan[j];
          }
      }
    return;
}
int main()
{
    scanf("%d",&n);
    ou();
    for(int i=1;i<=n;i++)
      sum[i]+=sum[i-1];
    for(int i=1;i<=tot;i++)
      ans+=2*sum[n/zhan[i]]-1;
    printf("%lld
",ans);
    return 0;
}

枚举每个质数，转化后求欧拉函数。